Incorporation of spectral data into multivariate geostatistical models to map soil phosphorus variability in a Florida wetland

Hybrid geostatistical prediction methods incorporate (i) spatially-explicit soil observations and exhaustive grids of ancillary environmental variables (e.g. derived from remote sensing), (ii) spatial autocorrelation, (iii) spatial covariation, and/or (iv) combinations of the above. In numerous studies of terrestrial soils it has been shown that hybrid geostatistical methods outperform univariate spatial and regression (aspatial) methods. However, hybrid methods have rarely been employed to predict soil properties in wetlands. In this study we used spectral data and derived indices from two remote sensors (Landsat ETM+ and ASTER), with different spatial resolutions, from different seasons, but with similar spectral range, ancillary environmental data, as well as floc and soil total phosphorus (TP) observations from 111 sites. The specific objective of our study was to evaluate the performance of aspatial methods (multivariate regressions — REG), univariate spatial (Ordinary Kriging — OK) and hybrid/multivariate geostatistical methods (Regression Kriging — RK and Co-kriging — CK) in predicting the spatial variability and distribution of floc and soil TP in a subtropical wetland, WCA-2A, in the Florida Everglades. Measured floc TP ranged from 194 to 1865 mg kg^− 1 with a median of 751 mg kg^− 1 and standard deviation (SD) of 381 mg kg^− 1. According to cross-validation, predictions for floc TP based on the root mean square prediction error (RMSE) were best in the following order: RK_quadratic (134.9) > RK_multivariate (201.1) > OK (206.1) > CK (212.1) > REG_multivariate (218.3) > REG_quadratic (220.3) > REG_linear (264.4); and based on the mean prediction error (ME) followed the order RK_multivariate (0.9)  RK_quadratic (1.1) > CK (− 6.7) > REG_multivariate (18.2) > REG_linear (25.1) > OK (− 27.3) > REG_quadratic (27.3). The Normalized Difference Vegetation Index (NDVI)-green derived from Landsat ETM+ showed the largest predictive power for floc TP. Measured soil TP ranged from 155 to 1702 mg kg^− 1 with a median of 433 mg kg^− 1 and standard deviation of 316 mg kg^− 1. Predictions for soil TP based on RMSE were best in the following order: RK_ASTER (200.1) > CK_ASTER (238.2)  CK_ETM (239.0) > OK (258.0) > RK_ETM (279.2) > REG_ASTER (281.8) > REG_ETM (356.1); and based on ME followed the order: CK_ASTER (0.1)  CK_ETM (0.2) > RK_ASTER (− 5.2) > RK_ETM (− 31.5) > OK (− 41.8) > REG_ASTER (94.4) > REG_ETM (133.7). The NDVI showed the largest predictive power for soil TP. This comparative study in a subtropical wetland demonstrated the benefits of incorporating remote sensing data into floc and soil TP prediction models. Overall, hybrid geostatistical methods (CK and RK) performed better than regressions and spatial univariate models (OK) in the prediction of floc and soil TP. Depending on the strength of the spatial covariance between primary and secondary variables (CK) and the ability of the regression model in RK to explain the variability of a target variable (e.g., floc or soil TP), either CK or RK performed best. Our findings in this wetland confirmed results from earlier studies on terrestrial soils indicating the superior performance of hybrid geostatistical methods in predicting soil properties.     

Eight years of carbon dioxide exchange above a mixed forest at Borden, Ontario

This paper summarizes results from 8 years (1996–2003) of eddy covariance-based ecosystem CO₂ exchange measurements at the Borden Forest Research Station (44°19′N, 79°56′W). The site represents a mid-latitude, 100-year-old, mixed deciduous and coniferous forest dominated by red maple, aspen and white pine. The years 1996 and 1997 were relatively cold, had a late spring and received below average photosynthetic photon flux density (PPFD). This contrasts with an early spring, warmer soil and air temperatures during 1998–1999, and with distinctly wet year of 2000 and dry years of 2001–2003. The combination of early spring, warmer air and soil temperature and relatively high level of PPFD was associated with higher net ecosystem productivity (NEP) that peaked during 1999. Photosynthetic capacity was reduced and NEP showed a mid-growing season depression during the dry years of 2001–2003. Annual average ecosystem respiration (R) determined from a light response model was 30% less than R derived from a logistic respiration equation, relating night time CO₂ flux and soil temperature. However these independently determined R values were well correlated indicating that the site is unaffected by fetch and spatial heterogeneity problems. Based on the combined 8 years of growing season daytime data, an air temperature of 20–25 °C and a vapor pressure deficit (VPD) of 1.3 kPa were found to be the optimal conditions for CO₂ uptake by the canopy. Over the 1996–2003 period, the forest sequestered carbon at an average rate of 140 ± 111 gC m⁻² y⁻¹. The corresponding gross ecosystem photosynthesis (GEP) and R over this period were 1116 ± 93 gC m⁻² y⁻¹ and 976 ± 68 gC m⁻² y⁻¹, respectively. The annual carbon sequestration ranged from 19 gC m⁻² in 1996 to 281 gC m⁻² in 1999. However, these estimates were sensitive to frictional velocity threshold () used for screening data associated with poor turbulent mixing at night. Increasing from 0.2 m s⁻¹ (based on the inflection point in the nighttime CO₂ flux vs. u_* relationship) to 0.35 m s⁻¹ (determined using a selection algorithm based on change-point detection) modified the 8-year mean NEP estimate from 140 ± 111 gC m⁻² y⁻¹ to 65 ± 120 gC m⁻² y⁻¹. Both approaches show that the Borden forest was a low to moderate sink of carbon over the 8-year period.     

Interpreting the dependence of soil respiration on soil temperature and water content in a boreal aspen stand

Continuous half-hourly measurements of soil CO₂ efflux made between January and December 2001 in a mature trembling aspen stand located at the southern edge of the boreal forest in Canada were used to investigate the seasonal and diurnal dependence of soil respiration (R_s) on soil temperature (T_s) and water content (θ). Daily mean R_s varied from a minimum of 0.1 μmol m⁻² s⁻¹ in February to a maximum of 9.2 μmol m⁻² s⁻¹ in mid-July. Daily mean T_s at the 2-cm depth was the primary variable accounting for the temporal variation of R_s and no differences between Arrhenius and Q₁₀ response functions were found to describe the seasonal relationship. R_s at 10 °C (R_s10) and the temperature sensitivity of R_s (Q_10Rs) calculated at the seasonal time scale were 3.8 μmol m⁻² s⁻¹ and 3.8, respectively. Temperature normalization of daily mean R_s (R_sN) revealed that θ in the 0–15 cm soil layer was the secondary variable accounting for the temporal variation of R_s during the growing season. Daily R_sN showed two distinctive phases with respect to soil water field capacity in the 0–15 cm layer (θ_fc, 0.30 m³ m⁻³): (1) R_sN was strongly reduced when θ decreased below θ_fc, which reflected a reduction in microbial decomposition, and (2) R_sN slightly decreased when θ increased above θ_fc, which reflected a restriction of CO₂ or O₂ transport in the soil profile.Diurnal variations of half-hourly R_s were usually out of phase with T_s at the 2-cm depth, which resulted in strong diurnal hysteresis between the two variables. Daily nighttime R_s10 and Q_10Rs parameters calculated from half-hourly nighttime measurements of R_s and T_s at the 2-cm depth (when there was steady cooling of the soil) varied greatly during the growing season and ranged from 6.8 to 1.6 μmol m⁻² s⁻¹ and 5.5 to 1.3, respectively. On average, daily nighttime R_s10 (4.5 μmol m⁻² s⁻¹) and Q_10Rs (2.8) were higher and lower, respectively, than the values obtained from the seasonal relationship. Seasonal variations of these daily parameters were highly correlated with variations of θ in the 0–15 cm soil layer, with a tendency of low R_s10 and Q_10Rs values at low θ. Overall, the use of seasonal R_s10 and Q_10Rs parameters led to an overestimation of daily ranges of half-hourly R_s (ΔR_s) during drought conditions, which supported findings that the short-term temperature sensitivity of R_s was lower during periods of low θ. The use of daily nighttime R_s10 and Q_10Rs parameters greatly helped at simulating ΔR_s during these periods but did not improve the estimation of half-hourly R_s throughout the year as it could not account for the diurnal hysteresis effect.     

Friction velocity and aerodynamic roughness of conventional and undercutter tillage within the Columbia Plateau, USA

Friction velocity (u_*) and aerodynamic roughness (z_o) at the soil–plant–atmosphere interface affect wind erosion, but no attempts have been made to quantify these parameters as affected by tillage systems within the Columbia Plateau region of the Pacific Northwest United States. Wind velocity profiles above adjacent field plots (>2 ha), with plots subject to conventional or undercutter tillage during the summer fallow phase of a winter wheat–summer fallow rotation, were measured over 50 high wind events (wind velocities in excess of 6.4 m s⁻¹ at a height of 3 m) during 2005 and 2006 near Lind, Washington to determine u_* and z_o of tillage treatments. Wheat stubble plots were subject to either conventional (disks) or undercutter (wide V-shaped blades) tillage in spring and then periodically rodweeded prior to sowing winter wheat in late summer. Prior to sowing, u_* for conventional and undercutter tillage respectively averaged 0.36 and 0.46 m s⁻¹ in 2005 and 0.38 and 0.40 m s⁻¹ in 2006 while z_o for conventional and undercutter tillage respectively averaged 2 and 7 mm in 2005 and 2 and 4 mm in 2006. The aerodynamically rougher surface of undercutter tillage was predicted to suppress vertical dust flux; this was collaborated with observations in the field where undercutter tillage reduced dust flux as compared with conventional disk tillage. Undercutter tillage, therefore, appears to be an effective management practice to roughen the surface and thereby suppress dust emissions from agricultural land subject to summer fallow within the Columbia Plateau.     

N2O and NO fluxes between a Norway spruce forest soil and atmosphere as affected by prolonged summer drought

Global change scenarios predict an increasing frequency and duration of summer drought periods in Central Europe especially for higher elevation areas. Our current knowledge about the effects of soil drought on nitrogen trace gas fluxes from temperate forest soils is scarce. In this study, the effects of experimentally induced drought on soil N₂O and NO emissions were investigated in a mature Norway spruce forest in the Fichtelgebirge (northeastern Bavaria, Germany) in two consecutive years. Drought was induced by roof constructions over a period of 46 days. The experiment was run in three replicates and three non-manipulated plots served as controls. Additionally to the N₂O and NO flux measurements in weekly to monthly intervals, soil gas samples from six different soil depths were analysed in time series for N₂O concentration as well as isotope abundances to investigate N₂O dynamics within the soil. N₂O fluxes from soil to the atmosphere at the experimental plots decreased gradually during the drought period from 0.2 to −0.0 μmol m⁻² h⁻¹, respectively, and mean cumulative N₂O emissions from the manipulated plots were reduced by 43% during experimental drought compared to the controls in 2007. N₂O concentration as well as isotope abundance analysis along the soil profiles revealed that a major part of the soil acted as a net sink for N₂O, even during drought. This N₂O sink, together with diminished N₂O production in the organic layers, resulted in successively decreased N₂O fluxes during drought, and may even turn this forest soil into a net sink of atmospheric N₂O as observed in the first year of the experiment. Enhanced N₂O fluxes observed after rewetting up to 0.1 μmol m⁻² h⁻¹ were not able to compensate for the preceding drought effect. During the experiment in 2006, with soil matric potentials in 20 cm depth down to −630 hPa, cumulative NO emissions from the throughfall exclusion plots were reduced by 69% compared to the controls, whereas cumulative NO emissions from the experimental plots in 2007, with minimum soil matric potentials of −210 hPa, were 180% of those of the controls. Following wetting, the soil of the throughfall exclusion plots showed significantly larger NO fluxes compared to the controls (up to 9 μmol m⁻² h⁻¹ versus 2 μmol m⁻² h⁻¹). These fluxes were responsible for 44% of the total emission of NO throughout the whole course of the experiment. NO emissions from this forest soil usually exceeded N₂O emissions by one order of magnitude or more except during wintertime.     

The role of hydropedologic vegetation zones in greenhouse gas emissions for agricultural wetland landscapes

Net greenhouse gas (GHG) source strength for agricultural wetland ecosystems in the Prairie Pothole Region (PPR) is currently unknown. In particular, information is lacking to constrain spatial variability associated with GHG emissions (CH₄, CO₂, and N₂O). GHG fluxes typically vary with edaphic, hydrologic, biologic, and climatic factors. In the PPR, characteristic wetland plant communities integrate hydropedologic factors and may explain some variability associated with trace gas fluxes at ecosystem and landscape scales. We addressed this question for replicate wetland basins located in central North Dakota stratified by hydropedologic vegetation zone on Jul 12 and Aug 3, 2003. Data were collected at the soil-atmosphere interface for six plant zones: deep marsh, shallow marsh, wet meadow, low prairie, pasture, and cropland. Controlling for soil moisture and temperature, CH₄ fluxes varied significantly with zone (p < 0.05). Highest CH₄ emissions were found near the water in the deep marsh (277,800 μg m^− 2 d^− 1 CH₄), which declined with distance from water to − 730 μg m^− 2 d^− 1 CH₄ in the pasture. Carbon dioxide fluxes also varied significantly with zone. Nitrous oxide variability was greater within zones than between zones, with no significant effects of zone, moisture, or temperature. Data were extrapolated for a 205.6 km² landscape using a previously developed synoptic classification for PPR plant communities. For this landscape, we found croplands contributed the greatest proportion to the net GHG source strength on Jul 12 (45,700 kg d^− 1 GHG-C equivalents) while deep marsh zones contributed the greatest proportion on Aug 3 (26,145 kg d^− 1 GHG-C equivalents). This was driven by a 30-fold reduction in cropland N₂O–N emissions between dates. The overall landscape average for each date, weighted by zone, was 462.4 kg km^− 2 d^− 1 GHG-C equivalents on Jul 12 and 314.3 kg km^− 2 d^− 1 GHG-C equivalents on Aug 3. Results suggest GHG fluxes vary with hydropedologic soil zone, particularly for CH₄, and provide initial estimates of net GHG emissions for heterogeneous agricultural wetland landscapes.     

Magnetic properties of the urban soils in Shanghai and their environmental implications

To study the feasibility of using magnetic techniques for monitoring soil pollution in Shanghai, magnetic properties and heavy metals in the topsoils in an urban site (Songnan Town) and a less-urbanized agricultural site (Luojing Town) in Baoshan District, Shanghai, were studied. Compared with the background, magnetic signals of the urban topsoils are extremely enhanced with magnetic susceptibility (χ_lf) from 127.3–1959 × 10^− 8 m³ kg^− 1; while those of the agricultural topsoils are only slightly increased. However, both the urban and agricultural topsoils contain few pedogenic SP grains, as indicated by their low χ_fd% (< 3.6%). Ratios of χ_arm/SIRM, χ_arm/χ_lf and SIRM/χ_lf indicate that the grain size of magnetic minerals in the urban topsoils is significantly coarser than that in the background and the agricultural topsoils. Furthermore, the urban topsoils show low coercivity and magnetic soft behaviors, as indicated by higher SOFT%, lower HARD%, higher IRM_300 mT/SIRM (close to 1) and lower IRM_− 200 mT/SIRM (close to − 1). It suggests that the urban topsoils have received some coarse ferrimagnetic particles. Heavy metals are highly enriched in the magnetic fractions of the topsoils. Geochemical properties of the magnetic fraction of the urban topsoils are significantly different from those of the agricultural topsoils, further indicating that the extra magnetic minerals accumulated in the urban topsoils are neither inherited from soil parent materials nor from pedogenic processes, but originate from anthropogenic activities. The significant correlations between heavy metals and χ_lf, χ_arm, SIRM, SOFT and HIRM of the topsoils in the district indicate that the magnetic techniques can be used for monitoring soil pollution in Shanghai. The soils with χ_lf from 39–50 × 10^− 8 m³ kg^− 1 in the district are tentatively defined as “slightly polluted soils”; those with χ_lf > 50 × 10^− 8 m³ kg^− 1 are defined as “polluted soils”.     

Prediction of potentially mineralizable N from amidohydrolase activities in a manure-applied, corn residue-amended soil

Prediction of potentially mineralizable N as an important N pool from soil amidohydrolases was investigated. Composite soil samples were collected from plots of a field experiment in which 0, 50 and 100 Mg cow manure ha⁻¹ year⁻¹ had been applied for five consecutive years. The soils were treated with corn shoots or roots or remained untreated in a factorial combination with the manure treatments, with three replications. The mineralized inorganic N was measured periodically in 20-week incubations and potentially mineralizable N (N₀) was calculated based on a first-order kinetic model. Urease, l-glutaminase and l-asparaginase activities were measured before and after incubation. The values of N₀ ranged from 208.6 in the controls to 388.4 in soils that had received 50 Mg ha⁻¹ year⁻¹ of cow manure and were amended with corn shoots. Corn residue amendment in the manure treated soils, increased the values of N₀ or changed the N mineralization kinetic pattern from a first-order to a zero-order model. According to a relative sensitivity index, l-asparaginase was the most sensitive enzyme to the treatments. Multiple regression analysis revealed that 92% of N₀ variations can be described by the activities of urease and l-asparaginase and therefore the soil amidohydrolase activities have the potential to evaluate potentially mineralizable N.     

Effect of di-(2-ethylhexyl) phthalate (DEHP) on microbial biomass C and enzymatic activities in soil

The effects of di-(2-ethylhexyl) phthalate (DEHP) at five different doses from 10 to 1000 mg kg⁻¹ soil on biological properties were investigated over a period of 56 days. Meanwhile, the dissipation of DEHP was also monitored. The results indicated that the microbial biomass C (C_mic) fluctuated at around 70 mg kg⁻¹ soil for the control, whereas the C_mic varied significantly for the soil samples contaminated by DEHP. The catalase activities in all five treatments were stimulated at most time, and the activities of phosphatase in the soils treated by DEHP with 500 mg kg⁻¹ or 1000 mg kg⁻¹ were significantly higher than the other treatments from the 20th day. Urease was more sensitive and inhibited significantly during the initial period of incubation. Additionally, the dose–response relationship of invertase was presented in the later phase of incubation. The activities of urease and invertase might indicate soil perturbations caused by the introduction of DEHP. The dissipation of DEHP was found to follow the pseudo first-order kinetics behavior.     

Temporal and spatial variability of soil respiration in a boreal mixedwood forest

Temporal and spatial variability of soil respiration (R_s) was measured and analyzed in a 74-year-old, mixedwood, boreal forest in Ontario, Canada, over a period of 2 years (August 2003–July 2005). The ranges of R_s measured during the two study years were 0.5–6.9 μmol CO₂ m⁻² s⁻¹ for 2003–2004 (Year 1) and 0.4–6.8 μmol CO₂ m⁻² s⁻¹ for 2004–2005 (Year 2). Mean annual R_s for the stand was the same for both years, 2.7 μmol CO₂ m⁻² s⁻¹. Temporal variability of R_s was controlled mainly by soil temperature (T_s), but soil moisture had a confounding effect on T_s. Annual estimates of total soil CO₂ emissions at the site, calculated using a simple empirical R_s–T_s relationship, showed that R_s can account for about 88 ± 27% of total annual ecosystem respiration at the site. The majority of soil CO₂ emissions came from the upper 12 to 20 cm organic LFH (litter–fibric–humic) soil layer. The degree of spatial variability in R_s, along the measured transect, was seasonal and followed the seasonal trend of mean R_s: increasing through the growing season and converging to a minimum in winter (coefficient of variation (CV) ranged from 4 to 74% in Year 1 and 4 to 62% in Year 2). Spatial variability in R_s was found to be negatively related to spatial variability in the C:N ratio of the LHF layer at the site. Spatial variability in R_s was also found to depend on forest tree species composition within the stand. R_s was about 15% higher in a broadleaf deciduous tree patch compared to evergreen coniferous area. However, the difference was not always significant (at 95% CI). In general, R_s in the mixedwood patch, having both deciduous and coniferous species, was dominated by broadleaf trees, reflecting changing physiological controls on R_s with seasons. Our results highlight the importance of discerning soil CO₂ emissions at a variety of spatial and temporal scales. They also suggest including the LFH soil layer and allowing for seasonal variability in CO₂ production within that layer, when modeling soil respiration in forest ecosystems.     

Influence of temperature and drought on seasonal and interannual variations of soil, bole and ecosystem respiration in a boreal aspen stand

Continuous half-hourly measurements of soil (R_s) and bole respiration (R_b), as well as whole-ecosystem CO₂ exchange, were made with a non steady-state automated chamber system and with the eddy covariance (EC) technique, respectively, in a mature trembling aspen stand between January 2001 and December 2003. Our main objective was to investigate the influence of long-term variations of environmental and biological variables on component-specific and whole-ecosystem respiration (R_e) processes. During the study period, the stand was exposed to severe drought conditions that affected much of the western plains of North America. Over the 3 years, daily mean R_s varied from a minimum of 0.1 μmol m⁻² s⁻¹ during winter to a maximum of 9.2 μmol m⁻² s⁻¹ in mid-summer. Seasonal variations of R_s were highly correlated with variations of soil temperature (T_s) and water content (θ) in the surface soil layers. Both variables explained 96, 95 and 90% of the variance in daily mean R_s from 2001 to 2003. Aspen daily mean R_b varied from negligible during winter to a maximum of 2.5 μmol m⁻² bark s⁻¹ (2.2 μmol m⁻² ground s⁻¹) during the growing season. Maximum R_b occurred at the end of the aspen radial growth increment and leaf emergence period during each year. This was 2 months before the peak in bole temperature (T_b) in 2001 and 2003. Nonetheless, R_b was highly correlated with T_b and this variable explained 77, 87 and 62% of the variance in R_b in the respective years. Partitioning of R_b between its maintenance (R_bm) and growth (R_bg) components using the mature tissue method showed that daily mean R_bg occurred at the same time as aspen radial growth increment during each growing season. This method led, however, to systematic over- and underestimations of R_bm and R_bg, respectively, during each year. Annual totals of R_s, R_b and estimated foliage respiration (R_f) from hazelnut and aspen trees were, on average, 829, 159 and 202 g C m⁻² year⁻¹, respectively, over the 3 years. These totals corresponded to 70, 14 and 16%, respectively, of scaled-up respiration estimates of R_e from chamber measurements. Scaled R_e estimates were 25% higher (1190 g C m⁻² year⁻¹) than the annual totals of R_e obtained from EC (949 g C m⁻² year⁻¹). The independent effects of temperature and drought on annual totals of R_e and its components were difficult to separate because the two variables co-varied during the 3 years. However, recalculation of annual totals of R_s to remove the limitations imposed by low θ, suggests that drought played a more important role than temperature in explaining interannual variations of R_s and R_e.     

Estimated Holocene soil carbon-soil degradation in Nevada and Western Utah, USA

There is much interest in predicting future carbon-soil degradation and that occurring today. We have National Soil Survey Laboratory data to assess some of the soil carbon degradation in the Great Basin and western Utah. For this we included data on 32 Nevada and Utah soils on Pleistocene geomorphic surfaces at elevations of 973 to 3172 m. Their mean annual precipitation (MAP) ranges from 20 to 55 cm and mean annual soil temperatures (MAST) from 5 to 12 °C. The MAP and MAST closely correlate with elevation (E) (r = 0.96 and − 0.97 respectively). Mountain big sagebrush (Artemisia tridentata Nutt. ssp. vaseyana (Rydb.) Beetle) dominates vegetation at the higher, colder elevations. Wyoming big sagebrush (Artemisia tridentata Nutt. ssp. wyomingensis Beetle and Young) and juniper (Juniperus L.) dominate at intermediate elevations. Little sagebrush (Artemisia arbuscula Nutt.) and related desert species dominate at the lower, warmer elevations. We used acid dichromate digestion and FeSO₄ titration to analyze for soil organic carbon (SOC) and bulk density and coarse fragments in the soils to put the data on a volume basis. The soils are well drained and uncultivated. Accumulation of organic carbon in each pedon (OCp) is correlated to MAP and MAST (r = 0.81, and − 0.78 respectively). We predicted OCp from the relationship,

(1)

r² = 0.64, S.E. = 1.30, n = 32. The soil OC degradation that may have occurred through the Holocene ranges from 35% at sites of the present Aridisols and Vertisols to 22% for the sites of the Mollisols and Alfisols. Eq. (1) shows that today, MAST rises of 1 to 3 °C would produce further OCp degradation from today's levels of 1% to 13% in Aridisols and Vertisols and 12% to 25% in the Mollisols and Alfisols respectively. It also shows that if the MAST drop of 6 °C predicted for the Pleistocene occurred, many of the Aridisols and Vertisols likely would have been Mollisols or Alfisols during the Pleistocene. A temperature rise of 1 °C in a century would likely move the Mollisol–Aridisol boundary from its present 2300 m elevation to an elevation of about 2900 m. A temperature rise of 3 °C in a century would likely move the Mollisol–Aridisol boundary from its present 2300 m elevation to elevations of the highest elevations in Nevada and to the middle of Idaho. Increasing the temperature by 3 °C will likely also increase the area affected by severe desertification in the southern Great Basin north by about 20%.     

The utility of process-based models for simulating N2O emissions from soils: A case study based on Costa Rican coffee plantations

Soil moisture and gaseous N-flux (N₂O, N₂) dynamics in Costa Rican coffee plantations were successively simulated using a mechanistic model (PASTIS) and two process-based models (NGAS and NOE). Two fertilized (250 kg N ha⁻¹ y⁻¹) coffee plantations were considered, namely a monoculture and a system shaded by the N₂ fixing legume species Inga densiflora. In situ N₂O fluxes were previously measured in these plantations. NGAS and NOE used specific microbial activities for the soils. To parameterize NGAS, we estimated N mineralization via in situ incubations and the contribution of heterotrophic soil respiration to total soil respiration. Potential denitrification rates and the proportion of denitrified N emitted as N₂O were measured in the laboratory to define the values of NOE parameters, as well as nitrification rates and related N₂O production rates for parameterizing both models. Soil moisture and both NGAS and NOE N₂O fluxes were best modelled on an hourly time step. Soil moisture dynamics were satisfactorily simulated by PASTIS. Simulated N₂O fluxes by both NGAS and NOE (3.2 and 2.1 kg N ha⁻¹ y⁻¹ for NGAS; 7.1 and 3.7 kg N ha⁻¹ y⁻¹ for NOE, for the monoculture and shaded plantations respectively) were within a factor of about 2 of the observed annual fluxes (4.3 and 5.8 kg N ha⁻¹ y⁻¹, for the monoculture and shaded plantations respectively). Statistical indicators of association and coincidence between simulated and measured values were satisfactory for both models. Nevertheless, the two models differed greatly in describing the nitrification and denitrification processes. Some of the algorithms in the model NGAS were apparently not applicable to these tropical acidic Andosols. Therefore, more detailed information about microbial processes in different agroecosystems would be needed, notably if process-oriented models were to be used for testing strategies for mitigating N₂O emissions.     

Alteration of secondary minerals along a time series in young alkaline soils derived from carbonatic wastes of soda production

Industrial activities result in increasing amounts of technical substrates being deposited in landfills. These substrates are subject to weathering and pedogenic processes. We studied the chemical and mineralogical transformations on naturally weathered waste deposits of soda industry. Four sites differing in age (15, 19, 57, and 70 years) and derived from carbonatic slurry (mainly CaO^.H₂O, CaCO₃, NaCl) were selected. The formed soils, calcareous spolic Regosols, are weakly to strongly alkaline with pH values ranging from 8 to 12. Within 15 years, the substrate's initial pH of 12 drops rapidly in the topsoil due to the reaction of dissolved Ca either with CO₂ from the atmosphere or evolved by microbial respiration and finally stabilizes at around 8.1. All soils show high electrical conductivity, up to 12.3 mS cm^− 1 at the youngest site. The electrical conductivity strongly decreases within 70 years of weathering due to leaching processes and the formation of less soluble secondary minerals. The content of organic C in the studied soils ranges from 2.4 to 70.8 g kg^− 1 and stocks increase with site age. Soil structure and soil color change distinctly. The binding of CO₂ results in large amounts of carbonate, increasing with time. Seventy years after deposition, calcite [CaCO₃] dominates the topsoil (0–30 cm depth), comprising about 80% of the soil material. The mineral composition was characterized by X-ray diffraction. Besides calcite, we found different quantities and different distributions of the less common minerals ettringite [Ca₆(Al(OH)₆)₂(SO₄)₃  26H₂O], thaumasite [Ca₆(Si(OH)₆)₂(CO₃)₂(SO₄)₂  24H₂O], hydrocalumite [Ca₂Al(OH)₇  2H₂O] and hydrotalcite [Mg₆Al₂(CO₃)(OH)₁₆  4H₂O]. Formation and alteration of these minerals are basically influenced by changes in the soil pH. With progressing weathering neither thaumasite nor ettringite are stable due to the non-favorable soil reaction (pH  8.1). In contrast, hydrocalumite and hydrotalcite exist in all investigated soils. They are stable also under weakly alkaline conditions and thus may exist in all carbonatic soils. Results indicate a surprisingly rapid soil development driven by the highly dynamic formation and alteration of minerals in carbonatic substrates under alkaline conditions.     

Land-use impacts on surface runoff and soil detachment within agricultural sloping lands in Northern Vietnam

Two consecutive years of investigation on soil surface features, surface runoff and soil detachment within 1-m² microplots on 40% slope highlighted the effects of land-use change, vegetation cover and biological activity on the water pathways in Northern Vietnam. Three replicate plots were set up on each of five land-uses: cassava (CAS), grass fodder of Bracharia ruziziensis (BRA), a 3-year old fallow (FAL), tree stands of Acacia mangium and Venicia montana (FOR), and a fallow with regrowth of Eucalyptus regularly cut (EUC). The second year, two of the microplots under FAL and EUC were treated with herbicide (FALh, EUCh), one of them was burnt (FALh+b, EUCh+b). The highest yearly surface runoff coefficient of 16%, and soil detachment rate of 700 g m^− 2 yr^− 1 in average with a maximum of 1305 g m^− 2 yr^− 1 have been recorded under CAS. On FALh and FALh+b, runoff ratios were 8.7 and 13.5%, respectively and detachment rates were 86 and 389 g m^− 2. On FAL and BRA the yearly runoff ratio varied from 5.9 to 9.8% but the detachment rate was limited at 24 to 35 g m^− 2. FOR and EUC annual runoff was ≤ 3.1% and annual soil detachment ≤ 71 g m^− 2. These values were very low compared to the values reported on steep slopes in Laos within similar climate and vegetation cover.The runoff and detachment rates underlined the importance of rainfall intensities, soil physical properties, soil surface features, soil vegetation cover and biological activity. The annual surface runoff was highly correlated to the soil surface crusting. CAS and BRA plots were prone to crusting especially after weeding at the onset of the rainy season, when the soil surface was still uncovered. Soil bioturbation (earthworm casting activity) was the second factor that explains local variation of surface runoff and soil detachment. The continuous production of earthworms casts on soil surface, especially on FOR and EUC microplots, induced a marked surface roughness and reduced the surface runoff. The production of casts was very limited in FAL and completely absent in CAS microplots. So it is evident that our results confirm the deleterious effects of cassava on soil and water conservation.     

Environmental control of net ecosystem CO2 exchange in a treed, moderately rich fen in northern Alberta

Peatlands cover about 21% of the landscape and contain about 80% of the soil carbon stock in western Canada. However, the current rates of carbon accumulation and the environmental controls on ecosystem photosynthesis and respiration in peatland ecosystems are poorly understood. As part of Fluxnet-Canada, we continuously measured net ecosystem carbon dioxide exchange (NEE) using the eddy covariance technique in a treed fen dominated by stunted Picea mariana and Larix laricina trees during August 2003–December 2004. The total carbon stock in the ecosystem was approximately 51,000 g C m⁻², with only 540 g C m⁻² contributed by live above ground vegetation. The NEE measurements were used to parameterize simple physiological models to assess temporal variation in maximum ecosystem photosynthesis (A_max) and ecosystem respiration rate at 10 °C (R₁₀). During mid-summer the ecosystem had a relatively high A_max (approx. 30 μmol m⁻² s⁻¹) with relatively low R₁₀ (approx. 4 μmol m⁻² s⁻¹). The peak mid-day NEE uptake rate during July and August was 10 μmol m⁻² s⁻¹. The ecosystem showed large seasonal variation in photosynthetic and respiratory activity that was correlated with shifts in temperature, with both spring increases and fall decreases in A_max well predicted by the mean daily air temperature averaged over the preceding 21 days. Leaf-level gas exchange and spectral reflectance measurements also suggested that seasonal changes in photosynthetic activity were primarily controlled by shifts in temperature. Ecosystem respiration was strongly correlated with changes in ecosystem photosynthesis during the growing season, suggesting important links between plant activity and mycorrhizae and microbial activity in the shallow layers of the peat. Only very low rates of respiration were observed during the winter months. During 2004, the peatland recorded a net annual gain of 144 g C m⁻² year⁻¹, the result of a difference between gross photosynthesis of 713 and total ecosystem respiration of 569 g C m⁻² year⁻¹.     

Soil structure and organic carbon relationships following 10 years of wheat straw management in no-till

Crop residue retention is important for sequestering soil organic carbon (SOC), controlling soil erosion, and improving soil quality. Magnitude of residue management impacts on soil structural properties and SOC sequestration is, however, site specific. This study assessed long-term (10 year) impacts of three levels (0, 8, and 16 Mg ha⁻¹ on a dry matter basis) of wheat (Triticum aestivum L.) straw applied annually on SOC concentration and physical properties of the bulk soil and individual 5- to 8-mm aggregates for the 0- to 50-cm soil depth under no-till (NT) on a Crosby silt loam (fine, mixed, active, mesic Aeric Epiaqualfs) in central Ohio. This study also quantified relationships between soil properties and straw-induced changes in SOC concentration. Changes in soil properties due to straw mulching were mostly confined to the upper 5 cm of the soil. Mulching increased SOC concentration, but it did not significantly change cone index (CI) and shear strength (SHEAR). Within the upper 0–5-cm soil depth, mulching decreased bulk density (ρ_b) by 40–50%, aggregate density (ρ_agg) by 30–40%, and particle density (ρ_s) by 10–15%, and increased tensile strength (TS) of aggregates by up to 14 times as compared to unmulched soil. At the same depth, soil with mulch retained >30% more water than soil without mulch from 0 to −1500 kPa potentials. The SOC amount was 16.0 Mg ha⁻¹ under no straw, 25.3 Mg ha⁻¹ under 8 Mg ha⁻¹ straw, and 33.5 Mg ha⁻¹ under 16 Mg ha⁻¹ straw in the 0- to 10-cm depth. Below 10 cm, differences in SOC pool between mulched and unmulched soil were not significant. Overall, SOC from 0- to 50-cm depth was 82.5 Mg ha⁻¹ for unmulched soil, 94.1 Mg ha⁻¹ for 8 Mg ha⁻¹ mulch, and 104.9 Mg ha⁻¹ for 16 Mg ha⁻¹. About 33% of C added with straw over the 10-year period was sequestered in soil. This means that 2/3 of the wheat straw applied was not converted to SOC and most probably was lost as emissions of CO₂ and CH₄. The annual rate of total C accrual was 1.2 Mg ha⁻¹ in soil mulched with 8 Mg ha⁻¹ and 2.2 Mg ha⁻¹ in soil mulched with 16 Mg ha⁻¹ of straw in the 0- to 50-cm depth. The percentage of macroaggregates (>5-mm) was six times higher under 8 Mg ha⁻¹ of straw and 12 times higher under 16 Mg ha⁻¹ compared to unmulched treatments. Macroaggregates contained greater SOC than microaggregates in mulched soil. The SOC concentration explained the variability in aggregate properties by as much as 96%. Overall, long-term straw mulching increased SOC concentration and improved near-surface aggregate properties.     

Assessing spatial variability of soil enzyme activities in pasture topsoils using geostatistics

The purpose of this study was to assess the spatial variability of the activity of three hydrolytic enzymes, i.e. urease activity (UAc), alkaline phosphatase activity (APAc), and arylsulfatase activity (ASAc), in pasture topsoils using geostatistics. Enzyme activities along a transect in a 1.35-ha pasture were determined using 77 soil samples from the upper 20 cm of soil. UAc varied from 101.0 to 182.7 μg N g⁻¹ soil h⁻¹; APAc varied from 1.56 to 3.62 μg p-nitrophenol g⁻¹ soil h⁻¹; and ASAc varied from 1.50 to 3.26 μg p-nitrophenol g⁻¹ soil h⁻¹. The linear models fit the best semivariogram models for UAc, APAc, and ASAc. Semivariograms for enzyme activities exhibited spatial dependence with ranges of influence of approximately 124.7 m.     

Modeling daily soil water dynamics during vertical drainage using the incoming flow concept

Crop management models require simulation of daily soil water dynamics. The objective of this study was to develop a model to simulate the daily soil water dynamics during vertical drainage with reasonable accuracy using the incoming flow concept. The execution of this model, which has been developed based on the conservation of mass law, consists of two steps. First, calculating the potential daily change of soil water content (Δθp) for each soil layer in the profile assuming each one receives no water from the above layer. Then, calculating the actual daily change of soil water (Δθa) for each soil layer in the profile by adjusting Δθp using the incoming water flow, which can be defined as the amount of drainage water that reaches a layer in a soil profile from the above layer. The model was compared with the Suleiman and Ritchie [Suleiman, A.A., Ritchie, J.T., 2004. Modifications to the DSSAT vertical drainage model for more accurate soil water dynamics estimation. Soil Sci. 169 (11), 745–757] vertical drainage model (SRVDM) and HYDRUS-1D for diverse soils and was tested using drainage experimental data of a Eutric Regosol in Bekkevoort, Belgium and a sandy soil in Georgia, U.S. The difference in Δθp between the new model and HYDRUS-1D for diverse soils ranged from − 0.01 to 0.016 m³ m^− 3 for the first day and from − 0.005 to − 0.025 m³ m^− 3 for the second day while the difference in Δθp between the SRVDM and HYDRUS-1D for these soils ranged from 0.014 to 0.062 m³ m^− 3 for the first day and from − 0.01 to 0.026 m³ m^− 3 for the second day. The relative maximum absolute errors in Δθa between the new model and HYDRUS-1D was 10% while the relative maximum absolute errors in Δθa between the SRVDM and HYDRUS-1D was 112%. In the experiments, the root mean square difference of the soil water content for the new model was lower than that for the SRVDM at the different soil depths. These results indicated that the new model outperformed the SRVDM in simulating Δθp and Δθa for diverse soil. It can be concluded that the new model was robust and reasonably accurate for diverse soils at different soil depths. The implementation of such model will improve the accuracy and applicability of regional soil water dynamics simulation and will reduce considerably the computational time and the required inputs.     

Cover crop and tillage effects on soil enzyme activities following tomato

Increasing numbers of vegetable growers are adopting conservation tillage practices and including cover crops into crop rotations. The practice helps to increase or maintain an adequate level of soil organic matter and improves vegetable yields. The effects of the practices, however, on enzyme activities in southeastern soils of the United States have not been well documented. Thus, the objectives of the study were to investigate the effects of cover crops and two tillage systems on soil enzyme activity profiles following tomato and to establish relationships between enzyme activities and soil organic carbon (C) and nitrogen (N). The cover crops planted late in fall 2005 included black oat (Avena strigosa), crimson clover (Trifolium incarnatum L.), or crimson clover–black oat mixed. A weed control (no cover crop) was also included. Early in spring 2006, the plots were disk plowed and incorporated into soil (conventional tillage) or mowed and left on the soil surface (no-till). Broiler litter as source of N fertilizer was applied at a rate of 4.6 Mg ha⁻¹, triple super phosphate at 79.0 kg P ha⁻¹, and potassium chloride at 100 kg K ha⁻¹ were also applied according to soil testing recommendations. Tomato seedlings were transplanted and grown for 60 days on a Marvyn sandy loam soil (fine-loamy, kaolinitic, thermic Typic Kanhapludults). Ninety-six core soil samples were collected at incremental depths (0–5, 5–10, and 10–15 cm) and passed through a 2-mm sieve and kept moist to study arylamidase (EC 3.4.11.2), l-asparaginase (EC 3.5.1.1), l-glutaminase (EC 3.5.1.2), and urease (EC 3.5.1.5) activities. Tillage systems affected only l-glutaminase activity in soil while cover crops affected activities of all the enzymes studied with the exception of urease. The research clearly demonstrated that in till and no-till systems, l-asparaginase activity is greater (P ≤ 0.05) in plots preceded by crimson clover than in those preceded by black oat or their mixture. Activity of the enzyme decreased from 11.7 mg NH₄⁺–N kg⁻¹ 2 h⁻¹ at 0–5 cm depth to 8.73 mg NH₄⁺–N kg⁻¹ 2 h⁻¹ at 5–10 cm and 10–15 cm depths in the no-till crimson clover plots. Arylamidase activity significantly correlated with soil organic C (r = 0.699**) and soil organic N (r = 0.764***). Amidohydrolases activities significantly correlated with soil organic N but only urease significantly correlated with soil organic C (r = 0.481*). These results indicated that incorporation of cover crops into rotations may increase enzyme activities in soils.