Data preparation - soil parameters

Hanjiang_soil.txt：

0.1 ° and GIS grid resampling:

VIC requires coordinates are latitude and longitude coordinates, so Arcgis processing data should be unified into geographic coordinates without projection coordinates.

1: SRTM90m effects include the study area [deg.] Is resampled to 0.1;
2: Get vector range larger study area,Cartography Tools -> Masking Tools -> Feature Outline Masks

3: greater than with the actual study shp mask rasters extraction zone; to ensure almost complete coverage of the study area;
4: After cutting raster to points; (since acquired by the DEM, so there attribute list for the latitude and longitude information and altitude information. empty information table may be: open the attribute table of points.shp, add filed recoding 'long' and 'lat' ( latitude and longitude), field selection right computational geometry, respectively the latitude and longitude information of each site)

5. .dbf file into txt or xlxs
6. Matlab

Soil properties (HWSD dataset + SPAW small software acquisition)

You must fill all the columns of the input file, but only need to define certain parameters upon the activation of the heat balance model or frozen. Only when FULL_ENERGY is true, red variables only need to have a meaningful value, otherwise they will be ignored ( but placeholder -9999 must still exist in this field, unless the document made clear that this field is optional ) . You can add a header row by ##.

Bulk_density

The first way:
soil data HWSD process flow:

1) Open the raster image (directly open, you can not view the attribute table) - Right - Properties - stretched by the percentage change without truncation - unique values - back raster image - Right - Open Attribute Table to see value, count property
2) open the attribute data --HWSD.dbf - HWSD-DATA.dbf - after conversion tool arcgis open --toolbox-- - excel-- turn table referred to as a roll-out excel hhhh - arcgis open hhhh
3) raster image - right - connection - select raster images hhhh value-- connection of mu_global-- test connections may waring reserved word (it does not matter) - connection - open the gate grid image attribute table, can be seen in all fields hhhh
after symbology - 4) show image - raster image right - properties - symbology - --value set to unique values like pH or attribute field sym90 field - primary display field - pH and the like, or selected attribute field sym90

NOTE: HWSD raster data is divided into two parts, and the attribute data, wherein the attribute value is displayed raster data mu_global attribute field.

Topsoil（0-30cm)
Subsoil（30-100cm）

Top soil bulk density and a second thin layer uses Topsoil data, the third data layer uses Subsoil.

Water is determined by HWSD: 519/577/689/690/730/731/763/769/770/807

The second way: based on their own research projects South interpolation :( here to select IDW)

Topsoil (survey sites around 0-20cm)
subsoil (survey point 20)

Introducing xy coordinate point - - --Field selection mode selected at the IDW bulk_density data - so that RSM by adjusting the parameters of the minimum - the first double-click showing the same extent and range of choice FIG web DJK_kuoda - and right data - raster (resolution 0.1): FIG environment third selection range DJK_kuoda consistent web - extraction mask - raster point - point data output

expt

exp index parameter, from Campbell exponential equation, = 3 + 2 / λ (or 3 + 2b), the value of not less than 3.

Wcr_FRACT

Soil The soil moisture as a parameter (the largest proportion of soil moisture, soil moisture = maximum porosity soil thickness *) at the critical point, it is the evaporation of soil, soil moisture saturated hydraulic conductivity when water begins to drop. Depending on the value of soil texture, it represents 70% of field capacity . (Field capacity moisture content here is under the -33KPa state) - net definition Officer

Wpwp_FRACT

The moisture in the soil as a soil parameter wilting point (the largest proportion of soil moisture, soil moisture = maximum porosity soil thickness *), wilting point is represented as soil moisture at 1500KPa, approximations of each soil texture. - the definition of official website

Ksat

Saturated hydraulic conductivity of soil each.
(Found in the study area: the subsurface of the whole is greater than the value of the surface, the reason?)

From the above variables SPAW - obtained swc_excel formula:

init_moistInitial

init_moistInitial moisture content of each layer can be set at any reasonable value. One approach is to use fractional soil moisture content (expressed as a fraction of the maximum soil moisture; max. soil moisture = porosity * layer depth) at the critical point, Wcr, which can be computed for each layer as a depth in millimeters by multiplying Wcr by the thickness of the layer in meters, and then multiplying by 1000.

avg_T

Soil heat damping depth (dp temperature at that depth remains constant) temperature, the average annual temperature is usually assumed to be the temperature calculated for all soil column bottom heat flux boundary.
(Meteorological data required average temperature is obtained in each grid)