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

Thomas E. Parsons

doi:10.1029/2020av000277

Home/Archive/USGS

USGS · 70218755

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

Thomas E. Parsons

View original source ↗Export metadata ↓

Abstract

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.

Explore related subjects

AGU Advances Article

Keep this discovery

Explore connections, maps & timelines

BibTeX RIS

Thomas E. Parsons. 2021-01-14. The weight of cities: Urbanization effects on Earth’s subsurface. https://doi.org/10.1029/2020av000277

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related USGS reports

USGS2025

A U.S. scientific community vision for sustained earth observations of greenhouse gases to support local to global action

Managing carbon stocks in the land, ocean, and atmosphere under changing climate requires a globally-integrated view of carbon cycle processes at local and regional scales. The growing Earth Observation (EO) record is the backbone of this multi-scale system, providing local information with discrete coverage from surface measurements and regional information at global scale from satellites. Carbon flux information, anchored by inverse estimates from spaceborne Greenhouse Gas (GHG) concentrations, provides an important top-down view of carbon emissions and sinks, but currently lacks global continuity at assessment and management scales (< 100 km). Partial-column data can help separate signals in the boundary layer from the overlying atmosphere, providing an opportunity to enhance surface sensitivity and bring flux resolution down from that of column-integrated data (100–500 km). Based on a workshop held in September 2024, the carbon cycle community envisions a carbon observation system leveraging GHG partial columns in the lower and upper troposphere to weave together information across scales from surface and satellite EO data, and integration of top- down/bottom-up analyses to link process understanding to global assessment.

AGU Advances↗

USGS2025

Toward co-designed Earth System Models: Reflecting end-user priorities in local applications from a modeler's perspective

Earth System Models (ESM) are crucial for quantifying climate impacts across Earth's interconnected systems and supporting science-based adaptation and mitigation. However, not including end-users, especially decision-makers representing communities vulnerable to climate change, can limit model utility, increase epistemic risks, and lead to information misuse in decision-making. While the ESM community increasingly values broad community engagement, end-users may not initially perceive models as useful for local planning. Co-designing models with end-users fosters two-way learning: users better understand models and their outputs, while modelers gain insights into fine-scale local processes like monitoring practices and management priorities. Higher-level co-design can lead to more customized, priority-driven, and useful modeling products. Despite these benefits, modelers often struggle to initiate meaningful partnerships with local communities. Therefore, this paper explores model co-design from the perspective of modelers. This study presents two case studies where modelers and social scientists collaborated with Indigenous communities' decision-makers to reflect their priorities in model design and application. In the Arctic Rivers Project, high-resolution climate and hydrology data sets for Alaska were developed with guidance from an Indigenous Advisory Council, using optimized, coupled land-atmosphere models. In the Mid-Klamath Project, we partnered with the Karuk Tribe's Department of Natural Resources to assess climate change and prescribed burning impacts on terrestrial hydrology in the Klamath River Basin. Drawing from these studies, we introduce a four-level framework: (a) Co-design Configuration; (b) Model Tuning; (c) Incorporate Contextual Knowledge; (d) Co-develop New Model Functions. We aim to help researchers consider and compare co-design across diverse modeling projects systematically and coherently.

Alaska, California, Oregon↗

USGS2025

Structural controls on splay fault rupture dynamics during Cascadia megathrust earthquakes

Great subduction earthquakes ( M w ≥ 8.0) can generate devastating tsunamis by rapidly displacing the seafloor and overlying water column. These potentially tsunamigenic seafloor offsets result from coseismic fault slip and deformation beneath or within the accretionary wedge. The mechanics of these shallow rupture phenomena and their dependence on subduction zone properties remain unresolved, partly due to the sparsity of offshore observations of shallow megathrust earthquake deformation. Here, we analyze how offshore structure influences shallow rupture mechanics and slip partitioning using 3D dynamic earthquake simulations of the Cascadia subduction zone (CSZ) megathrust with and without variably dipping seaward- or landward-vergent splay faults in the wedge that sole into the megathrust. Resulting tradeoffs between splay and megathrust slip reveal structural controls on rupture partitioning, with greater splay slip leading to less shallow megathrust slip updip. Gently dipping and seaward-vergent splays host more slip than those with steeper, landward-vergent splays. To isolate the underlying mechanisms, we compare models with Andersonian and plunging principal stresses. Results suggest distinct static and dynamic processes control the dip- and vergence-dependence of splay rupture: static (mis)alignment relative to far-field tectonic loading favors slip on more optimally oriented, shallowly dipping splay faults. In contrast, dynamic stress interactions of an updip-propagating megathrust rupture front with the free surface and potential branch faults favor forward branching onto seaward-vergent splays and inhibit backward branching onto landward-vergent splays. Resulting seafloor displacements suggest splay fault structure may influence coseismic tsunami source processes, highlighting the importance of dynamically viable rupture scenarios in subduction hazard assessments.

Cascadia subduction zone↗