Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Reservoir simulation”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Optimizing district energy systems by integrating Borehole Thermal Energy Storage Using a Mixed-Integer Linear Programming g-function framework with a Multi-Timescale Rolling Horizon method

Shallow geothermal has gained increasing attention in recent years; however, a reliable framework for its accurate incorporation into large-scale energy system optimization remains lacking. This study proposes a Mixed-Integer Linear Programming (MILP) framework combined with the g-function approach to integrate Borehole Thermal Energy Storage (BTES) technology into energy system optimization. Validation against a Modelica-based reservoir network simulation demonstrates that the proposed framework effectively captures the ground thermal response under varying energy loads and accurately estimates the borefield energy supply. To enhance scalability, a Rolling Horizon with Multi-Timescale (RH-MTS) method is further introduced, reducing computational time by 73 % for the 1-year optimization model with only minor loss of optimality. The framework is demonstrated through the case study of the UC Berkeley campus. Results indicate that BTES is a cost-effective and low-carbon solution: two borefields comprising 382 boreholes can meet 8.0 % and 6.6 % of the total campus heating and cooling demand, respectively, at an average energy rate of 0.70–0.77 USD/kWh and carbon intensity of 0.54 kg-CO2/kWh. Short-term analysis reveals a 35%–65% decline in BTES energy flow after 3–6 months of continuous heating/cooling operation, while long-term simulation shows that annual energy production of BTES can vary by up to 12.0 % after four years before stabilizing. Overall, this study develops a novel optimization framework that couples physics-based g-function method with MILP optimization framework, thereby advancing methodological development for shallow-geothermal integration and providing actionable guidance for BTES deployment in district-energy systems.

Yang, Jiahui↗
DOE OSTI2026

Impact of Rock Wettability and Mineralization on CO 2 Storage Efficiency in Basalt Reservoirs

Basalt formations have emerged as highly promising targets for geological CO 2 storage due to their abundance of reactive silicate minerals that can rapidly convert dissolved CO 2 into stable carbonate minerals. The efficiency of CO 2 trapping in basalts, however, depends not only on their geochemical reactivity but also on reservoir wettability, which governs how CO 2 partitions among structural, residual, solubility, and mineral trapping mechanisms. However, wettability and mineral kinetics have largely been examined separately, wettability for plume migration and kinetics in terms of mineralization, leaving their coupled impact unquantified. Here, we systematically integrate th ese processes in a reservoir-scale reactive transport model to quantify the impact of three distinct wettability conditions (water-wet, mixed-wet, and CO 2 -wet) on plume migration, mineral precipitation, pore-scale blockage, and long-term immobilization. Results show that under ideal conditions (no mineralization), water-wet systems exhibited plume aspect ratios up to 28% higher than CO 2 -wet systems and promoted at least 60% greater residual trapping. When mineralization was incorporated, mineral trapping accounted for 57%, 53%, and 46% of the total immobilized CO 2 in water-wet, mixed-wet, and CO 2 -wet systems, respectively. Dissolved CO 2 contributed 43-46%, while residual trapping remained below 5%, and mobile CO 2 was substantially reduced due to mineralization, demonstrating that as the system becomes more water-wet, CO 2 is more effectively converted into stable carbonate minerals, although extensive mineral precipitation may lead to pore clogging. A global Sobol-Morris sensitivity analysis further reveals that the kinetic parameters of reactive minerals are the most influential factor, followed by wettability, while salinity has only a minor influence. Furthermore, these findings highlight that wettability is not merely a flow property but a fundamental control on CO 2 mineralization in basalts, with direct implications for site selection for storage projects.

Basalts↗
DOE OSTI2026

WBM-Hydro-ReEDS

WBM-hydro-ReEDS is a tool for converting daily reservoir release simulations from the WMB model to monthly, plant-level hydropower generation estimates, and then to inputs for the ReEDS model.

Turner, Sean [Oak Ridge National Laboratory (ORNL)↗
DOE OSTI2024

Impacts of Alternative Operations and Renewable Energy Deployment on Columbia River Hydropower

The Columbia River Treaty Tribes in the Pacific Northwest - the Nez Perce, Umatilla, Warm Springs, and Yakama - hold treaty-reserved fishing rights for the Columbia River, one of the world's most productive salmon rivers and a critical resource for these tribes. However, the tribes have expressed that current operating regimes of hydropower dams throughout the Columbia River Basin (CRB) do not fully account for tribal fishing rights and have negatively impacted fish populations. Four tribes acting together through the Columbia River Inter-Tribal Fish Commission (CRITFC)'s 2022 Energy Vision stated that current power system models often do not fully account for the many constraints faced by hydropower facilities in the CRB. As the deployment of variable renewable energy (VRE) like solar and wind continues to increase, power system flexibility will become increasingly important. As a result, there is a growing need to better understand the true capabilities of the hydropower generation fleet while accurately accounting for ecological constraints (Northwest Power and Conservation Council 2022). Understanding hydropower's role in a future grid with a higher share of VRE can inform water resources planning to address ecological needs. The goal of this study is to examine the impacts of alternative hydropower operation rules and weather variability on hydropower generation and grid operations in VRE and transmission infrastructure deployment scenarios. The study evaluates how hydropower scheduling practices can reduce ecological impacts to the region's salmon populations. More specifically, the study examines the impact of today's dam water release rules (called "adjusted base water rules" in this study) and new, ecologically informed water release rules (called "ecosystem water rules" in this study) on two VRE scenarios (current renewable energy levels and higher renewable energy levels) using six weather years for each scenario (2008-2013). CRITFC developed the ecosystem water rules, which capture a portion of the changes they recommend for CRB hydropower operations. These analyses use a water resource planning model (OASIS) and a production cost model (PLEXOS) to simulate CRB reservoir cascade and Western Interconnection power grid operation.

13 HYDRO ENERGY↗
DOE OSTI2026

Oil in Alaska North Slope Gas Hydrate Reservoir: Micro-CT and Flow Simulation Insights into Permeability

Gas hydrate-bearing sands on the Alaska North Slope (ANS) host minor volumes of crude oil whose impact on formation permeability has never been quantified relative to a gas hydrate reservoir system. Here, we combine in situ pressure-core microcomputed-tomography (μ-CT), thermogravimetric analysis (TGA), gas-chromatography–mass-spectrometry (GC–MS), scanning-electron microscopy coupled with energy-dispersive spectroscopy (SEM-EDS), and Stokes-flow simulation to (i) quantify oil saturation, (ii) infer its source, and (iii) evaluate its influence on permeability after gas-hydrate dissociation. μ-CT and SEM-EDS imaging identify the host sediment as silty, with some grain-coating clay present. Eight representative 3003-voxel subvolumes extracted from a preserved pressure core (Core 15P-3c, from the B1 sand (unit B), HYDRATE 02 Geo Data Well (GDW), 886.96–887.29 m measured depth, MD) exhibit porosities of 0.34–0.37 and an average oil saturation of 0.09 ± 0.03, in comparison to an independent TGA-based estimate of 0.16. GC-MS fingerprints obtained from an adjacent depressurized core (Core 17P-1, 891.37–891.44 m MD), together with oil pore habit revealed by μ-CT, suggest that the oil is partly native but also includes contributions from mineral oil-based drilling fluid contamination. Further analysis of oil–porewater interactions and flushing tests confirms that native oil saturation in the analyzed samples ranges from 0.04 to 0.08 with an average of approximately 0.06, and the oil is highly immobile. Flow simulations demonstrate that the native oil saturation of ≈0.06 reduces permeability by approximately 50% due to both pore blockage and increased flow-path tortuosity. These findings reveal the potential presence of native oil within the ANS gas hydrate reservoir and suggest that future simulation models may need to account for its impact on permeability to improve long-term performance predictions of gas and water production.

02 PETROLEUM↗
DOE OSTI2025

High-Resolution Simulations of Geological CO 2 Injection: Application to the SPE11 Benchmark

Geological carbon sequestration (GCS) will play a critical role in decarbonization and in facilitating the transition to clean energy systems. Because CO 2 is highly mobile, ensuring its safe and permanent injection into subsurface geological formations involves monitoring over larger spatial domains and longer time periods than is typical for hydrocarbon reservoirs. This can benefit from simulation tools capable of modeling key CO 2 trapping mechanisms, particularly those optimized for speed and scalability on high-performance computing systems. Using isothermal versions of the SPE11B and SPE11C benchmark cases, we conduct a mesh refinement study simulating CO 2 injection into kilometer-scale rock formations at centimeter resolution with the GEOS open-source simulation framework. We focus on how mesh refinement improves the accuracy of convective mixing in both 2D and 3D simulations. The computational costs associated with achieving a converged solution highlight the need for predictive upscaling techniques. A systematic performance scaling analysis—including both central processing unit (CPU) and graphics processing unit (GPU) architectures—complements the “Results” section.

Geosciences↗
DOE OSTI2026

ELM‐MOSART‐DOC: A Large‐Scale Riverine Dissolved Organic Carbon Model and Its Application Over the United States

Riverine dissolved organic carbon (DOC), primarily sourced from soil organic carbon (SOC), plays a crucial role in regional and global carbon cycles. However, the complexities of the underlying mechanisms and limited observations present significant challenges for predictive understanding of DOC at regional or larger scales. Recently, we developed a machine learning‐based (ML) map of DOC transformation rates, bridging the gap between SOC and DOC leaching flux and simplifying terrestrial DOC representation. Building on this advancement, we introduce ELM‐MOSART‐DOC, a DOC module integrated into the riverine component of the Energy Exascale Earth System Model (E3SM)—the Model for Scale Adaptive River Transport (MOSART). ELM‐MOSART‐DOC simulates DOC transport and transformation across both headwater streams and river networks, including those managed. Model validation demonstrates the ability of ELM‐MOSART‐DOC to accurately capture long‐term average DOC concentrations, with Kling‐Gupta Efficiency (KGE) scores of 0.58 and 0.76 at large and local stations, respectively. We further assess the impact of reservoirs through different simulation schemes, revealing that reservoirs significantly alter DOC fluxes by regulating streamflow patterns and promoting DOC mineralization. Model simulations indicate that reservoirs reduce total DOC flux from the Mississippi River into the ocean by 7.5%, with the long‐term average annual export decreasing from 3.34 to 3.14 teragrams (Tg) per year. ELM‐MOSART‐DOC integrates process‐based modeling with ML parameterization to enhance the predictive understanding of riverine biogeochemical processes. This approach reduces uncertainties in modeling regional and global carbon cycle ESMs and provides new insights into carbon cycling and its implications for global environmental change.

Li, Lingbo [Univ. of Houston, TX (United States); ↗
DOE OSTI2025

Reservoir Dynamics in Proposed Operational Scenarios Where CO 2 -EOR Fields are Transitioned From CO 2 -Flood Enhanced Oil Recovery to Dedicated Carbon Storage: A Field Case Study

This study models the transition of CO 2 -enhanced oil recovery (CO 2 -EOR) fields to dedicated carbon storage using a generalized reservoir model informed by real field data, SACROC. Through scenario-based simulations, the work evaluates how reservoir depletion levels, boundary conditions, fluid properties, and domain size affect pressure buildup, CO 2 plume migration, and long-term containment performance. The results inform practices for repurposing oil fields into reliable CO 2 storage sites.

02 PETROLEUM↗
DOE OSTI2026

WellPINN: Accurate Well Representation for Transient Fluid Pressure Diffusion in Subsurface Reservoirs With Physics‐Informed Neural Networks

Accurate representation of pumping wells is essential for reliable reservoir characterization and simulation of operational scenarios in subsurface flow models. Physics-informed neural networks (PINNs) are emerging as a promising alternative to numerical models for reservoir modeling, offering seamless integration of monitoring data and governing physical equations. However, existing PINN-based studies face major challenges in capturing fluid pressure near wells when using a source/sink term, particularly during the early stages after pumping begins. We address this problem by introducing WellPINN, a workflow in which an initially trained PINN infers fluid pressure across the entire reservoir domain using a large equivalent well radius. This initial PINN solution is then locally refined around the well by a set of subdomain PINNs that are trained for smaller equivalent well radii. Continuity across these subdomain interfaces as well as at the initial condition is ensured by hard-constraining each PINN on its subdomain boundary. Our results demonstrate WellPINN as the first workflow of its kind to focus on accurate inference of fluid pressure from pumping rates throughout the entire injection period, significantly advancing the potential of PINNs for inverse modeling and operational scenario simulations. All data and code for this paper are openly available at https://doi.org/10.20350/DIGITALCSIC/17260.

58 GEOSCIENCES↗
DOE OSTI2026

Energy to Communities Program Short Introduction: Evaluation of Enhanced Geothermal System for Cornell District Heating [Slides]

Through the U.S. Department of Energy's Energy to Communities (E2C) program, NLR, other national laboratory experts, and select organizations provide Expert Match - free, short-term technical assistance to address near-term energy challenges and questions. Expert Match is for community stakeholders who have decision-making power or influence in their community but need access to additional energy expertise to inform key upcoming decisions. This Expert Match request supported Cornell University (Itaca, NY) community efforts in developing Enhanced Geothermal System (EGS) reservoir performance, by simulating fracture flow to estimate pressure drop, flow distribution, and heat extraction efficiency in fractured rock. Thermal drawdown modeling was used to assess long-term resource sustainability, while system optimization studies refined key parameters such as well configuration and flow rate to improve overall reservoir efficiency.

15 GEOTHERMAL ENERGY↗
DOE OSTI2024

EGS Reservoir Modeling for Developing Geothermal District Heating at Cornell University: Preprint

Cornell University is pursuing development of an enhanced geothermal system (EGS) for providing heating to its main campus in Upstate New York. A ~10,000 ft (~3 km) deep vertical observation well ("CUBO") was drilled in 2022 to characterize the subsurface using wellbore logging, borehole imaging, fluid sampling, mini-frac tests, coring and drill cutting analysis. Down-hole temperatures measured at 3 km depth are about 80 degrees C, sufficiently high for direct-use heating. The well drilled through generally low porosity and low permeability Paleozoic sedimentary formations and into metamorphic basement rock, encountered at about 9,400 ft depth. Leveraging subsurface data obtained through CUBO, we investigated technical feasibility and design requirements of a doublet well system with horizontal laterals connected to a fracture network created through hydraulic fracturing. The EGS reservoir is sized to provide a nominal heat output in the range 5 to 10 MWth of continuous heating over a 15 year-lifetime with limited thermal drawdown. We applied the Gringarten Multiple parallel fracture model, the Cornell Discrete Fracture Simulator FOXFEM and the commercial simulator ResFrac to estimate required heat transfer area and design a potential hydraulic stimulation treatment. Reservoir simulations indicate that, depending on fluid flow rate and injection temperature, 2 to 3 km2 of effective fracture heat transfer area is required to supply the target heat output of 5 to 10 MWth over 15 years.

district heating↗
DOE OSTI2025

EGS Reservoir Modeling for Developing Geothermal District Heating at Cornell University

Cornell University is pursuing development of an enhanced geothermal system (EGS) for providing heating to its main campus in Upstate New York. A ~10,000 ft (~3 km) deep vertical observation well ('CUBO') was drilled in 2022 to characterize the subsurface using wellbore logging, borehole imaging, fluid sampling, mini-frac tests, coring and drill cutting analysis. Down-hole temperatures measured at 3 km depth are about 80 degrees C, sufficiently high for direct-use heating. The well drilled through generally low porosity and low permeability Paleozoic sedimentary formations and into metamorphic basement rock, encountered at about 9,400 ft depth. Leveraging subsurface data obtained through CUBO, we investigated technical feasibility and design requirements of a doublet well system with horizontal laterals connected to a fracture network created through hydraulic fracturing. The EGS reservoir is sized to provide a nominal heat output in the range 5 to 10 MWth of continuous heating over a 15 year-lifetime with limited thermal drawdown. We applied the Gringarten Multiple parallel fracture model, the Cornell Discrete Fracture Simulator FOXFEM and the commercial simulator ResFrac to estimate required heat transfer area and design a potential hydraulic stimulation treatment. Reservoir simulations indicate that, depending on fluid flow rate and injection temperature, 2 to 3 km2 of effective fracture heat transfer area is required to supply the target heat output of 5 to 10 MWth over 15 years.

district heating↗
DOE OSTI2025

Foundational Dataset for Developing Large-Sample Stream Temperature Models in the Conterminous United States

This dataset provides inputs, evaluation results, and trained weights from a large-sample Long Short-Term Memory (LSTM) model designed to predict daily stream temperatures across unregulated river reaches in the conterminous United States (CONUS). It includes dynamic meteorological and hydrologic forcings, static physiographic attributes, and model outputs from cross-validation experiments spanning 300 basins. It supports reproducible modeling, direct application for new basins, and provides data suitable for integration with reservoir and river simulations under current and future climates. It contains two .zip files described below · RQ-AI_runs.zip: Model outputs from 10-fold cross-validation experiments, including observed and predicted daily stream temperatures, along with test performance metrics for water years 2017–2019. Two versions are included: 1. Model trained and validated using subbasin-area weighted dynamic features. 2. Model trained and validated using whole-basin area weighted dynamic features. · RQ-AI_inputs.zip: Collection of all formatted dynamic and static predictor datasets (meteorological, hydrologic, and physiographic features) used in model training and analysis. Detailed instructions and data structure is held at the following GitLab repository: https://code.ornl.gov/tempwise/training.

Gomez-Velez, Jesus [Oak Ridge National Laboratory ↗
NASA NTRS1984

Vectorization of a penalty function algorithm for well scheduling

In petroleum engineering, the oil production profiles of a reservoir can be simulated by using a finite gridded model. This profile is affected by the number and choice of wells which in turn is a result of various production limits and constraints including, for example, the economic minimum well spacing, the number of drilling rigs available and the time required to drill and complete a well. After a well is available it may be shut in because of excessive water or gas productions. In order to optimize the field performance a penalty function algorithm was developed for scheduling wells. For an example with some 343 wells and 15 different constraints, the scheduling routine vectorized for the CYBER 205 averaged 560 times faster performance than the scalar version.

Absar, I.↗
NASA NTRS2016

Validation of a Cephalad Fluid Shift Countermeasure

INTRODUCTION: This project will provide critical data required to objectively determine how an optimized thigh cuff could be incorporated into the NASA integrated physiological countermeasure suite. This project will determine if thigh cuffs used during simulated spaceflight impact intracranial pressure (ICP), ocular structure and function, and intraocular pressure (IOP) using state of-the-art techniques. Additionally, some of the same methods, hardware, and protocols will be employed in the present investigation to enable direct comparisons to the International Space Station (ISS) "Fluid Shifts" experiment with Chibis-Lower Body Negative Pressure (LBNP). This study will determine the temporal physiological responses of thigh cuff application and removal on ocular and cerebral variables (including invasive ICP) in a microgravity analog. Furthermore, this proposed study will determine tissue pressure distribution applied by thigh cuffs in order to improve comfort, mobility, and efficacy of the countermeasure. Our specific aim is to determine the efficacy of a novel thigh cuff device to mitigate cephalad fluid shifts. We hypothesize that a thigh cuff countermeasure employed in a microgravity analog will temporarily reverse or attenuate ocular and cerebral-volume-pressure variables, approaching normal Earth-based seated posture, the most frequent posture assumed in daily life. In addition, we hypothesize that the magnitude of fluid and pressure redistribution using a thigh cuff countermeasure may require a longer exposure time than that of Chibis-LBNP (using ground-based data from our "Fluid Shifts" project). This project directly addresses Critical Path Roadmap Risks and Questions regarding "Risk of Spaceflight-Induced Intracranial Hypertension/Vision Alterations," and IRP Gap VIIP13: We need to identify preventative and treatment countermeasures to mitigate changes in ocular structure and function and intracranial pressure during spaceflight. METHODS: Noninvasive measures and tissue pressure distributions beneath thigh cuffs The objectives of this study are to: 1) determine the distribution of skin surface pressures beneath the advanced thigh cuff in ten subjects, 2) calibrate the built-in pressure measurement system of the advanced thigh cuff using an industry standard device, and 3) collect subjective feedback and data on the new cuff design to allow for further adjustments prior to invasive studies. A Tekscan Industrial Sensing (I-Scan) system will measure the pressure distribution of the advanced thigh cuff against the skin. In addition, we will measure blood pooling in the thigh and record the circumference of the thigh using Hokanson strain gauge plethysmography. The advanced thigh cuff will be adjusted to obtain a skin contact pressure of 30-50 mmHg as visualized on the Tekscan system. The built-in advanced thigh cuff pressure monitor will be recorded simultaneously to allow direct comparison to the Tekscan measurements. The volunteer will then remove the thigh cuff and remain at rest for five minutes with no legging applied. The thigh cuff will be donned again and pressure measurements will be taken in the same manner for up to 10 repetitions to show reproducibility of pressure after donning. At the conclusion of the study, subjects will be asked to flex their knee, stand, walk, and sit with the thigh cuff activated. During each of these maneuvers the subject will rate their pain/comfort using a modified Borg scale. Effect of thigh cuffs on ICP during simulated microgravity Ommaya reservoir patients will be recruited from the John Wayne Cancer Institute. Ommaya reservoirs provide safe and direct access for the measurement of ICP. Subjects will be instrumented for continuous blood pressure, ECG, and invasive ICP measures. The subjects will be positioned in the upright sitting posture for a 10-minute stabilization period. After the 10-minute stabilization period, imaging measures [ICP, Optical Coherence Tomography, IOP, ocular and vascular ultrasound] will be performed. Following baseline seated measures, the subject will be positioned randomly in the supine, 15deg head-down-tilt, and 15deg head-down-tilt with thigh cuffs and measures repeated. DISCUSSION: Tests to down-select thigh cuff designs will occur in early 2016. Invasive ICP and noninvasive eye imaging tests will begin in spring 2016. Supported by NSBRI through NCC 9-58.

Macias, B.↗
NASA NTRS2019

Measuring the D/H Ratios of Exoplanets and Brown Dwarfs

The relative abundance of deuterium and hydrogen is a potent tracer of planet formation and evolution. Jupiter and Saturn have protosolar atmospheric D/H ratios, a relic of substantial gas accretion from the nebula, while the atmospheres of Neptune and Uranus are enhanced in D by accretion of ices into their envelopes. For terrestrial planets, D/H ratios are used to determine the mechanisms of volatile delivery and subsequent atmosphere loss over the lifetime of the planet. Planets and brown dwarfs more massive than ~13 M(J) quickly fuse their initial D reservoir. Here, we simulate spectra for giant exoplanets and brown dwarfs (2 M(Neptune) to ~10 M(Jupiter)) from T(eff) = 200–1800 K including both CH3D and HDO to determine the observability of these dominant deuterium isotopologues in mid-infrared thermal emission spectra. Colder objects have stronger molecular features in their spectra, due to the temperature dependence of molecular cross sections. CH3D is easier to observe than HDO at all temperatures considered, due to the strength of its absorption bands and locations of features at wavelengths with few other strong absorption features. We predict that for nearby cool brown dwarfs, the CH3D abundance can be measured at high signal to noise with the James Webb Space Telescope (JWST); for objects from 200 to 800 K closer than 10 pc, a protosolar D/H ratio would be readily observable in 2.5 hr. Moderately young Jupiter-mass planets (100–300 Myr) and young Neptunes (10 Myr) may be discovered with JWST and provide the best targets for detecting deuterium on an exoplanet in the coming decade. Future telescope designs should consider the importance of isotopes for understanding the formation and evolution of planetary atmospheres.

Caroline V. Morley↗
NASA NTRS2012

The Connection between Model Performance on the CCMVal Transport Diagnostics and Simulated Sensitivity of Ozone to Chlorine Change

The chemistry climate models that contributed simulations for past and future ozone evolution to the 2010 Scientific Assessment of Ozone Depletion were subject to extensive evaluation by the SPARC (Stratospheric Processes and their Role in Climate) CCMVal (Chemistry-Climate Model Validation) activity. The sensitivity of ozone to changes in composition and climate varies among the models, but the relationship between these variations and the model evaluations of CCMVal is not obvious. We have learned that the transport evaluation can be used to interpret the comparisons between observed and simulated columns of chlorine reservoirs, hydrochloric acid (HCl) and chlorine nitrate (ClONO2); these comparisons were part of the CCMVal evaluation of chemistry. The simulations with best performance on the transport diagnostics most faithfully reproduce the evolution and seasonal variation of the chlorine reservoirs as observed at NDACC (Network for Detection of Atmospheric Composition Change) stations (NyAlesund 78.9N, Kiruna 67.8N, Harestua 60.2N, Jungfraujoch 46.6N, Toronto 43.6N, Kitt Peak 31.9N, Izana 28.3N, Mauna Loa 19.5N, Lauder 45S and Arrival Heights 77.8S). In the simulations, the HCl in the lower stratosphere depends on total inorganic chlorine (Cly) and partitioning between HCl and ClON02. Total inorganic chlorine depends on the fractional release of chlorine from source gases, and ratio of ClON02 to HCl is inversely dependent on methane and varies quadratically with ozone. Simulated HCl from various models may agree with observations even though Cly is in error, partitioning is in error, or both. Simulated ozone sensitivity to chlorine is shown to be greater for models that produce larger values of chlorine nitrate for background chlorine levels, and vice versa. Comparisons with the NDACC data show why the models with 'best' transport have similar sensitivity to chlorine change. The realistic evolution of the simulated HCl and ClONO2 columns suggests realistic levels of Cly in the lower atmosphere. In addition, the wide range values for the sensitivity of ozone to chlorine obtained from the CCMVal simulations is explained by the wide range in lower atmospheric columns of ClONO2 and the concomitant wide range of levels for chlorine monoxide.

Douglass, Anne↗