USGS Science⌕ Search

SEARCH · USGS Science

Results for “SOIL”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 541 records · Page 30Linked to original sources

Aeolian controls of soil geochemistry and weathering fluxes in high-elevation ecosystems of the Rocky Mountains, Colorado

When dust inputs are large or have persisted for long periods of time, the signature of dust additions are often apparent in soils. The of dust will be greatest where the geochemical composition of dust is distinct from local sources of soil parent material. In this study the influence of dust accretion on soil geochemistry is quantified for two different soils from the San Juan Mountains of southwestern Colorado, USA. At both study sites, dust is enriched in several trace elements relative to local rock, especially Cd, Cu, Pb, and Zn. Mass-balance calculations that do not explicitly account for dust inputs indicate the accumulation of some elements in soil beyond what can be explained by weathering of local rock. Most observed elemental enrichments are explained by accounting for the long-term accretion of dust, based on modern isotopic and geochemical estimates. One notable exception is Pb, which based on mass-balance calculations and isotopic measurements may have an additional source at one of the study sites. These results suggest that dust is a major factor influencing the development of soil in these settings and is also an important control of soil weathering fluxes. After accounting for dust inputs in mass-balance calculations, Si weathering fluxes from San Juan Mountain soils are within the range observed for other temperate systems. Comparing dust inputs with mass-balanced based flux estimates suggests dust could account for as much as 50–80% of total long-term chemical weathering fluxes. These results support the notion that dust inputs may sustain chemical weathering fluxes even in relatively young continental settings. Given the widespread input of far-traveled dust, the weathering of dust is likely and important and underappreciated aspect of the global weathering engine.

Colorado↗

Development of soil radiocarbon profiles in a reactive transport framework

Today, there is a greater appreciation for the importance of the physical protection of carbon (C) through interactions with mineral surfaces, isolation from microbes, and the important role of transport in shaping soil properties and controlling moisture limitations on decomposition. As our paradigm for soil organic carbon (SOC) preservation changes, so too should our representation of the underlying processes in soil models. Reactive transport models (RTMs) provide a framework capable of assessing the interactive influence of soil chemistry and transport processes on the accumulation and turnover of SOC. In this study, we present new developments in the isotopically enabled RTM “CrunchTope,” which is capable of explicitly tracking the three isotopes of carbon ( 12 C, 13 C, and 14 C) and their fractionation between multiple coexisting and interacting solid, liquid and gas phases. This modeling framework opens the door to new applications of depth-resolved RTMs models in application to SOC and deeper subsurface carbon reservoirs. Here, we demonstrate SOC accumulation and radiocarbon aging for long-timescale models of soil development in CrunchTope. Our goal is to assess advantages and limitations of such an approach and to identify the type and complexity of reaction networks that are required to adequately apply this model to SOC dynamics. We assess the behavior of this model relative to a high-resolution dataset of SOC content, stable isotope composition, and radiocarbon ages as well as physical and hydrologic data measured from a chronosequence of soils located near Santa Cruz, California. Starting from a previously published model using a simplified reaction network with a single class of carbon, we sequentially incorporate multiple C reservoirs subject to both reactivity and transport pathways. Our results indicate that multiple SOC pools with different mean ages of C do not inherently emerge as a result of including reactions which are conventionally expected to provide a diversity of transit times, i.e., sorption and complexation of SOC on mineral surfaces. Instead, transit times emerge as a result of the timescales of the reactions represented in the reaction network. For mineral associated C, the RTM framework imposes dynamic equilibrium with the fluid phase dissolved organic C, such that no distinction in radiocarbon ages is achieved between these pools. Aged C can be produced by including a solid-phase C reservoir, with a rate-limited solubilization coefficient. Aging of SOC in this way is more akin to selective preservation than to mineral protection and, while such a mechanism may be at play in many soils, mineral protection is thought to be at least as important. As such, our results indicate that additional parameterization is required to reproduce the heterogeneity of carbon transit times that result from organo-mineral interactions. These efforts show the promise of a modeling approach where the varied transit time of soil C emerges from the dynamic physical and hydrologic properties of the model rather than from the a priori assignment of operationally defined pools.

Geochimica et Cosmochimica Acta↗

Changes in soil erosion caused by wildfire: A conceptual biogeographic model

Soil erosion rates after wildfire are strongly controlled by intrinsic properties such as topography, weather, climate, soil, and vegetation. These landscape and hydroclimatic properties are important in determining post-fire erosion rates; however, their influence on post-fire erosion and their interaction with the intensity of a wildfire remains uncertain. A key limitation in resolving this uncertainty is the lack of conceptual models and frameworks for organising data related to the geomorphic sensitivity of landscapes to wildfire. Our aim is to develop a framework for consolidating understanding of post-fire erosion in the context of hydroclimatic conditions which contribute to system states, for example soil and vegetation properties, and wildfire regime. The framework is developed around a simple conceptual model where the change in erosion due to wildfire is a product of change in runoff generation and sediment supply, which is strongly related to landscape net primary productivity (NPP). We hypothesised that geomorphic sensitivity to wildfire should vary as a unimodal humped relationship across a gradient of NPP, peaking at an intermediate level. To develop this framework and to test the hypothesis, we first review intrinsic soil and vegetation properties related to the supply and transport of sediment from burned and unburned hillslopes. Net primary productivity is systematically related to these intrinsic properties because it integrates many processes involved in soil and vegetation development. Empirical data indicate a trend in the change in surface runoff generation with NPP after wildfire, peaking at an NPP of approximately 15 Mg C ha −1 y −1 . A simple model of fuel availability and soil heating are correlated with a similar “humped” trend in sediment supply. These results are consistent with our conceptual model, which indicates that sediment supply and runoff contribute towards a distinct peak in wildfire effects on erosion at an intermediate level of NPP. We propose that landscapes of intermediate NPP typically have the highest quantity of fuel available to burn, which cause large changes to the soil surface properties. Landscapes at intermediate NPP also tend to produce intrinsic soil and vegetation properties that promote erosion after wildfire. The interplay between these short and long-term landscape characteristics is strongest at intermediate levels of NPP. Our proposed biogeographic model of geomorphic sensitivity to wildfire was supported by erosion data from burned hillslope and zero-order catchments studies from a range fire-prone landscapes in Australia and North America. Our proposed conceptual model will help identify areas most vulnerable to post-fire erosion changes.

Geomorphology↗

Abundance and distribution of selected elements in soils, stream sediments, and selected forage plants from desert tortoise habitats in the Mojave and Colorado deserts, USA

A baseline and background chemical survey was conducted in southeastern California, USA, to identify potential sources of toxicants in natural and anthropogenically-altered habitats of the threatened desert tortoise (Gopherus agassizii). Soil, stream sediment, and plant samples were collected from six tortoise habitat study areas in the Mojave and Colorado deserts and analysed for up to 66 different elements. The chemical analyses provided new information on the abundances and distributions of selected elements in this region. Soil, stream-sediment, and plant analyses showed distinct variations in bulk chemistries from locality to locality. Variations were, in general, consistent with the many types of exposed rock units in the region, their highly variable bulk mineralogies, and chemical contents. Of elements in soils that might have been toxic to tortoises, only As seemed to be anomalous region-wide. Some soil and plant anomalies were clearly anthropogenic. In the Rand and Atolia mining districts, soil anomalies for As, Au, Cd, Hg, Sb, and(or) W and plant anomalies for As, Sb, and(or) W extend as far as ???15 km outward from the present area of mining; soils containing anomalous Hg were found at least 6 km away from old piles of tailings. The anomalous concentrations of As and Hg may have been the source of elevated levels of these elements found in ill tortoises from the region. In the Goldstone mining district, soil anomalies extended several km from the mining area. These areas probably represented anthropogenic surface contamination of dust redistributed by wind, vehicles, and rainfall. One of two study areas transected by a paved road (Chemehuevi Valley) showed weakly elevated levels of Pb, which extended as far as ???22 m from the pavement edge and were probably related to vehicle exhaust. No soil or plant samples from historically used military areas (Goldstone, Goffs, Chemehuevi Valley, Chuckwalla Bench) contained anomalous concentrations of the elements As, Cu, Cr, Fe, Pb, or Zn that could be ascribed to military maneuvers, vehicles, or ordnance. For future studies, the distribution and abundance of elements in the tortoise forage plants need to be evaluated for the respective roles of dust and systemic uptake. Additional chemical data from tortoise necropsies and nutritional studies are needed to determine the effects of potentially toxic elements in tortoise habitats on their health. ?? 2006 Elsevier Ltd. All rights reserved.

Journal of Arid Environments↗

Disruption rates for one vulnerable soil in Organ Pipe Cactus National Monument, Arizona, USA

Rates of soil disruption from hikers and vehicle traffic are poorly known, particularly for arid landscapes. We conducted an experiment in Organ Pipe Cactus National Monument (ORPI) in western Arizona, USA, on an air-dry very fine sandy loam that is considered to be vulnerable to disruption. We created variable-pass tracks using hikers, an all-terrain vehicle (ATV), and a four-wheel drive vehicle (4WD) and measured changes in cross-track topography, penetration depth, and bulk density. Hikers (one pass = 5 hikers) increased bulk density and altered penetration depth but caused minimal surface disruption up to 100 passes; a minimum of 10 passes were required to overcome surface strength of this dry soil. Both ATV and 4WD traffic significantly disrupted the soil with one pass, creating deep ruts with increasing passes that rendered the 4WD trail impassable after 20 passes. Despite considerable soil loosening (dilation), bulk density increased in the vehicle trails, and lateral displacement created berms of loosened soil. This soil type, when dry, can sustain up to 10 passes of hikers but only one vehicle pass before significant soil disruption occurs; greater disruption is expected when soils are wet. Bulk density increased logarithmically with applied pressure from hikers, ATV, and 4WD.

Arizona↗

The influence of use-related, environmental, and managerial factors on soil loss from recreational trails

Recreational uses of unsurfaced trails inevitably result in their degradation, with the type and extent of resource impact influenced by factors such as soil texture, topography, climate, trail design and maintenance, and type and amount of use. Of particular concern, the loss of soil through erosion is generally considered a significant and irreversible form of trail impact. This research investigated the influence of several use-related, environmental, and managerial factors on soil loss on recreational trails and roads at Big South Fork National River and Recreation Area, a unit of the U.S. National Park Service. Regression modeling revealed that trail position, trail slope alignment angle, grade, water drainage, and type of use are significant determinants of soil loss. The introduction of individual and groups of variables into a series of regression models provides improved understanding and insights regarding the relative influence of these variables, informing the selection of more effective trail management actions. Study results suggest that trail erosion can be minimized by avoiding “fall-line” alignments, steep grades, and valley-bottom alignments near streams, installing and maintaining adequate densities of tread drainage features, applying gravel to harden treads, and reducing horse and all-terrain vehicle use or restricting them to more resistant routes. This research also sought to develop a more efficient Variable Cross-Sectional Area method for assessing soil loss on trails. This method permitted incorporation of CSA measures in a representative sampling scheme applied to a large (24%) sample of the park's 526 km trail system. The variety of soil loss measures derived from the Variable CSA method, including extrapolated trail-wide soil loss estimates, permit an objective quantification of soil erosion on recreational trails and roads. Such data support relational analyses to increase understanding of trail degradation, and long-term monitoring of the natural and recreational integrity of the trail system infrastructure.

Journal of Environmental Management↗

NEHRP soil classifications for estimating site-dependent seismic coefficients in the Upper Mississippi Embayment

Local soil conditions have a profound influence on the characteristics of ground shaking during an earthquake. Exceptionally deep soil deposits, on the order of 100-1000 m deep, are found in the Upper Mississippi Embayment of the central United States. Shear waves (SH) from earthquakes in the New Madrid seismic zone are expected to be strongly affected by the sharp impedance contrasts at the bedrock/sediment interface, attenuation of seismic waves in the soil column, and the SH-wave velocities of the more poorly consolidated near-surface (???50 m) soils. SH-wave velocities of the near-surface soils at nearly 400 sites in the Upper Mississippi Embayment were determined using conventional seismic SH-wave refraction and reflection techniques. Based on the average SH-wave velocities of the upper 30 m of the soils, sites in the Mississippi River floodplain portion of the study area are predominantly classified as Site Class D (180-360 m/s) in accordance with the 1997 NEHRP provisions. Sites away from the active floodplains in western Kentucky and western Tennessee, the SH-wave velocities of the upper 30 m of soils typically ranged from mid-200 to mid-300 m/s. Several sites in western Kentucky had averaged SH-wave velocities greater than 360 m/s, thereby qualifying them as Site Class C (360-760 m/s) in accordance with the 1997 NEHRP provisions. One dimensional site effects, including amplification and dynamic site period, were calculated for a representative suite of sites across the Upper Mississippi Embayment at latitude ?? 38.5??. Although seismic attenuation is greater in the Mississippi River floodplain (i.e. thicker, lower velocity material), the site effects tend to be greater than in the uplands of western Tennessee because of larger impedance contrasts within the near-surface soils. ?? 2001 Elsevier Science B.V. All rights reserved.

Engineering Geology↗

Soil N and 15N variation with time in a California annual grassland ecosystem

The %N and ??15N values of soils and plants were measured along a chronosequence spanning 3 to 3000 Ky in a California annual grassland. Total soil N decreased with increasing soil age (1.1 to 0.4 kg N m-2) while the mean ?? 15N values of the soil N increased by several ??? from the youngest to oldest sites (+3.5 to +6.2 ???). The ?? 15N values of plants varied along the gradient, reflecting changing soil N pools and differences in the form of N uptake. The decline in total N storage with time is hypothesized to be due to a shift from N to P limitation with increasing soil age. The general increase in ?? 15N values with time is interpreted using a N mass balance model, and appears to reflect a shift toward an increasing proportional losses of inorganic mineral forms of N (vs. organic forms) with increasing soil age. We develop a quantitative index of this trend (mineral vs. organic forms of N loss) using mass balance considerations and parameters. The %N and ?? 15N values along the California age gradient were compared to the published data for a comparably aged chronosequence in Hawaii. Most striking in this comparison is the observation that the California soil and plant ?? 15N values are several ??? greater than those on comparably aged Hawaiian sites. Multiple explanations are plausible, but assuming the sites have a similar range in ?? 15N values of atmospheric inputs, the isotopic differences suggest that N may be, at least seasonally, in greater excess in the strongly seasonal, semi-arid, California grassland. Copyright ?? 2001 Elsevier Science Ltd.

Geochimica et Cosmochimica Acta↗

Chemical weathering of a soil chronosequence on granitoid alluvium: II. Mineralogic and isotopic constraints on the behavior of strontium

The use of strontium isotopes to evaluate mineral weathering and identify sources of base cations in catchment waters requires an understanding of the behavior of Sr in the soil environment as a function of time. Our approach is to model the temporal evolution of 87 Sr/ 86 Sr of the cation exchange pool in a soil chronosequence developed on alluvium derived from central Sierra Nevada granitoids during the past 3 Ma. With increasing soil age, 87 Sr/ 86 Sr of ammonium-acetate extractable Sr initially decreases from values typical of K-feldspar to those of plagioclase and hornblende and then remains constant, even though plagioclase and hornblende are absent from the soils after approximately 1 Ma of weathering. The temporal variation of 87 Sr/ 86 Sr of exchangeable Sr is modeled by progressively equilibrating Sr derived from mineral weathering and atmospheric deposition with Sr on exchange sites as waters infiltrate a soil column. Observed decreases in quartz-normalized modal abundances of plagioclase, hornblende, and K-feldspar with time, and the distinct 87 Sr/ 86 Sr values of these minerals can be used to calculate Sr flux from weathering reactions. Hydrobiotites in the soils have nearly constant modal abundances, chemistry, and 87 Sr/ 86 Sr over the chronosequence and provide negligible Sr input to weathering solutions. The model requires time and soil horizon-dependent changes in the amount of exchangeable Sr and the efficiency of Sr exchange, as well as a biologic cycling term. The model predicts that exchangeable Sr initially has 87 Sr/ 86 Sr identical to that of K-feldspar, and thus could be dominated by Sr leached from K-feldspar following deposition of the alluvium. The maximum value of 87 Sr/ 86 Sr observed in dilute stream waters associated with granitoids of the Yosemite region is likewise similar to that of the K-feldspars, suggesting that K-feldspar and not biotite may be the dominant source of radiogenic Sr in the streams. This study reveals that, when attempting to use Strontium isotopes to identify sources of base cations in catchment waters and biomass, both preferential leaching of Sr from minerals during incipient soil development and changing Sr exchange efficiency must be considered along with chemical contributions due to mineral dissolution.

Geochimica et Cosmochimica Acta↗

The surface area of soil organic matter

The previously reported surface area for soil organic matter (SOM) of 560-800 m2/g as determined by the ethylene glycol (EG) retention method was reexamined by the standard BET method based on nitrogen adsorption at liquid nitrogen temperature. Test samples consisted of two high organic content soils, a freeze-dried soil humic acid, and an oven-dried soil humic acid. The measured BET areas for these samples were less than 1 m2/g, except for the freeze-dried humic acid. The results suggest that surface adsorption of nonionic organic compounds by SOM is practically insignificant in comparison to uptake by partition. The discrepancy between the surface areas of SOM obtained by BET and EG methods was explained in terms of the 'free surface area' and the 'apparent surface area' associated with these measurements.The previously reported surface area for soil organic matter (SOM) of 560-800 m2/g as determined by the ethylene glycol (EG) retention method was reexamined by the standard BET method based on nitrogen adsorption at liquid nitrogen temperature. Test samples consisted of two high organic content soils, a freeze-dried soil humic acid, and an oven-dried soil humic acid. The measured BET areas for these samples were less than 1 m2/g, except for the freeze-dried humic acid. The results suggest that surface adsorption of nonionic organic compounds by SOM is practically insignificant in comparison to uptake by partition. The discrepancy between the surface areas of SOM obtained by BET and EG methods was explained in terms of the 'free surface area' and the 'apparent surface area' associated with these measurements.

Environmental Science & Technology↗

Partition characteristics of polycyclic aromatic hydrocarbons on soils and sediments

The partition behavior was determined for three polycyclic aromatic hydrocarbons (PAHs), i.e., naphthalene, phenanthrene, and pyrene, from water to a range of soil and sediment samples. The measured partition coefficients of the individual PAHs between soil/sediment organic matter (SOM) and water (i.e., K(oc) values) are relatively invariant either for the 'clean' (uncontaminated) soils or for the clean sediments; however, the mean K(oc) values on the sediments are about twice the values on the soils. This disparity is similar to the earlier observation for other nonpolar solutes and reflects the compositional differences between soil and sediment organic matters. No significant differences in K(oc) are observed between a clean coastal marine sediment and freshwater sediments. The coastal sediments that are significantly impacted by organic contaminants exhibit higher K(oc) values. At given K(ow) values (octanol-water), the PAHs exhibit much higher K(oc) values than other relatively nonpolar solutes (e.g., chlorinated hydrocarbons). This effect is shown to result from the enhanced partition of PAHs to SOM rather than from lower K(ow) values of PAHs at given supercooled liquid solute solubilities in water. The enhanced partition of PAHs over other nonpolar solutes in SOM provides an account of the markedly different correlations between log K(oc) and log K(ow) for PAHs and for other nonpolar solutes. The improved partition of PAHs in SOM stems apparently from the enhanced compatibility of their cohesive energy densities with those of the aromatic components in SOM. The approximate aromatic fraction in soil/sediment organic matter has been assessed by solid-state 13C-NMR spectroscopy.The partition behavior was determined for three polycyclic aromatic hydrocarbons (PAHs), i.e., naphthalene, phenanthrene, and pyrene, from water to a range of soil and sediment samples. The measured partition coefficients of the individual PAHs between soil/sediment organic matter (SOM) and water (i.e., Koc values) are relatively invariant either for the `clean' (uncontaminated) soils or for the clean sediments; however, the mean Koc values on the sediments are about twice the values on the soils. This disparity is similar to the earlier observation for other nonpolar solutes and reflects the compositional differences between soil and sediment organic matters. No significant differences in Koc are observed between a clean coastal marine sediment and freshwater sediments. The coastal sediments that are significantly impacted by organic contaminants exhibit higher Koc values. At given Kow values (octanol-water), the PAHs exhibit much higher Koc values than other relatively nonpolar solutes (e.g., chlorinated hydrocarbons). This effect is shown to result from the enhanced partition of PAHs to SOM rather than from lower Kow values of PAHs at given supercooled liquid solute solubilities in water. The enhanced partition of PAHs over other nonpolar solutes in SOM provides an account of the markedly different correlations between log Koc and log Kow for PAHs and for other nonpolar solutes. The improved partition of PAHs in SOM stems apparently from the enhanced compatibility of their cohesive energy densities with those of the aromatic components in SOM. The approximate aromatic fraction in soil/sediment organic matter has been assessed by solid-state 13C-NMR spectroscopy.

Environmental Science & Technology↗

Modeling soil thermal and carbon dynamics of a fire chronosequence in interior Alaska

In this study, the dynamics of soil thermal, hydrologic, and ecosystem processes were coupled to project how the carbon budgets of boreal forests will respond to changes in atmospheric CO 2 , climate, and fire disturbance. The ability of the model to simulate gross primary production and ecosystem respiration was verified for a mature black spruce ecosystem in Canada, the age-dependent pattern of the simulated vegetation carbon was verified with inventory data on aboveground growth of Alaskan black spruce forests, and the model was applied to a postfire chronosequence in interior Alaska. The comparison between the simulated soil temperature and field-based estimates during the growing season (May to September) of 1997 revealed that the model was able to accurately simulate monthly temperatures at 10 cm ( R > 0.93) for control and burned stands of the fire chronosequence. Similarly, the simulated and field-based estimates of soil respiration for control and burned stands were correlated ( R = 0.84 and 0.74 for control and burned stands, respectively). The simulated and observed decadal to century-scale dynamics of soil temperature and carbon dynamics, which are represented by mean monthly values of these variables during the growing season, were correlated among stands ( R = 0.93 and 0.71 for soil temperature at 20- and 10-cm depths, R = 0.95 and 0.91 for soil respiration and soil carbon, respectively). Sensitivity analyses indicate that along with differences in fire and climate history a number of other factors influence the response of carbon dynamics to fire disturbance. These factors include nitrogen fixation, the growth of moss, changes in the depth of the organic layer, soil drainage, and fire severity.

Alaska↗

Short-term response of methane fluxes and methanogen activity to water table and soil warming manipulations in an Alaskan peatland

Growing season CH 4 fluxes were monitored over a two year period following the start of ecosystem-scale manipulations of water table position and surface soil temperatures in a moderate rich fen in interior Alaska. The largest CH 4 fluxes occurred in plots that received both flooding (raised water table position) and soil warming, while the lowest fluxes occurred in unwarmed plots in the lowered water table treatment. A combination of treatment and soil hydroclimate variables explained more than 70% of the variation in ln-transformed CH 4 fluxes, with mean daily water table position representing the strongest predictor. We used quantitative PCR of the α -subunit of mcr operon to explore the influence of soil climate manipulations on methanogen abundances. Methanogen abundances were greatest in warmed plots, and showed a positive relationship with mean daily CH 4 fluxes. Our results show that water table manipulations that led to soil inundation (flooding) had a stronger effect on CH 4 fluxes than water table drawdown. Seasonal CH 4 fluxes increased by 80–300% under the combined wetter and warmer soil climate treatments. Thus, while warming is expected to increase CH 4 emissions from Alaskan wetlands, higher water table positions caused by increases in precipitation or disturbances such as permafrost thaw that lead to thermokarst and flooding in wetlands will stimulate CH 4 emissions beyond the effects of soil warming alone. Consequently, we argue that modeling the effects of climate change on Alaskan wetland CH 4 emissions needs to consider the interactive effects of soil warming and water table position on CH 4 production and transport.

Alaska↗

Evaluation of MODIS NDVI and NDWI for vegetation drought monitoring using Oklahoma Mesonet soil moisture data

The evaluation of the relationship between satellite-derived vegetation indices (normalized difference vegetation index and normalized difference water index) and soil moisture improves our understanding of how these indices respond to soil moisture fluctuations. Soil moisture deficits are ultimately tied to drought stress on plants. The diverse terrain and climate of Oklahoma, the extensive soil moisture network of the Oklahoma Mesonet, and satellite-derived indices from the Moderate Resolution Imaging Spectroradiometer (MODIS) provided an opportunity to study correlations between soil moisture and vegetation indices over the 2002-2006 growing seasons. Results showed that the correlation between both indices and the fractional water index (FWI) was highly dependent on land cover heterogeneity and soil type. Sites surrounded by relatively homogeneous vegetation cover with silt loam soils had the highest correlation between the FWI and both vegetation-related indices (r???0.73), while sites with heterogeneous vegetation cover and loam soils had the lowest correlation (r???0.22). Copyright 2008 by the American Geophysical Union.

Geophysical Research Letters↗

Estimating soil respiration in a subalpine landscape using point, terrain, climate and greenness data

Landscape carbon (C) flux estimates are necessary for assessing the ability of terrestrial ecosystems to buffer further increases in anthropogenic carbon dioxide (CO2) emissions. Advances in remote sensing have allowed for coarse-scale estimates of gross primary productivity (GPP) (e.g., MODIS 17), yet efforts to assess spatial patterns in respiration lag behind those of GPP. Here, we demonstrate a method to predict growing season soil respiration at a regional scale in a forested ecosystem. We related field measurements (n=144) of growing season soil respiration across subalpine forests in the Southern Rocky Mountains ecoregion to a suite of biophysical predictors with a Random Forest model (30 m pixel size). We found that Landsat Enhanced Vegetation Index (EVI), growing season AI, temperature, precipitation, elevation, and slope aspect explained spatiotemporal variability in soil respiration. Our model had a psuedo-r2 of 0.45 and root mean squared error (RMSE) of roughly one-quarter of the mean value of respiration. Predicted growing season soil respiration across the region was remarkably consistent across 2004, 2005 and 2006 (150-d averages of 542.8, 544.3, and 536.5 g C m-2, respectively). Yet, we observed substantial variability in spatial patterns of soil respiration predictions that varied between years, suggesting that our method is sensitive to changes in respiration drivers. We compared our estimates to MODIS GPP and nocturnal net ecosystem exchange (NEE) derived from eddy covariance towers as a proxy for ecosystem respiration. Averaged across the predictive region, mean predicted growing season soil respiration was 73% of MODIS GPP, while predicted soil respiration was generally within 20% of nocturnal NEE from eddy covariance towers. This study demonstrated that geospatial and remotely-sensed datasets can be used in a statistical modeling framework to estimate soil respiration at landscape scales.

Colorado, Wyoming↗

Multi-stage soil-hydraulic recovery and limited ravel accumulations following the 2017 Nuns and Tubbs wildfires in Northern California

Wildfire can impact soil-hydraulic properties by reducing saturated hydraulic conductivity and sorptivity, making recently burned landscapes prone to debris flows and flash floods. The post-fire hazard window can range from years to decades. In Northern California, where wildfire frequency is steadily increasing, the impact and soil-hydraulic recovery from wildfires is unknown. Following the October 2017 Nuns and Tubbs fires in the Northern Bay Area of California, we established 41 monitoring sites for repeat tension-disc infiltrometer measurements of field-saturated hydraulic conductivity ( K fs ) over 3.5 years. Our site arrays, which encompass grasslands, chaparral, and oak and conifer forests across a range in lithology, show a marked decrease in K fs following the wildfires and a swift partial recovery following the initial post-fire rainy season. Our time series reveals a complex path to soil-hydraulic recovery marked by distinct seasonal stages. Analysis of changing K fs , sorptivity, and infiltration model residuals collectively suggests that these stages are related to transitions between soil-hydraulic processes like structural soil sealing from rainsplash, thermal cracking of bare soil, and vegetation regrowth. While soil infiltration rates were strongly impacted by the 2017 fires, dry ravel estimates are an order of magnitude less for similar slopes than the 2009 Station fire in the San Gabriel mountains of Southern California, suggesting that limited ravel flux may insufficiently load channels for debris flows that initiate from within-channel failure. Our analysis suggests that burned landscapes in the Northern Bay Area of California may experience rapid soil-hydraulic recovery and limited pathways toward post-fire debris flow initiation.

California↗

A model of the spatiotemporal dynamics of soil carbon following coastal wetland loss applied to a Louisiana salt marsh in the Mississippi River Deltaic Plain

The potential for carbon sequestration in coastal wetlands is high due to protection of carbon (C) in flooded soils. However, excessive flooding can result in the conversion of the vegetated wetland to open water. This transition results in the loss of wetland habitat in addition to the potential loss of soil carbon. Thus, in areas experiencing rapid wetland submergence, such as the Mississippi River Delta, coastal wetlands could become a significant source of carbon emissions if land loss is not mitigated. To accurately assess the capacity of wetlands to store (or emit) carbon in dynamic environments, it is critical to understand the fate of soil carbon following the transition from vegetated wetland to open water. We developed a simple soil carbon model representing soil depths to 1 m using the data collected from a Louisiana coastal salt marsh in the Mississippi River Deltaic Plain to predict soil carbon density and stock following the transition from a vegetated salt marsh to an open water pond. While immediate effects of ponding on the distribution of carbon within the 1-m soil profile were apparent, there were no effects of ponding on the overall, integrated, carbon stocks 14 years, following wetland submergence. Rather, the model predicts that soil carbon losses in the first meter will be realized over long periods of time (∼200 years) due to changes in the source of carbon (biomass vs. mineral sediment) with minimal losses through mineralization.

Louisiana↗

The mobilization of aluminum in a natural soil system: Effects of hydrologic pathways

A two-component soil water flow model was used in conjunction with an equilibrium speciation model WATEQF to study aluminum mobility in soils of a forested watershed, White Oak Run, in the Shenandoah National Park, Virginia. Soil solution samples, taken from the O, E, B, C 1 , and C 2 horizons, were collected from zero-tension lysimeters designed to collect faster gravitational macropore flow and tension lysimeters designed to collect slower capillary micropore flow. Dissolved aluminum was fractionated into acid-soluble, inorganic monomeric, and organic monomeric aluminum. Soil water aluminum concentrations decreased with depth indicating that the deep soil is a sink for aluminum. All waters contained significant concentrations of acid-soluble aluminum and exhibited a negative correlation between p H and the inorganic monomeric aluminum concentrations. Water in the shallow soil showed distinctly different chemical compositions for the two flow types, while C horizon micropore and macropore waters were more similar. Because of its shorter residence time, water flowing in deep soil macropores underwent less extensive neutralization and immobilization of aqueous aluminum than micropore water. The O horizon macropore waters were undersaturated for all hydroxide, silicate, and sulfate mineral phases considered. The C horizon samples from both flow types were near equilibrium with respect to kaolinite and synthetic gibbsite, indicating that mineral solubility controls water chemistry in the deep soil, while organic substances are the key control in the shallow macropore waters.

Water Resources Research↗