Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Subsidence”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Increased flood exposure in the Pacific Northwest following earthquake-driven subsidence and sea-level rise

Climate-driven sea-level rise is increasing the frequency of coastal flooding worldwide, exacerbated locally by factors like land subsidence from groundwater and resource extraction. However, a process rarely considered in future sea-level rise scenarios is sudden (over minutes) land subsidence associated with great (>M8) earthquakes, which can exceed 1 m. Along the Washington, Oregon, and northern California coasts, the next great Cascadia subduction zone earthquake could cause up to 2 m of sudden coastal subsidence, dramatically raising sea level, expanding floodplains, and increasing the flood risk to local communities. Here, we quantify the potential expansion of the 1% floodplain (i.e., the area with an annual flood risk of 1%) under low (~0.5 m), medium (~1 m), and high (~2 m) earthquake-driven subsidence scenarios at 24 Cascadia estuaries. If a great earthquake occurred today, floodplains could expand by 90 km 2 (low), 160 km 2 (medium), or 300 km 2 (high subsidence), more than doubling the flooding exposure of residents, structures, and roads under the high subsidence scenario. By 2100, when climate-driven sea-level rise will compound the hazard, a great earthquake could expand floodplains by 170 km 2 (low), 240 km 2 (medium), or 370 km 2 (high subsidence), more than tripling the flooding exposure of residents, structures, and roads under the high subsidence scenario compared to the 2023 floodplain. Our findings can support decision-makers and coastal communities along the Cascadia subduction zone as they prepare for compound hazards from the earthquake cycle and climate-driven sea-level rise and provide critical insights for tectonically active coastlines globally.

Science & Technology - Other Topics↗
DOE OSTI2023

Drying of tundra landscapes will limit subsidence-induced acceleration of permafrost thaw

We used a model for permafrost hydrology informed by detailed measurements of soil ice content to better understand the potential risk of abrupt permafrost thaw triggered by melting ground ice, a key open question associated with permafrost response to a warming Arctic. Our spatially resolved simulations of a well-characterized site in polygonal tundra near Utqiaġvik, Alaska, agree well with multiple types of observations in the current climate. Projections indicate 63 cm of bulk subsidence from 2006 to 2100 in the strong-warming Representative Concentration Pathway 8.5 climate. Permafrost thaw as measured by the increase in active layer thickness (ALT)—the thickness of the soil layer that thaws each summer—is accelerated by subsidence, but the effect is relatively small. The ALT increases from the current-day value of approximately 50 cm to approximately 180 cm by 2100 when subsidence is included compared to about 160 cm when it is neglected. In these simulations, previously identified positive feedbacks between subsidence and thaw are self-limiting on decadal time frames because landscape runoff and increasing evapotranspiration result in drier tundra with weaker surface/atmosphere coupling. These results for a tundra site that is representative of large swathes of the Alaska North Slope suggest that subsidence is unlikely to lead to abrupt thaw over large areas. However, subsidence does have significant effects on the hydrology of polygonal tundra. Specifically, subsidence increases landscape runoff, which helps maintain streamflow in the face of increased evapotranspiration but also causes drier tundra conditions that could have deleterious effects on sensitive Arctic wetland ecosystems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Subsidence Measurements at the Teller 27 Field Site, Seward Peninsula, Alaska, 2021

In September 2021, differential GPS (dGPS) points were collected in areas of visual subsidence in the Teller watershed near mile marker 27 of the Bob Blodgett Nome-Teller Memorial Highway on the Seward Peninsula of Alaska. This study aimed to understand how structural permafrost loss is affecting meter-scale ground elevation. As ice thaws, heaving and slumping of the landscape creates new microtopographical features on the landscape that may alter surface hydrology and soil moisture and thus plant community composition and biogeochemical cycling. The Teller 27 field site is underlain with discontinuous permafrost, so the landscape features a variety of permafrost features. We collected ground elevations from eight subsidence areas in two forms: 1) Transects were sampled laterally through areas of ground subsidence to capture the high edges and low spots where the earth slumped. 2) The circumference of the subsidence area was sampled in order to measure the ground area affected by subsidence. These data were collected with an Emlid Reach RS2 dGPS and a base station. This dataset includes one *.csv of dGPS locations of the eight subsidence locations and one *.kml of measurement locations.The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic) was a 15-year research effort (2012-2027) to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy’s Office of Biological and Environmental Research.The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska.Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy’s Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Permafrost Thaw, Uneven Subsidence and Projected Drying of Ice-wedge Polygon Tundra: Modeling Archive

This dataset is a model archive of the paper Permafrost Thaw, Uneven Subsidence and Projected Drying of Ice-wedge Polygon Tundra (in prep) to support a modeling study investigating how projected increases in Arctic temperature and precipitation will jointly influence hydrologic conditions in ice-rich tundra landscapes. With this dataset, this study is to address the research question: Will Arctic tundra landscapes become wetter or drier with increasing precipitation and temperature in the future when thaw-induced ground subsidence and associated microtopographic evolution are represented? The simulations focus on ice-wedge polygon tundra, a widespread form of ice-rich permafrost terrain that is highly sensitive to thaw-driven landscape change. This dataset contains model input and output data for four study watersheds in Alaska: Anaktuvuk, Utqiagvik (formerly Barrow), Brooks Foothills, and Prudhoe Bay. Simulations were performed using the Advanced Terrestrial Simulator (ATS, v1.5), a physics-rich integrated surface–subsurface hydrologic model. For each watershed, ten modeling cases were performed representing two landscape evolution conditions (with subsidence and without subsidence) combined with five climate forcing scenarios derived from Shared Socioeconomic Pathways (SSP5, SSP5 with precipitation trend, SSP2, SSP2 with precipitation trend, and SSP2 with double precipitation trend). Particularly, for each watershed under the forcing SSP2 with precipitation trend, there are two additional simulations considering spatially heterogeneous subsidence distributions: one assumes randomly distributed scaling and the other includes elevation dependent distribution scaling. These simulations span 1980–2099 and include spin-up runs (1980–2009) followed by transient projections (2010–2099). To facilitate reproducibility of simulations, all datasets are organized by watershed. For each study watershed, the dataset contains: (1) Pre-partitioned mesh files for 32-core modeling (.par.32.XX), located in EACH_WATERSHED/mesh/basin; and also a non-partitioned mesh file (.exo) located in EACH_WATERSHED/mesh; (2) Climate forcings corresponding to the five SSP scenarios (.h5), located in EACH_WATERSHED/data; (3) Final states (.h5) from column spin-up modeling used to initialize historical watershed-scale spin-up runs from 1980 to 2009, located in EACH_WATERSHED/PreSpinupHistorical; (4) Final states (.h5) of historical watershed-scale spin-up runs from 1980 to 2009 used to initialize projection runs, located in EACH_WATERSHED/Spinup_daymetERA5; (5) ATS modeling input files (.xml), located in EACH_WATERSHED/EACH_SIMULATION_SCENARIO/inputfiles; (6) ATS modeling output files (.dat), located in in EACH_WATERSHED/EACH_SIMULATION_SCENARIO/combined_obs; (7) For the Brooks Foothills watershed, additional spatial model outputs are provided (.h5) for selected years (2033 and 2093) used to generate spatial figures in this study, located in Brooksfoothills/EACH_SIMULATION_SCENARIO/results-WITH/WITHOUT_SUBSIDENCE-year2033/2093. All data files with suffix .h5 can be accessible through Python h5py, and all data files with suffix of .dat can be imported by Python pandas. Mesh file with .exo can be visualized through Paraview or read by Python netCDF. The Next-Generation Ecosystem Experiments in the Arctic (NGEE Arctic) project is a research effort to reduce uncertainty in the Department of Energy’s Energy Exascale Earth System Model (E3SM) by developing a predictive understanding of Arctic tundra ecosystems underlain by permafrost and to quantify feedbacks from the Arctic tundra to the Earth system. NGEE Arctic is supported by the Department of Energy's Office of Biological and Environmental Research. Over Phases 1–3, observations made by the NGEE Arctic team across a gradient of permafrost landscapes in Arctic Alaska improved the representation of tundra processes in the land surface component of E3SM (the E3SM Land Model, ELM). Model improvements emphasized unique aspects of permafrost environments and explored reductions in model complexity while retaining predictive power. The Arctic-informed ELM developed by NGEE Arctic has been used to make novel predictions on processes ranging from permafrost thaw to soil biogeochemical cycling to Earth system feedbacks associated with the unique characteristics of tundra plants. In Phase 4, the NGEE Arctic team is evaluating our new predictive understanding under novel conditions across the Arctic domain. In collaboration with partners at long-term pan-Arctic research sites we are examining whether an Arctic-informed ELM can faithfully simulate interactions among surface and subsurface processes at site, regional, and pan-Arctic scales. In turn, we are using variety of tools to dynamically extend and evaluate ELM inference, with an emphasis on data synthesis and pan-Arctic model evaluation, reintegration of code with an evolving E3SM, scaling across heterogeneous Arctic landscapes, and the appropriate representation of the impacts of increasingly frequent Arctic disturbances.

EARTH SCIENCE > ATMOSPHERE > PRECIPITATION↗
DOE OSTI2025

Determination of Ground Subsidence Around Snow Fences in the Arctic Region

In this study, we analyzed the effects of snow cover changes caused by snow fences (SFs) installed in 2017 in the Alaskan tundra to examine ground subsidence. Digital surface model data obtained through LiDAR-based remote sensing in 2019 and 2022, combined with a field survey in 2021, revealed approximately 0.2 m of ground subsidence around the SF. To investigate the relationship between SF-induced snow cover changes and ground subsidence, geophysical methods, electrical resistivity tomography (ERT) and ground-penetrating radar (GPR), were applied in 2023 to analyze subsurface characteristics. The increased snow cover due to the SF-enhanced insulation, delaying the penetration of winter cold into the subsurface. This delay caused subsurface temperatures to decrease more slowly, melting the upper permafrost and increasing the thickness of the active layer. ERT and GPR surveys well delineated the boundary between the active layer and permafrost, confirming that the increased snow cover thickened the active layer. This thickening led to the melting of pore ice, causing water runoff and ground compaction, which resulted in subsidence. The runoff also formed channels flowing eastward over the SF. This study highlights how changes in snow cover can influence active layer properties, leading to localized environmental changes and ground subsidence.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Drying of tundra landscapes will limit subsidence-induced acceleration of permafrost thaw: Modeling Archive

This Modeling Archive is in support of a NGEE Arctic publication in review "Drying of tundra landscapes will limit subsidence-induced acceleration of permafrost thaw". The study used a cryohydrology model to assess the potential risk of abrupt permafrost thaw triggered by melting ground ice, a key open question associated with permafrost response to a warming Arctic. The spatially resolved simulations are for a small catchment 465 ice-wedge polygons in polygonal tundra near Utqiagvik, Alaska in the high-emissions RCP8.5 climate scenario. The simulations are compared to runoff, evapotranspiration and subsidence in the current climate and agree well. The study used the ATS code configured as an intermediate-scale cryohydrology model (Advanced Terrestrial Simulator, Version v2). The archive includes input files for spinup (1985 to 2005) and projections (2006 to 2100 - the manuscript reports 2006 to 2098). Three of the projections include the effects of subsidence and microtopography change with different depth profiles of ice content corresponding to the median, 20th percentile and 80th percentile. The fourth projection has subsidence turnoff and uses the reference case (median) ice content. Spatially averaged or aggregated output variables are included for the projections. Selected checkpoint files for the projections and postprocessing scripts are also included. Files included are *.xml; *.h5; *.exo; *.py; *.sh, *. nb, mesh files, and one file in the original Excel format plus *.csv and *.pdf to conserve formatting contained in the original file. The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research. The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska. Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Rising concerns of climate extremes and land subsidence impacts

A recent article in Reviews of Geophysics explores land subsidence drivers, rates, and impacts across the globe. It also discusses the need for improved process representations and the inclusion of the interplay among land subsidence and climatic extremes, including their effects in models and risk assessments. Here, we asked the lead author to explain the concept of land subsidence, its impacts, and future directions needed for improved mitigation.

Earth science↗
DOE OSTI2024

Sensitivity of modelled passive margin stratigraphy to variations in sea level, sediment supply and subsidence

Abstract We produced a 10 Myr synthetic stratigraphic section using a forward stratigraphic model that generates marine deltaic stratigraphy over geological timescales. We recursively fit the model using a Bayesian inversion algorithm to test: (1) if it could be accurately reconstructed; (2) if the parameters used to create it could be recovered; and (3) the sensitivity of the model output to given model parameters and the attendant physical processes. The original synthetic stratigraphic section was produced with cyclical sea‐level variations of 40 and 30 m with 2.4 and 10 Myr periods respectively. Sediment was also supplied cyclically, in 2.4 and 10 Myr cycles with amplitudes of 30 and 80 tons/100 kyr, respectively, varying from a mean of 232 tons/100 kyr. Parameter values were sampled to fit the model using a Markov chain Monte Carlo algorithm, resulting in a ±5 m (1σ) variation between the experimental output and the original. Sea level varied by ±7 m (1σ) within the posterior distribution of parameters. As a result, both the 10 Myr and 2.4 Myr sea‐level cycles could be extracted from the original output. The variation in sediment supply was approximately ±38 tons/100 kyr (1σ) and, as a result, only the larger long‐term supply variations could be accurately recovered in refitting the model. The variation in thermal, flexural and total subsidence across those parameter sets is less than ±10 m (1σ). The original section experienced 150 m of total subsidence at the depocentre. Our results demonstrate the distinct and interpretable imprint of sea level and subsidence on continental margin stratigraphy can be quantified. Moreover, we conclude that sea‐level change produces a defined effect on the geometries of stratigraphic architecture, and that techniques applied for the purpose of delineating sea‐level variation from continental margin strata have a well‐founded conceptual basis.

Schmelz, William J.↗
DOE OSTI2024

A Path to Clean Energy: Cross-Subsidization Concerns From Local Solar Development in Frankfort, Kentucky That Can Apply to Other Communities

Municipal electricity utilities and communities with public power are increasingly interested in the opportunities afforded by locally-sited solar projects. Understanding the allocation of costs and benefits of a solar project across affected entities and stakeholders such as customer types can be critical for community buy-in. This case study looked specifically at an analysis done for the City of Frankfort and how plans to meet 100% of city government electricity loads with renewable energy by 2023 could affect an all requirements contract the city's municipal utility has with a regional electricity provider. While the results are specific to Frankfort and the electricity contracts in place at the time of the study, it is an example of how the question of cross-subsidization can be addressed using a quantitative study, with the goal of increasing transparency and buy-in across multiple stakeholders.

community solar↗
DOE OSTI2023

Abrupt permafrost thaw drives spatially heterogeneous soil moisture and carbon dioxide fluxes in upland tundra

Permafrost thaw causes the seasonally thawed active layer to deepen, causing the Arctic to shift toward carbon release as soil organic matter becomes susceptible to decomposition. Ground subsidence initiated by ice loss can cause these soils to collapse abruptly, rapidly shifting soil moisture as microtopography changes and also accelerating carbon and nutrient mobilization. The uncertainty of soil moisture trajectories during thaw makes it difficult to predict the role of abrupt thaw in suppressing or exacerbating carbon losses. Here, in this study, we investigated the role of shifting soil moisture conditions on carbon dioxide fluxes during a 13-year permafrost warming experiment that exhibited abrupt thaw. Warming deepened the active layer differentially across treatments, leading to variable rates of subsidence and formation of thermokarst depressions. In turn, differential subsidence caused a gradient of moisture conditions, with some plots becoming consistently inundated with water within thermokarst depressions and others exhibiting generally dry, but more variable soil moisture conditions outside of thermokarst depressions. Experimentally induced permafrost thaw initially drove increasing rates of growing season gross primary productivity (GPP), ecosystem respiration (R eco ), and net ecosystem exchange (NEE) (higher carbon uptake), but the formation of thermokarst depressions began to reverse this trend with a high level of spatial heterogeneity. Plots that subsided at the slowest rate stayed relatively dry and supported higher CO 2 fluxes throughout the 13-year experiment, while plots that subsided very rapidly into the center of a thermokarst feature became consistently wet and experienced a rapid decline in growing season GPP, R eco , and NEE (lower carbon uptake or carbon release). These findings indicate that Earth system models, which do not simulate subsidence and often predict drier active layer conditions, likely overestimate net growing season carbon uptake in abruptly thawing landscapes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Shallow-Water Habitat in the Lower Columbia River Estuary: A Highly Altered System

Abstract Decreases in shallow-water habitat area (SWHA) in the Lower Columbia River and Estuary (LCRE) have adversely affected salmonid populations. We investigate the causes by hindcasting SWHA from 1928 to 2004, system-wide, based on daily higher high water (HHW) and system hypsometry. Physics-based regression models are used to represent HHW along the system as a function of river inflow, tides, and coastal processes, and hypsometry is used to estimate the associated SWHA. Scenario modeling is employed to attribute SWHA losses to levees, flow regulation, diversion, navigational development, and climate-induced hydrologic change, for subsidence scenarios of up to 2 m, and for 0.5 m fill. For zero subsidence, the system-wide annual-average loss of SWHA is 55 ± 5%, or 51 × 10 5 ha/year; levees have caused the largest decrease ( $${54}_{-14}^{+5}$$ 54 - 14 + 5 %, or ~ 50 × 10 5 ha/year). The loss in SWHA due to operation of the hydropower system is small, but spatially and seasonally variable. During the spring freshet critical to juvenile salmonids, the total SWHA loss was $${63}_{-3}^{+2}$$ 63 - 3 + 2 %, with the hydropower system causing losses of 5–16% (depending on subsidence). Climate change and navigation have caused SWHA losses of $${5}_{-5}^{+16}$$ 5 - 5 + 16 % and $${4}_{-6}^{+14}$$ 4 - 6 + 14 %, respectively, but with high spatial variability; irrigation impacts have been small. Uncertain subsidence causes most of the uncertainty in estimates; the sum of the individual factors exceeds the total loss, because factors interact. Any factor that reduces mean or peak flows (reservoirs, diversion, and climate change) or alters tides and along-channel slope (navigation) becomes more impactful as assumed historical elevations are increased to account for subsidence, while levees matter less.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Impacts of a Cascadia Subduction Zone Earthquake on Water Levels and Wetlands of the Lower Columbia River and Estuary

Abstract Subsidence after a subduction zone earthquake can cause major changes in estuarine bathymetry. Here, we quantify the impacts of earthquake‐induced subsidence on hydrodynamics and habitat distributions in a major system, the lower Columbia River Estuary, using a hydrodynamic and habitat model. Model results indicate that coseismic subsidence increases tidal range, with the smallest changes at the coast and a maximum increase of ∼10% in a region of topographic convergence. All modeled scenarios reduce intertidal habitat by 24%–25% and shifts ∼93% of estuarine wetlands to lower‐elevation habitat bands. Incorporating dynamic effects of tidal change from subsidence yields higher estimates of remaining habitat by multiples of 0–3.7, dependent on the habitat type. The persistent tidal change and chronic habitat disturbance after an earthquake poses strong challenges for estuarine management and wetland restoration planning, particularly when coupled with future sea‐level rise effects.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Model spread in tropical low cloud feedback tied to overturning circulation response to warming

Among models participating in the Coupled Model Intercomparison Project phase 6 (CMIP6), here we show that the magnitude of the tropical low cloud feedback, which contributes considerably to uncertainty in estimates of climate sensitivity, is intimately linked to tropical deep convection and its effects on the tropical atmospheric overturning circulation. First, a reduction in tropical ascent area and an increased frequency of heavy precipitation result in high cloud reduction and upper-tropospheric drying, which increases longwave cooling and reduces subsidence weakening, favoring low cloud reduction (Radiation-Subsidence Pathway). Second, increased longwave cooling decreases tropospheric stability, which also reduces subsidence weakening and low cloudiness (Stability-Subsidence Pathway). In summary, greater high cloud reduction and upper-tropospheric drying (negative longwave feedback) lead to a more positive cloud feedback among CMIP6 models by contributing to a greater reduction in low cloudiness (positive shortwave feedback). Varying strengths of the two pathways contribute considerably to the intermodel spread in climate sensitivity.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Surface Deformation and Seismicity Linked to Fluid Injection in the Raton Basin

Abstract It is suggested that in addition to seismicity deep fluid injection may cause surface uplift and subsidence in oil and gas‐producing regions. This study uses the Raton Basin as an example to investigate the hydromechanical processes of surface uplift and subsidence during wastewater injection. The Raton Basin, in southern central Colorado and northern central New Mexico, has experienced wastewater injection related to coalbed methane and gas production starting in 1994. In this study, we estimate the extent and magnitude of total vertical deformation in the Raton Basin from 1994 to 2020 and incremental deformation between the years 2017 to 2020. Results indicate a modeled uplift as much as 15 cm occurring between 1994 and 2020. Between 2017 and 2020, up to 3 cm of uplift occurred, largely near the northwestern injection wells. Most modeled uplift between 1994 and 2020 occurred near the southern wells, where the greatest cumulative volume of wastewater was injected. However, modeled subsidence occurred around the southern and eastern wells between 2017 and 2020, after the rate of injection decreased. Modeling indicates that while the magnitude of modeled uplift corresponds to cumulative injection volume and maximum rate in the long‐term, short‐term incremental deformation (uplift or subsidence) is controlled by changes in the rate of injection. The number of yearly earthquake events follows periods of rapid modeled uplifting throughout the Basin, suggesting that measurable surface deformation may be caused by the same injection‐induced pore pressure perturbations that initiate seismicity.

Geology↗
DOE OSTI2023

Groundwater Volume Loss in Mexico City Constrained by InSAR and GRACE Observations and Mechanical Models

Abstract Groundwater withdrawal can cause localized and rapid poroelastic subsidence, spatially broad elastic uplift of low amplitude, and changes in the gravity field. Constraining groundwater loss in Mexico City, we analyze data from the Gravity Recovery and Climate Experiment and its follow‐on mission (GRACE/FO) and Synthetic Aperture Radar (SAR) Sentinel‐1A/B images between 2014 and 2021. GRACE/FO observations yield a groundwater loss of 0.85–3.87 km 3 /yr for a region of ∼300 × 600 km surrounding Mexico City. Using the high‐resolution interferometric SAR data set, we measure >35 cm/yr subsidence within the city and up to 2 cm/yr of uplift in nearby areas. Attributing the long‐term subsidence to poroelastic aquifer compaction and the long‐term uplift to elastic unloading, we apply respective models informed by local geology, yielding groundwater loss of 0.86–12.57 km 3 /yr. Our results suggest Mexico City aquifers have been depleting at faster rates since 2015, exacerbating the socioeconomic and health impacts of long‐term groundwater overdrafts.

58 GEOSCIENCES↗
DOE OSTI2022

An Energetic Diagnostic of Tropical Cyclone Size in f –Plane Simulations

As a major feature of tropical cyclones (TCs), controlling factors of TC outer size or size scaling remains a fundamental scientific question. The Rossby deformation radius and a natural extent associated with potential intensity have been proposed as two scalings of TC size. But neither of them satisfactorily captures the sensitivity of TC size to sea surface temperature (SST) in idealized f-plane simulations. Inspired by the studies of the Hadley circulation, here we proposed a new TC scaling based on an energetic diagnostic scaling. TC size is primarily a ratio of the secondary circulation strength to subsidence velocity, further determined by the total atmospheric heating in the ascending area, the gross moist stability, the diabatic cooling, and the dry static stability. The former two is based on the moist energetic budget applied to the whole storm structure, while the latter two is based on the dry thermodynamic budget applied to the subsidence areas. The new scaling well captured the sensitivity of TC size to SST in idealized f-plane simulations, partly resulted from expanded ascending area, increased surface moisture deficit, and weakened subsidence with increased SST.

58 GEOSCIENCES↗
DOE OSTI2022

Cloud adjustments from large-scale smoke–circulation interactions strongly modulate the southeastern Atlantic stratocumulus-to-cumulus transition

Smoke from southern Africa blankets the southeastern Atlantic Ocean from June to October, producing strong and competing aerosol radiative effects. Smoke effects on the transition between overcast stratocumulus and scattered cumulus clouds are investigated along a Lagrangian (air-mass-following) trajectory in regional climate and large eddy simulation models. Results are compared with observations from three recent field campaigns that took place in August 2017: ObseRvations of Aerosols above CLouds and their intEractionS (ORACLES), CLouds and Aerosol Radiative Impacts and Forcing: Year 2017 (CLARIFY), and Layered Atlantic Smoke Interactions with Clouds (LASIC). The case study is set up around the joint ORACLES–CLARIFY flight that took place near Ascension Island on 18 August 2017. Smoke sampled upstream on an ORACLES flight on 15 August 2017 likely entrained into the marine boundary layer later sampled during the joint flight. The case is first simulated with the WRF-CAM5 regional climate model in three distinct setups: (1) FireOn, in which smoke emissions and any resulting smoke–cloud–radiation interactions are included; (2) FireOff, in which no smoke emissions are included; (3) RadOff, in which smoke emissions and their microphysical effects are included but aerosol does not interact directly with radiation. Over the course of the Lagrangian trajectory, differences in free tropospheric thermodynamic properties between FireOn and FireOff are nearly identical to those between FireOn and RadOff, showing that aerosol–radiation interactions are primarily responsible for the free tropospheric effects. These effects are non-intuitive: in addition to the expected heating within the core of the smoke plume, there is also a “banding” effect of cooler temperature (~1–2 K) and greatly enhanced moisture (>2 g kg -1 ) at the plume top. This banding effect is caused by a vertical displacement of the former continental boundary layer in the free troposphere in the FireOn simulation resulting from anomalous diabatic heating due to smoke absorption of sunlight that manifests primarily as a few hundred meters per day reduction in large-scale subsidence over the ocean. A large eddy simulation (LES) is then forced with free tropospheric fields taken from the outputs for the WRF-CAM5 FireOn and FireOff runs. Cases are run by selectively perturbing one variable (e.g., aerosol number concentration, temperature, moisture, vertical velocity) at a time to better understand the contributions from different indirect (microphysical), “large-scale” semi-direct (above-cloud thermodynamic and subsidence changes), and “local” semi-direct (below-cloud smoke absorption) effects. Despite a more than 5-fold increase in cloud droplet number concentration when including smoke aerosol concentrations, minimal differences in cloud fraction evolution are simulated by the LES when comparing the base case with a perturbed aerosol case with identical thermodynamic and dynamic forcings. A factor of 2 decrease in background free tropospheric aerosol concentrations from the FireOff simulation shifts the cloud evolution from a classical entrainment-driven “deepening–warming” transition to trade cumulus to a precipitation-driven “drizzle-depletion” transition to open cells, however. The thermodynamic and dynamic changes caused by the WRF-simulated large-scale adjustments to smoke diabatic heating strongly influence cloud evolution in terms of both the rate of deepening (especially for changes in the inversion temperature jump and in subsidence) and in cloud fraction on the final day of the simulation (especially for the moisture “banding” effect). Such large-scale semi-direct effects would not have been possible to simulate using a small-domain LES model alone.

54 ENVIRONMENTAL SCIENCES↗