Results for “flux vertical gradients” — Page 4

NASA NTRS1992

The coupling of the troposphere and stratosphere

A linear, time-dependent, primitive-equation, 3D numerical model is used to investigate the dynamical coupling between the troposphere and stratosphere through the vertical propagation of planetry waves. Attention is given to the roles of the tropopause in affecting the propagation of planetary waves and to the extent to which the stratosphere can determine its own wave driving, i.e., the sensitivity of the generation of the wave activity in the troposphere to the change of zonal wind in the stratosphere. Two maxima in Eliassen-Palm (EP) flux convergence are found for wave 1 in high latitudes, one in the upper stratosphere and the other in the upper troposphere. For wave 2 there is a third maximum north of the zero wind line in the tropical stratosphere. The formation of the maximum in EP-flux convergence in the upper troposphere is related to the dissipation of non-WKBJ waves near a local maximum of the refractive index. The existence of the local maximum in refractive index at the tropopause is a result of the great gradient of the buoyancy frequency as well as the vertical shear of the zonal flow there.

Chen, Ping↗

DOE OSTI2025

Topography-Induced TKE Budget Behavior Over an Amazon Forest

Data from three different heights (35, 50 and 81 m) of one of the ATTO Project towers in the Amazon forest were used to calculate the TKE (turbulence kinetic energy) budget and some other statistics within the RSL (roughness sublayer). The statistical analyses were carried out for unstable and stable cases. The vertical transport and the vertical advection terms do not explain the imbalances found in the TKE budget, highlighting the impor- tance of horizontal terms in complex terrain. Indeed, patterns independent of the time of day were found for vertical transport, horizontal and vertical velocity skewness, the mean vertical component velocity, and momentum and vertical TKE flux. Here, the highest values of mean vertical velocity were observed from the direction where there is a valley near the tower, which may indicate a topographical effect. Monin-Obukhov Similarity Theory (MOST) predictions for the mean velocity gradient, strictly applicable only in the inertial sublayer (ISL), gave reasonable predictions at 35 and 50 m, but dimensionless dissipation rates displayed large scatter. At 81 m, there were clear deviations from MOST for the two functions, disclosing the effect of topography. In addition, a neutral LES (large-eddy simu- lation) study showed good agreement with tower data for the same general wind direction.

ATTO project↗

DOE OSTI2023

Estimating the Meridional Extent of Adiabatic Mixing in the Stratosphere Using Age-Of-Air

Wave-induced adiabatic mixing in the winter midlatitudes is one of the key processes impacting stratospheric transport. Understanding its strength and structure is vital to understanding the distribution of trace gases and their modulation under a changing climate. Age-of-air is often used to understand stratospheric transport, and this study proposes refinements to the vertical age gradient theory of Linz et al. (2021), https://doi.org/10.1029/2021JD035199. The theory assumes exchange of air between a well-mixed tropics and a well-mixed extratropics, separated by a transport barrier, quantifying the adiabatic mixing flux across the interface using age-based measures. These assumptions are re-evaluated and a refined framework that includes the effects of meridional tracer gradients is established to quantify the mixing flux. This is achieved, in part, by computing a circulation streamfunction in age-potential temperature coordinates to generate a complete distribution of parcel ages being mixed in the midlatitudes. The streamfunction quantifies the “true” age of parcels mixed between the tropics and the extratropics. Applying the revised theory to an idealized and a comprehensive climate model reveals that ignoring the meridional gradients in age leads to an underestimation of the wave-driven mixing flux. Stronger, and qualitatively similar fluxes are obtained in both models, especially in the lower-to-middle stratosphere. While the meridional span of adiabatic mixing in the two models exhibits some differences, they show that the deep tropical pipe, that is, latitudes equatorward of 15° barely mix with older midlatitude air. The novel age-potential temperature circulation can be used to quantify additional aspects of stratospheric transport.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2026

Hydrogeological assessment of CO2 containment assurance and wellbore integrity at a Gulf Coast storage site

Abstract A large-scale carbon capture and storage (CCS) initiative on the Texas Gulf Coast serves as a premier demonstration of the U.S. Department of Energy’s CarbonSAFE program. Targeting deep saline formations, specifically Oligo-Miocene deltaic sequences, the project aims to establish technical and commercial viability for geologic CO2 storage within a major industrial corridor. This study provides a rigorous hydrogeological assessment to support Class VI permitting by quantifying the high degree of containment security. Utilizing a compositional reservoir simulator, we developed a suite of 27 distinct simulation cases to evaluate vertical plume dynamics near both planned injection wells and proximal legacy infrastructure. To ensure numerical accuracy near wellbores, we implemented a refined mesh strategy, determining that a 5.6 ft × 5.6 ft grid refinement offered the optimal balance between computational efficiency and descriptive precision. The modeling framework utilized a systematic sensitivity-based approach to evaluate the mechanical redundancy of the subsurface system by performing a bounding analysis of wellbore interfaces against hypothetical high-permeability microannuli. By systematically isolating competing physical drivers, including permeability, porosity, gas hysteresis, thermal gradients, salinity, and solubility trapping (quantified via Henry’s law with dynamically adjusted coefficients), this work moves beyond binary assessments to establish a nuanced hierarchy of containment factors. The results confirm that primary trapping mechanisms (e.g., gas hysteresis and solubility), combined with the site's unique geomechanical stratigraphy, significantly restrict vertical mobility and reinforce the robust containment security of the reservoir. Baseline results demonstrate substantial vertical separation between the CO2 plume and the upper confining system, ensuring robust containment. Sensitivity analysis reveals that even under highly conservative bounding scenarios—assuming theoretical 10-Darcy pathways at specific wellbore locations—the 2,900-ft thick multi-layered confining zone remains a reliable barrier. In these hypothetical upper-bound cases, peak upward fluxes of CO2 and saltwater after 15 years of injection remain localized and dissipate rapidly within the lower sections of the confining interval, leaving the integrity of the seal uncompromised. Furthermore, the study identifies that while localized wellbore pathways define theoretical upper bounds of vertical migration, the Area of Review (AoR) is primarily sensitive to regional thermal gradients and hysteresis, which can influence the AoR by over 3,000 acres in pessimistic configurations. Also, primary trapping mechanisms, specifically gas hysteresis and solubility, work in tandem with the Gulf Coast’s unique geomechanical stratigraphy to significantly restrict vertical mobility. Ductile, smectite-rich mudstones facilitate natural borehole convergence and the self-healing of potential conduits, creating a natural geomechanical bridge that effectively mitigates migration potential at both current injection points and legacy-well locations. This comprehensive modeling effort demonstrates that the integration of high-resolution wellbore simulations and regional geomechanical observations confirms the long-term storage security of the studied site, providing a physics-based foundation for industrial-scale CCS deployments. This modeling framework establishes a baseline for future research into coupled geomechanical effects, such as time-dependent borehole convergence, to further refine long-term containment projections. Acknowledgements We thank the Gulf Coast Carbon Center (GCCC) at the Bureau of Economic Geology for foundational research support. We appreciate Alex Bump for technical guidance and David Hoffman for model mesh generation. This work used TACC’s Frontera cluster for simulations and CMG Ltd. software licenses provided to UT-Austin. This material is based upon work supported by the Department of Energy under Award Number DE-FE0032338. Disclaimer This material is based upon work supported by the U.S. Department of Energy’s Fossil Energy and Carbon Management Office under the CarbonSAFE program, award Number DE-FE0032338. The views expressed herein do not necessarily represent the views of the U.S. Department of Energy or the United States Government.

58 GEOSCIENCES↗

DOE OSTI2026

Hydrogeological assessment of CO2 containment assurance and wellbore integrity at a Gulf Coast storage site

Abstract A large-scale carbon capture and storage (CCS) initiative on the Texas Gulf Coast serves as a premier demonstration of the U.S. Department of Energy’s CarbonSAFE program. Targeting deep saline formations, specifically Oligo-Miocene deltaic sequences, the project aims to establish technical and commercial viability for geologic CO2 storage within a major industrial corridor. This study provides a rigorous hydrogeological assessment to support Class VI permitting by quantifying the high degree of containment security. Utilizing a compositional reservoir simulator, we developed a suite of 27 distinct simulation cases to evaluate vertical plume dynamics near both planned injection wells and proximal legacy infrastructure. To ensure numerical accuracy near wellbores, we implemented a refined mesh strategy, determining that a 5.6 ft × 5.6 ft grid refinement offered the optimal balance between computational efficiency and descriptive precision. The modeling framework utilized a systematic sensitivity-based approach to evaluate the mechanical redundancy of the subsurface system by performing a bounding analysis of wellbore interfaces against hypothetical high-permeability microannuli. By systematically isolating competing physical drivers, including permeability, porosity, gas hysteresis, thermal gradients, salinity, and solubility trapping (quantified via Henry’s law with dynamically adjusted coefficients), this work moves beyond binary assessments to establish a nuanced hierarchy of containment factors. The results confirm that primary trapping mechanisms (e.g., gas hysteresis and solubility), combined with the site's unique geomechanical stratigraphy, significantly restrict vertical mobility and reinforce the robust containment security of the reservoir. Baseline results demonstrate substantial vertical separation between the CO2 plume and the upper confining system, ensuring robust containment. Sensitivity analysis reveals that even under highly conservative bounding scenarios—assuming theoretical 10-Darcy pathways at specific wellbore locations—the 2,900-ft thick multi-layered confining zone remains a reliable barrier. In these hypothetical upper-bound cases, peak upward fluxes of CO2 and saltwater after 15 years of injection remain localized and dissipate rapidly within the lower sections of the confining interval, leaving the integrity of the seal uncompromised. Furthermore, the study identifies that while localized wellbore pathways define theoretical upper bounds of vertical migration, the Area of Review (AoR) is primarily sensitive to regional thermal gradients and hysteresis, which can influence the AoR by over 3,000 acres in pessimistic configurations. Also, primary trapping mechanisms, specifically gas hysteresis and solubility, work in tandem with the Gulf Coast’s unique geomechanical stratigraphy to significantly restrict vertical mobility. Ductile, smectite-rich mudstones facilitate natural borehole convergence and the self-healing of potential conduits, creating a natural geomechanical bridge that effectively mitigates migration potential at both current injection points and legacy-well locations. This comprehensive modeling effort demonstrates that the integration of high-resolution wellbore simulations and regional geomechanical observations confirms the long-term storage security of the studied site, providing a physics-based foundation for industrial-scale CCS deployments. This modeling framework establishes a baseline for future research into coupled geomechanical effects, such as time-dependent borehole convergence, to further refine long-term containment projections. Acknowledgements We thank the Gulf Coast Carbon Center (GCCC) at the Bureau of Economic Geology for foundational research support. We appreciate Alex Bump for technical guidance and David Hoffman for model mesh generation. This work used TACC’s Frontera cluster for simulations and CMG Ltd. software licenses provided to UT-Austin. This material is based upon work supported by the Department of Energy under Award Number DE-FE0032338. Disclaimer This material is based upon work supported by the U.S. Department of Energy’s Fossil Energy and Carbon Management Office under the CarbonSAFE program, award Number DE-FE0032338. The views expressed herein do not necessarily represent the views of the U.S. Department of Energy or the United States Government.

58 GEOSCIENCES↗

NASA NTRS1978

Baroclinic adjustment

A detailed comparison is presented of the actual shear in the atmosphere with the critical shear given by the two-layer model of Phillips (1954), in which there is a critical temperature gradient separating stable conditions from baroclinically unstable ones. A very simple parameterization of the effect of eddy fluxes on atmospheric temperature is suggested, where the parameterization includes beta effects. The parameterization is illustrated by applying it in a one-dimensional heat-balance climate model. Enhancement of the eddy flux in a continuous atmosphere under supercritical conditions is stressed. This enhancement leads to a negative feedback between the meridional eddy flux of heat and the meridional temperature gradient. The feedback restricts gradients to values near the critical value, a process referred to as baroclinic adjustment. This should facilitate the development of simple climate models involving feedbacks associated with both the meridional vertical temperature structure.

Stone, P. H.↗

NASA NTRS1986

A study of the response of deep tropical clouds to mesoscale processes - Three-dimensional numerical experiments

A three-dimensional numerical cloud model is utilized to study the collective feedback effects of cloud systems on the large-scale environment. The characteristic features of the cloud ensemble model are described. The model is applied to a case of a well-defined ITCZ rainband over the eastern tropical Atlantic ocean. The location, number, and configuration of the clouds that formed in the model during the simulations are examined. The cloud heating and moistening, environmental wind shear, and vertical transports of momentum are analyzed. The collective feedbacks of cloud systems on the large-scale temperature, moisture, and horizontal momentum fields are computed after each simulation. The effect of the pressure gradient force on the v-momentum flux is investigated. It is observed that there is good correlation between the model simulations and the observed data.

Tao, Wei-Kuo↗

NASA NTRS1976

Geothermal flux through palagonitized tephra, Surtsey, Iceland - The Surtsey temperature-data-relay experiment via Landsat-1

The net geothermal flux through palagonitized basaltic tephra rims of the Surtur I and Surtur II craters at Surtsey, Iceland, in 1972, is estimated at 780 plus or minus 325 microcal/sq cm/s, indicating a decline since 1969 when a flux of 1500 microcal/sq cm/s was estimated. Heat flux in this range characterizes the postvolcanic environment on Surtsey in which the subaerial polagonitization of basaltic tephra is associated with mass transfer of hydrothermal vapor, either of meteoric or sea-water origin, only a few years after cessation of eruptive activity. The flux estimation is the result of the Surtsey data-relay experiment via Landsat-1 which was carried out in several phases. Temperature data were transmitted for a 38-day period in November and December 1972. A near-surface vertical gradient of 69.4 C/m was obtained, suggesting a mixed mechanism of heat transfer, partitioned between conduction and convection.

Friedman, J. D.↗

NASA NTRS1981

Transport of infrared radiation in cuboidal clouds

The transport of infrared radiation in a single cuboidal cloud using a vertical two steam approximation was modeled. The emittance of the top face of the model cloud is always less than that for a plane parallel cloud of the same optical depth. The hemisphere flux escaping from the cloud top has a gradient from the center to the edges which brighten when the cloud is over warmer ground. Cooling rate calculations in the 8 to 13.6 micrometer region show that there is cooling from the sides of the cloud at all levels even when there is heating of the core from the ground below. The radiances exiting from model cuboidal clouds were computed by path integration over the source function obtained with the two stream approximation. It is suggested that the brightness temperature measured from finite clouds will overestimate the cloud top temperature.

HARSHVARDHAN↗

NASA NTRS1989

Flux divergence of thermal radiation within stratiform clouds

A longwave radiative transfer model based on the narrow-band transmission functions of Lowtran 6 has been developed to investigate the profile of radiative flux divergence at the tops of stratiform clouds. The model has a vertical resolution of 1 m and a spectral resolution of 20/cm. These features are shown to permit modeling of broadband radiative transfer in the presence of a cloud top inversion and finite liquid water gradient. Radiative cooling rates in excess of 40 K/hr are shown to be possible for extreme conditions, with values of 10 K/hr being representative of more typical conditions. The values obtained are shown to be strongly influenced by the magnitude of liquid water concentration, by the gradient of liquid water in the top 50 m, and by the strength of the temperature inversion.

Davies, Roger↗

NASA NTRS1995

Adaptively Refined Euler and Navier-Stokes Solutions with a Cartesian-Cell Based Scheme

A Cartesian-cell based scheme with adaptive mesh refinement for solving the Euler and Navier-Stokes equations in two dimensions has been developed and tested. Grids about geometrically complicated bodies were generated automatically, by recursive subdivision of a single Cartesian cell encompassing the entire flow domain. Where the resulting cells intersect bodies, N-sided 'cut' cells were created using polygon-clipping algorithms. The grid was stored in a binary-tree data structure which provided a natural means of obtaining cell-to-cell connectivity and of carrying out solution-adaptive mesh refinement. The Euler and Navier-Stokes equations were solved on the resulting grids using an upwind, finite-volume formulation. The inviscid fluxes were found in an upwinded manner using a linear reconstruction of the cell primitives, providing the input states to an approximate Riemann solver. The viscous fluxes were formed using a Green-Gauss type of reconstruction upon a co-volume surrounding the cell interface. Data at the vertices of this co-volume were found in a linearly K-exact manner, which ensured linear K-exactness of the gradients. Adaptively-refined solutions for the inviscid flow about a four-element airfoil (test case 3) were compared to theory. Laminar, adaptively-refined solutions were compared to accepted computational, experimental and theoretical results.

Coirier, William J.↗

NASA NTRS2008

Photosynthetic Control of Atmospheric Carbonyl Sulfide during the Growing Season

Climate models incorporate photosynthesis-climate feedbacks, yet we lack robust tools for large-scale assessments of these processes. Recent work suggests that carbonyl sulfide (COS), a trace gas consumed by plants, could provide a valuable constraint on photosynthesis. Here we analyze airborne observations of COS and carbon dioxide concentrations during the growing season over North America with a three-dimensional atmospheric transport model. We successfully modeled the persistent vertical drawdown of atmospheric COS using the quantitative relation between COS and photosynthesis that has been measured in plant chamber experiments. Furthermore, this drawdown is driven by plant uptake rather than other continental and oceanic fluxes in the model. These results provide quantitative evidence that COS gradients in the continental growing season may have broad use as a measurement-based photosynthesis tracer.

Campbell, J. Elliott↗

DOE OSTI2020

Summertime Post-Cold-Frontal Marine Stratocumulus Transition Processes over the Eastern North Atlantic

Abstract Marine boundary layer (MBL) cloud morphology associated with two summertime cold fronts over the eastern North Atlantic (ENA) is investigated using high-resolution simulations from the Weather Research and Forecasting (WRF) Model and observations from the Atmospheric Radiation Measurement (ARM) ENA Climate Research Facility. Lagrangian trajectories are used to study the evolution of post-cold-frontal MBL clouds from solid stratocumulus to broken cumulus. Clouds within specified domains in the vicinity of transitions are classified according to their degree of decoupling, and cloud-base and cloud-top breakup processes are evaluated. The Lagrangian derivative of the surface latent heat flux is found to be strongly correlated with that of the cloud fraction at cloud base in the simulations. Cloud-top entrainment instability (CTEI) is shown to operate only in the decoupled MBL. A new indicator of inversion strength at cloud top that employs the vertical gradients of equivalent potential temperature and saturation equivalent potential temperature, which can be computed directly from soundings, is proposed as an alternative to CTEI. Overall, results suggest that the deepening–warming hypothesis suggested by Bretherton and Wyant explains many of the characteristics of the summertime postfrontal MBL evolution of cloud structure over the ENA, thereby widening the phase space over which the hypothesis may be applied. A subset of the deepening–warming hypothesis involving warming initially dominating over moistening is proposed. It is postulated that changes in climate change–induced modifications in cold-frontal structure over the ENA may be accompanied by coincident changes in the location and timing of MBL cloud transitions in the post-cold-frontal environment.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1991

Cold-air outbreak during GALE - Lidar observations and modeling of boundary layer dynamics

Two cold-air outbreaks were studied during the Genesis of Atlantic Lows Experiment. A lidar system was operated to observe the boundary layer evolution and the development of clouds. On the first day (January 30, 1986) boundary layer rise was less than 50 percent of the value for the second day (March 2, 1986). On the first day only a thin broken cloud cover formed, while on the second day a thick solid cloud deck formed - although the average moisture content was 60 percent of that on the first day. A trajectory slab model was employed to simulate the evolution of the layer over the ocean near the east Atlantic shore. The model allows for vertical gradients in conservative variables under neutrally buoyant conditions. The primary effect of these assumptions, which are based on observed thermodynamic profiles, is to reduce cloudiness to be more in line with observations. Boundary-layer depth was reasonably well predicted as was sensible and latent heat flux.

Boers, Reinout↗

NASA NTRS1998

A Parameterization for the Triggering of Landscape Generated Moist Convection

A set of relatively high resolution three-dimensional (3D) simulations were produced to investigate the triggering of moist convection by landscape generated mesoscale circulations. The local accumulated rainfall varied monotonically (linearly) with the size of individual landscape patches, demonstrating the need to develop a trigger function that is sensitive to the size of individual patches. A new triggering function that includes the effect of landscapes generated mesoscale circulations over patches of different sizes consists of a parcel's perturbation in vertical velocity (nu(sub 0)), temperature (theta(sub 0)), and moisture (q(sub 0)). Each variable in the triggering function was also sensitive to soil moisture gradients, atmospheric initial conditions, and moist processes. The parcel's vertical velocity, temperature, and moisture perturbation were partitioned into mesoscale and turbulent components. Budget equations were derived for theta(sub 0) and q(sub 0). Of the many terms in this set of budget equations, the turbulent, vertical flux of the mesoscale temperature and moisture contributed most to the triggering of moist convection through the impact of these fluxes on the parcel's temperature and moisture profile. These fluxes needed to be parameterized to obtain theta(sub 0) and q(sub 0). The mesoscale vertical velocity also affected the profile of nu(sub 0). We used similarity theory to parameterize these fluxes as well as the parcel's mesoscale vertical velocity.

Lynn, Barry H.↗

NASA NTRS2012

Upper-Ocean Heat Balance Processes and the Walker Circulation in CMIP5 Model Projections

Considerable uncertainty remains as to the importance of mechanisms governing decadal and longer variability of the Walker Circulation, its connection to the tropical climate system, and prospects for tropical climate change in the face of anthropogenic forcing. Most contemporary climate models suggest that in response to elevated CO2 and a warmer but more stratified atmosphere, the required upward mass flux in tropical convection will diminish along with the Walker component of the tropical mean circulation as well. Alternatively, there is also evidence to suggest that the shoaling and increased vertical stratification of the thermocline in the eastern Pacific will enable a muted SST increase there-- preserving or even enhancing some of the dynamical forcing for the Walker cell flow. Over the past decade there have been observational indications of an acceleration in near-surface easterlies, a strengthened Pacific zonal SST gradient, and globally-teleconnected dislocations in precipitation. But is this evidence in support of an ocean dynamical thermostat process posited to accompany anthropogenic forcing, or just residual decadal fluctuations associated with variations in warm and cold ENSO events and other stochastic forcing? From a modeling perspective we try to make headway on this question by examining zonal variations in surface energy fluxes and dynamics governing tropical upper ocean heat content evolution in the WCRP CMIP5 model projections. There is some diversity among model simulations; for example, the CCSM4 indicates net ocean warming over the IndoPacific region while the CSIRO model concentrates separate warming responses over the central Pacific and Indian Ocean regions. The models, as with observations, demonstrate strong local coupling between variations in column water vapor, downward surface longwave radiation and SST; but the spatial patterns of changes in the sign of this relationship differ among models and, for models as a whole, with observations. Our analysis focuses initially on probing the inter-model differences in energy fluxes / transports and Walker Circulation response to forcing. We then attempt to identify statistically the El Nino- / La Nina-related ocean heat content variability unique to each model and regress out the associated energy flux, ocean heat transport and Walker response on these shorter time scales for comparison to that of the anthropogenic signals.

Robertson, F. R.↗

NASA NTRS1978

A two-dimensional solar model

A two-dimensional, nearly flux constant solar atmosphere which includes a physical model of granulation is presented. If the inhomogeneous character of the solar photosphere is neglected in the construction of empirical models, the result is an underestimate of the temperature gradient in the layers below a Rosseland optical depth of 0.5 and an overestimate of the convective flux penetration into the observable layers. The two-dimensional model adequately reproduces the mean limb darkening, the magnitude and center-to-limb variation of the intensity fluctuations, and the rms vertical and horizontal convective velocities.

Nelson, G. D.↗

NASA NTRS1984

A dynamical basis for the parameterization of organized deep convection in large-scale numerical models

A hierarchy of steady, nonlinear, semianalytic models of different types of convection were produced. These provide a theoretical framework for determining cloud outflow fluxes of both dynamic and thermodynamic quantities, which can be used to formulate dynamical transports in parameterization schemes. This was achieved by exploiting certain Lagrangian conservation properties of steady flow, from which an equation for the vertical displacement of particles can be obtained and the outflow entropy, energy and momentum fluxes and the infow/outflow mass fluxes can be determined from solution to the equation. These fluxes are determined in terms of grid scale parameters such as convective available potential energy (CAPE), cloud layer shear, and horizontal pressure gradients. Five main types of system models are identified, respectively representing archtypes of convection in zero shear, large shear, midlatitude squall lines, tropical squall lines and cellular convection. The downdraught is an important aspect in the first four of these and the cloud scale transport of momentum is very distinctive.

Moncrieff, M. W.↗