SEARCH · Search USGS

Results for “Earth”

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

SourceResearch connectionFrom yearThrough year

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

At least 343 records · Page 19Linked to original sources

USGS2020

Mapping croplands of Europe, Middle East, Russia, and Central Asia using Landsat 30-m data, machine learning algorithms and Google Earth Engine

Accurate and timely information on croplands is important for environmental, food security, and policy studies. Spatially explicit cropland datasets are also required to derive information on crop type, crop yield, cropping intensity, as well as irrigated areas. Large area defined as continental to global cropland mapping is challenging due to differential manifestation of croplands, wide range of cultivation practices and limited reference data availability. This study presents the results of a cropland extent mapping of 64 Countriescovering large parts of Europe, Middle East, Russia and Central Asia. To cover such a vast area, roughly 160,000 Landsat scenes from 3,351 footprints between 2014 and 2016 were processed within the Google Earth Engine (GEE) cloud-platform. We used the pixel-based supervised Random Forest (RF) machine learning algorithm with a set of satellite data inputs capturing diverse spectral, temporal and topographical characteristics across twelve agroecological zones (AEZs). The reference data to train the classification model were collected from very high spatial resolution imagery (VHRI) and ancillary datasets. The result is a binary map showing cultivated/non-cultivated areas ca. 2015. The map produced an overall accuracy of 94 percent with roughly 14 percent omission and commission errors for the cropland class based on a large set of independent validation samples. The map suggests the entire study area has a total 546 million hectares (Mha) of croplands occupying 18 percent of the land area. Comparison between national cropland area estimates from United Nations Food and Agricultural Organizations (FAO) and those derived from this work also showed an R-square value of 0.95. For the entire Landsat-derived 30-m product the overall accuracy was 93.8% with cropland class providing producers accuracy of 86.5% (errors of omissions = 13.5%) and users accuracy of 85.7% (errors of commissions = 14.3%). This Landsat-derived 30-m cropland product (GFSAD30) provided 10-30% greater cropland areas compared to UN FAO in the 64 Countries. Finally, the map-to-map comparison between GFSAD30 with several other cropland products revealed that the best similarity matrix was with the 30m global land cover (GLC30) product providing an overall accuracy of 88.8 percent (Kappa 0.7) with producers cropland similarity of 89.2 percent (errors of omissions = 10.8%) and users cropland similarity of 81.8 percent (errors of commissions = 8.1%). GFSAD30 captured the missing croplands in GLC30 product around significantly irrigated agricultural areas in Germany and Belgium and rainfed agriculture in Italy. This study also established that the real strength of GFSAD30 product, compared to other products, were in: 1. Identifying precise location of croplands, and 2. Capturing fragmented croplands. The cropland extent map dataset is available through NASAs Land Processes Distributed Active Archive Center (LP DAAC) at https://doi.org/10.5067/MEaSUREs/GFSAD/GFSAD30EUCEARUMECE.001, while the training and reference data as well as visualization are available at the Global Croplands website.

ISPRS Journal of Photogrammetry and Remote Sensing↗

USGS2006

The behavior of rare earth elements in naturally and anthropogenically acidified waters

In this paper, the behavior of rare earth elements (REE) in a watershed impacted by acid-mine drainage (Fisher Creek, Montana) is compared to that in a volcanically acidified watershed (Rio Agrio and Lake Caviahue, Argentina). The REE behave conservatively in acidic waters with pH values less than approximately 5.5. However, above pH 5.5, REE concentrations are controlled by adsorption onto or co-precipitation with a variety of Fe or Al oxyhydroxides. The heavy REE partition to a greater extent into the solid phase than the light REE as pH rises above 6. Concentrations of REE exhibit diel (24-h) cycling in waters that were initially acidic, but have become neutralized downstream. In Fisher Creek, at the most downstream sampling station investigated (pH 6.8), concentrations of dissolved REE were 190–840% higher in the early morning versus the late afternoon. This cycling can be related to temperature-dependent, cyclic adsorption–desorption of REE onto hydrous ferric or aluminum oxide or both. Similar but gentler diel cycling of the REE was found at Rio Agrio. The existence of such cycling has important ramifications for the study of REE in natural waters.

Journal of Alloys and Compounds↗

USGS2025

Alteration mapping in granitic gneiss using handheld geophysical and geochemical instruments: Implications for iron oxide-apatite and rare earth elements exploration

The Adirondack Mountains of New York, U.S.A. contain iron oxide-apatite (IOA) mineral deposits with variable concentrations of rare earth elements (REE). The IOA mineral deposits are typically hosted in the Lyon Mountain Granite Gneiss and are spatially correlated with extensive Na metasomatism (albitization) of the surrounding country rocks, although some mineral deposits also occur in metagabbro, paragneiss, and anorthosite hosts. The location of albitization zones is key to finding new IOA mineral deposits and better understanding the mineralization processes associated with their genesis. However, the Na alteration zones are generally not visibly identifiable in outcrop or hand specimens because the color and textural changes are subtle and are thus difficult to map. Here, we discuss the results from testing two spectroscopic techniques (handheld, field gamma-ray spectroscopy, and portable X-ray fluorescence (pXRF)) to characterize albitization zones surrounding IOA mineral deposits. The gamma-ray and pXRF spectrometers successfully distinguished the relative intensity of alteration in the Lyon Mountain Granite Gneiss based on the K abundance. The measured K content decreases towards the IOA mineral deposits, and the magnitude and width of the gradient are similar along the strike of each mineral deposit. Elevated Th and U values are present in host rocks adjacent to REE-bearing IOA mineral deposits. The pXRF and gamma-ray spectrometer K values are remarkably consistent with laboratory-based, whole-rock XRF compositional data and, therefore, useful for semi-quantitative analysis. Notably, albitization aureoles are consistent around REE-rich and REE-poor mineral deposits, suggesting that saline fluids are essential for the petrogenesis of IOA mineral deposits but may not be critical to REE mineralization. Ultimately, this study demonstrates the utility of handheld gamma-ray and pXRF spectrometry for identifying otherwise cryptic albitization gradients associated with IOA mineral deposits in granitic gneiss.

New York↗

USGS2019

Assessing plant production responses to climate across water-limited regions using Google Earth Engine

(Munson) Climate variability and change acting at broad scales can lead to divergent changes in plant production at local scales. Quantifying how production responds to variation in climate at local scales is essential to understand underlying ecological processes and inform land management decision-making, but has historically been limited in spatiotemporal scale based on the use of discrete ground-based measurements or coarse resolution satellite observations. With the advent of cloud-based computing through Google Earth Engine (GEE), production responses to climate can be evaluated across broad landscapes though time at a resolution useful for ecological and land management applications. Here, GEE was employed to synthesize a multi-platform Landsat time series (1988 – 2014) and evaluate relationships between the soil-adjusted vegetation index (a proxy for plant production) and climate across deserts and plant communities of the southwestern U.S. A “climate pivot point” approach was adopted in GEE to assess the trade-off between production responses to increasing wetness and resistances to drought at 30-m resolution. Consistent with a long-term seasonal climate gradient, production was most related to climate variance during the cool-season in the western deserts, during the warm-season in the eastern deserts, and equally related to both seasons within several desert areas. Communities dominated by grasses and deciduous trees displayed large production responses to an increase in wetness and low resistances to water deficit, while shrublands and evergreen woodlands had variable responses and high drought resistances. Production in plant communities that spanned multiple deserts responded differently to seasonal climate variability in each desert. Defining these plant production sensitivities to climate at 30-m resolution in GEE advances forecasts of how long-term climate trajectories may affect carbon storage, wildlife habitat, and the vulnerability of water-limited ecosystems.

Arizona, California, Colorado, Nevada, New Mexico,↗

USGS2022

Mapping actual evapotranspiration using Landsat for the conterminous United States: Google Earth Engine implementation and assessment of the SSEBop model

The estimation and mapping of actual evapotranspiration (ETa) is an active area of applied research in the fields of agriculture and water resources. Thermal remote sensing-based methods, using coarse resolution satellites, have been successful at estimating ETa over the conterminous United States (CONUS) and other regions of the world. In this study, we present CONUS-wide ETa from Landsat thermal imagery-using the Operational Simplified Surface Energy Balance (SSEBop) model in the Google Earth Engine (GEE) cloud computing platform. Over 150,000 Landsat satellite images were used to produce 10 years of annual ETa (2010–2019) at unprecedented scale. The accuracy assessment of the SSEBop results included point-based evaluation using monthly Eddy Covariance (EC) data from 25 AmeriFlux stations as well as basin-scale comparison with annual Water Balance ETa (WBET) for more than 1000 sub-basins. Evaluations using EC data showed generally mixed performance with weaker (R 2 < 0.6) correlation on sparsely vegetated surfaces such as grasslands or woody savanna and stronger correlation (R 2 > 0.7) over well-vegetated surfaces such as croplands and forests, but location-specific conditions rather than cover type were attributed to the variability in accuracy. Croplands performed best with R 2 of 0.82, root mean square error of 29 mm/month, and average bias of 12%. The WBET evaluation indicated that the SSEBop model is strong in explaining the spatial variability (up to R 2 > 0.90) of ETa across large basins, but it also identified broad hydro-climatic regions where the SSEBop ETa showed directional biases, requiring region-specific model parameter improvement and/or bias correction with an overall 7% bias nationwide. Annual ETa anomalies over the 10-year period captured widely reported drought-affected regions, for the most part, in different parts of the CONUS, indicating their potential applications for mapping regional- and field-scale drought and fire effects. Due to the coverage of the Landsat Path/Row system, the availability of cloud-free image pixels ranged from less than 12 (mountainous cloud-prone regions and U.S. Northeast) to more than 60 (U.S. Southwest) per year. However, this study reinforces a promising application of Landsat satellite data with cloud-computing for quick and efficient mapping of ETa for agricultural and water resources assessments at the field scale.

Remote Sensing of Environment↗

USGS2019

Reimagining the potential of Earth observations for ecosystem service assessments

The benefits nature provides to people, called ecosystem services, are increasingly recognized and accounted for in assessments of infrastructure development, agricultural management , conservation prioritization, and sustainable sourcing. These assessments are often limited by data, however, a gap with tremendous potential to be filled through Earth observations (EO), which produce a variety of data across spatial and temporal extents and resolutions. Despite widespread recognition of this potential, in practice few ecosystem service studies use EO. Here, we identify challenges and opportunities to using EO in ecosystem service modeling and assessment. Some challenges are technical, related to data awareness, processing, and access. These challenges require systematic investment in model platforms and data management. Other challenges are more conceptual but still systemic; they are byproducts of the structure of existing ecosystem service models and addressing them requires scientific investment in solutions and tools applicable to a wide range of models and approaches. We also highlight new ways in which EO can be leveraged for ecosystem service assessments, identifying promising new areas of research. More widespread use of EO for ecosystem service assessment will only be achieved if all of these types of challenges are addressed. This will require non-traditional funding and partnering opportunities from private and public agencies to promote data exploration, sharing, and archiving. Investing in this integration will be reflected in better and more accurate ecosystem service assessments worldwide.

Science of the Total Environment↗

USGS2002

Rare earth element variations resulting from inversion of pigeonite and subsolidus reequilibration in lunar ferroan anorthosites

We present results of a secondary ion mass spectrometry study of the rare earth elements (REEs) in the minerals of two samples of lunar ferroan anorthosite, and the results are applicable to studies of REEs in all igneous rocks, no matter what their planet of origin. Our pyroxene analyses are used to determine solid-solid REE distribution coefficients (D = C REE in low-Ca pyroxene/C REE in augite) in orthopyroxene-augite pairs derived by inversion of pigeonite. Our data and predictions from crystal-chemical considerations indicate that as primary pigeonite inverts to orthopyroxene plus augite and subsolidus reequilibration proceeds, the solid-solid Ds for orthopyroxene-augite pairs progressively decrease for all REEs; the decrease is greatest for the LREEs. The REE pattern of solid-solid Ds for inversion-derived pyroxene pairs is close to a straight line for Sm-Lu and turns upward for REEs lighter than Sm; the shape of this pattern is predicted by the shapes of the REE patterns for the individual minerals. Equilibrium liquids calculated for one sample from the compositions of primary phases, using measured or experimentally determined solid-liquid Ds, have chondrite-normalized REE patterns that are very slightly enriched in LREEs. The plagioclase equilibrium liquid is overall less rich in REEs than pyroxene equilibrium liquids, and the discrepancy probably arises because the calculated plagioclase equilibrium liquid represents a liquid earlier in the fractionationsequence than the pyroxene equilibrium liquids. “Equilibrium” liquids calculated from the compositions of inversion-derived pyroxenes or orthopyroxene derived by reaction of olivine are LREE depleted (in some cases substantially) in comparison with equilibrium liquids calculated from the compositions of primary phases. These discrepancies arise because the inversion-derived and reaction-derived pyroxenes did not crystallize directly from liquid, and the use of solid-liquid Ds is inappropriate. The LREE depletion of the calculated liquids is a relic of formation of these phases from primary LREE-depleted minerals. Thus, if one attempts to calculate the compositions of equilibrium liquids from pyroxene compositions, it is important to establish that the pyroxenes are primary. In addition, our data suggest that experimental studies have underestimated solid-liquid Ds for REEs in pigeonite and that REE contents of liquids calculated using these Ds are overestimates. Our results have implications for Sm-Nd age studies. Our work shows that if pigeonite inversion and/or subsolidus reequilibration between augite and orthopyroxene occurred significantly after crystallization, and if pyroxene separates isolated for Sm-Nd studies do not have the bulk composition of the primary pyroxenes, then the Sm-Nd isochron age and ε Nd will be in error.

Geochimica et Cosmochimica Acta↗

USGS1998

Rare earth element metasomatism in hydrothermal systems: The Willsboro-Lewis wollastonite ores, New York, USA

Wollastonite ores and garnet-pyroxene skarns in the Willsboro-Lewis district, New York, USA were formed in a complex hydrothermal system associated with the emplacement of a large anorthosite pluton. Contact-metamorphic marbles were replaced by wollastonite, garnet, and clinopyroxene during infiltration metasomatism involving large volumes of water of chiefly meteoric origin. Rare earth elements (REE) in these rocks show large departures from the protolith REE distribution, indicative of substantial REE mobility. Three types of chondrite-normalized REE distribution patterns are present. The most common, found in ores and skarns containing andradite-rich garnet, is convex-up in the light REE (LREE) with a maximum at Pr and a positive Eu anomaly. Europium anomalies and Pr/Yb ratios are correlated with X(Ad) in garnet. This pattern (type C) results from uptake of REE from hydrothermal fluids by growing crystals of calcsilicate minerals, principally andradite, with amounts of LREE controlled by the difference in ionic radius between Ca++ and REE3+ in garnet X sites. The Eu anomaly results either from prior interaction of the fluids with plagioclase-rich, Eu-positive anorthositic rocks in and near the ore zone, or by enrichment of divalent Eu on growth surfaces of garnet followed by entrapment, or both. Relative enrichment in heavy REE (type H) occurs in ores and skarn where calcsilicates, including grossularitic garnet, in contact-metamorphic marble have been concentrated by dissolution of calcite. In most cases a negative Eu anomaly is inherited from the marble protolith. Skarns containing titanite and apatite exhibit high total REE, relative light REE enrichment, and negative Eu anomalies (type L). These appear to be intrusive igneous rocks (ferrodiorites or anorthositic gabbros) that have been converted to skarn by Ca metasomatism. REE, sequestered in titanite, apatite, and garnet, preserve the approximate REE distribution pattern of the igneous protolith. Post-ore granulite facies metamorphism homogenized zoned mineral grains without causing complete intergranular reequilibration and does not appear to have significantly affected the whole-rock REE distributions. These results demonstrate that extensive REE metasomatism can occur in hydrothermal systems at shallow to intermediate depths and that REE geochemistry may be useful in discerning the origin of skarns and skarn-related ore deposits.

Geochimica et Cosmochimica Acta↗

USGS1999

A topological system for delineation and codification of the Earth's river basins

A comprehensive reference system for the Earth's river basins is proposed as a support to fiver basin management, global change research, and the pursuit of sustainable development. A natural system for delineation and codification of basins is presented which is based upon topographic control and the topology of the fiver network. These characteristics make the system well suited for implementation and use with digital elevation models (DEMs) and geographic information systems. A demonstration of these traits is made with the 30-arcsecond GTOPO30 DEM for North America. The system has additional appeal owing to its economy of digits and the topological information that they carry. This is illustrated through presentation of comparisons with USGS hydrologic unit codes and demonstration of the use of code numbers to reveal dependence or independence of water use activities within a basin.

Journal of Hydrology↗

USGS2002

34 - State of stress in the Earth's lithosphere

The state of stress in the lithosphere is the result of the forces acting upon and within it. Knowledge of the magnitude and distribution of these forces can be combined with mechanical, thermal and rheological constraints to examine a broad range of lithospheric deformational processes. For example, such knowledge contributes to a better understanding of the processes that both drive and inhibit lithospheric plate motions as well as the forces responsible for the occurrence of crustal earthquakes − both along plate boundaries and in intraplate regions. Although the topic of this chapter is the state of stress in the Earth's lithosphere, the comments below come primarily from the perspective of the state of stress in the brittle upper crust. As defined by the depth of shallow earthquakes, the brittle crust extends to ˜15–20 km depth at most continental locations around the world. This perspective is adopted here because nearly all the data available on lithospheric stress come from the upper crust of continents. Furthermore, in the sections that follow, we argue that, to first order, the state of stress in the brittle crust results from relatively large-scale lithospheric processes so that knowledge of crustal stress can be used to constrain the forces involved in these processes.

International Geophysics↗

USGS2001

Standard reference water samples for rare earth element determinations

Standard reference water samples (SRWS) were collected from two mine sites, one near Ophir, CO, USA and the other near Redding, CA, USA. The samples were filtered, preserved, and analyzed for rare earth element (REE) concentrations (La, Ce, Pr, Nd, Sm, Eu, Gd, Tb, Dy, Ho, Er, Tm, Yb, and Lu) by inductively coupled plasma-mass spectrometry (ICP-MS). These two samples were acid mine waters with elevated concentrations of REEs (0.45–161 μg/l). Seventeen international laboratories participated in a ‘round-robin’ chemical analysis program, which made it possible to evaluate the data by robust statistical procedures that are insensitive to outliers. The resulting most probable values are reported. Ten to 15 of the participants also reported values for Ba, Y, and Sc. Field parameters, major ion, and other trace element concentrations, not subject to statistical evaluation, are provided.

Applied Geochemistry↗

USGS2010

Evaluating the behavior of gadolinium and other rare earth elements through large metropolitan sewage treatment plants

A primary pathway for emerging contaminants (pharmaceuticals, personal care products, steroids, and hormones) to enter aquatic ecosystems is effluent from sewage treatment plants (STP), and identifying technologies to minimize the amount of these contaminants released is important. Quantifying the flux of these contaminants through STPs is difficult. This study evaluates the behavior of gadolinium, a rare earth element (REE) utilized as a contrasting agent in magnetic resonance imaging (MRI), through four full-scale metropolitan STPs that utilize several biosolids thickening, conditioning, stabilization, and dewatering processing technologies. The organically complexed Gd from MRIs has been shown to be stable in aquatic systems and has the potential to be utilized as a conservative tracer in STP operations to compare to an emerging contaminant of interest. Influent and effluent waters display large enrichments in Gd compared to other REEs. In contrast, most sludge samples from the STPs do not display Gd enrichments, including primary sludges and end-product sludges. The excess Gd appears to remain in the liquid phase throughout the STP operations, but detailed quantification of the input Gd load and residence times of various STP operations is needed to utilize Gd as a conservative tracer.

Environmental Science & Technology↗

USGS2004

Terrestrial microorganisms at an altitude of 20,000 m in Earth's atmosphere

A joint effort between the U.S. Geological Survey's (USGS) Global Desert Dust and NASA's Stratospheric and Cosmic Dust Programs identified culturable microbes from an air sample collected at an altitude of 20,000 m. A total of 4 fungal (Penicillium sp.) and 71 bacteria colonyforming units (70 colonies of Bacillus luciferensis believed to have originated from a single cell collected at altitude and one colony of Bacillus sphaericus) were enumerated, isolated and identified using a morphological key and 16S rDNA sequencing respectively. All of the isolates identified were sporeforming pigmented fungi or bacteria of terrestrial origin and demonstrate that the presence of viable microorganisms in Earth's upper atmosphere may not be uncommon.

Aerobiologia↗

USGS2009

EMAG2: A 2-arc min resolution Earth Magnetic Anomaly Grid compiled from satellite, airborne, and marine magnetic measurements

A global Earth Magnetic Anomaly Grid (EMAG2) has been compiled from satellite, ship, and airborne magnetic measurements. EMAG2 is a significant update of our previous candidate grid for the World Digital Magnetic Anomaly Map. The resolution has been improved from 3 arc min to 2 arc min, and the altitude has been reduced from 5 km to 4 km above the geoid. Additional grid and track line data have been included, both over land and the oceans. Wherever available, the original shipborne and airborne data were used instead of precompiled oceanic magnetic grids. Interpolation between sparse track lines in the oceans was improved by directional gridding and extrapolation, based on an oceanic crustal age model. The longest wavelengths (>330 km) were replaced with the latest CHAMP satellite magnetic field model MF6. EMAG2 is available at http://geomag.org/models/EMAG2 and for permanent archive at http://earthref.org/ cgi-bin/er.cgi?s=erda.cgi?n=970. ?? 2009 by the American Geophysical Union.

Geochemistry, Geophysics, Geosystems↗

USGS2009

The Mauna Loa carbon dioxide record: lessons for long-term Earth observations

The Mauna Loa carbon dioxide record is an iconic symbol of the human capacity to alter the planet. Yet this record would not have been possible without the remarkable work of one man, Charles David Keeling. We describe three emergent themes that characterized his work: (1) his desire to study and understand the processes that control atmospheric CO 2 and the global carbon cycle, (2) his campaign to identify and minimize systematic measurement error, and (3) his tenacious efforts to maintain continuous funding despite changing government priorities and institutions. In many ways, the story of the Mauna Loa record demonstrates that distinctions between research and “routine” measurements are not very useful in long-term monitoring of Earth properties and processes.

Book chapter↗

USGS2019

Upper mantle earth structure in Africa from full-wave ambient noise tomography

Our understanding of the tectonic development of the African continent and the interplay between its geological provinces is hindered by unevenly distributed seismic instrumentation. In order to better understand the continent, we used long-period ambient noise full waveform tomography on data collected from 186 broadband seismic stations throughout Africa and surrounding regions to better image the upper mantle structure. We extracted empirical Green’s functions from ambient seismic noise using a frequency-time normalization method and retrieved coherent signal at periods of 7-340 seconds. We simulated wave propagation through a heterogeneous Earth using a spherical finite-difference approach to obtain synthetic waveforms, measured the misfit as phase delay between the data and synthetics, calculated numerical sensitivity kernels using the scattering integral approach, and iteratively inverted for structure. The resulting images of isotropic, shear wavespeed for the continent reveal segmented, low-velocity upper mantle beneath the highly magmatic northern and eastern sections of the East African Rift System (EARS). In the southern and western sections, high-velocity upper mantle dominates, and distinct, low-velocity anomalies are restricted to regions of current volcanism. At deeper depths, the southern and western EARS transitions to low-velocities. In addition to the EARS, several low-velocity anomalies are scattered through the shallow upper mantle beneath Angola and North Africa, and some of these low-velocity anomalies may be connected to a deeper feature. Distinct upper mantle high-velocity anomalies are imaged throughout the continent and suggest multiple cratonic roots within the Congo region, and possible cratonic roots within the Sahara Metacraton.

Geochemistry, Geophysics, Geosystems↗

USGS2019

Hillslope hydrology in global change research and earth system modeling

Earth System Models (ESMs) are essential tools for understanding and predicting global change, but they cannot explicitly resolve hillslope‐scale terrain structures that fundamentally organize water, energy, and biogeochemical stores and fluxes at subgrid scales. Here we bring together hydrologists, Critical Zone scientists, and ESM developers, to explore how hillslope structures may modulate ESM grid‐level water, energy, and biogeochemical fluxes. In contrast to the one‐dimensional (1‐D), 2‐ to 3‐m deep, and free‐draining soil hydrology in most ESM land models, we hypothesize that 3‐D, lateral ridge‐to‐valley flow through shallow and deep paths and insolation contrasts between sunny and shady slopes are the top two globally quantifiable organizers of water and energy (and vegetation) within an ESM grid cell. We hypothesize that these two processes are likely to impact ESM predictions where (and when) water and/or energy are limiting. We further hypothesize that, if implemented in ESM land models, these processes will increase simulated continental water storage and residence time, buffering terrestrial ecosystems against seasonal and interannual droughts. We explore efficient ways to capture these mechanisms in ESMs and identify critical knowledge gaps preventing us from scaling up hillslope to global processes. One such gap is our extremely limited knowledge of the subsurface, where water is stored (supporting vegetation) and released to stream baseflow (supporting aquatic ecosystems). We conclude with a set of organizing hypotheses and a call for global syntheses activities and model experiments to assess the impact of hillslope hydrology on global change predictions.

Water Resources Research↗

USGS2021

The weight of cities: Urbanization effects on Earth’s subsurface

Across the world, people increasingly choose to live in cities. By 2050, 70% of Earth's population will live in large urban areas. Upon considering a large city, questions arise such as, how much does that weigh? What are its effects on the landscape? Does it cause measurable subsidence? Here I calculate the weight of San Francisco Bay region urbanization, where 7.75 million people live at, or near the coast. It is difficult to account for everything that is in a city. I assume that most of the weight is buildings and their contents, which allows the use of base outline and height data to approximate their mass, which is cumulatively 1.6·10 12 kg. I build a series of finite element models to study effects of pressure exerted by the weight distribution. Within the elastic realm, I look at compression, flexure, isostatic compensation, stress change, dilatation, and fluid flow changes. Within the nonlinear realm I show example calculations of primary and secondary settlement of soils under load. The combined modeled subsidence from building loads is at least 5–80 mm, with the largest contributions coming from nonlinear settlement and creep in soils. A general result is closing of pore space and redirection of pore fluids. While the calculated subsidence of the Bay Area is relatively small compared with other sources of elevation change such as pumping and recharge of aquifers, all sources of subsidence are concerning given an expected 200–300 mm sea level rise at San Francisco by the year 2050.

AGU Advances↗