Engineering topics

Lee, Y. H.

Publications and source records attributed to Lee, Y. H..

At least 19 records

DOE OSTI2024

Overview of the KSTAR experiments toward fusion reactor

The Korean Superconducting Tokamak Advanced Research has been focused on exploring the key physics and engineering issues for future fusion reactors by demonstrating the long pulse operation of high beta steady-state discharge. Advanced scenarios are being developed with the goal for steady-state operation, and significant progress has been made in high ℓ i , hybrid and high beta scenarios with β N of 3. In the new operation scenario called fast ion regulated enhanced (FIRE), fast ions play an essential role in confinement enhancement. GK simulations show a significant reduction of the thermal energy flux when the thermal ion fraction decreases and the main ion density gradient is reversed by the fast ions in FIRE mode. Optimization of 3D magnetic field techniques, including adaptive control and real-time machine learning control algorithm, enabled long-pulse operation and high-performance ELM-suppressed discharge. Symmetric multiple shattered pellet injections (SPIs) and real-time disruption event characterization and forecasting are being performed to mitigate and avoid the disruptions associated with high-performance, long-pulse ITER-like scenarios. Finally, the near-term research plan will be addressed with the actively cooled tungsten divertor, a major upgrade of the NBI and helicon current drive heating, and transition to a full metallic wall.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Investigation of performance enhancement by balanced double-null shaping in KSTAR

We report experimental observations on the effect of plasma boundary shaping towards balanced double-null (DN) configuration on the plasma performance in KSTAR. The transition from a single-null to a DN configuration resulted in improved plasma performance, manifested through changes in the pedestal region, decreased density, and core MHD activity variation. Specifically, the DN transition led to a wider and higher pedestal structure, accompanied by grassy edge-localized modes (ELMs) characteristics. The density decrease was a prerequisite for performance enhancement during DN shaping, increasing fast ion confinement. Optimizing the plasma near the core region was associated with the suppression of sawtooth instabilities and the occurrence of fishbone modes during the DN transition. Integrated modeling demonstrated that secondary effects of the DN shaping could increase core thermal energy confinement.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

A sustained high-temperature fusion plasma regime facilitated by fast ions

Nuclear fusion is one of the most attractive alternatives to carbon-dependent energy sources. Harnessing energy from nuclear fusion in a large reactor scale, however, still presents many scientific challenges despite the many years of research and steady advances in magnetic confinement approaches. State-of-the-art magnetic fusion devices cannot yet achieve a sustainable fusion performance, which requires a high temperature above 100 million kelvin and sufficient control of instabilities to ensure steady-state operation on the order of tens of seconds. Here, in this study, we report experiments at the Korea Superconducting Tokamak Advanced Research device producing a plasma fusion regime that satisfies most of the above requirements: thanks to abundant fast ions stabilizing the core plasma turbulence, we generate plasmas at a temperature of 100 million kelvin lasting up to 20 seconds without plasma edge instabilities or impurity accumulation. A low plasma density combined with a moderate input power for operation is key to establishing this regime by preserving a high fraction of fast ions. This regime is rarely subject to disruption and can be sustained reliably even without a sophisticated control, and thus represents a promising path towards commercial fusion reactors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Exploration of RMP ELM control on ITER-similar shape (ISS) in KSTAR

Recent control advances at KSTAR enabled us to not only establish but also to sustain the ITER-similar shape (ISS) in a stationary manner, matching the ITER-like dimensionless shape parameters in lower single null configuration with vastly contrasting upper/lower triangularities. Here, the kinetic parameters of the ISS are different from the typical KSTAR configurations seen in ELM-RMP suppression experiments. We experimentally demonstrate the ISS-compatible resonant magnetic perturbation (RMP) control in KSTAR for the first time, using the n = 2, +90 degrees phasing RMP.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

On hybrid scenarios in KSTAR

We report the status of hybrid scenario experiments in Korea Superconducting Tokamak Advanced Research (KSTAR). The hybrid scenario is defined as stationary discharges with ${\beta _{\rm{N}}}$ $ \geqslant $ 2.4 and ${H_{89}}$ $ \geqslant $ 2.0 at ${q_{95}}$ < 6.5 without or with very mild sawtooth activities in KSTAR. It is being developed towards reactor-relevant conditions. High performance of ${\beta _{\rm{N}}} \lesssim 3.0$, ${H_{89}} \lesssim 2.4$ and G-factor ($ \equiv {\beta _{\rm{N}}}{H_{89}}/q_{95}^2$) $ \lesssim 0.4$6 has been achieved and sustained for $ \gtrsim 40{{\tau }_{\text{E}}}$ at ne/nGW ~ 0.7 with heating power of $ \lesssim $5 MW. Some KSTAR hybrid discharges exhibit a unique feature of a slow transition from conventional H-mode to hybrid mode after the third neutral beam injection. The reason for the confinement enhancement is extensively studied in this transition period of a representative discharge exhibiting a common feature of KSTAR hybrid scenarios. 0D performance analysis with magnetohydrodynamic activities, 1D kinetic profile dynamics, power balance analysis, linear gyro-kinetic analysis and edge pedestal stability analysis were conducted. The enhancement is thought to be from both the core and the pedestal. The improvement in the core region of the ion energy channel is observed from the linear gyro-kinetic analysis considering the electromagnetic, the fast ion, the Shafranov shift, ${\it{\omega }_{E \times B}}$, and the magnetic shear effect. The electromagnetic finite $\beta $ stabilisation plays a role in the inner core region at ${\rho _{\rm{tor}}}\sim0.35$ together with the fast ion effect. The alpha stabilisation effect is also found at ${\rho _{\rm{tor}}}\sim0.5$. ${\omega _{E \times B}}$, which could reduce the linear growth of the ion temperature gradient mode in the outer core region at ${\rho _{\rm{tor}}}\sim0.5 - 0.7$ with the highest contribution from the toroidal rotation. Furthermore, regarding the improvement in the pedestal, Shafranov shift broadens the stability boundary of the pedestal in support of the diamagnetic effect. The pedestal height and width could be reproduced by the EPED model, while a realistic current profile is used to calculate the internal inductance for Shafranov shift. Based on these findings, a comprehensive confinement enhancement mechanism has been proposed by considering the core-edge interplay.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS2016

Potential Impact of a US Climate Policy and Air Quality Regulations on Future Air Quality and Climate Change

We have investigated how future air quality and climate change are influenced by the US air quality regulations that existed or were proposed in 2013 and a hypothetical climate mitigation policy that aims to reduce 2050 CO2 emissions to be 50% below 2005 emissions. Using the NASA GISS ModelE2 general circulation model, we look at the impacts for year 2030 and 2055. The US energy-sector emissions are from the GLIMPSE project (GEOS-Chem LIDORT Integrated with MARKAL (MARKet ALlocation) for the Purpose of Scenario Exploration), and other US emissions data sets and the rest of the world emissions data sets are based on the RCP4.5 scenario. The US air quality regulations are projected to have a strong beneficial impact on US air quality and public health in year 2030 and 2055 but result in positive radiative forcing. Under this scenario, no more emission constraints are added after 2020, and the impacts on air quality and climate change are similar between year 2030 and 2055. Surface particulate matter with a diameter smaller than 2.5 micron PM(sub 2:5) is reduced by 2 approximately μg/m(sup -3) on average over the USA, and surface ozone by approximately 8 ppbv. The improved air quality prevents about 91 400 premature deaths in the USA, mainly due to the PM(sub 2:5) reduction approximately (74 200 lives saved). The air quality regulations reduce the light-reflecting aerosols (i.e., sulfate and organic matter) more than the light-absorbing species (i.e., black carbon and ozone), leading to a strong positive radiative forcing (RF) over the USA by both aerosols' direct and indirect forcing: the total RF is approximately 0.04 W m(sup -2) over the globe, and approximately 0.8 W m(sup -2) over the USA. Under the hypothetical climate policy, a future CO2 emissions cut is achieved in part by relying less on coal, and thus SO2 emissions are noticeably reduced. This provides air quality co-benefits, but it could lead to potential climate disbenefits over the USA. In 2055, the US mean total RF is +C02 W m(sup -2) due to positive aerosol direct and indirect forcing, while the global mean total RF is -0.06 W m(sup -2) due to the dominant negative CO2 RF (instantaneous RF). To achieve a regional-scale climate benefit via a climate policy, it is critical (1) to have multinational efforts to reduce greenhouse gas (GHG) emissions and (2) to simultaneously target emission reduction of light-absorbing species (e.g., BC and O3) on top of long-lived species. The latter is very desirable as the resulting climate benefit occurs faster and provides cobenefits to air quality and public health.

Air quality↗

NASA NTRS2016

Evaluation of Observed and Modelled Aerosol Lifetimes Using Radioactive Tracers of Opportunity and an Ensemble of 19 Global Models

Aerosols have important impacts on air quality and climate, but the processes affecting their removal from the atmosphere are not fully understood and are poorly constrained by observations. This makes modelled aerosol lifetimes uncertain. In this study, we make use of an observational constraint on aerosol lifetimes provided by radionuclide measurements and investigate the causes of differences within a set of global models. During the Fukushima Dai-Ichi nuclear power plant accident of March 2011, the radioactive isotopes cesium-137 (Cs-137) and xenon-133 (Xe-133) were released in large quantities. Cesium attached to particles in the ambient air, approximately according to their available aerosol surface area. Cs-137 size distribution measurements taken close to the power plant suggested that accumulation mode (AM) sulfate aerosols were the main carriers of cesium. Hence, Cs-137 can be used as a proxy tracer for the AM sulfate aerosol's fate in the atmosphere. In contrast, the noble gas Xe-133 behaves almost like a passive transport tracer. Global surface measurements of the two radioactive isotopes taken over several months after the release allow the derivation of a lifetime of the carrier aerosol. We compare this to the lifetimes simulated by 19 different atmospheric transport models initialized with identical emissions of Cs-137that were assigned to an aerosol tracer with each model's default properties of AM sulfate, and Xe-133 emissions that were assigned to a passive tracer. We investigate to what extent the modelled sulfate tracer can reproduce the measurements, especially with respect to the observed loss of aerosol mass with time. Modelled Cs-137and Xe-133 concentrations sampled at the same location and times as station measurements allow a direct comparison between measured and modelled aerosol lifetime. The e-folding lifetime e, calculated from station measurement data taken between 2 and 9 weeks after the start of the emissions, is 14.3 days (95% confidence interval 13.1-15.7 days). The equivalent modelled e lifetimes have a large spread, varying between 4.8 and 26.7 days with a model median of 9.42.3 days, indicating too fast a removal in most models. Because sufficient measurement data were only available from about 2 weeks after the release, the estimated lifetimes apply to aerosols that have undergone long-range transport, i.e. not for freshly emitted aerosol. However, modelled instantaneous lifetimes show that the initial removal in the first 2 weeks was quicker (lifetimes between 1 and 5 days) due to the emissions occurring at low altitudes and co-location of the fresh plume with strong precipitation. Deviations between measured and modelled aerosol lifetimes are largest for the northernmost stations and at later time periods, suggesting that models do not transport enough of the aerosol towards the Arctic. The models underestimate passive tracer (Xe-133) concentrations in the Arctic as well but to a smaller extent than for the aerosol (Cs-137) tracer. This indicates that in addition to too fast an aerosol removal in the models, errors in simulated atmospheric transport towards the Arctic in most models also contribute to the underestimation of the Arctic aerosol concentrations.

Radioactive isotopes↗

NASA NTRS2015

Evaluation of the Global Aerosol Microphysical ModelE2-TOMAS Model Against Satellite and Ground-Based Observations

The TwO-Moment Aerosol Sectional (TOMAS) microphysics model has been integrated into the state-of-the- art general circulation model, GISS ModelE2. This paper provides a detailed description of the ModelE2-TOMAS model and evaluates the model against various observations including aerosol precursor gas concentrations, aerosol mass and number concentrations, and aerosol optical depths. Additionally, global budgets in ModelE2-TOMAS are compared with those of other global aerosol models, and the ModelE2-TOMAS model is compared to the default aerosol model in ModelE2, which is a one-moment aerosol (OMA) model (i.e. no aerosol microphysics). Overall, the ModelE2- TOMAS predictions are within the range of other global aerosol model predictions, and the model has a reasonable agreement (mostly within a factor of 2) with observations of sulfur species and other aerosol components as well as aerosol optical depth. However, ModelE2-TOMAS (as well as ModelE2-OMA) cannot capture the observed vertical distribution of sulfur dioxide over the Pacific Ocean, possibly due to overly strong convective transport and over-predicted precipitation. The ModelE2-TOMAS model simulates observed aerosol number concentrations and cloud condensation nuclei concentrations roughly within a factor of 2. Anthropogenic aerosol burdens in ModelE2-OMA differ from ModelE2-TOMAS by a few percent to a factor of 2 regionally, mainly due to differences in aerosol processes including deposition, cloud processing, and emission parameterizations. We observed larger differences for naturally emitted aerosols such as sea salt and mineral dust, as those emission rates are quite different due to different upper size cutoff assumptions.

Lee, Y. H.↗

NASA NTRS2013

Multi-Model Mean Nitrogen and Sulfur Deposition from the Atmospheric Chemistry and Climate Model Intercomparison Project (ACCMIP): Evaluation of Historical and Projected Future Changes

We present multi-model global datasets of nitrogen and sulfate deposition covering time periods from 1850 to 2100, calculated within the Atmospheric Chemistry and Climate Model Intercomparison Project (ACCMIP). The computed deposition fluxes are compared to surface wet deposition and ice core measurements.We use a new dataset of wet deposition for 2000-2002 based on critical assessment of the quality of existing regional network data. We show that for present day (year 2000 ACCMIP time slice), the ACCMIP results perform similarly to previously published multi-model assessments. For this time slice, we find a multimodel mean deposition of approximately 50 Tg(N) yr−1 from nitrogen oxide emissions, 60 Tg(N) yr−1 from ammonia emissions, and 83 Tg(S) yr−1 from sulfur emissions. The analysis of changes between 1980 and 2000 indicates significant differences between model and measurements over the United States but less so over Europe. This difference points towards a potential misrepresentation of 1980 NH3 emissions over North America. Based on ice core records, the 1850 deposition fluxes agree well with Greenland ice cores, but the change between 1850 and 2000 seems to be overestimated in the Northern Hemisphere for both nitrogen and sulfur species. Using the Representative Concentration Pathways (RCPs) to define the projected climate and atmospheric chemistry related emissions and concentrations, we find large regional nitrogen deposition increases in 2100 in Latin America, Africa and parts of Asia under some of the scenarios considered. Increases in South Asia are especially large, and are seen in all scenarios, with 2100 values more than double their 2000 counterpart in some scenarios and reaching >1300 mg(N)m−2 yr−1 averaged over regional to continental-scale regions in RCP 2.6 and 8.5, 30-50% larger than the values in any region currently (circa 2000). However, sulfur deposition rates in 2100 are in all regions lower than in 2000 in all the RCPs. The new ACCMIP multimodel deposition dataset provides state-of-the-science, consistent and evaluated time slice (spanning 1850-2100) global gridded deposition fields for use in a wide range of climate and ecological studies.

Lamarque, J.-F.↗

NASA NTRS2013

Representation of Nucleation Mode Microphysics in a Global Aerosol Model with Sectional Microphysics

In models, nucleation mode (1 nm<Dp <10 nm) particle microphysics can be represented explicitly with aerosol microphysical processes or can be parameterized to obtain the growth and survival of nuclei to the model's lower size boundary. This study investigates how the representation of nucleation mode microphysics impacts aerosol number predictions in the TwO-Moment Aerosol Sectional (TOMAS) aerosol microphysics model running with the GISS GCM II-prime by varying its lowest diameter boundary: 1 nm, 3 nm, and 10 nm. The model with the 1 nm boundary simulates the nucleation mode particles with fully resolved microphysical processes, while the model with the 10 nm and 3 nm boundaries uses a nucleation mode dynamics parameterization to account for the growth of nucleated particles to 10 nm and 3 nm, respectively.We also investigate the impact of the time step for aerosol microphysical processes (a 10 min versus a 1 h time step) to aerosol number predictions in the TOMAS models with explicit dynamics for the nucleation mode particles (i.e., 3 nm and 1 nm boundary). The model with the explicit microphysics (i.e., 1 nm boundary) with the 10 min time step is used as a numerical benchmark simulation to estimate biases caused by varying the lower size cutoff and the time step. Different representations of the nucleation mode have a significant effect on the formation rate of particles larger than 10 nm from nucleated particles (J10) and the burdens and lifetimes of ultrafinemode (10 nm≤Dp ≤70 nm) particles but have less impact on the burdens and lifetimes of CCN-sized particles. The models using parameterized microphysics (i.e., 10 nm and 3 nm boundaries) result in higher J10 and shorter coagulation lifetimes of ultrafine-mode particles than the model with explicit dynamics (i.e., 1 nm boundary). The spatial distributions of CN10 (Dp ≥10 nm) and CCN(0.2 %) (i.e., CCN concentrations at 0.2%supersaturation) are moderately affected, especially CN10 predictions above ∼700 hPa where nucleation contributes most strongly to CN10 concentrations. The lowermost-layer CN10 is substantially improved with the 3 nm boundary (compared to 10 nm) in most areas. The overprediction in CN10 with the 3 nm and 10 nm boundaries can be explained by the overprediction of J10 or J3 with the parameterized microphysics, possibly due to the instantaneous growth rate assumption in the survival and growth parameterization. The errors in CN10 predictions are sensitive to the choice of the lower size boundary but not to the choice of the time step applied to the microphysical processes. The spatial distribution of CCN(0.2 %) with the 3 nm boundary is almost identical to that with the 1 nm boundary, but that with the 10 nm boundary can differ more than 10-40% in some areas. We found that the deviation in the 10 nm simulations is partly due to the longer time step (i.e., 1 h time step used in the 10 nm simulations compared to 10 min time step used in the benchmark simulations), but, even with the same time step, the 10 nm cutoff showed noticeably higher errors than the 3 nm cutoff. In conclusion, we generally recommend using a lower diameter boundary of 3 nm for studies focused on aerosol indirect effects but down to 1 nm boundary for studies focused on CN10 predictions or nucleation.

microphysics↗

NASA NTRS2013

Preindustrial to Present-Day Changes in Tropospheric Hydroxyl Radical and Methane Lifetime from the Atmospheric Chemistry and Climate Model Intercomparison Project (ACCMIP)

We have analysed time-slice simulations from 17 global models, participating in the Atmospheric Chemistry and Climate Model Intercomparison Project (ACCMIP), to explore changes in present-day (2000) hydroxyl radical (OH) concentration and methane (CH4) lifetime relative to preindustrial times (1850) and to 1980. A comparison of modeled and observation-derived methane and methyl chloroform lifetimes suggests that the present-day global multi-model mean OH concentration is overestimated by 5 to 10% but is within the range of uncertainties. The models consistently simulate higher OH concentrations in the Northern Hemisphere (NH) compared with the Southern Hemisphere (SH) for the present-day (2000; inter-hemispheric ratios of 1.13 to 1.42), in contrast to observation-based approaches which generally indicate higher OH in the SH although uncertainties are large. Evaluation of simulated carbon monoxide (CO) concentrations, the primary sink for OH, against ground-based and satellite observations suggests low biases in the NH that may contribute to the high north–south OH asymmetry in the models. The models vary widely in their regional distribution of present-day OH concentrations (up to 34%). Despite large regional changes, the multi-model global mean (mass-weighted) OH concentration changes little over the past 150 yr, due to concurrent increases in factors that enhance OH (humidity, tropospheric ozone, nitrogen oxide (NOx) emissions, and UV radiation due to decreases in stratospheric ozone), compensated by increases in OH sinks (methane abundance, carbon monoxide and non-methane volatile organic carbon (NMVOC) emissions). The large inter-model diversity in the sign and magnitude of preindustrial to present-day OH changes (ranging from a decrease of 12.7% to an increase of 14.6%) indicate that uncertainty remains in our understanding of the long-term trends in OH and methane lifetime. We show that this diversity is largely explained by the different ratio of the change in global mean tropospheric CO and NOx burdens (Delta CO/Delta NOx, approximately represents changes in OH sinks versus changes in OH sources) in the models, pointing to a need for better constraints on natural precursor emissions and on the chemical mechanisms in the current generation of chemistry-climate models. For the 1980 to 2000 period, we find that climate warming and a slight increase in mean OH (3.5 +/- 2.2%) leads to a 4.3 +/- 1.9% decrease in the methane lifetime. Analysing sensitivity simulations performed by 10 models, we find that preindustrial to present-day climate change decreased the methane lifetime by about four months, representing a negative feedback on the climate system. Further, we analysed attribution experiments performed by a subset of models relative to 2000 conditions with only one precursor at a time set to 1860 levels. We find that global mean OH increased by 46.4 +/- 12.2% in response to preindustrial to present-day anthropogenic NOx emission increases, and decreased by 17.3 +/-2.3%, 7.6 +/- 1.5%, and 3.1 +/- 3.0% due to methane burden, and anthropogenic CO, and NMVOC emissions increases, respectively.

nitrogen oxides↗

NASA NTRS2013

A 4-D Climatology (1979-2009) of the Monthly Tropospheric Aerosol Optical Depth Distribution over the Mediterranean Region from a Comparative Evaluation and Blending of Remote Sensing and Model Products

Since the 1980s several spaceborne sensors have been used to retrieve the aerosol optical depth (AOD) over the Mediterranean region. In parallel, AOD climatologies coming from different numerical model simulations are now also available, permitting to distinguish the contribution of several aerosol types to the total AOD. In this work, we perform a comparative analysis of this unique multiyear database in terms of total AOD and of its apportionment by the five main aerosol types (soil dust, seasalt, sulfate, black and organic carbon). We use 9 different satellite-derived monthly AOD products: NOAA/AVHRR, SeaWiFS (2 products), TERRA/MISR, TERRA/MODIS, AQUA/MODIS, ENVISAT/MERIS, PARASOL/POLDER and MSG/SEVIRI, as well as 3 more historical datasets: NIMBUS7/CZCS, TOMS (onboard NIMBUS7 and Earth- Probe) and METEOSAT/MVIRI. Monthly model datasets include the aerosol climatology from Tegen et al. (1997), the climate-chemistry models LMDz-OR-INCA and RegCM-4, the multi-model mean coming from the ACCMIP exercise, and the reanalyses GEMS and MACC. Ground-based Level- 2 AERONET AOD observations from 47 stations around the basin are used here to evaluate the model and satellite data. The sensor MODIS (on AQUA and TERRA) has the best average AOD scores over this region, showing a relevant spatio-temporal variability and highlighting high dust loads over Northern Africa and the sea (spring and summer), and sulfate aerosols over continental Europe (summer). The comparison also shows limitations of certain datasets (especially MERIS and SeaWiFS standard products). Models reproduce the main patterns of the AOD variability over the basin. The MACC reanalysis is the closest to AERONET data, but appears to underestimate dust over Northern Africa, where RegCM-4 is found closer to MODIS thanks to its interactive scheme for dust emissions. The vertical dimension is also investigated using the CALIOP instrument. This study confirms differences of vertical distribution between dust aerosols showing a large vertical spread, and other continental and marine aerosols which are confined in the boundary layer. From this compilation, we propose a 4-D blended product from model and satellite data, consisting in monthly time series of 3-D aerosol distribution at a 50 km horizontal resolution over the Euro-Mediterranean marine and continental region for the 2003-2009 period. The product is based on the total AOD from AQUA/MODIS, apportioned into sulfates, black and organic carbon from the MACC reanalysis, and into dust and sea-salt aerosols from RegCM-4 simulations, which are distributed vertically based on CALIOP climatology.We extend the 2003-2009 reconstruction to the past up to 1979 using the 2003-2009 average and applying the decreasing trend in sulfate aerosols from LMDz-OR-INCA, whose AOD trends over Europe and the Mediterranean are median among the ACCMIP models. Finally optical properties of the different aerosol types in this region are proposed from Mie calculations so that this reconstruction can be included in regional climate models for aerosol radiative forcing and aerosol-climate studies.

troposphere↗

NASA NTRS2013

Evaluation of ACCMIP Outgoing Longwave Radiation from Tropospheric Ozone Using TES Satellite Observations.

We use simultaneous observations of tropospheric ozone and outgoing longwave radiation (OLR) sensitivity to tropospheric ozone from the Tropospheric Emission Spectrometer (TES) to evaluate model tropospheric ozone and its effect on OLR simulated by a suite of chemistry-climate models that participated in the Atmospheric Chemistry and Climate Model Intercomparison Project (ACCMIP). The ensemble mean of ACCMIP models show a persistent but modest tropospheric ozone low bias (5-20 ppb) in the Southern Hemisphere (SH) and modest high bias (5-10 ppb) in the Northern Hemisphere (NH) relative to TES ozone for 2005-2010. These ozone biases have a significant impact on the OLR. Using TES instantaneous radiative kernels (IRK), we show that the ACCMIP ensemble mean tropospheric ozone low bias leads up to 120mW/ sq. m OLR high bias locally but zonally compensating errors reduce the global OLR high bias to 39+/- 41mW/ sq. m relative to TES data. We show that there is a correlation (Sq. R = 0.59) between the magnitude of the ACCMIP OLR bias and the deviation of the ACCMIP preindustrial to present day (1750-2010) ozone radiative forcing (RF) from the ensemble ozone RF mean. However, this correlation is driven primarily by models whose absolute OLR bias from tropospheric ozone exceeds 100mW/ sq. m. Removing these models leads to a mean ozone radiative forcing of 394+/- 42mW/ sq. m. The mean is about the same and the standard deviation is about 30% lower than an ensemble ozone RF of 384 +/- 60mW/ sq. m derived from 14 of the 16 ACCMIP models reported in a companion ACCMIP study. These results point towards a profitable direction of combining satellite observations and chemistry-climate model simulations to reduce uncertainty in ozone radiative forcing.

troposphere↗

NASA NTRS2013

Radiative Forcing in the ACCMIP Historical and Future Climate Simulations

A primary goal of the Atmospheric Chemistry and Climate Model IntercomparisonProject (ACCMIP) was to characterize the short-lived drivers of preindustrial to 2100climate change in the current generation of climate models. Here we evaluate historicaland 5 future radiative forcing in the 10 ACCMIP models that included aerosols, 8 of whichalso participated in the Coupled Model Intercomparison Project phase 5 (CMIP5).The models generally reproduce present-day climatological total aerosol opticaldepth (AOD) relatively well. components to this total, however, and most appear to underestimate AOD over East10 Asia. The models generally capture 1980-2000 AOD trends fairly well, though theyunderpredict AOD increases over the YellowEastern Sea. They appear to strongly underestimate absorbing AOD, especially in East Asia, South and Southeast Asia, SouthAmerica and Southern Hemisphere Africa.We examined both the conventional direct radiative forcing at the tropopause (RF) and the forcing including rapid adjustments (adjusted forcing AF, including direct andindirect effects). The models calculated all aerosol all-sky 1850 to 2000 global meanannual average RF ranges from 0.06 to 0.49 W m(sup -2), with a mean of 0.26 W m(sup -2) and a median of 0.27 W m(sup -2. Adjusting for missing aerosol components in some modelsbrings the range to 0.12 to 0.62W m(sup -2), with a mean of 0.39W m(sup -2). Screen20ing the models based on their ability to capture spatial patterns and magnitudes ofAOD and AOD trends yields a quality-controlled mean of 0.42W m(sup -2) and range of0.33 to 0.50 W m(sup -2) (accounting for missing components). The CMIP5 subset of ACCMIPmodels spans 0.06 to 0.49W m(sup -2), suggesting some CMIP5 simulations likelyhave too little aerosol RF. A substantial, but not well quantified, contribution to histori25cal aerosol RF may come from climate feedbacks (35 to 58). The mean aerosol AF during this period is 1.12W m(sup -2) (median value 1.16W m(sup -2), range 0.72 to1.44W m(sup -2), indicating that adjustments to aerosols, which include cloud, water vaporand temperature, lead to stronger forcing than the aerosol direct RF.

Climate Models↗

NASA NTRS2013

Tropospheric Ozone Changes, Radiative Forcing and Attribution to Emissions in the Atmospheric Chemistry and Climate Model Intercomparison Project (ACCMIP)

Ozone (O3) from 17 atmospheric chemistry models taking part in the Atmospheric Chemistry and Climate Model Intercomparison Project (ACCMIP) has been used to calculate tropospheric ozone radiative forcings (RFs). All models applied a common set of anthropogenic emissions, which are better constrained for the present-day than the past. Future anthropogenic emissions follow the four Representative Concentration Pathway (RCP) scenarios, which define a relatively narrow range of possible air pollution emissions. We calculate a value for the pre-industrial (1750) to present-day (2010) tropospheric ozone RF of 410 mW m−2. The model range of pre-industrial to present-day changes in O3 produces a spread (+/-1 standard deviation) in RFs of +/-17%. Three different radiation schemes were used - we find differences in RFs between schemes (for the same ozone fields) of +/-10 percent. Applying two different tropopause definitions gives differences in RFs of +/-3 percent. Given additional (unquantified) uncertainties associated with emissions, climate-chemistry interactions and land-use change, we estimate an overall uncertainty of +/-30 percent for the tropospheric ozone RF. Experiments carried out by a subset of six models attribute tropospheric ozone RF to increased emissions of methane (44+/-12 percent), nitrogen oxides (31 +/- 9 percent), carbon monoxide (15 +/- 3 percent) and non-methane volatile organic compounds (9 +/- 2 percent); earlier studies attributed more of the tropospheric ozone RF to methane and less to nitrogen oxides. Normalising RFs to changes in tropospheric column ozone, we find a global mean normalised RF of 42 mW m(−2) DU(−1), a value similar to previous work. Using normalised RFs and future tropospheric column ozone projections we calculate future tropospheric ozone RFs (mW m(−2); relative to 1750) for the four future scenarios (RCP2.6, RCP4.5, RCP6.0 and RCP8.5) of 350, 420, 370 and 460 (in 2030), and 200, 300, 280 and 600 (in 2100). Models show some coherent responses of ozone to climate change: decreases in the tropical lower troposphere, associated with increases in water vapour; and increases in the sub-tropical to mid-latitude upper troposphere, associated with increases in lightning and stratosphere-to-troposphere transport. Climate change has relatively small impacts on global mean tropospheric ozone RF.

climate change↗

NASA NTRS2013

Evaluation of preindustrial to present-day black carbon and its albedo forcing from Atmospheric Chemistry and Climate Model Intercomparison Project (ACCMIP)

As part of the Atmospheric Chemistry and Climate Model Intercomparison Project (ACCMIP), we evaluate the historical black carbon (BC) aerosols simulated by 8 ACCMIP models against observations including 12 ice core records, long-term surface mass concentrations, and recent Arctic BC snowpack measurements. We also estimate BC albedo forcing by performing additional simulations using offline models with prescribed meteorology from 1996-2000. We evaluate the vertical profile of BC snow concentrations from these offline simulations using the recent BC snowpack measurements. Despite using the same BC emissions, the global BC burden differs by approximately a factor of 3 among models due to differences in aerosol removal parameterizations and simulated meteorology: 34 Gg to 103 Gg in 1850 and 82 Gg to 315 Gg in 2000. However, the global BC burden from preindustrial to present-day increases by 2.5-3 times with little variation among models, roughly matching the 2.5-fold increase in total BC emissions during the same period.We find a large divergence among models at both Northern Hemisphere (NH) and Southern Hemisphere (SH) high latitude regions for BC burden and at SH high latitude regions for deposition fluxes. The ACCMIP simulations match the observed BC surface mass concentrations well in Europe and North America except at Ispra. However, the models fail to predict the Arctic BC seasonality due to severe underestimations during winter and spring. The simulated vertically resolved BC snow concentrations are, on average, within a factor of 2-3 of the BC snowpack measurements except for Greenland and the Arctic Ocean. For the ice core evaluation, models tend to adequately capture both the observed temporal trends and the magnitudes at Greenland sites. However, models fail to predict the decreasing trend of BC depositions/ice core concentrations from the 1950s to the 1970s in most Tibetan Plateau ice cores. The distinct temporal trend at the Tibetan Plateau ice cores indicates a strong influence from Western Europe, but the modeled BC increases in that period are consistent with the emission changes in Eastern Europe, the Middle East, South and East Asia. At the Alps site, the simulated BC suggests a strong influence from Europe, which agrees with the Alps ice core observations. At Zuoqiupu on the Tibetan Plateau, models successfully simulate the higher BC concentrations observed during the non-monsoon season compared to the monsoon season but overpredict BC in both seasons. Despite a large divergence in BC deposition at two Antarctic ice core sites, some models with a BC lifetime of less than 7 days are able to capture the observed concentrations. In 2000 relative to 1850, globally and annually averaged BC surface albedo forcing from the offline simulations ranges from 0.014 to 0.019Wm−2 among the ACCMIP models. Comparing offline and online BC albedo forcings computed by some of the same models, we find that the global annual mean can vary by up to a factor of two because of different aerosol models or different BC-snow parameterizations and snow cover. The spatial distributions of the offline BC albedo forcing in 2000 show especially high BC forcing (i.e., over 0.1W/sq. m) over Manchuria, Karakoram, and most of the Former USSR. Models predict the highest global annual mean BC forcing in 1980 rather than 2000, mostly driven by the high fossil fuel and biofuel emissions in the Former USSR in 1980.

snow cover↗

NASA NTRS2013

Pre-industrial to End 21st Century Projections of Tropospheric Ozone from the Atmospheric Chemistry and Climate Model Intercomparison Project (ACCMIP)

Present day tropospheric ozone and its changes between 1850 and 2100 are considered, analysing 15 global models that participated in the Atmospheric Chemistry and Climate Model Intercomparison Project (ACCMIP). The ensemble mean compares well against present day observations. The seasonal cycle correlates well, except for some locations in the tropical upper troposphere. Most (75 %) of the models are encompassed with a range of global mean tropospheric ozone column estimates from satellite data, but there is a suggestion of a high bias in the Northern Hemisphere and a low bias in the Southern Hemisphere, which could indicate deficiencies with the ozone precursor emissions. Compared to the present day ensemble mean tropospheric ozone burden of 337+/-23 Tg, the ensemble mean burden for 1850 time slice is approx. 30% lower. Future changes were modelled using emissions and climate projections from four Representative Concentration Pathways (RCPs). Compared to 2000, the relative changes in the ensemble mean tropospheric ozone burden in 2030 (2100) for the different RCPs are: −4% (−16 %) for RCP2.6, 2% (−7%) for RCP4.5, 1% (−9%) for RCP6.0, and 7% (18 %) for RCP8.5. Model agreement on the magnitude of the change is greatest for larger changes. Reductions in most precursor emissions are common across the RCPs and drive ozone decreases in all but RCP8.5, where doubled methane and a 40-150% greater stratospheric influx (estimated from a subset of models) increase ozone. While models with a high ozone burden for the present day also have high ozone burdens for the other time slices, no model consistently predicts large or small ozone changes; i.e. the magnitudes of the burdens and burden changes do not appear to be related simply, and the models are sensitive to emissions and climate changes in different ways. Spatial patterns of ozone changes are well correlated across most models, but are notably different for models without time evolving stratospheric ozone concentrations. A unified approach to ozone budget specifications and a rigorous investigation of the factors that drive tropospheric ozone is recommended to help future studies attribute ozone changes and inter-model differences more clearly.

Troposphere↗

NASA NTRS2013

The Atmospheric Chemistry and Climate Model Intercomparison Project (ACCMIP): Overview and Description of Models, Simulations and Climate Diagnostics

The Atmospheric Chemistry and Climate Model Intercomparison Project (ACCMIP) consists of a series of time slice experiments targeting the long-term changes in atmospheric composition between 1850 and 2100, with the goal of documenting composition changes and the associated radiative forcing. In this overview paper, we introduce the ACCMIP activity, the various simulations performed (with a requested set of 14) and the associated model output. The 16 ACCMIP models have a wide range of horizontal and vertical resolutions, vertical extent, chemistry schemes and interaction with radiation and clouds. While anthropogenic and biomass burning emissions were specified for all time slices in the ACCMIP protocol, it is found that the natural emissions are responsible for a significant range across models, mostly in the case of ozone precursors. The analysis of selected present-day climate diagnostics (precipitation, temperature, specific humidity and zonal wind) reveals biases consistent with state-of-the-art climate models. The model-to- model comparison of changes in temperature, specific humidity and zonal wind between 1850 and 2000 and between 2000 and 2100 indicates mostly consistent results. However, models that are clear outliers are different enough from the other models to significantly affect their simulation of atmospheric chemistry.

humidity↗