Search USGSSearch

SEARCH · Search USGS

Results for “Soil Systems”

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 145 records · Page 8Linked to original sources

Relative prediction intervals reveal larger uncertainty in 3D approaches to predictive digital soil mapping of soil properties with legacy data

Fine scale maps of soil properties enable efficient land management and inform earth system models. Recent efforts to create soil property maps from field observations tend to use similar tree-based machine learning interpolation approaches, but often deal with depth of predictions, validation, and uncertainty differently. One of the main differences in approaches is whether to model individual depths of interest separately as ‘2D’ models, or to create models that incorporate depth as a predictor variable creating a ‘3D’ model that can make pre-dictions for all depths. It is unclear how choice of 2D or 3D approach influences model accuracy and uncertainty due to lack of direct comparison and inconsistent presentation of results in past studies. This study compares 2D and 3D methods for mapping soil electrical conductivity (salinity), pH, sum of fine and very fine sands, and organic carbon at 30 m resolution for the upper 432,000 km 2 of the Colorado River Watershed of the United States of America. A new, simple, model-agnostic relative prediction interval (RPI) approach to report un-certainty is presented that scales prediction interval width to the 95% interquantile width of the original training sample distribution. The RPI approach enables direct comparison of uncertainty between properties and depths and is easily interpretable by end users. Results indicate that 3D mapping of soil properties with strong variation with depth can result in substantial areas with much higher uncertainty that coincide with unrealistic predictions relative to 2D models, even though 3D models had slightly better global cross-validation scores. Maps and global model summaries of RPI proved helpful in identifying these issues with 3D models. These results suggest that the use of RPI or similar approaches to evaluate models can identify accuracy problems not evident in global va-lidation diagnostics.

Geoderma

In situ soil moisture sensors in undisturbed soils

Soil moisture directly affects operational hydrology, food security, ecosystem services, and the climate system. However, the adoption of soil moisture data has been slow due to inconsistent data collection, poor standardization, and typically short record duration. Soil moisture, or quantitatively volumetric soil water content (SWC), is measured using buried, in situ sensors that infer SWC from an electromagnetic response. This signal can vary considerably with local site conditions such as clay content and mineralogy, soil salinity or bulk electrical conductivity, and soil temperature; each of these can have varying impacts depending on the sensor technology., Furthermore, poor soil contact and sensor degradation can affect the quality of these readings over time. Unlike more traditional environmental sensors, there are no accepted standards, maintenance practices, or quality controls for SWC data. As such, SWC is a challenging measurement for many environmental monitoring networks to implement. Here, we attempt to establish a community-based standard of practice for in situ SWC sensors so that future research and applications have consistent guidance on site selection, sensor installation, data interpretation, and long-term maintenance of monitoring stations., The videography focuses on a multi-agency consensus of best-practices and recommendations for the installation of in situ SWC sensors. This paper presents an overview of this protocol along with the various steps essential for high-quality and long-term SWC data collection. This protocol will be of use to scientists and engineers hoping to deploy a single station or an entire network.

Journal of Visualized Experiments

Karstification without carbonic acid: Bedrock dissolution by gypsum- driven dedolomitization

Aggressive karstification can take place where dolomite and gypsum are in contact with the same aquifer. Gypsum dissolution drives the precipitation of calcite, thus consuming carbonate ions released by dolomite. Lake Banyoles, in northeastern Spain, is a karst lake supplied by sublacustrine springs, and karstic collapse is occurring in the immediate vicinity of the lake. Lake water is dominated by Mg-Ca and SO 4 -HCO 3 , and is supersaturated with calcite that is actively accumulating in lake sediments. Water chemistry, sulfur isotope composition, local stratigraphy, and mass-balance modelling suggest that the primary karst-forming process at at Lake Banyoles is dedolomitization of basement rocks driven by gypsum dissolution. Karstification takes place along the subsurface contact between the gypsiferous Beuda Formation and the dolomitic Perafita Formation. This process is here recognized for the first time to cause karstification on a large scale; this is significant because it proceeds without the addition of soil-generated carbonic acid. Gypsum-driven dedolomitization may be responsible for other karstic systems heretofore attributed to soil-generated carbonic acid.

Geology

Understanding the role of initial soil moisture and precipitation magnitude in flood forecast using a hydrometeorological modelling system

We adapted the WRF-Hydro modelling system to Hurricane Florence (2018) and performed a series of diagnostic experiments to assess the influence of initial soil moisture and precipitation magnitude on flood simulation over the Cape Fear River basin in the United States. Model results suggest that: (1) The modulation effect of initial soil moisture on the flood peak is non-linear and weakens as precipitation magnitude increases. There is a threshold value of the soil saturation, below and above which the sensitivity of flood peak to the soil moisture differentiates substantially; (2) For model spin-up, streamflow needs longer time to reach the ‘practical’ equilibrium (10%) than the soil moisture and latent heat flux. The model uncertainty from spin-up can propagate through the hydrometeorological modelling chain and get amplified into the flood peak; (3) For ensemble flood modelling with a hydrometeorological system, modelling uncertainty is dominated by the precipitation forecast. Spin-up induced uncertainty can be minimized once the model reaches the ‘practical’ equilibrium.

North Carolina

Preliminary investigation of soil and ground-water contamination at a U.S. Army Petroleum Training Facility, Fort Lee, Virginia, September-October 1989

Fuel-oil constituents in the soil and groundwater at the Fort Lee Petroleum Training Facility near Petersburg, Virginia, were studied by the U.S. Geological Survey (USGS) in cooperation with the Department of Defense, U.S. Army. The study included installation of 25 groundwater monitoring wells and description of groundwater flow patterns of the shallow-aquifer system underlying the facility. Soil and groundwater samples were collected to determine the concentrations of fuel-oil constituents and to determine the potential for off-site migration of the constituents. Total petroleum hydrocarbon concentrations up to 18,400 mg/km were reported in soil samples. Concentrations of benzene in water from wells at the facility were up to 130 micrograms per liter (ug/L), and concentrations of ethylbenzene and xylene were up to 54 and 120 ug/L, respectively. Potential exists for off-site migration of the contaminants and migration of contaminants downward to deeper aquifers. Further investigations of these potential contamination-migration pathways are warranted. Risk identification at the Petroleum Training Facility cannot be properly addressed because the distribution of the fuel-oil constituents has not been fully characterized. Preliminary identification of risk, however is presented by an examination of toxicity data for the chemical constituents reported in the groundwater at the facility. Concentrations of constituents were compared to the maximum contaminant levels (MCLs) for drinking water established by the U.S. Environmental Protection Agency (USEPA). Concentrations of benzene in water from wells at the facility exceed the USEPA 's 5 ug/L MCL by as much as 26 times. Sufficient data are not available to fully design the remedial-action plan for the facility; however, general responses to contamination of the type associated with the facility include no-action, monitoring, institutional controls, removal, and treatment.

Virginia

Effects of flood inundation, invasion by Phalaris arundinacea, and nitrogen enrichment on extracellular enzyme activity in an Upper Mississippi River floodplain forest

The community structures and ecosystem functions of floodplains are primarily driven by variation in flood inundation. However, global changes, such as invasive species and nutrient enrichment, may alter the effects of flooding in these systems. We added nitrogen (N) to correspond with twice the annual atmospheric deposition rate of the south-west Wisconsin, USA region within mature floodplain forest plots and patches of an invasive grass (reed canarygrass, Phalaris arundinacea ) along a floodplain elevation gradient in an Upper Mississippi River floodplain forest. We measured soil physicochemical properties and the activity of six extracellular enzymes during 3 months that varied in flooding conditions. Multivariate analyses (distance-based redundancy analysis) revealed that floodplain elevation, month of sampling, and vegetation type were all significant predictors of variation in soil physicochemical properties, while elevation and month were significant predictors of multivariate extracellular enzyme activity (EEA). The best model for predicting EEA consisted of nitrogen availability, soil porosity, and water filled pore space. Although the categorical fertilization and invasion treatments were not significant predictors of EEA, our results suggest that their effects depend on the degree to which they modify N availability and soil moisture. In this system, spatial and temporal patterns in flooding appear to be the main driver of these properties, but N enrichment and invasion may have the potential to further modify them.

Wisconsin

Binational digital soils map of the Ambos Nogales watershed, southern Arizona and northern Sonora, Mexico

We have prepared a digital map of soil parameters for the international Ambos Nogales watershed to use as input for selected soils-erosion models. The Ambos Nogales watershed in southern Arizona and northern Sonora, Mexico, contains the Nogales wash, a tributary of the Upper Santa Cruz River. The watershed covers an area of 235 km2, just under half of which is in Mexico. Preliminary investigations of potential erosion revealed a discrepancy in soils data and mapping across the United States-Mexican border due to issues including different mapping resolutions, incompatible formatting, and varying nomenclature and classification systems. To prepare a digital soils map appropriate for input to a soils-erosion model, the historical analog soils maps for Nogales, Ariz., were scanned and merged with the larger-scale digital soils data available for Nogales, Sonora, Mexico using a geographic information system.

Arizona, Sonora

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

Stratigraphy of the Younger Dryas Chronozone and paleoenvironmental implications: Central and Southern Great Plains

The Great Plains of the United States was the setting for some of the earliest research in North America into patterns and changes in the character of late Pleistocene environments and their effects on contemporary human populations. Many localities in the region have well-stratified records of terminal Pleistocene and early Holocene human (Paleoindian) activity and past environments. These have proven important in debates over the character of the Younger Dryas Chronozone (YDC; 11,000–10,000 14 C BP; 12,900–11,700 cal BP) in the continental interior. This paper reviews the lithostratigraphic record of the YDC on the Central and Southern Great Plains and summarizes paleobiological records (largely isotopic). The goal is to determine if there is any uniformity in the timing, character, direction and/or magnitude of changes in depositional environments or broader geomorphic systems before, during or after the YDC in order to address the question of the character of environments through this time. The stratigraphic records of the late Pleistocene to early Holocene transition, and in particular, the stratigraphic records of the YDC vary through time and space. The data clearly show that a host of geomorphic processes produced the terminal Pleistocene and early Holocene stratigraphic records of the Great Plains. Moreover, the YDC is not necessarily manifest as a distinct lithostratigraphic or biostratigraphic entity in these different types of deposits and soils. The various geomorphic systems of the Great Plains did not behave synchronously in response to any common climate driver. These stratigraphic records reflect local environmental conditions and probably a complex response to the reorganization of mid-latitude climates in the terminal Pleistocene and early Holocene.

Montana, Wyoming, Colorado, New Mexico, Texas, Okl

Emission of pesticides into the air

During and after the application of a pesticide in agriculture, a substantial fraction of the dosage may enter the atmosphere and be transported over varying distances downwind of the target. The rate and extent of the emission during application, predominantly as spray particle drift, depends primarily on the application method (equipment and technique), the formulation and environmental conditions, whereas the emission after application depends primarily on the properties of the pesticide, soils, crops and environmental conditions. The fraction of the dosage that misses the target area may be high in some cases and more experimental data on this loss term are needed for various application types and weather conditions. Such data are necessary to test spray drift models, and for further model development and verification as well. Following application, the emission of soil fumigants and soil incorporated pesticides into the air can be measured and computed with reasonable accuracy, but further model development is needed to improve the reliability of the model predictions. For soil surface applied pesticides reliable measurement methods are available, but there is not yet a reliable model. Further model development is required which must be verified by field experiments. Few data are available on pesticide volatilization from plants and more field experiments are also needed to study the fate processes on the plants. Once this information is available, a model needs to be developed to predict the volatilization of pesticides from plants, which, again, should be verified with field measurements. For regional emission estimates, a link between data on the temporal and spatial pesticide use and a geographical information system for crops and soils with their characteristics is needed.

Conference Paper

Can grazing by elk and bison stimulate herbaceous plant productivity in semiarid ecosystems?

Plant communities in rangeland ecosystems vary widely in the degree to which they can compensate for losses to herbivores. Ecosystem-level factors have been proposed to affect this compensatory capacity, including timing and intensity of grazing, and availability of soil moisture and nutrients. Arid ecosystems are particularly challenging to predict because of their high degree of temporal variability in moisture inputs. We used a replicated herbivore exclusion experiment to evaluate herbaceous plant responses to grazing by large ungulates to test current theory and identify constraints on plant compensation in a dryland ecosystem. We measured nitrogen (N) yield and herbaceous production in three plant communities: meadows, willow-associated herbaceous communities, and riparian communities. We implemented grazing exclusion treatments from 2005 to 2008 in areas with elk and bison and areas with only elk. Grazing by large ungulates increased herbaceous production and N yield in herbaceous riparian communities. In willow communities, herbaceous plants displayed equal compensation in response to grazing in total aboveground production and N yield. Our results support the idea that plant compensation in this semiarid system is contingent on soil moisture availability, wherein the most productive sites (that received substantial moisture inputs from subsurface flow) exhibited overcompensation. Although the herbaceous riparian communities we studied are isolated patches of productive grassland in an otherwise shrub-dominated and minimally productive semiarid landscape, grazing by a combination of bison and elk removed only 44%–53% of aboveground net primary productivity (ANPP) during the growing season, and 25%–38% of production over winter. Consumption by ungulates was a positive linear function of herbaceous production, similar to reported patterns from other temperate and tropical grazing ecosystems. The slope of this relationship was affected by the analytical method used to calculate ANPP and consumption rates, but, regardless of the method, was lower or similar to reported slopes for other intensively grazed systems (Yellowstone, Serengeti, Laikipia) that have sustained high ungulate densities for decades to centuries. Given that the vegetation communities exhibited equal or overcompensation in terms of total herbaceous ANPP in both years, elk and bison population levels during our study period did not appear to occur at densities leading to degradation of herbaceous communities.

Colorado

Assessing nonpoint-source uranium pollution in an irrigated stream-aquifer system

Uranium (U) in rocks and soils of arid and semi-arid environments can be mobilized by irrigation and fertilization, posing environmental and health risks. Elevated U, along with selenium (Se) and nitrate (NO 3 ) co-constituents, necessitates careful monitoring and management. We developed a distributed-parameter numerical model to assess U pollution in an irrigated stream-aquifer system, applying it to a 552 km 2 region in Colorado's Lower Arkansas River Valley (LARV) over 14 years. A MODFLOW model, describing groundwater and stream flow, was coupled with an RT3D-OTIS model to portray reactive U transport. Calibration using the PESTPP-iES iterative ensemble smoother (iES) software indicated good agreement with observed U concentrations. The model revealed substantial and variable U levels across the LARV, highlighting potential hotspots and possible contributing factors, such as geological composition of the bedrock and near-surface shale and aquifer sediments derived from them, irrigation practices, and riparian landscape. U levels exceed the chronic standard (85th percentile = 30 μg/L, set by the US Environmental Protection Agency), which is the permissible regulatory threshold, in groundwater across 44 % of the region and along the river by an average factor of 2.9. Simulated average U concentrations in the non-riparian aquifer and river are 124 μg/L and 60 μg/L, respectively, compared with 112 μg/L and 62 μg/L for measured values. The average 85th percentile U concentration is 222 μg/L in the aquifer and 82 μg/L in the river. Average simulated U mass loading to the river is 0.17 kg/day per km, compared to an estimated 0.23 kg/day per km. Findings provide a baseline for comparing future simulated outcomes of alternative best management practices (BMPs) for U pollution mitigation and offer a methodology applicable to other irrigated regions.

Colorado

Hydrogeomorphology influences soil nitrogen and phosphorus mineralization in floodplain wetlands

Conceptual models of river–floodplain systems and biogeochemical theory predict that floodplain soil nitrogen (N) and phosphorus (P) mineralization should increase with hydrologic connectivity to the river and thus increase with distance downstream (longitudinal dimension) and in lower geomorphic units within the floodplain (lateral dimension). We measured rates of in situ soil net ammonification, nitrification, N, and P mineralization using monthly incubations of modified resin cores for a year in the forested floodplain wetlands of Difficult Run, a fifth order urban Piedmont river in Virginia, USA. Mineralization rates were then related to potentially controlling ecosystem attributes associated with hydrologic connectivity, soil characteristics, and vegetative inputs. Ammonification and P mineralization were greatest in the wet backswamps, nitrification was greatest in the dry levees, and net N mineralization was greatest in the intermediately wet toe-slopes. Nitrification also was greater in the headwater sites than downstream sites, whereas ammonification was greater in downstream sites. Annual net N mineralization increased with spatial gradients of greater ammonium loading to the soil surface associated with flooding, soil organic and nutrient content, and herbaceous nutrient inputs. Annual net P mineralization was associated negatively with soil pH and coarser soil texture, and positively with ammonium and phosphate loading to the soil surface associated with flooding. Within an intensively sampled low elevation flowpath at one site, sediment deposition during individual incubations stimulated mineralization of N and P. However, the amount of N and P mineralized in soil was substantially less than the amount deposited with sedimentation. In summary, greater inputs of nutrients and water and storage of soil nutrients along gradients of river–floodplain hydrologic connectivity increased floodplain soil nutrient mineralization rates.

Maryl;Virginia

Hydrologic study of green infrastructure in poorly drained urbanized soils at RecoveryPark, Detroit, Michigan, 2014–21

Uncontrolled stormwater runoff volume is a legacy stressor on sewer-system capacity that is further compromised by the effects of aging infrastructure. Green stormwater infrastructure (GSI) has been used in a variety of designs and configurations (for example, bioretention) with the goal of increasing evapotranspiration and infiltration in the local water cycle. In practice, GSIs have variable effectiveness in reducing runoff volume. An urban residential site near Detroit, Michigan, called RecoveryPark was monitored for 8 years before and after GSI construction to evaluate how effectively the GSI reduced volumes of stormwater flowing to Detroit’s Water Resource Recovery Facility through combined sewer systems. In addition to the GSI, the study site included an urban farm where salad crops were grown in hoop houses. The monitoring approach was to characterize the urban water cycle through high-frequency measurements of inflows and outflows. Datasets included meteorological data, soils and sediment characteristics, groundwater levels, flows within the combined sewer system, and soils and water chemistry with specific focus on the disposition of road salt. Although land cover within the RecoveryPark sewershed was high-density residential in the 1950s, the sewershed included only one residence within the 8.74-acre sewershed during this study. Measurements of annual precipitation at the site exceeded long-term annual averages by more than 10 inches during 3 of the 8 years of study. Potential evapotranspiration was often greater than the measured precipitation that averaged 28–34 inches per year. As compared to underlying clay-rich sediments, soils data indicated relatively permeable sediments near land surface with estimated hydraulic conductivity of 0.75 inches per hour; however, these values decreased with increasing depth. Groundwater-level data revealed increases in groundwater storage as indicated by increases in seasonal groundwater levels and development of a groundwater mound adjacent to the GSI. These increases in groundwater levels were directly adjacent to swales designed to infiltrate stormwater and only became evident after installing the GSI. Flows within the combined sewer system included rainwater, septic effluent, groundwater infiltration, leakage from water-supply lines, and release of water stored in abandoned foundations. Dry-weather flows (no rain fell within the prior 3 days) averaged 7–10 gallons per minute, which were much greater than flows estimated by septic outflow alone. A set of estimated water budgets were compiled, and results showed large discrepancies in unaccounted flows. To further examine these discrepancies, dye-tracing within the combined sewer system helped examine the sources of water by relating flow volumes to drainage area. For one of the monitoring sites within the combined sewer system along the southeast side of the study area, flows estimated by dye concentrations were more than 10 percent greater than those measured by standard methods. Through peak-flow-regression analysis, a minimum of 2.4 million gallons of water per year were infiltrated or lost to evapotranspiration because of GSI construction. After site modifications were made by excavating gravel drains to improve drainage characteristics, estimated stormwater volumes within the combined sewer system returned to near preconstruction levels. The GSI was effectively bypassed to address slow infiltration rates and standing water; the bypass all but eliminated the potential benefits of volume reduction. Late in the project, a water-quality study was added to examine the transport of road salt and associated chloride within the GSI and the combined sewer system. Continuous specific conductance was used as a surrogate for chloride concentrations to estimate that 2,790 pounds of dissolved chloride passed through the sewershed during the winter months of late 2020 through early 2021. These data were collected after GSI modification, therefore most, if not all, of the chloride was transported directly to Detroit’s Water Resource Recovery Facility via the combined sewer system. Mixing diagrams using chloride and bromide concentrations of road salt, potable water, rainwater, groundwater, and water from the combined sewer system confirmed that water within the combined sewer system is a mix of these sources. The poor condition of the combined sewer system pipes and resulting unaccounted inflows added to the challenge of accurately monitoring and identifying sources and sinks of water within the RecoveryPark sewershed. Our research results suggest that—along with clear and quantifiable objectives—the catchment and site conditions should be well-characterized before determining the GSI design. In addition, the work presented in this report provides implications and lessons learned for effectiveness and future studies of GSI in urban settings. These efforts can be improved through increased communication between stakeholders, use of high-quality soils in GSI that have suitable hydraulic characteristics, redundant data-collection networks for critical data streams, and focusing meteorological-data collection within the GSI to obtain relevant evapotranspiration data.

Michigan

Archaeal diversity and CO2 fixers in carbonate-/siliciclastic-rock groundwater ecosystems

Groundwater environments provide habitats for diverse microbial communities, and although Archaea usually represent a minor fraction of communities, they are involved in key biogeochemical cycles. We analysed the archaeal diversity within a mixed carbonate-rock/siliciclastic-rock aquifer system, vertically from surface soils to subsurface groundwater including aquifer and aquitard rocks. Archaeal diversity was also characterized along a monitoring well transect that spanned surface land uses from forest/woodland to grassland and cropland. Sequencing of 16S rRNA genes showed that only a few surface soil-inhabiting Archaea were present in the groundwater suggesting a restricted input from the surface. Dominant groups in the groundwater belonged to the marine group I (MG-I) Thaumarchaeota and the Woesearchaeota. Most of the groups detected in the aquitard and aquifer rock samples belonged to either cultured or predicted lithoautotrophs (e.g., Thaumarchaeota or Hadesarchaea). Furthermore, to target autotrophs, a series of 13 CO 2 stable isotope-probing experiments were conducted using filter pieces obtained after filtration of 10,000 L of groundwater to concentrate cells. These incubations identified the SAGMCG Thaumarchaeota and Bathyarchaeota as groundwater autotrophs. Overall, the results suggest that the majority of Archaea on rocks are fixing CO 2 , while archaeal autotrophy seems to be limited in the groundwater.

Archaea

Soil phosphorus cycling in tropical soils: An ultisol and oxisol perspective

Phosphorus (P) is essential for life. It is the backbone of our DNA, provides energy for biological reactions, and is an integral component of cell membranes. As such, it is no surprise that P availability plays a strong role in regulating ecosystem structure and function (Wassen et al. 2005, Elser et al. 2007, Condit et al. 2013), and in determining our capacity to grow food for a burgeoning human population (Sharpley et al. 1997, Sims and Sharpley 2005, Lal 2009). Concerns that P supplies are insufficient to meet our species’ growing demands are on the rise (Richardson and Simpson 2011) and scientific and media outlets increasingly discuss P as an element worthy of our attention and concern (e.g., Cordell et al. 2009, Lougheed 2011, Edixhoven et al. 2013, Ulrich et al. 2013). Indeed, a number of groups are calling for the explicit stewardship of our planet’s P stocks (Schipper 2014, Withers et al. 2015). Yet a focus on P as a vital and limited resource is not new in the tropics, where an abundance of soils characterized by low P has resulted in a substantial, longstanding reliance on P inputs for tropical ecosystem function in both unmanaged and agriculture settings (Table 1, Figure 2; Sanchez 1976, Swap et al. 1992, Chadwick et al. 1999, Okin et al. 2004, Lal 2009). Indeed, there is a long history of cultivation in the tropics, where for thousands of years land management practices have included methods that effectively modify P availability for plant growth (e.g., Giardina et al. 2000, Lawrence and Schlesinger 2001, Vitousek et al. 2004, Lewis et al. 2015). Nevertheless, low soil fertility in tropical systems where fertilizer is scarce has enduringly been recognized as a major source of hunger and starvation (Sanchez and Buol 1975, Sanchez 2002, Sanchez and Swaminathan 2005).

Book chapter