Modeling long-term variability and change of soil moisture and groundwater level - from catchment to global scale

The water stored in and flowing through the subsurface is fundamental for sustaining human activities and needs, feeding water and its constituents to surface water bodies and supporting the functioning of their ecosystems. Quantifying the changes that affect the subsurface water is crucial for our understanding of its dynamics and changes driven by climate change and other changes in the landscape, such as in land-use and water-use. It is inherently difficult to directly measure soil moisture and groundwater levels over large spatial scales and long times. Models are therefore needed to capture the soil moisture and groundwater level dynamics over such large spatiotemporal scales. This thesis develops a modeling framework that allows for long-term catchment-scale screening of soil moisture and groundwater level changes. The novelty in this development resides in an explicit link drawn between catchment-scale hydroclimatic and soil hydraulics conditions, using observed runoff data as an approximation of soil water flux and accounting for the effects of snow storage-melting dynamics on that flux. Both past and future relative changes can be assessed by use of this modeling framework, with future change projections based on common climate model outputs. By direct model-observation comparison, the thesis shows that the developed modeling framework can reproduce the temporal variability of large-scale changes in soil water storage, as obtained from the GRACE satellite product, for most of 25 large study catchments around the world. Also compared with locally measured soil water content and groundwater level in 10 U.S. catchments, the modeling approach can reasonably well reproduce relative seasonal fluctuations around long-term average values. The developed modeling framework is further used to project soil moisture changes due to expected future climate change for 81 catchments around the world. The future soil moisture changes depend on the considered radiative forcing scenario (RCP) but are overall large for the occurrence frequency of dry and wet events and the inter-annual variability of seasonal soil moisture. These changes tend to be higher for the dry events and the dry season, respectively, than for the corresponding wet quantities, indicating increased drought risk for some parts of the world.

Spatial variability of soil properties and yield of a grazed alfalfa pasture in Brazil.

2016

Spatial variability of satured soil hydraulic conductivity in the region of Araguaia River - Brazil.

This study evaluates the spatial variability of saturated hydraulic conductivity in the soil in an area of 51,850 ha at the headwaters of the Araguaia River MT/GO. This area is highly vulnerable because it is a location of recharging through natural water infiltration of the Guarani Aquifer System and an area of intense increases in agriculture since its adoption by growers in the last 30 years. Soil samples were collected at 383 points, geographically located by GPS. The samples were collected from depths of 0 - 20 cm and 60 - 80 cm. Exploratory statistics and box-plot were used in the descriptive analysis and semivariogram were constructed to determine the spatial model. The exploratory analysis showed that the mean hydraulic conductivity in the superficial layer was less than at the level of 60-80 cm; however, the greatest variability evaluated with a coefficient of variation also was from this layer. Data tended towards a normal distribution. These results can be explained by the greater soil compaction in the superficial layer. The semivariogram models, adjusted for the two layers, were exponential and demonstrated moderate and strong dependence, with ranges of 5000 and 3000 utm respectively. It was concluded that soil use is influencing the spatial distribution model of the hydraulic conductivity in the region.

Spatial variability of satured soil hydraulic conductivity in the region of Araguaia River - Brazil.

2008

Do plant species influence soil CO(2) and N(2)O fluxes in a diverse tropical forest?

To test whether plant species influence greenhouse gas production in diverse ecosystems, we measured wet season soil CO(2) and N(2)O fluxes close to similar to 300 large (>35 cm in diameter at breast height (DBH)) trees of 15 species at three clay-rich forest sites in central Amazonia. We found that soil CO(2) fluxes were 38% higher near large trees than at control sites >10 m away from any tree (P < 0.0001). After adjusting for large tree presence, a multiple linear regression of soil temperature, bulk density, and liana DBH explained 19% of remaining CO(2) flux variability. Soil N(2)O fluxes adjacent to Caryocar villosum, Lecythis lurida, Schefflera morototoni, and Manilkara huberi were 84%-196% greater than Erisma uncinatum and Vochysia maxima, both Vochysiaceae. Tree species identity was the most important explanatory factor for N(2)O fluxes, accounting for more than twice the N(2)O flux variability as all other factors combined. Two observations suggest a mechanism for this finding: (1) sugar addition increased N(2)O fluxes near C. villosum twice as much (P < 0.05) as near Vochysiaceae and (2) species mean N(2)O fluxes were strongly negatively correlated with tree growth rate (P = 0.002). These observations imply that through enhanced belowground carbon allocation liana and tree species can stimulate soil CO(2) and N(2)O fluxes (by enhancing denitrification when carbon limits microbial metabolism). Alternatively, low N(2)O fluxes potentially result from strong competition of tree species with microbes for nutrients. Species-specific patterns in CO(2) and N(2)O fluxes demonstrate that plant species can influence soil biogeochemical processes in a diverse tropical forest.

Spatio-temporal patterns of throughfall and solute deposition in an open tropical rain forest

The brief interaction of precipitation with a forest canopy can create a high spatial variability of both throughfall and solute deposition. We hypothesized that (i) the variability in natural forest systems is high but depends on system-inherent stability, (ii) the spatial variability of solute deposition shows seasonal dynamics depending on the increase in rainfall frequency, and (iii) spatial patterns persist only in the short-term. The study area in the north-western Brazilian state of Rondonia is subject to a climate with a distinct wet and dry season. We collected rain and throughfall on an event basis during the early wet season (n = 14) and peak of the wet season (n = 14) and analyzed the samples for pH and concentrations of NH4+, Na+, K+, Ca2+ Mg2+,, Cl-, NO3-, SO42- and DOC. The coefficient 3 4 cient of variation for throughfall based on both sampling intervals was 29%, which is at the lower end of values reported from other tropical forest sites, but which is higher than in most temperate forests. Coefficients of variation of solute deposition ranged from 29% to 52%. This heterogeneity of solute deposition is neither particularly high nor particularly tow compared with a range of tropical and temperate forest ecosystems. We observed an increase in solute deposition variability with the progressing wet season, which was explained by a negative correlation between heterogeneity of solute deposition and antecedent dry period. The temporal stability of throughfall. patterns was Low during the early wet season, but gained in stability as the wet season progressed. We suggest that rapid plant growth at the beginning of the rainy season is responsible for the lower stability, whereas less vegetative activity during the later rainy season might favor the higher persistence of ""hot"" and ""cold"" spots of throughfall. quantities. The relatively high stability of throughfall patterns during later stages of the wet season may influence processes at the forest floor and in the soil. Solute deposition patterns showed less clear trends but all patterns displayed a short-term stability only. The weak stability of those patterns is apt to impede the formation of solute deposition -induced biochemical microhabitats in the soil. (C) 2008 Elsevier B.V. All rights reserved.

Influence of the microstructure in the collapse of a residual clayey tropical soil

A typical residual clayey soil originating from basalt in southern Brazil has been analyzed in order to assess the influence of wetting-induced deformation and microstructural features on the collapse behavior. Single and double oedometer tests were undertaken on a soil profile to 9 m depth. The results indicated collapsible behaviour at all profile depths. The influence of pre-consolidation stress and pedogenetic factors in the variability of the physical characteristics of the soil and in the magnitude of the collapse was noted. The collapse coefficient has been shown to be related to the both the microaggregate plasma and the varying nature of the pores and their interconnectivity.

Changes in optical properties caused by UV-irradiation of aquatic humic substances from the amazon river basin: Seasonal variability evaluation

Aquatic humic substances (AHS) isolated from two characteristic seasons of the Negro river, winter and summer corresponding to floody and dry periods, were structurally characterized by (13)C nuclear magnetic ressonance. Subsequently, AHS aqueous solutions were irradiated with a polychromatic lamp (290-475 nm) and monitored by its total organic carbon (TOC) content, ultraviolet-visible (UV-vis) absorbance, fluorescence and Fourier transformed infrared spectroscopy (FTIR). As a result, a photobleaching upto 80% after irradiation of 48 h was observed. Conformational rearrangements and formation of low molecular complexity structures were formed during the irradiation, as deduced from the pH decrement and the fluorescence shifting to lower wavelengths. Additionally a significant mineralization with the formation Of CO(2), CO, and inorganic carbon compounds was registered, as assumed by TOC losses of up to 70%. The differences in photodegradation between samples expressed by photobleaching efficiency were enhanced in the summer sample and related to its elevated aromatic content. Aromatic structures are assumed to have high autosensitization capacity effects mediated by the free radical generation from quinone and phenolic moieties.

Watershed scale temporal stability of soil water content

The recognition of temporally stable locations with respect to soil water content is of importance for soil water management decisions, especially in sloping land of watersheds. Neutron probe soil water content (0 to 0.8 m), evaluated at 20 dates during a year in the Loess Plateau of China, in a 20 ha watershed dominated by Ust-Sandiic Entisols and Aeolian sandy soils, were used to define their temporal stability through two indices: the standard deviation of relative difference (SDRD) and the mean absolute bias error (MABE). Specific concerns were (a) the relationship of temporal stability with soil depth, (b) the effects of soil texture and land use on temporal stability, and (c) the spatial pattern of the temporal stability. Results showed that temporal stability of soil water content at 0.2 m was significantly weaker than those at the soil depths of 0.6 and 0.8 m. Soil texture can significantly (P<0.05) affect the stability of soil water content except for the existence of an insignificant difference between sandy loam and silt loam textures, while temporal stability of areas covered by bunge needlegrass land was not significantly different from those covered by korshinsk peashrub. Geostatistical analysis showed that the temporal stability was spatially variable in an organized way as inferred by the degree of spatial dependence index. With increasing soil depth, the range of both temporal stability indices showed an increasing trend, being 65.8-120.5 m for SDRD and 148.8-214.1 m for MABE, respectively. This study provides a valuable support for soil water content measurements for soil water management and hydrological applications on sloping land areas. (C) 2010 Elsevier B.V. All rights reserved.

Using a New Criterion to Identify Sites for Mean Soil Water Storage Evaluation

Establishing a few sites in which measurements of soil water storage (SWS) are time stable significantly reduces the efforts involved in determining average values of SWS. This study aimed to apply a new criterion the mean absolute bias error (MABE)-to identify temporally stable sites for mean SWS evaluation. The performance of MABE was compared with that of the commonly used criterion, the standard deviation of relative difference (SDRD). From October 2004 to October 2008, SWS of four soil layers (0-1.0, 1.0-2.0,2.0-3.0, and 3.0-4.0 m) was measured, using a neutron probe, at 28 sites on a hillslope of the Loess Plateau, China. A total of 37 SWS data sets taken over time were divided into two subsets, the first consisting of 22 dates collected during the calibration period from October 2004 to September 2006, and the second with 15 dates collected during the validation period from October 2006 to October 2008. The results showed that if a critical value of 5% for MABE was defined, more than half the sites were temporally stable for both periods, and the number of temporally stable sires generally increased with soil depth. Compared with SDRD, MABE was more suitable for the identification of time-stable sites for mean SS prediction. Since the absolute prediction error of drier sites is more sensitive to changes in relative difference in terms of mean SWS prediction, the sites of wet sectors should be preferable for mean SWS prediction for the same changes in relative difference.

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.

Characterization of fungal soil communities by F-RISA and arbuscular mycorrhizal fungi from Araucaria angustifolia forest soils after replanting and wildfire disturbances

The Fungal Ribosomal Intergenic Spacer Analysis (F-RISA) was used to characterize soil fungal communities from three ecosystems of Araucaria angustifolia from Brazil: a native forest and two replanted forest ecosystems, one of them with a past history of wildfire. The arbuscular mycorrhizal fungi (AMF) infection was evaluated in Araucaria roots of 18-month-old axenic plants previously inoculated with soils collected from those areas in a greenhouse experiment. The principal component analysis of F-RISA profiles showed different soil fungal community between the three studied areas. Sixty three percent of F-RISA fragments amplified in the soil and the substrate samples presented lengths between 500 and 700 bp. The number of Operational Taxonomic Units (OTUs) was 34 for soil and 38 for substrate, however, more fragments were detected in soil (214) than in substrate (163). An in silico F-RISA analysis to compare our data with ITS1-5.8S-ITS2 sequences from NCBI database showed the presence of Ascomycota, Basidiomycota and Glomeromycota among the soil and substrate fungal communities. AMF infection was higher in plants inoculated with soil from the native forest and the replanted forest with wildfire, both presenting similar chemical characteristics but with different disturbance levels. These results indicate that soil chemical composition may influence the soil fungal community structures rather than the anthropogenic or fire disturbances.

Soil CO(2) emission estimated by different interpolation techniques

Soil CO(2) emissions are highly variable, both spatially and across time, with significant changes even during a one-day period. The objective of this study was to compare predictions of the diurnal soil CO(2) emissions in an agricultural field when estimated by ordinary kriging and sequential Gaussian simulation. The dataset consisted of 64 measurements taken in the morning and in the afternoon on bare soil in southern Brazil. The mean soil CO(2) emissions were significantly different between the morning (4.54 mu mol m(-2) s(-1)) and afternoon (6.24 mu mol m(-2) s(-1)) measurements. However, the spatial variability structures were similar, as the models were spherical and had close range values of 40.1 and 40.0 m for the morning and afternoon semivariograms. In both periods, the sequential Gaussian simulation maps were more efficient for the estimations of emission than ordinary kriging. We believe that sequential Gaussian simulation can improve estimations of soil CO(2) emissions in the field, as this property is usually highly non-Gaussian distributed.

Contribution of soil organic carbon to the ion exchange capacity of tropical soils

Highly weathered soils represent about 3 billion ha of the tropical region. Oxisols represent about 60% of the Brazilian territory (more than 5 million km 2), in areas of great agricultural importance. Soil organic carbon (SOC) can be responsible for more than 80% of the cation exchange capacity (CEC) of highly weathered soils, such as Oxisols and Ultisols. The objective of this study was to estimate the contribution of the SOC to the CEC of Brazilian soils from different orders. Surface samples (0.0 to 0.2 m) of 30 uncultivated soils (13 Oxisols, 6 Ultisols, 5 Alfisols, 3 Entisols, I Histosol, 1 Inceptisol. and I Molisol), under native forests and from reforestation sites from Sao Paulo State, Brazil, were collected in order to obtain a large variation of (electro)chemical, physical, and mineralogical soil attributes. Total content of SOC was quantified by titulometric and colorimetric methods. Effective cation exchange capacity (ECEC) was obtained by two methods: the indirect method-summation-estimated the ECECi from the sum of basic cations (Ca+ Mg+ K+ Na) and exchangeable Al; and the direct ECECd obtained by the compulsive exchange method, using unbuffered BaCl2 solution. The contribution of SOC to the soil CEC was estimated by the Bennema statistical method. The amount of SOC var ied from 6.6 g kg(-1) to 213.4 g kg(-1). while clay contents varied from 40 g kg(-1) to 716 g kg(-1). Soil organic carbon contents were strongly associated to the clay contents, suggesting that clay content was the primary variable in controling the variability of SOC contents in the samples. Cation exchange capacity varied from 7.0 mmol(c) kg(-1) to 137.8 mmol(c) kg(-1) and had a positive Correlation with SOC. The mean contribution (per grain) of the SOC (1.64 mmol(c)) for the soil CEC was more than 44 times higher than the contribution of the clay fraction (0.04 mmol(c),). A regression model that considered the SOC content as the only significant variable explained 60% of the variation in the soil total CEC. The importance of SOC was related to soil pedogenetic process, since its contribution to the soil CEC was more evident in Oxisols with predominance of Fe and Al (oxihydr)oxides in the mineral fraction or in Ultisols, that presented illuviated clay. The influence of SOC in the sign and in the magnitude of the net charge of soils reinforce the importance of agricultural management systems that preserve high levels of SOC, in order to improve their sustainability.

Spatial patterns of soil attributes and components in a mangrove system in Southeast Brazil (Sao Paulo)

Purpose Among environmental factors governing innumerous processes that are active in estuarine environments, those of edaphic character have received special attention in recent studies. With the objectives of determining the spatial patterns of soil attributes and components across different mangrove forest landscapes and obtaining additional information on the cause-effect relationships between these variables and position within the estuary, we analyzed several soil attributes in 31 mangrove soil profiles from the state of So Paulo (Guaruja, Brazil). Materials and methods Soil samples were collected at low tide along two transects within the CrumahA(0) mangrove forest. Samples were analyzed to determine pH, Eh, salinity, and the percentages of sand, silt, clay, total organic carbon (TOC), and total S. Mineralogy of the clay fraction (< 2 mm) was also studied by X-ray diffraction analysis, and partitioning of solid-phase Fe was performed by sequential extraction. Results and discussion The results obtained indicate important differences in soil composition at different depths and landscape positions, causing variations in physicochemical parameters, clay mineralogy, TOC contents, and iron geochemistry. The results also indicate that physicochemical conditions may vary in terms of different local microtopographies. Soil salinity was determined by relative position in relation to flood tide and transition areas with highlands. The proportions of TOC and total S are conditioned by the sedimentation of organic matter derived from vegetation and by the prevailing redox conditions, which clearly favored intense sulfate reduction in the soils (similar to 80% of the total Fe is Fe-pyrite). Particle-size distribution is conditioned by erosive/deposition processes (present and past) and probably by the positioning of ancient and reworked sandy ridges. The existing physicochemical conditions appear to contribute to the synthesis (smectite) and transformation (kaolinite) of clay minerals. Conclusions The results demonstrate that the position of soils in the estuary greatly affects soil attributes. Differences occur even at small scales (meters), indicating that both edaphic (soil classification, soil mineralogy, and soil genesis) and environmental (contamination and carbon stock) studies should take such variability into account.