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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Carbon Thaw Rate Doubles When Accounting for Subsidence in a Permafrost Warming Experiment

Permafrost soils, which are perennially frozen soils found throughout cold regions, contain vast quantities of carbon and ice. When permafrost thaws, carbon can be lost to the atmosphere, contributing to climate change. This means it is important to track permafrost thaw, which is often done using active layer thickness, or the depth of the seasonally thawed surface layer of soil. However, ice volume can be lost from thawing permafrost, causing the soil surface to drop. Conventional measurements do not account for this surface drop, and the rate of thaw could therefore be underestimated. We found that experimentally warmed soils dropped at a rate of 6 cm year-1, mostly due to loss of ice volume and also due to the loss of soil mass. When accounting for the change in soil surface height over time, the full depth of permafrost thaw was 49% greater. The increased depth of thaw resulted in more than twice as much carbon being thawed as was estimated with standard methods that did not account for subsidence. These findings suggest that permafrost is thawing more quickly than long-term records indicate and that this could result in additional carbon release contributing to climate change.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Pan-Arctic Permafrost Subsidence [Slides]

The Arctic is warming twice as fast as the rest of the globe, resulting in dramatic loss of ice across the Arctic system. Models project 7% to 90% of permafrost could be lost by the end of the century. Permafrost thaw is driving landscape degradation and infrastructure disruption through thermokarst and thermal erosion. Previous pan-arctic modeling shows potential for significant subsidence (Lee et al. 2014). CLM4.5 pan-Arctic simulations from PCNMIP were used to calculate change in projected ALT from 2010-2030.

54 ENVIRONMENTAL SCIENCES↗
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 OSTI2021

Displacement Analysis of Geothermal Field Based on PSInSAR And SOM Clustering Algorithms A Case Study of Brady Field, Nevada—USA

The availability of free and high temporal resolution satellite data and advanced SAR techniques allows us to analyze ground displacement cost-effectively. Our aim was to properly define subsidence and uplift areas to delineate a geothermal field and perform time-series analysis to identify temporal trends. A Persistent Scatterer Interferometry (PSI) algorithm was used to estimate vertical displacement in the Brady geothermal field located in Nevada by analyzing 70 Sentinel-1A Synthetic-Aperture Radar (SAR) images, between January 2017 and December 2019. To classify zones affected by displacement, an unsupervised Self-Organizing Map (SOM) algorithm was applied to classify points based on their behavior in time, and those clusters were used to determine subsidence, uplift, and stable regions automatically. Finally, time-series analysis was applied to the clustered data to understand the inflection dates. The maximum subsidence is –19 mm/yr with an average value of –6 mm/yr within the geothermal field. The maximum uplift is 14 mm/yr with an average value of 4 mm/yr within the geothermal field. The uplift occurred on the NE of the field, where the injection wells are located. On the other hand, subsidence is concentrated on the SW of the field where the production wells are located. The coupling of the PSInSAR and the SOM algorithms was shown to be effective in analyzing the direction and pattern of the displacements observed in the field.

54 ENVIRONMENTAL SCIENCES↗
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 OSTI2020

Training dataset and results for geothermal exploration artificial intelligence, applied to Brady Hot Springs and Desert Peak

The submission includes the labeled datasets, as ESRI Grid files (.gri, .grd) used for training and classification results for our machine leaning model: - brady_som_output.gri, brady_som_output.grd, brady_som_output.* - desert_som_output.gri, desert_som_output.grd, desert_som_output.* The data corresponds to two sites: Brady Hot Springs and Desert Peak, both located near Fallon, NV. Input layers include: - Geothermal: Labeled data (0: Non-geothermal; 1: Geothermal) - Minerals: Hydrothermal mineral alterations, as a result of spectral analysis using Chalcedony, Kaolinite, Gypsum, Hematite and Epsomite - Temperature: Land surface temperature (% of times a pixel was classified as "Hot" by K-Means) - Faults: Fault density with a 300mradius - Subsidence: PSInSAR results showing subsidence displacement of more than 5mm - Uplift: PSInSAR results showing subsidence displacement of more than 5mm Also, the results of the classification using Brady and Desert Peak to build 2 Convolutional Neural Networks. These were applied to the training site as well as the other site, the results are in GeoTiff format. - brady_classification: Results of classification of the Brady-trained model - desert_classification: Results of classification of the Desert Peak-trained model - b2d_classification: Results of classification of Desert Peak using the Brady-trained model - d2b_classification: Results of classification of Brady using the Desert Peak-trained model

15 GEOTHERMAL ENERGY↗
DOE OSTI2020

Time–series Analysis of Volume Change at Brady Hot Springs, Nevada, USA using Geodetic Data from 2003 – 2018

Brady Hot Springs geothermal field has exhibited subsidence, as measured by interferometric synthetic aperture radar (InSAR). Previous studies have examined both the temporal evolution of the deformation from 2004 through 2016 and the spatial extent of the deformation, directly relating the observed subsidence to volumetric changes below the surface. We extend the modeling at Brady to analyze a data set of interferometric pairs spanning from the end of 2003 through 2018. We examine spatial and temporal trends in the observed deformation by time–series analysis of each of the 1656 cubic voxels in a parameterized elastic dislocation model to identify areas where the subsurface volume changes as a function of time. Joint time–series analysis of Global Positioning System and InSAR pairs confirm significant changes in rates of volume change during time intervals when well operations were varied. Here, the rate of subsidence increases with increased injection, consistent with the identification of thermal contraction of the rock matrix as the dominant driving mechanism. Conversely, the modeled volume increases when pumping ceases, suggesting thermal expansion of the rock matrix.

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↗
DOE OSTI2023

Predicting Short-Term Deformation in the Central Valley Using Machine Learning

Land subsidence caused by excessive groundwater pumping in Central Valley, California, is a major issue that has several negative impacts such as reduced aquifer storage and damaged infrastructures which, in turn, produce an economic loss due to the high reliance on crop production. This is why it is of utmost importance to routinely monitor and assess the surface deformation occurring. Two main goals that this paper attempts to accomplish are deformation characterization and deformation prediction. The first goal is realized through the use of Principal Component Analysis (PCA) applied to a series of Interferomtric Synthetic Aperture Radar (InSAR) images that produces eigenimages displaying the key characteristics of the subsidence. Water storage changes are also directly analyzed by the use of data from the Gravity Recovery and Climate Experiment (GRACE) twin satellites and the Global Land Data Assimilation System (GLDAS). The second goal is accomplished by building a Long Short-Term Memory (LSTM) model to predict short-term deformation after developing an InSAR time series using LiCSBAS, an open-source InSAR time series package. The model is applied to the city of Madera and produces better results than a baseline averaging model and a one dimensional convolutional neural network (CNN) based on a mean squared error metric showing the effectiveness of machine learning in deformation prediction as well as the potential for incorporation in hazard mitigation models. The model results can directly aid policy makers in determining the appropriate rate of groundwater withdrawal while maintaining the safety and well-being of the population as well as the aquifers’ integrity.

58 GEOSCIENCES↗
DOE OSTI2022

Analysis of short‐ and long‐term system response during gas production from a gas hydrate deposit at the UBGH2‐6 site of the Ulleung Basin in the Korean East Sea

Abstract This study is a continuation of an investigation into the feasibility of long‐term production from a marine hydrate accumulation that has the properties and conditions of the UBGH2‐6 (UBGH2, Ulleung Basin Gas Hydrates 2; 2 is the number of the scientific expedition) site at the Ulleung Basin in the Korean East Sea. The 12.7 m‐thick system is in deep water (2157 m), but at 140 m below the seafloor. It is characterized by alternating hydrate‐free clays and muds and hydrate‐rich sand layers. The layered stratigraphy and the presence of mud layers favours the use of vertical wells rather than horizontal wells for production. The analysis indicates that production from such a hydrate accumulation is technically feasible, but the gas production rates are generally low. Water production accompanying gas production from this deposit appears manageable under all the scenarios investigated in this study, however, the water‐to‐gas ratio is high. Subsidence at the ocean floor at the end of a 14 day test is quite limited. However, there is significant uncertainty in the predictions of the geomechanical system's behaviour because they are not based on measured system properties but only on estimates/assumptions from analogues. The long‐term production potential of the reservoir at the site investigated here appears challenging because of the limited effectiveness of dissociation and large water production, in addition to substantial subsidence.

Moridis, George J.↗
DOE OSTI2021

Analysis of Coastal Fog from a Ship During the C-FOG Campaign

This work presents ship-based measurements of fog off St John’s Newfoundland on 13 September 2018 during the Coastal Fog (C-FOG) field campaign. The measurements included cloud-particle spectra, cloud-base height and aerosol backscatter, radiation, turbulence, visibility, and sea surface temperature. Radiosonde soundings were made at intervals less than 2 hours. Fog occurred in two episodes under the center of eastward-moving synoptic low-pressure system. The boundary layer structure during the first fog episode consisted of 3 layers, separated by 2 saturated air temperature inversions and capped by a subsidence inversion. The lowest layer was fog and the upper layers were cloud. The second fog episode consisted of one well-mixed fog layer capped by a subsidence inversion. Low wind speeds and stable stratification maintained low surface-layer turbulence during fog. Droplet size distributions had typical bimodal distributions. The visibility correlated with the droplet number concentration and liquid water content. The air temperature was higher than the sea surface temperature for the first 30 minutes of the first fog episode but was colder than the sea for the remainder of all fog. The sensible heat flux was upward, from sea to air, for the first 62 % of the first fog episode and then reversed to downward, from air to sea, for the remainder of the first fog episode and the second fog episode. Here, the counter-gradient heat fluxes observed (i.e., opposite to what is expected from the instantaneous air-sea temperature difference) appear to be related to turbulence, entrainment and stratification in the fog layer that overwhelmed the influence of air-sea temperature difference. While the synoptic-scale dynamics preconditioned the area for fog formation, the final step of fog appearance in this case was nuanced by stratification-turbulence interactions, local advective processes and microphysical environment.

54 ENVIRONMENTAL SCIENCES↗