This is the first part of a series of two articles analyzing the global thermal properties of atmosphere–ocean coupled general circulation models (AOGCMs) within the framework of a two-layer ...energy-balance model (EBM). In this part, the general analytical solution of the system is given and two idealized climate change scenarios, one with a step forcing and one with a linear forcing, are discussed. These solutions give a didactic description of the contributions from the equilibrium response and of the fast and slow transient responses during a climate transition. Based on these analytical solutions, a simple and physically based procedure to calibrate the two-layer model parameters using an AOGCM step-forcing experiment is introduced. Using this procedure, the global thermal properties of 16 AOGCMs participating in phase 5 of the Coupled Model Intercomparison Project (CMIP5) are determined. It is shown that, for a given AOGCM, the EBM tuned with only the abrupt 4×CO₂ experiment is able to reproduce with a very good accuracy the temperature evolution in both a step-forcing and a linear-forcing experiment. The role of the upper-ocean and deep-ocean heat uptakes in the fast and slow responses is also discussed. One of the main weaknesses of the simple EBM discussed in this part is its ability to represent the evolution of the top-of-the-atmosphere radiative imbalance in the transient regime. This issue is addressed in Part II by taking into account the efficacy factor of deep-ocean heat uptake.
Celotno besedilo
Dostopno za:
BFBNIB, DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
Policies to control air quality focus on mitigating emissions of aerosols and their precursors, and other short-lived climate pollutants (SLCPs). On a local scale, these policies will have beneficial ...impacts on health and crop yields, by reducing particulate matter (PM) and surface ozone concentrations; however, the climate impacts of reducing emissions of SLCPs are less straightforward to predict. In this paper we consider a set of idealized, extreme mitigation strategies, in which the total anthropogenic emissions of individual SLCP emissions species are removed. This provides an upper bound on the potential climate impacts of such air quality strategies. We focus on evaluating the climate responses to changes in anthropogenic emissions of aerosol precursor species: black carbon (BC), organic carbon (OC) and sulphur dioxide (SO2). We perform climate integrations with four fully coupled atmosphere-ocean global climate models (AOGCMs), and examine the effects on global and regional climate of removing the total land-based anthropogenic emissions of each of the three aerosol precursor species. We find that the SO2 emissions reductions lead to the strongest response, with all models showing an increase in surface temperature focussed in the Northern Hemisphere mid and (especially) high latitudes, and showing a corresponding increase in global mean precipitation. Changes in precipitation patterns are driven mostly by a northward shift in the ITCZ (Intertropical Convergence Zone), consistent with the hemispherically asymmetric warming pattern driven by the emissions changes. The BC and OC emissions reductions give a much weaker response, and there is some disagreement between models in the sign of the climate responses to these perturbations. These differences between models are due largely to natural variability in sea-ice extent, circulation patterns and cloud changes. This large natural variability component to the signal when the ocean circulation and sea-ice are free-running means that the BC and OC mitigation measures do not necessarily lead to a discernible climate response.
Within the framework of the AeroCom (Aerosol Comparisons between Observations and Models) initiative, the state-of-the-art modelling of aerosol optical properties is assessed from 14 global models ...participating in the phase III control experiment (AP3). The models are similar to CMIP6/AerChemMIP Earth System Models (ESMs) and provide a robust multi-model ensemble. Inter-model spread of aerosol species lifetimes and emissions appears to be similar to that of mass extinction coefficients (MECs), suggesting that aerosol optical depth (AOD) uncertainties are associated with a broad spectrum of parameterised aerosol processes.
Total AOD is approximately the same as in AeroCom phase I (AP1) simulations. However, we find a 50 % decrease in the optical depth (OD) of black carbon (BC), attributable to a combination of decreased emissions and lifetimes. Relative contributions from sea salt (SS) and dust (DU) have shifted from being approximately equal in AP1 to SS contributing about 2∕3 of the natural AOD in AP3. This shift is linked with a decrease in DU mass burden, a lower DU MEC, and a slight decrease in DU lifetime, suggesting coarser DU particle sizes in AP3 compared to AP1.
Relative to observations, the AP3 ensemble median and most of the participating models underestimate all aerosol optical properties investigated, that is, total AOD as well as fine and coarse AOD (AODf, AODc), Ångström exponent (AE), dry surface scattering (SCdry), and absorption (ACdry) coefficients. Compared to AERONET, the models underestimate total AOD by ca. 21 % ± 20 % (as inferred from the ensemble median and interquartile range). Against satellite data, the ensemble AOD biases range from −37 % (MODIS-Terra) to −16 % (MERGED-FMI, a multi-satellite AOD product), which we explain by differences between individual satellites and AERONET measurements themselves. Correlation coefficients (R) between model and observation AOD records are generally high (R>0.75), suggesting that the models are capable of capturing spatio-temporal variations in AOD. We find a much larger underestimate in coarse AODc (∼ −45 % ± 25 %) than in fine AODf (∼ −15 % ± 25 %) with slightly increased inter-model spread compared to total AOD. These results indicate problems in the modelling of DU and SS. The AODc bias is likely due to missing DU over continental land masses (particularly over the United States, SE Asia, and S. America), while marine AERONET sites and the AATSR SU satellite data suggest more moderate oceanic biases in AODc.
Column AEs are underestimated by about 10 % ± 16 %. For situations in which measurements show AE > 2, models underestimate AERONET AE by ca. 35 %. In contrast, all models (but one) exhibit large overestimates in AE when coarse aerosol dominates (bias ca. +140 % if observed AE < 0.5). Simulated AE does not span the observed AE variability. These results indicate that models overestimate particle size (or underestimate the fine-mode fraction) for fine-dominated aerosol and underestimate size (or overestimate the fine-mode fraction) for coarse-dominated aerosol. This must have implications for lifetime, water uptake, scattering enhancement, and the aerosol radiative effect, which we can not quantify at this moment.
Comparison against Global Atmosphere Watch (GAW) in situ data results in mean bias and inter-model variations of −35 % ± 25 % and −20 % ± 18 % for SCdry and ACdry, respectively. The larger underestimate of SCdry than ACdry suggests the models will simulate an aerosol single scattering albedo that is too low. The larger underestimate of SCdry than ambient air AOD is consistent with recent findings that models overestimate scattering enhancement due to hygroscopic growth. The broadly consistent negative bias in AOD and surface scattering suggests an underestimate of aerosol radiative effects in current global aerosol models.
Considerable inter-model diversity in the simulated optical properties is often found in regions that are, unfortunately, not or only sparsely covered by ground-based observations. This includes, for instance, the Sahara, Amazonia, central Australia, and the South Pacific. This highlights the need for a better site coverage in the observations, which would enable us to better assess the models, but also the performance of satellite products in these regions.
Using fine-mode AOD as a proxy for present-day aerosol forcing estimates, our results suggest that models underestimate aerosol forcing by ca. −15 %, however, with a considerably large interquartile range, suggesting a spread between −35 % and +10 %.
Characteristics and trends of surface ocean dimethylsulfide (DMS) concentrations and fluxes into the atmosphere of four Earth system models (ESMs:
CNRM-ESM2-1, MIROC-ES2L, NorESM2-LM, and ...UKESM1-0-LL) are analysed over the recent past (1980–2009) and into the future, using Coupled Model
Intercomparison Project 6 (CMIP6) simulations. The DMS concentrations in historical simulations systematically underestimate the most widely used
observed climatology but compare more favourably against two recent observation-based datasets. The models better reproduce observations in mid to
high latitudes, as well as in polar and westerlies marine biomes. The resulting multi-model estimate of contemporary global ocean DMS emissions is
16–24 Tg S yr−1, which is narrower than the observational-derived range of 16 to 28 Tg S yr−1. The four models
disagree on the sign of the trend of the global DMS flux from 1980 onwards, with two models showing an increase and two models a decrease. At the
global scale, these trends are dominated by changes in surface DMS concentrations in all models, irrespective of the air–sea flux parameterisation
used. In turn, three models consistently show that changes in DMS concentrations are correlated with changes in marine productivity; however, marine productivity is poorly constrained in the current generation of ESMs, thus limiting the predictive ability of this relationship. In contrast, a consensus
is found among all models over polar latitudes where an increasing trend is predominantly driven by the retreating sea-ice extent. However, the
magnitude of this trend between models differs by a factor of 3, from 2.9 to 9.2 Gg S decade−1 over the period 1980–2014, which is
at the low end of a recent satellite-derived analysis. Similar increasing trends are found in climate projections over the 21st century.
In this second part of a series of two articles analyzing the global thermal properties of atmosphere–ocean coupled general circulation models (AOGCMs) within the framework of a two-layer ...energy-balance model (EBM), the role of the efficacy of deep-ocean heat uptake is investigated. Taking into account such an efficacy factor is shown to amount to representing the effect of deep-ocean heat uptake on the local strength of the radiative feedback in the transient regime. It involves an additional term in the formulation of the radiative imbalance at the top of the atmosphere (TOA), which explains the nonlinearity between radiative imbalance and the mean surface temperature observed in some AOGCMs. An analytical solution of this system is given and this simple linear EBM is calibrated for the set of 16 AOGCMs of phase 5 of the Coupled Model Intercomparison Project (CMIP5) studied in Part I. It is shown that both the net radiative fluxes at TOA and the global surface temperature transient response are well represented by the simple EBM over the available period of simulations. Differences between this two-layer EBM and the previous version without an efficacy factor are analyzed and relationships between parameters are discussed. The simple model calibration applied to AOGCMs constitutes a new method for estimating their respective equilibrium climate sensitivity and adjusted radiative forcing amplitude from short-term step-forcing simulations and more generally a method to compute their global thermal properties.
Celotno besedilo
Dostopno za:
BFBNIB, DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
Aerosol-induced absorption of shortwave radiation can modify the climate through local atmospheric heating, which affects lapse rates, precipitation, and cloud formation. Presently, the total amount ...of aerosol absorption is poorly constrained, and the main absorbing aerosol species (black carbon (BC), organic aerosols (OA), and mineral dust) are diversely quantified in global climate models. As part of the third phase of the Aerosol Comparisons between Observations and Models (AeroCom) intercomparison initiative (AeroCom phase III), we here document the distribution and magnitude of aerosol absorption in current global aerosol models and quantify the sources of intermodel spread, highlighting the difficulties of attributing absorption to different species. In total, 15 models have provided total present-day absorption at 550 nm (using year 2010 emissions), 11 of which have provided absorption per absorbing species. The multi-model global annual mean total absorption aerosol optical depth (AAOD) is 0.0054 (0.0020 to 0.0098; 550 nm), with the range given as the minimum and maximum model values. This is 28 % higher compared to the 0.0042 (0.0021 to 0.0076) multi-model mean in AeroCom phase II (using year 2000 emissions), but the difference is within 1 standard deviation, which, in this study, is 0.0023 (0.0019 in Phase II). Of the summed component AAOD, 60 % (range 36 %–84 %) is estimated to be due to BC, 31 % (12 %–49 %) is due to dust, and 11 % (0 %–24 %) is due to OA; however, the components are not independent in terms of their absorbing efficiency. In models with internal mixtures of absorbing aerosols, a major challenge is the lack of a common and simple method to attribute absorption to the different absorbing species. Therefore, when possible, the models with internally mixed aerosols in the present study have performed simulations using the same method for estimating absorption due to BC, OA, and dust, namely by removing it and comparing runs with and without the absorbing species. We discuss the challenges of attributing absorption to different species; we compare burden, refractive indices, and density; and we contrast models with internal mixing to models with external mixing. The model mean BC mass absorption coefficient (MAC) value is 10.1 (3.1 to 17.7) m2 g−1 (550 nm), and the model mean BC AAOD is 0.0030 (0.0007 to 0.0077). The difference in lifetime (and burden) in the models explains as much of the BC AAOD spread as the difference in BC MAC values. The difference in the spectral dependency between the models is striking. Several models have an absorption Ångstrøm exponent (AAE) close to 1, which likely is too low given current knowledge of spectral aerosol optical properties. Most models do not account for brown carbon and underestimate the spectral dependency for OA.
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 (137Cs) and xenon-133 (133Xe) were released in large quantities. Cesium attached to particles in the ambient air, approximately according to their available aerosol surface area. 137Cs size distribution measurements taken close to the power plant suggested that accumulation-mode (AM) sulfate aerosols were the main carriers of cesium. Hence, 137Cs can be used as a proxy tracer for the AM sulfate aerosol's fate in the atmosphere. In contrast, the noble gas 133Xe 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 137Cs that were assigned to an aerosol tracer with each model's default properties of AM sulfate, and 133Xe 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 137Cs and 133Xe 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.4 ± 2.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 (133Xe) concentrations in the Arctic as well but to a smaller extent than for the aerosol (137Cs) 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.
For the period 1860–2100 (SRES scenario A1B for 2000–2100), the impact of road transport, maritime shipping and aviation on climate is studied using an Atmosphere Ocean General Circulation Model ...(AOGCM). In addition to carbon dioxide (CO2) emissions from these transport sectors, most of their non-CO2 emissions are also taken into account, i.e. the forcing from ozone, methane, black carbon, organic carbon, sulfate, CFC-12 and HFC-134a from air conditioning systems in cars, and contrails. For the year 2000, the CO2 emissions from all sectors together induce a global annual-mean surface air temperature increase of around 0.1 K. In 2100, the CO2 emissions from road transport induce a global mean warming of 0.3 K, while shipping and aviation each contribute 0.1 K. For road transport, the non-CO2 impact is largest between 2000 and 2050 (of the order of 0.1 K) becoming smaller at the end of the 21st century. The non-CO2 impact from shipping is negative, reaching −0.1 K between 2050 and 2100, while for aviation it is positive and its estimate varies between 0 and 0.15 K in 2100. The largest changes in sea-level from thermal expansion in 2000 are 1.6 mm for the CO2 emissions from road transport, and around −3 mm from the non-CO2 effects of shipping. In 2100, sea-level rises by 18 mm due to the CO2 emissions from road transport and by 4.6 mm due to shipping or aviation CO2 emissions. Non-CO2 changes are of the order of 1 mm for road transport, −6.6 mm for shipping, and the estimate for aviation varies between −1.2 and 4.3 mm. When focusing on the geographical distribution, the non-CO2 impact from road transport and shipping on the surface air temperature is only slightly stronger in northern than in southern mid-latitudes, while the impact from aviation can be a factor of 5 stronger in the northern than in the southern hemisphere. Further it is observed that most of the impacts are more pronounced at high latitudes, and that the non-CO2 emissions from aviation strongly impact the NAO index. The impacts on the oceanic meridional overturning circulation and the Niño3.4 index are also quantified.
Oceanic bromoform (CHBr3) is an important precursor of atmospheric bromine. Although highly relevant for the future halogen burden and ozone layer in the stratosphere, global CHBr3 production in the ...ocean and its emissions are still poorly constrained in observations and are mostly neglected in climate models. Here, we newly implement marine CHBr3 in the second version of the state-of-the-art Norwegian Earth System Model (NorESM2) with fully coupled interactions of ocean, sea ice, and atmosphere. Our results are validated using oceanic and atmospheric observations from the HalOcAt (Halocarbons in the Ocean and Atmosphere) database. The simulated mean oceanic concentrations (6.61 ± 3.43 pmol L−1) are in good agreement with observations from open-ocean regions (5.02 ± 4.50 pmol L−1), while the mean atmospheric mixing ratios (0.76 ± 0.39 ppt) are lower than observed but within the range of uncertainty (1.45 ± 1.11 ppt). The NorESM2 ocean emissions of CHBr3 (214 Gg yr−1) are within the range of or higher than previously published estimates from bottom-up approaches but lower than estimates from top-down approaches. Annual mean fluxes are mostly positive (sea-to-air fluxes); driven by oceanic concentrations, sea surface temperature, and wind speed; and dependent on season and location. During winter, model results imply that some oceanic regions in high latitudes act as sinks of atmospheric CHBr3 due to their elevated atmospheric mixing ratios. We further demonstrate that key drivers for oceanic and atmospheric CHBr3 variability are spatially heterogeneous. In the tropical West Pacific, which is a hot spot for oceanic bromine delivery to the stratosphere, wind speed is the main driver for CHBr3 fluxes on an annual basis. In the North Atlantic, as well as in the Southern Ocean region, atmospheric and oceanic CHBr3 variabilities interact during most of the seasons except for the winter months, when sea surface temperature is the main driver. Our study provides an improved process-based understanding of the biogeochemical cycling of CHBr3 and more reliable natural emission estimates, especially on seasonal and spatial scales, compared to previously published model estimates.
Differences between a 1.5 and 2.0 ∘C warmer climate than 1850 pre-industrial conditions are investigated using a suite of uncoupled (Atmospheric Model Intercomparison Project; AMIP), fully coupled, ...and slab-ocean experiments performed with Norwegian Earth System Model (NorESM1)-Happi, an upgraded version of NorESM1-M. The data from the AMIP-type runs with prescribed sea-surface temperatures (SSTs) and sea ice were provided to a model intercomparison project (HAPPI – Half a degree Additional warming, Prognosis and Projected Impacts; http://www.happimip.org/, last access date: 14 September 2019). This paper compares the AMIP results to those from the fully coupled version and the slab-ocean version of the model (NorESM1-HappiSO) in which SST and sea ice are allowed to respond to the warming, focusing on Arctic amplification of the global change signal. The fully coupled and the slab-ocean runs generally show stronger responses
than the AMIP runs in the warmer worlds. The Arctic polar amplification
factor is stronger in the fully coupled and slab-ocean runs than in the AMIP runs, both in the 1.5 ∘C warming run and with the additional 0.5 ∘C warming. The low-level Equator-to-pole temperature gradient consistently weakens more between the present-day climate and the 1.5 ∘C warmer climate in the experiments with an active ocean component. The magnitude of the upper-level Equator-to-pole temperature gradient increases in a warmer climate but is not systematically larger in the experiments with an active ocean component. Implications for storm tracks and blocking are investigated. We find considerable reductions in the Arctic sea-ice cover in the slab-ocean model runs; while ice-free summers are rare under 1.5 ∘C warming, they occur 18 % of the time in the 2.0 ∘C warming simulation. The fully coupled model does not, however, reach ice-free conditions as it is too cold and has too much ice in the present-day climate. Differences between the experiments with active ocean and sea-ice models and those with prescribed SSTs and sea ice can be partially due to ocean and
sea-ice feedbacks that are neglected in the latter case but can also in
part be due to differences in the experimental setup.