Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Spreading models”

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 217 records · Page 12

Integrated assessment model diagnostics: key indicators and model evolution

Integrated assessment models (IAMs) form a prime tool in informing climate mitigation strategies. Diagnostic indicators that allow to compare these models can help to describe and explain differences in model projections. This also increases transparency and comparability. Earlier, the IAM community has developed an approach to diagnose models (Kriegler et al., 2015). Here we build on this, by proposing a selected set of well-defined indicators as a community standard, similar to metrics used for other modeling communities such as climate models. These indicators are the relative abatement index (RAI), emission reduction type index (ERT), inertia timescale (IT), fossil fuel reduction (FFR), transformation index (TI) and cost per abatement value (CAV). We apply the approach to 17 IAMs, including both older version as well as their latest versions, as applied in the IPCC 6th Assessment Report (AR6). The study shows that the approach can be easily applied and allows for comparison of model versions in time. The indicators and their trends can often be explained in terms of model characteristics and changes. We show that together, the set of six indicators can provide an useful indication of the main traits of the model and can roughly indicate the general model behavior. The results also show that there is often a considerable spread across the models. Interestingly, the diagnostic values often change for different model versions, but there does not seem to be a distinct trend across the different models.

54 ENVIRONMENTAL SCIENCES↗

Subpolar North Atlantic Mean State Affects the Response of the Atlantic Meridional Overturning Circulation to the North Atlantic Oscillation in CMIP6 Models

Abstract The Atlantic meridional overturning circulation (AMOC) plays an important role in climate, transporting heat and salt to the subpolar North Atlantic. The AMOC’s variability is sensitive to atmospheric forcing, especially the North Atlantic Oscillation (NAO). Because AMOC observations are short, climate models are a valuable tool to study the AMOC’s variability. Yet, there are known issues with climate models, like uncertainties and systematic biases. To investigate this, preindustrial control experiments from models participating in the phase 6 of Coupled Model Intercomparison Project (CMIP6) are evaluated. There is a large, but correlated, spread in the models’ subpolar gyre mean surface temperature and salinity. By splitting models into groups of either a warm–salty or cold–fresh subpolar gyre, it is shown that warm–salty models have a lower sea ice cover in the Labrador Sea and, hence, enable a larger heat loss during a positive NAO. Stratification in the Labrador Sea is also weaker in warm–salty models, such that the larger NAO-related heat loss can also affect greater depths. As a result, subsurface density anomalies are much stronger in the warm–salty models than in those that tend to be cold and fresh. As these anomalies propagate southward along the western boundary, they establish a zonal density gradient anomaly that promotes a stronger delayed AMOC response to the NAO in the warm–salty models. These findings demonstrate how model mean state errors are linked across variables and affect variability, emphasizing the need for improvement of the subpolar North Atlantic mean states in models.

54 ENVIRONMENTAL SCIENCES↗

A semi-empirical model for heat transfer and flame spreading in slotted rocket motors

This paper presents analytical and experimental results for the cold-flow aerothermodynamics of the head-end star slot region of the Space Shuttle solid rocket motor (SRM). Particular attention is paid to the determination of the heat transfer from the developing igniter flow to the propellant grain (slot wall). Igniter geometries considered include a single port configuration with flow along the axis of the motor (present Shuttle configuration) and a multiport 45-deg canted nozzle configuration which is similar (but not identical) to the proposed ASRM igniter geometry.

Jenkins, R. M.↗

Estimating the Value of Nuclear Integrated Hydrogen Production and the Dependency of Electricity and Hydrogen Markets on Natural Gas

Producing low carbon Hydrogen at a competitive price is one of the challenges to hydrogen being part of the solution to reach net-zero emission targets set by the U.S. DOE by 2050. With projected near-term improvements in technology, hydrogen production via solid oxide electrolysis cell (SOEC) / high-temperature steam electrolysis (HTSE) integrated with existing light water reactor (LWR) Nuclear Power Plants (NPP-HTSE) can produce carbon-free hydrogen competitively. In the near term, a 10-year production tax credit (PTC) found in the Inflation Reduction Act (IRA) has been passed, which will catalyze the development and improvement of hydrogen production technology to be competitive. The “1-1-1” target set by the U.S. DOE is to reduce the cost of carbon-free hydrogen by 80% to $1 per kilogram in 1 decade. Several models are available to analyze the profitability, opportunity, and technical capability of NPP-HTSE systems. In order of complexity from most complex to least complex some of these models include: RAVEN/HERON, process models using Aspen HYSYS and capital expense estimations using Aspen Process Economic Analyzer (APEA) and levelized cost of hydrogen (LCOH) calculation using the H2A model (Hydrogen Analysis Model), and custom spread sheets built by the interested party. Though some of the more advanced existing models provide detailed analysis to complex grid integrated problems, they also can take considerable time to setup and run. These advanced models are well suited to complex grid integrated analysis and the consideration of flexibility and variability of regulated and de-regulated electricity price and advanced estimation of capital and operating expenses and heat and material balances.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Sensitivity Analysis and Effective Parametrization of PEM Fuel Cell Models

The cost of proton-exchange-membrane fuel cells (PEMFCs) remains a major hurdle in large-scale commercialization of this technology. To improve their performance and reduce cost, novel materials and electrode designs are continuously envisioned, e.g., non-PGM catalyst layers, ultra-thin Pt/Pt-Ni based catalyst layers, structured ionomer arrays or NSTF catalyst layers.1 Understanding the impact of these improvement strategies can be extremely time and cost intensive due to complex physical phenomena and large design space. We have previously developed a PEMFC modeling framework2 which has been a time and cost effective tool for understanding and optimizing the complex multi-physics phenomena within PEMFCs; however, several of the cell parameters used in the modeling have large spread in measured data.3 Furthermore, several transport parameters such as water adsorption kinetics have not been accurately measured and the approximations are spread over several orders of magnitude. These uncertainties cause problems in ascertaining accuracy of the modeling approach and they reduce the predictive power of the numerical models. The aim of this work is to identify the sensitivity of PEFC numerical model outputs to various input parameters. The previously in-house developed MEA modeling framework2 is used for PEMFC modeling. The sensitivity of the model outputs with respect to inputs parameters is obtained by analyzing the condition numbers for different output-input pairs at varying operating conditions. An example of the sensitivity analysis is shown in Figure 1. The condition numbers are obtained for the entire possible range of input parameters at varying operating conditions to identify the most crucial parameters of the PEMFC model. Based on our preliminary analysis, parameters related to kinetics (exchange current density and ECSA) and heat/water management in electrodes and ionomer (ionomer fraction, thermal conductivity) are most crucial. One of the major goals of this work is to identify the most crucial set of parameters towards which the model shows maximum sensitivity. This will guide future experimentalists to measure these properties with higher accuracy. Furthermore, the sensitivity analysis will also enable us to optimize the PEMFC performance by selectively targeting the most sensitive parameters and thereby making the largest impact. Acknowledgements The work is funded under the Fuel Cell Performance and Durability Consortium (FC-PAD), by the Fuel Cell Technologies Office (FCTO), Office of Energy Efficiency and Renewable Energy (EERE), of the U.S. Department of Energy under contract number DE-AC02-05CH11231. The authors would like to thank Nathan Craig at Robert Bosch LLC for his valuable input in designing the sensitivity analysis. The authors would also like to thank Giovanna Bucci and Matthias Hanauer at Robert Bosch for their valuable inputs and discussion. References P. K. Sinha, W. Gu, A. Kongkanand and E. Thompson, J. Electrochem. Soc., 158, B831 (2011). L. M. Pant, M. R. Gerhardt, N. Macauley, R. Mukundan, R. L. Borup and A. Z. Weber, Electrochim. Acta, 326, 134963 (2019). R. Vetter and J. O. Schumacher, ArXiv181110091 Phys. (2018). Figure 1

Pant, Lalit↗

Observationally constrained analysis on the distribution of fine- and coarse-mode nitrate in global models

Nitrate plays an important role in the Earth system and air quality. A key challenge in simulating the life cycle of nitrate aerosol in global models is to accurately represent mass size distribution of nitrate aerosol. In this study, we evaluate the performance of the Energy Exascale Earth System Model version 2 (E3SMv2) and the Community Earth System Model version 2 (CESM2), along with Aerosol Comparisons between Observations and Models (AeroCom) phase III models, in simulating spatial distribution of fine-mode nitrate, the mass size distribution of fine- and coarse-mode nitrate, and the gas–aerosol partitioning between nitric acid gas and nitrate, using long-term ground-based observations and measurements from multiple aircraft campaigns. We find that most models underestimate the annual mean PM 2.5 (particulate matter with diameter less than 2.5 µm) nitrate surface concentration averaged over all sites. The observed nitrate PM 2.5 / PM 10 and PM 1 / PM 4 ratios are influenced by the relative contribution of fine sulfate or organic particles and coarse dust or sea salt particles. Overall, the ground-based observations give an annual mean surface nitrate PM 2.5 / PM 10 ratio of 0.7. Most models underestimate the annual mean PM 2.5 / PM 10 ratio in all regions. There are large spreads in the modeled nitrate PM 1 / PM 4 ratios, which span the full range from 0 to 1. Most models underestimate the surface molar ratio of nitrate to total inorganic nitrate averaged across all sites. Our study indicates the importance of gas–aerosol partition parameterization and the simulation of dust and sea salt in correctly simulating the mass size distribution of nitrate.

Nitrate↗

Exploring ship track spreading rates with a physics-informed Langevin particle parameterization

Abstract. The rate at which aerosols spread from a point source injection, such as from a ship or other stationary pollution source, is critical for accurately representing subgrid plume spreading in a climate model. Such climate model results will guide future decisions regarding the feasibility and application of large-scale intentional marine cloud brightening (MCB). Prior modeling studies have shown that the rate at which ship plumes spread may be strongly dependent on meteorological conditions, such as precipitating versus non-precipitating boundary layers and shear. In this study, we apply a Lagrangian particle model (PM-ABL v1.0), governed by a Langevin stochastic differential equation, to create a simplified framework for predicting the rate of spreading from a ship-injected aerosol plume in sheared, precipitating, and non-precipitating boundary layers. The velocity and position of each stochastic particle is predicted with the acceleration of each particle being driven by the turbulent kinetic energy, dissipation rate, momentum variance, and mean wind. These inputs to the stochastic particle velocity equation are derived from high-fidelity large-eddy simulations (LES) equipped with a prognostic aerosol–cloud microphysics scheme (UW-SAM) to simulate an aerosol injection from a ship into a cloud-topped marine boundary layer. The resulting spreading rate from the reduced-order stochastic model is then compared to the spreading rate in the LES. The stochastic particle velocity representation is shown to reasonably reproduce spreading rates in sheared, precipitating, and non-precipitating cases using domain-averaged turbulent statistics from the LES.

54 ENVIRONMENTAL SCIENCES↗

Accurate Simulation of Both Sensitivity and Variability for Amazonian Photosynthesis: Is It Too Much to Ask?

Abstract Estimates of Amazon rainforest gross primary productivity (GPP) differ by a factor of 2 across a suite of three statistical and 18 process models. This wide spread contributes uncertainty to predictions of future climate. We compare the mean and variance of GPP from these models to that of GPP at six eddy covariance (EC) towers. Only one model's mean GPP across all sites falls within a 99% confidence interval for EC GPP, and only one model matches EC variance. The strength of model response to climate drivers is related to model ability to match the seasonal pattern of the EC GPP. Models with stronger seasonal swings in GPP have stronger responses to rain, light, and temperature than does EC GPP. The model to data comparison illustrates a trade‐off inherent to deterministic models between accurate simulation of a mean (average) and accurate responsiveness to drivers. The trade‐off exists because all deterministic models simplify processes and lack at least some consequential driver or interaction. If a model's sensitivities to included drivers and their interactions are accurate, then deterministically predicted outcomes have less variability than is realistic. If a GPP model has stronger responses to climate drivers than found in data, model predictions may match the observed variance and seasonal pattern but are likely to overpredict GPP response to climate change. High or realistic variability of model estimates relative to reference data indicate that the model is hypersensitive to one or more drivers.

54 ENVIRONMENTAL SCIENCES↗

Gravitational spreading of high terrain in Ishtar Terra, Venus

A plane-strain viscoelastic finite element algorithm derived by Melosh and Rafesky (1980, 1983) is used to model the gravitational spreading of the high topography of Ishtar Terra. Brittle failure and relaxation of topographic relief are predicted as functions of time for ranges of crustal thickness and thermal gradient and for observed ranges in topographic elevation and slope. For a crustal thickness greater than 10 km, a thermal gradient of 15 K/km or more, and average values of relief (3 km) and slope (3 deg), the topography relaxes to 25 percent of its original height in less than 10 m.y. Although initial failure, in the form of normal faulting on the highlands and margin and occasional shallow thrusting in the lowlands, is predicted to occur much earlier, values of horizontal surface strain that are likely to be observable (about 1 percent) do not accumulate until significant relaxation of relief begins.

Smrekar, Suzanne E.↗

Comparison of Flow-Dependent and Static Error Correlation Models in the DAO Ozone Data Assimilation System

In a data assimilation system the forecast error covariance matrix governs the way in which the data information is spread throughout the model grid. Implementation of a correct method of assigning covariances is expected to have an impact on the analysis results. The simplest models assume that correlations are constant in time and isotropic or nearly isotropic. In such models the analysis depends on the dynamics only through assumed error standard deviations. In applications to atmospheric tracer data assimilation this may lead to inaccuracies, especially in regions with strong wind shears or high gradient of potential vorticity, as well as in areas where no data are available. In order to overcome this problem we have developed a flow-dependent covariance model that is based on short term evolution of error correlations. The presentation compares performance of a static and a flow-dependent model applied to a global three- dimensional ozone data assimilation system developed at NASA s Data Assimilation Office. We will present some results of validation against WMO balloon-borne sondes and the Polar Ozone and Aerosol Measurement (POAM) III instrument. Experiments show that allowing forecast error correlations to evolve with the flow results in positive impact on assimilated ozone within the regions where data were not assimilated, particularly at high latitudes in both hemispheres and in the troposphere. We will also discuss statistical characteristics of both models; in particular we will argue that including evolution of error correlations leads to stronger internal consistency of a data assimilation ,

Wargan, K.↗

Spatiotemporal spread of perturbations in power-law models at low temperatures: Exact results for classical out-of-time-order correlators

Here, we present exact results for the classical version of the out-of-time-order commutator (OTOC) for a family of power-law models consisting of N particles in one dimension and confined by an external harmonic potential. These particles are interacting via power-law interaction of the form ∝ $Σ_{i,j=(i≠j)}^{N}$ |x i -x j | -k for every k>1 where x i is the position of the i th particle. We present numerical results for the OTOC for finite N at low temperatures and short enough times so that the system is well approximated by the linearized dynamics around the many-body ground state. In the large-N limit, we compute the ground-state dispersion relation in the absence of external harmonic potential exactly and use it to arrive at analytical results for OTOC. We find excellent agreement between our analytical results and the numerics. We further obtain analytical results in the limit where only linear and leading nonlinear (in momentum) terms in the dispersion relation are included. The resulting OTOC is in agreement with numerics in the vicinity of the edge of the “light cone.” We find remarkably distinct features in OTOC below and above k=3 in terms of going from non-Airy behavior (1 < k < 3) to an Airy universality class (k>3). We present certain additional rich features for the case k=2 that stem from the underlying integrability of the Calogero-Moser model. We present a field theory approach that also assists in understanding certain aspects of OTOC such as the sound speed. Our findings are a step forward towards a more general understanding of the spatiotemporal spread of perturbations in long-range interacting systems.

1-dimensional systems↗

Centered CW interference rejection using spread spectrum techniques

A mathematical model is developed for the rejection of CW interference by spread spectrum techniques. When this interference is known to be exactly at the carrier frequency of the spread spectrum signal, this information can be used to design optimal IF filtering prior to despreading. The application of this approach to the pilot beam receiver of the Solar Power Satellite is considered as an example.

Scholtz, R. A.↗

The role of geography in the complex diffusion of innovations

The urban–rural divide is increasing in modern societies calling for geographical extensions of social influence modelling. Improved understanding of innovation diffusion across locations and through social connections can provide us with new insights into the spread of information, technological progress and economic development. In this work, we analyze the spatial adoption dynamics of iWiW, an Online Social Network (OSN) in Hungary and uncover empirical features about the spatial adoption in social networks. During its entire life cycle from 2002 to 2012, iWiW reached up to 300 million friendship ties of 3 million users. We find that the number of adopters as a function of town population follows a scaling law that reveals a strongly concentrated early adoption in large towns and a less concentrated late adoption. We also discover a strengthening distance decay of spread over the life-cycle indicating high fraction of distant diffusion in early stages but the dominance of local diffusion in late stages. The spreading process is modelled within the Bass diffusion framework that enables us to compare the differential equation version with an agent-based version of the model run on the empirical network. Although both model versions can capture the macro trend of adoption, they have limited capacity to describe the observed trends of urban scaling and distance decay. We find, however that incorporating adoption thresholds, defined by the fraction of social connections that adopt a technology before the individual adopts, improves the network model fit to the urban scaling of early adopters. Controlling for the threshold distribution enables us to eliminate the bias induced by local network structure on predicting local adoption peaks. Finally, we show that geographical features such as distance from the innovation origin and town size influence prediction of adoption peak at local scales in all model specifications.

97 MATHEMATICS AND COMPUTING↗

Comprehensive stellar population models and the disentanglement of age and metallicity effects

The construction of detailed models for intermediate and old stellar populations is described. Input parameters include metallicity (-2 less than (Fe/H) less than 0.5), single-burst age (between 1.5 and 17 Gyr), and initial mass function (IMF) exponent. Quantities output include broadband magnitudes, spectral energy distributions, surface brightness fluctuation magnitudes, and a suite of 21 absorption feature indices. The models are checked against a wide variety of available observations. Examinations of model output yield the following conclusions. (1) If the percentage change delta age/delta Z approximately equals 3/2 for two populations, they will appear almost identical in most indices. A few indices break this degeneracy by being either more abundance sensitive (Fe4668, Fe5015, Fe5709, and Fe5782) or more age sensitive (G4300, H beta, and presumably higher order Balmer lines) than usual. (2) Present uncertainties in stellar evolution are of the same magnitude as the effects of IMF and Y in the indices studied. (3) Changes in abundance ratios (like (Mg/Fe)) are predicted to be readily apparent in the spectra of old stellar populations. (4) The I-band flux of a stellar population is predicted to be nearly independent of metallicity and only modestly sensitive to age. The I band is therefore recommended for standard candle work or studies of M/L in galaxies. Other conclusions stem from this work. (1) Intercomparison of models and observations of two TiO indices seem to indicate variation of the (V/Ti) ratio among galaxies, but it is not clear how this observation ties into the standard picture of chemical enrichment. (2) Current estimates of (Fe/H) for the most metal-rich globulars that are based on integrated indices are probably slightly too high. (3) Colors of population models from different authors exhibit a substantial range. At solar metallicity and 13 Gyr, this range corresponds to an age error of roughly +/- 7 Gyr. Model colors from different authors applied in a differential sense have smaller uncertainties. (4) In the present models the dominant error for colors is probably the transformation from stellar atmospheric parameters to stellar colors. (5) Stellar B - V is difficult to model, and current spreads among different authors can reach 0.2 mag. (6) If known defects in the stellar flux library are corrected, the population model colors of this work in passbands redder than U would be accurate to roughly 0.03 mag in an absolute sense. These corrections are not made in the tables of model output.

Worthey, Guy↗

Flame Spread and Extinction in Partial-Gravity Environments

Considerable progress has been made in understanding the mechanisms of spreading flames under certain conditions, nearly all under the influence of normal Earth gravity. Recently, several investigators have studied some aspects of flame spread in purely forced flows in microgravity. However, very few have considered (especially experimentally) purely-buoyant flow influences, using gravity as a variable. In addition to the scientific interest in understanding how variable gravity affects flame spread in purely-buoyant flow, prospective human exploration of the Moon and Mars provides an incentive to obtain practical knowledge for use in fire-safety related engineering and mission operations in those partial-gravity environments. The purpose of this research effort is to conduct a focused experimental effort to observe the behavior of flames spreading both upward (concurrent flow) and downward (opposed flow) over thin fuels in partial-gravity environments, and to extend an existing numerical model of flame spread to predict flammability and flame spread behavior in these two regimes. A significant aspect of the experimental effort is to use a special device to improve the simulated partial-gravity environment achievable aboard reduced-gravity aircraft facilities.

K R Sacksteder↗

Formation of ferrobasalt at east Pacific midocean spreading centers

A numerical modeling using a least squares petrologic mixing program demonstrates that ridge ferrobasalts from three east Pacific regions are generated from normal midocean ridge balsalt by crystal fractionation of plagioclase, clinopyroxene, and minor olivine at shallow mantle depths. Clinopyroxene is an important fractionating phase in this model in contrast to most previously proposed fractionation models. Calculations demonstrate that the ferrobasalts represent a much smaller percentage of residual liquid than has previously been proposed. Implications of extensive shallow crystal fractionation for the determination of depths of magma generation are discussed. The different fields on an anorthite-olivine-hypersthene ternary diagram for basalts from east Pacific fast and moderate spreading ridges can be explained by variable degrees of shallow crystal fractionation rather than by different depths of magma generation as suggested by Bass (1971) and Scheidegger (1973).

Clague, D. A.↗

The Response of the Ozone Layer Under Abrupt 4xCO2 in CMIP6

Previous studies indicate a possible role of stratospheric ozone chemistry feedbacks in the climate response to 4xCO2, either via a reduction in equilibrium climate sensitivity (ECS) (Nowack et al., 2015), or via changes in the tropospheric circulation (Chiodo and Polvani, 2017). However, these effects are subject to uncertainty. Part of the uncertainty may stem from the dependency of the feedback on the pattern of the ozone response, as the radiative efficiency of ozone largely depends on its vertical distribution (Lacis et al., 1990). Here, an analysis is presented of the ozone layer response to 4xCO2 in chemistry–climate models (CCMs) which participated to CMIP intercomparisons. In a previous study using CMIP5 models, it has been shown that under 4xCO2, ozone decreases in the tropical lower stratosphere, and increases over the high latitudes and throughout the upper stratosphere. It was also found that a substantial portion of the spread in the tropical column ozone is tied to inter-model spread in tropical upwelling, which is in turn partly tied to ECS (Chiodo et al., 2018). Here, we revisit this connection using4xCO2 data from CMIP6, thereby exploiting the larger number of CCMs than in CMIP5. In addition, we explore the linearity of the ozone response, by complementing the analysis with simulations using transient CO2 forcing (1pc/year). We show that the pattern of the ozone response in is similar to CMIP5. In some models (e.g. WACCM), we find larger ozone responses in CMIP6 than in CMIP5, partly because of the larger ECS and thus larger upwelling response in the tropical pipe. In this presentation, we will discuss the relationship between radiative forcing, transport and ozone, as well as further implications for CMIP6 models.

Ozone layer↗

The Response of the Ozone Layer Under Abrupt 4xCO2 in CMIP6

Previous studies indicate a possible role of stratospheric ozone chemistry feedbacks in the climate response to 4xCO2, either via a reduction in equilibrium climate sensitivity (ECS) (Nowack et al.,2015), or via changes in the tropospheric circulation (Chiodo and Polvani, 2017). However, these effects are subject to uncertainty. Part of the uncertainty may stem from the dependency of the feedback on the pattern of the ozone response, as the radiative efficiency of ozone largely depends on its vertical distribution (Lacis et al., 1990). Here, an analysis is presented of the ozone layer response to 4xCO2 in chemistry–climate models (CCMs) which participated to CMIP intercomparisons. In a previous study using CMIP5 models, it has been shown that under 4xCO2, ozone decreases in the tropical lower stratosphere, and increases over the high latitudes and throughout the upper stratosphere. It was also found that a substantial portion of the spread in the tropical column ozone is tied to inter-model spread in tropical upwelling, which is in turn tied to ECS (Chiodo et al., 2018). Here, we revisit this connection using 4xCO2 data from CMIP6, thereby exploiting the larger number of CCMs than in CMIP5. In addition, we explore the linearity of the ozone response, by complementing the analysis with simulations using lower CO 2 forcing levels (2xCO2). We show that the pattern of the ozone response in is similar to CMIP5. In some models (e.g. WACCM), we find larger ozone responses in CMIP6 than in CMIP5, partly because of the larger ECS and thus larger upwelling response in the tropical pipe. In this presentation, we will discuss the relationship between radiative forcing, transport and ozone, as well as further implications for CMIP6 models.

Ozone layer↗