Time stability of soil water storage measured by neutron probe and the effects of calibration procedures in a small watershed

The knowledge of soil water storage (SWS) of soil profiles is crucial for the adoption of vegetation restoration practices. With the aim of identifying representative sites to obtain the mean SWS of a watershed, a time stability analysis of neutron probe evaluations of SWS was performed by the means of relative differences and Spearman rank correlation coefficients. At the same time, the effects of different neutron probe calibration procedures were explored on time stability analysis. mean SWS estimation. and preservation of the spatial variability of SWS. The selected watershed, with deep gullies and undulating slopes which cover an area of 20 ha, is characterized by an Ust-Sandiic Entisol and an Aeolian sandy soil. The dominant vegetation species are bunge needlegrass (Stipa bungeana Trim) and korshinsk peashrub (Carugano Korshinskii kom.). From June 11, 2007 to July 23,2008, SWS of the top1 m soil layer was evaluated for 20 dates, based on neutron probe data of 12 sampling sites. Three calibration procedures were employed: type 1, most complete, with each site having its own linear calibration equation (TrE); type II. with TrE equations extended over the whole field: and type III, with one single linear calibration curve for the whole field (UnE) and also correcting its intercept based on site specific relative difference analysis (RdE) and on linear fitting of data (RcE), both maintaining the same slope. A strong time stability of SWS estimated by TrE equations was identified. Soil particle size and soil organic matter content were recognized as the influencing factors for spatial variability of SWS. Land use influenced neither the spatial variability nor the time stability of SWS. Time stability analysis identified one site to represent the mean SWS of the whole watershed with mean absolute percentage errors of less than 10%, therefore. this site can be used as a predictor for the mean SWS of the watershed. Some equations of type II were found to be unsatisfactory to yield reliable mean SWS values or in preserving the associated soil spatial variability. Hence, it is recommended to be cautious in extending calibration equations to other sites since they might not consider the field variability. For the equations with corrected intercept (type III), which consider the spatial variability of calibration in a different way in relation to TrE, it was found that they can yield satisfactory means and standard deviation of SWS, except for the RdE equations, which largely leveled off the SWS values in the watershed. Correlation analysis showed that the neutron probe calibration was linked to soil bulk density and to organic matter content. Therefore, spatial variability of soil properties should be taken into account during the process of neutron probe calibration. This study provides useful information on the mean SWS observation with a time stable site and on distinct neutron probe calibration procedures, and it should be extended to soil water management studies with neutron probes, e.g., the process of vegetation restoration in wider area and soil types of the Loess Plateau in China. (C) 2009 Elsevier B.V. All rights reserved.

New Analytic Solution of Boltzmann Transform for Horizontal Water Infiltration into Sand

The use of the Boltzmann transform function, lambda(theta), to solve the Richards equation when the diffusivity, D, is a function of only soil water content,., is now commonplace in the literature. Nevertheless, a new analytic solution of the Boltzmann transform lambda(h) as a function of matric potential for horizontal water infiltration into a sand was derived without invoking the concept or use of D(theta). The derivation assumes that a similarity exists between the soil water retention function and the Boltzmann transform lambda(theta). The solution successfully described soil water content profiles experimentally measured for different infiltration times into a homogeneous sand and agrees with those presented by Philip in 1955 and 1957. The applicability of this solution for all soils remains open, but it is anticipated to hold for soils whose air-filled pore-size distribution before wetting is sufficiently narrow to yield a sharp increase of water content at the wetting front during infiltration. It also improves and provides a versatile alternative to the well-known analysis pioneered by Green and Ampt in 1911.

Spatio-temporal variability behavior of land surface soil water content in shrub- and grass-land

The knowledge of the relationship between spatial variability of the surface soil water content (theta) and its mean across a spatial domain (theta(m)) is crucial for hydrological modeling and understanding soil water dynamics at different scales. With the aim to compare the soil moisture dynamics and variability between the two land uses and to explore the relationship between the spatial variability of theta and theta(m), this study analyzed sets of surface theta measurements performed with an impedance soil moisture probe, collected 136 times during a period of one year in two transects covering different land uses, i.e., korshinsk peashrub transect (KPT) and bunge needlegrass transect (BNT), in a watershed of the Loess Plateau, China. Results showed that the temporal pattern of theta behaved similarly for the two land uses, with both relative wetter soils during wet period and relative drier soils during dry period recognized in BNT. Soil moisture tended to be temporally stable among different dates, and more stable patterns could be observed for dates with more similar soil water conditions. The magnitude of the spatial variation of theta in KPT was greater than that in ENT. For both land uses, the standard deviation (SD) of theta in general increased as theta(m) increased, a behavior that could be well described with a natural logarithmic function. Convex relationship of CV and theta(m) and the maximum CV for both land uses (43.5% in KPT and 41.0% in BNT) can, therefore, be ascertained. Geostatistical analysis showed that the range in KPT (9.1 m) was shorter than that in BNT (15.1 m). The nugget effects, the structured variability, hence the total variability increased as theta(m) increased. For both land uses, the spatial dependency in general increased with increasing theta(m). 2011 Elsevier B.V. All rights reserved.

Changes in wood-water relationship due to heat treatment assessed on micro-samples of three Eucalyptus species

Wood-water relationship of untreated and heat-treated wood was studied. Specimens of Eucalyptus grandis, E saligna, and E citriodora were submitted to five conditions of heat treatment: 180 degrees C and 220 degrees C with air; 220 degrees C, 250 degrees C, and 280 degrees C with N(2). The wood-water relationships were accurately studied in a special device, in which the moisture content (MC) of the sample was measured with a highly sensitive electronic microbalance placed in a climatic chamber. The dimensions of the sample were collected continuously without contact by means of two high-speed laser scan micrometers. Sorption curves and shrinkage-MC relationships were observed. To study the effects of heat treatment, the following parameters were also determined: fiber saturation point (FPS), wood anisotropy (T/R ratio), shrinkage slope, reduction in hygroscopicity, and anti-shrink efficiency (ASE). The physical properties were significantly affected only at 220 degrees C and above. At heat temperature levels higher than 220 degrees C, the reduction in hygroscopicity and ASE are higher than 40% and continue to be reduced with increasing temperature level. This work also demonstrates that heat treatment does not change the slope of the curves shrinkage vs. MC, proving that heat treatment affects the domain of alterations in wood properties, but not the behavior within this domain.

Production and carbon allocation in a clonal Eucalyptus plantation with water and nutrient manipulations

We examined resource limitations on growth and carbon allocation in a fast-growing, clonal plantation of Eucalyptus grandis x urophylla in Brazil by characterizing responses to annual rainfall, and response to irrigation and fertililization for 2 years. Productivity measures included gross primary production (GPP), total belowground carbon allocation (TBCA), bole growth, and net ecosystem production (NEP). Replicate plots within a single plantation were established at the midpoint of the rotation (end of year 3), with treatments of no additional fertilization or irrigation, heavy fertilization (to remove any nutrient limitation), irrigation (to remove any water limitation), and irrigation plus fertilization. Rainfall was unusually high in the first year (1769mm) of the experiment, and control plots had high rates of GPP (6.64 kg C m(-2) year(-1)), TBCA (2.14 kg C m(-2) year(-1)), and bole growth (1.81 kg C m(-2) year). Irrigation increased each of these rates by 15-17%. The second year of the experiment had average rainfall (1210 mm), and lower rainfall decreased production in control plots by 46% (GPP), 52% (TBCA), and 40% (bole growth). Fertilization treatments had neglible effects. The response to irrigation was much greater in the drier year, with irrigated plots exceeding the production in control plots by 83% (GPP), 239% (TBCA), and 24% (bole growth). Even though the rate of irrigation ensured no water limitation to tree growth, the high rainfall year showed higher production in irrigated plots for both GPP (38% greater than in drier year) and bole growth (23% greater). Varying humidity and supplies of water led to a range in NEP of 0.8-2.7 kg C m(-2) year. This difference between control and irrigated treatments, combined with differences between drier and wetter years, indicated a strong response of these Eucalyptus trees to both water supply and atmospheric humidity during the dry season. The efficiency of converting light energy into fixed carbon ranged from a low of 0.027 mol C to a high of 0.060 mol C per mol of absorbed photosynthetically active radiation (APAR), and the efficiency of bolewood production ranged from 0.78 to 1.98 g wood per MJ of APAR. Irrigation increased the efficiency of wood production per unit of water used from 2.55 kg wood m(-3) in the rainfed plot to 3.51 kg m(-3) in irrigated plots. Detailed information on the response of C budgets to environmental conditions and resource supplies will be necessary for accurate predictions of plantation yields across years and landscapes. (V) 2007 Elsevier B.V. All rights reserved.

PCJ River Basins` Water Availability Caused by Water Diversion Scenarios to Supply Metropolitan Areas of So Paulo

A dynamic systems simulation model of water resources was developed as a tool to help analyze alternatives to water resources management for the Piracicaba, Capivari and Jundiai River Water Basins (RB-PCJ), and used to run six 50-year simulations from 2004 to 2054. The model estimates water supply and demand, as well as contamination load by several consumers. Six runs were performed using a constant mean precipitation value, changing water supply and demand and different volumes diverted from RB-PCJ to RB-Alto Tiet. For the Business as Usual scenario, the Sustainability Index went from 0.44 in 2004 to 0.20 by 2054. The Water Sustainability Index changed from 74% in 2004 to 131% by 2054. The Falkenmark Index changed from 1,403 m(3) person (-aEuro parts per thousand 1) year (-aEuro parts per thousand 1) in 2004 to 734 m(3) person (-aEuro parts per thousand 1) year (-aEuro parts per thousand 1) by 2054. We concluded that sanitation is one of the major problems for the PCJ River Basins.

Water Resources Assessment at Piracicaba, Capivari and Jundiai River Basins: A Dynamic Systems Approach

The Piracicaba, Capivari, and Jundiai River Basins (RB-PCJ) are mainly located in the State of So Paulo, Brazil. Using a dynamics systems simulation model (WRM-PCJ) to assess water resources sustainability, five 50-year simulations were run. WRM-PCJ was developed as a tool to aid decision and policy makers on the RB-PCJ Watershed Committee. The model has 254 variables. The model was calibrated and validated using available information from the 80s. Falkenmark Water Stress Index went from 1,403 m(3) person (-aEuro parts per thousand 1) year (-aEuro parts per thousand 1) in 2004 to 734 m(3) P (-aEuro parts per thousand 1) year (-aEuro parts per thousand 1) in 2054, and Xu Sustainability Index from 0.44 to 0.20. In 2004, the Keller River Basin Development Phase was Conservation, and by 2054 was Augmentation. The three criteria used to evaluate water resources showed that the watershed is at crucial water resources management turning point. The WRM-PCJ performed well, and it proved to be an excellent tool for decision and policy makers at RB-PCJ.

WATER RESOURCES SITUATION AT PIRACICABA, CAPIVARI AND JUNDIAI WATERSHEDS USING A DYNAMIC SYSTEMS MODEL

Using a dynamic systems model specifically developed for Piracicaba, Capivari and Jundia River Water Basins (BH-PCJ) as a tool to help to analyze water resources management alternatives for policy makers and decision takers, five simulations for 50 years timeframe were performed. The model estimates water supply and demand, as well as wastewater generation from the consumers at BH-PCJ. A run was performed using mean precipitation value constant, and keeping the actual water supply and demand rates, the business as usual scenario. Under these considerations, it is expected an increment of about similar to 76% on water demand, that similar to 39% of available water volume will come from wastewater reuse, and that waste load increases to similar to 91%. Falkenmark Index will change from 1,403 m(3) person(-1) year(-1) in 2004, to 734 m(3) P(-1) year(-1) by 2054, and the Sustainability Index from 0.44 to 0.20. Another four simulations were performed by affecting the annual precipitation by 90 and 110%; considering an ecological flow equal to 30% of the mean daily flow; and keeping the same rates for all other factors except for ecological flow and household water consumption. All of them showed a tendency to a water crisis in the near future at BH-PCJ.

Root Water Extraction under Combined Water and Osmotic Stress

Using a numerical implicit model for root water extraction by a single root in a symmetric radial flow problem, based on the Richards equation and the combined convection-dispersion equation, we investigated some aspects of the response of root water uptake to combined water and osmotic stress. The model implicitly incorporates the effect of simultaneous pressure head and osmotic head on root water uptake, and does not require additional assumptions (additive or multiplicative) to derive the combined effect of water and salt stress. Simulation results showed that relative transpiration equals relative matric flux potential, which is defined as the matric flux potential calculated with an osmotic pressure head-dependent lower bound of integration, divided by the matric flux potential at the onset of limiting hydraulic conditions. In the falling rate phase, the osmotic head near the root surface was shown to increase in time due to decreasing root water extraction rates, causing a more gradual decline of relative transpiration than with water stress alone. Results furthermore show that osmotic stress effects on uptake depend on pressure head or water content, allowing a refinement of the approach in which fixed reduction factors based on the electrical conductivity of the saturated soil solution extract are used. One of the consequences is that osmotic stress is predicted to occur in situations not predicted by the saturation extract analysis approach. It is also shown that this way of combining salinity and water as stressors yields results that are different from a purely multiplicative approach. An analytical steady state solution is presented to calculate the solute content at the root surface, and compared with the outputs of the numerical model. Using the analytical solution, a method has been developed to estimate relative transpiration as a function of system parameters, which are often already used in vadose zone models: potential transpiration rate, root length density, minimum root surface pressure head, and soil theta-h and K-h functions.

Macroscopic root water uptake distribution using a matric flux potential approach

Hydrological models featuring root water uptake usually do not include compensation mechanisms such that reductions in uptake from dry layers are compensated by an increase in uptake from wetter layers. We developed a physically based root water uptake model with an implicit compensation mechanism. Based on an expression for the matric flux potential (M) as a function of the distance to the root, and assuming a depth-independent value of M at the root surface, uptake per layer is shown to be a function of layer bulk M, root surface M, and a weighting factor that depends on root length density and root radius. Actual transpiration can be calculated from the sum of layer uptake rates. The proposed reduction function (PRF) was built into the SWAP model, and predictions were compared to those made with the Feddes reduction function (FRF). Simulation results were tested against data from Canada (continuous spring wheat [(Triticum aestivum L.]) and Germany (spring wheat, winter barley [Hordeum vulgare L.], sugarbeet [Beta vulgaris L.], winter wheat rotation). For the Canadian data, the root mean square error of prediction (RMSEP) for water content in the upper soil layers was very similar for FRF and PRF; for the deeper layers, RMSEP was smaller for PRF. For the German data, RMSEP was lower for PRF in the upper layers and was similar for both models in the deeper layers. In conclusion, but dependent on the properties of the data sets available for testing,the incorporation of the new reduction function into SWAP was successful, providing new capabilities for simulating compensated root water uptake without increasing the number of input parameters or degrading model performance.

Optimal irrigation depth for bean crops considering water constraints and farmer`s risk aversion level

Nowadays, the rising competition for the use of water and environmental resources with consequent restrictions for farmers should change the paradigm in terms of irrigation concepts, or rather, in order to attain economical efficiency other than to supply water requirement for the crop. Therefore, taking into account the social and economical role of bean activity in Brazil, as well as the risk inherent to crop due to its high sensibility to both deficit and excessive water, the optimization methods regarding to irrigation management have become more interesting and essential. This study intends to present a way to determine the optimal water supply, considering different combinations between desired bean yield and level of risk, bringing as a result a graph with the former associated with the latter, depending on different water depths.

Simulation of irrigated Thanzania grass growth based on photothermal units, nitrogen fertilization and water availability

Simulation of irrigated Thanzania grass growth based on photothermal units, nitrogen fertilization and water availability. The mathematical model to predict the forage yield using photothennal units was utilized with success in Elephant grass, Thanzania and Brachiaria niziziensis in the absence of water stress and nitrogen stress. The aim of this study was to propose models to estimate the forage yield of Thanzania grass under different irrigation (25, 50,75, 100 e 125% of ETc) and nitrogen level in various regions of Brazil. As such, models were developed to estimate the dry matter production of Panicum maximum Jacq. frass cv Thanzania in different irrigation and nitrogen levels, using photothermal units. The models were adjusted to doses of 0, 30, 60, 110 and 270 kg of N ha(-1), doses were divided in applications after each evaluation, with a rest cycle of 35 days. The adjusted model presented good performance in predicting dry matter production of Thanzania grass, with r(2) = 0.9999. The results made it possible to verify that the proposed model can be used to predict forage production in different regions of Brazil. It can be estimated, with good precision. The production of Thanzania grass dry matter can be accurately estimated in specific places (in function of latitude and time of year), with the maximum and minimum temperature values.

Ecophysiological responses of water hyacinth exposed to Cr(3+) and Cr(6+)

Due to its wide industrial use, chromium (Cr) is considered a serious environmental pollutant of aquatic bodies. in order to investigate the ecophysiological responses of water hyacinth [Eichhornia crassipes (Mart.) Solms] to Cr treatment, plants were exposed to 1 and 10 mM Cr(2)O(3) (Cr(3+)) and K(2)Cr(2)O(7) (Cr(6+)) concentrations for two or 4 days in a hydroponic system. Plants exposed to the higher concentration of Cr(6+) for 4 days did not survive, whereas a 2 days treatment with 1 mM Cr(3+) apparently stimulated growth. Analysis of Cr uptake indicated that most of the Cr accumulated in the roots, but some was also translocated and accumulated in the leaves. However, in plants exposed to Cr(6+) (1 mM), a higher translocation of Cr from roots to shoots was observed. it is possible that the conversion from Cr(6+) to Cr(3+), which immobilizes Cr in roots, was not total due to the presence of Cr(6+), causing deleterious effects on gas exchange, chlorophyll a fluorescence and photosynthetic pigment contents. Chlorophyll a was more sensitive to Cr than chlorophyll b. Cr(3+) was shown to be less toxic than Cr(6+) and, in some cases even increased photosynthesis and chlorophyll content. This result indicated that the F(v)/F(0) ratio was more effective than the F(v)/F(m) ratio in monitoring the development of stress by Cr(6+). There was a linear relationship between qP and F(v)/F(m). No statistical differences were observed in NPQ and chlorophyll a/b ratio, but there was a tendency to decrease these values with Cr exposure. This suggests that there were alterations in thylakoid stacking, which might explain the data obtained for gas exchanges and other chlorophyll a fluorescence parameters. (C) 2008 Elsevier B.V. All rights reserved.

Chlorophyll a fluorescence and ultrastructural changes in chloroplast of water hyacinth as indicators of environmental stress

Chlorophyll a fluorescence parameters and transmission electron microscopy (TEM) were used to assess the stress conditions in water hyacinth along the Paraiba do Sul River (PSR), an important River in southeastern Brazil. The data were obtained at the end of the dry season of 2005 and at the end of the wet season of 2006. Changes in F-o and F-m parameters were observed as differentiated responses, depending on the season. Non-photochemical dissipation (qN and NPQ) from plants was greater in the most industrialized region of the PSR in both seasons. However, F-v/F-m for all samples ranged between 0.77 and 0.81, showing that high maximum quantum yield was maintained. Although the F-v/F-m suggests that the plants were exhibiting normal photochemical activities, ultrastructural changes in chloroplasts showed thylakoids disorganization. Plants from the most industrialized region showed non-stacking grana thylakoids disposition. In spite of these alterations, the membrane integrity was maintained, suggesting an adaptation to adjustment to adverse environmental conditions. (C) 2008 Elsevier B.V. All rights reserved.

Validation of a root water uptake model to estimate transpiration constraints

Experimental results obtained from a greenhouse trial with common bean (Phaseolus vulgaris L) plants performed to test model hypotheses regarding the onset of limiting hydraulic conditions and the shape of the transpiration reduction curve in the falling rate phase are presented. According to these hypotheses based on simulations with an upscaled single-root model, the matric flux potential at the onset of limiting hydraulic conditions is as a function of root length density and potential transpiration rate, while the relative transpiration in the falling rate phase equals the relative matric flux potential. Transpiration of bean plants in water stressed pots with four different soils was determined daily by weighing and compared to values obtained from non-stressed pots. This procedure allowed determining the onset of the falling rate phase and corresponding soil hydraulic conditions. At the onset of the falling rate phase, the value of matric flux potential M(I) showed to differ in order of magnitude from the model predicted value for three out of four soils. This difference between model and experiment can be explained by the heterogeneity of the root distribution which is not considered by the model. An empirical factor to deal with this heterogeneity should be included in the model to improve predictions. Comparing the predictions of relative transpiration in the falling rate phase using a linear shape with water content, pressure head or matric flux potential, the matric flux potential based reduction function, in agreement with the hypothesis, showed the best performance, while the pressure head based equation resulted in the highest deviations between observed and predicted values of relative transpiration rates. (C) 2010 Elsevier B.V. All rights reserved.