Automated determination of auroral breakup during the substorm expansion phase using all sky imager data

This technique paper describes a novel method for quantitatively and routinely identifying auroral breakup following substorm onset using the Time History of Events and Macroscale Interactions During Substorms (THEMIS) all-sky imagers (ASIs). Substorm onset is characterised by a brightening of the aurora that is followed by auroral poleward expansion and auroral breakup. This breakup can be identified by a sharp increase in the auroral intensity i(t) and the time derivative of auroral intensity i'(t). Utilising both i(t) and i'(t) we have developed an algorithm for identifying the time interval and spatial location of auroral breakup during the substorm expansion phase within the field of view of ASI data based solely on quantifiable characteristics of the optical auroral emissions. We compare the time interval determined by the algorithm to independently identified auroral onset times from three previously published studies. In each case the time interval determined by the algorithm is within error of the onset independently identified by the prior studies. We further show the utility of the algorithm by comparing the breakup intervals determined using the automated algorithm to an independent list of substorm onset times. We demonstrate that up to 50% of the breakup intervals characterised by the algorithm are within the uncertainty of the times identified in the independent list. The quantitative description and routine identification of an interval of auroral brightening during the substorm expansion phase provides a foundation for unbiased statistical analysis of the aurora to probe the physics of the auroral substorm as a new scientific tool for aiding the identification of the processes leading to auroral substorm onset.

Simulation of the Earth's radiation budget by the European Centre for Medium-Range Weather Forecasts 40-year reanalysis (ERA40)

The radiation budget simulated by the European Centre for Medium-Range Weather Forecasts (ECMWF) 40-year reanalysis (ERA40) is evaluated for the period 1979–2001 using independent satellite data and additional model data. This provides information on the quality of the radiation products and indirect evaluation of other aspects of the climate produced by ERA40. The climatology of clear-sky outgoing longwave radiation (OLR) is well captured by ERA40. Underestimations of about 10 W m−2 in clear-sky OLR over tropical convective regions by ERA40 compared to satellite data are substantially reduced when the satellite sampling is taken into account. The climatology of column-integrated water vapor is well simulated by ERA40 compared to satellite data over the ocean, indicating that the simulation of downward clear-sky longwave fluxes at the surface is likely to be good. Clear-sky absorbed solar radiation (ASR) and clear-sky OLR are overestimated by ERA40 over north Africa and high-latitude land regions. The observed interannual changes in low-latitude means are not well reproduced. Using ERA40 to analyze trends and climate feedbacks globally is therefore not recommended. The all-sky radiation budget is poorly simulated by ERA40. OLR is overestimated by around 10 W m−2 over much of the globe. ASR is underestimated by around 30 W m−2 over tropical ocean regions. Away from marine stratocumulus regions, where cloud fraction is underestimated by ERA40, the poor radiation simulation by ERA40 appears to be related to inaccurate radiative properties of cloud rather than inaccurate cloud distributions.

Simultaneous observations of reconnection pulses at Cluster and their effects on the cusp aurora seen at Chinese Yellow River Station

While the Cluster spacecraft were located near the high-latitude magnetopause, between 10:10 and 10:40 UT on 16 January 2004, three typical flux transfer event (FTE) signatures were observed. During this interval, simultaneous and conjugated all-sky camera measurements, recorded at Yellow River Station, Svalbard, are available at 630.0 and 557.7nm that show poleward-moving auroral forms (PMAFs), consistent with magnetic reconnection at dayside magnetopause. Simultaneous FTEs seen at the magnetopause mainly move northward, but having duskward (eastward) and tailward velocity components, roughly consistent with the observed direction of motion of the PMAFs in all-sky images. Between the PMAFs meridional keograms, extracted from the all-sky images, show intervals of lower intensity aurora which migrate equatorward just before the PMAFs intensify. This is strong evidence for an equatorward eroding and poleward moving open-closed boundary (OCB) associated with a variable magnetopause reconnection rate under variable IMF conditions. From the durations of the PMAFs we infer that the evolution time of FTEs is 5-11 minutes from its origin on magnetopause to its addition to the polar cap.

Simultaneous observations of reconnection pulses at Cluster and their effects on the cusp aurora observed at the Chinese Yellow River station

While the Cluster spacecraft were located near the high-latitude magnetopause, between 1010 and 1040 UT on 16 January 2004, three typical flux transfer event (FTE) signatures were observed. During this interval, simultaneous and conjugated all‐sky camera measurements, recorded at Yellow River Station, Svalbard, are available at 630.0 and 557.7 nm that show poleward‐moving auroral forms (PMAFs), consistent with magnetic reconnection at the dayside magnetopause. Simultaneous FTEs seen at the magnetopause mainly move northward, but having duskward (eastward) and tailward velocity components, roughly consistent with the observed direction of motion of the PMAFs in all‐sky images. Between the PMAFs meridional keograms, extracted from the all‐sky images, show intervals of lower intensity aurora which migrate equatorward just before the PMAFs intensify. This is strong evidence for an equatorward eroding and poleward moving open‐closed boundary associated with a variable magnetopause reconnection rate under variable IMF conditions. From the durations of the PMAFs, we infer that the evolution time of FTEs is 5–11 minutes from its origin on the magnetopause to its addition to the polar cap.

The roles of aerosol, water vapor and cloud in future global dimming/brightening

Observational evidence indicates significant regional trends in solar radiation at the surface in both all-sky and cloud-free conditions. Negative trends in the downwelling solar surface irradiance (SSI) have become known as ‘dimming’ while positive trends have become known as ‘brightening’. We use the Met Office Hadley Centre HadGEM2 climate model to model trends in cloud-free and total SSI from the pre-industrial to the present-day and compare these against observations. Simulations driven by CMIP5 emissions are used to model the future trends in dimming/brightening up to the year 2100. The modeled trends are reasonably consistent with observed regional trends in dimming and brightening which are due to changes in concentrations in anthropogenic aerosols and, potentially, changes in cloud cover owing to the aerosol indirect effects and/or cloud feedback mechanisms. The future dimming/brightening in cloud-free SSI is not only caused by changes in anthropogenic aerosols: aerosol impacts are overwhelmed by a large dimming caused by increases in water vapor. There is little trend in the total SSI as cloud cover decreases in the climate model used here, and compensates the effect of the change in water vapor. In terms of the surface energy balance, these trends in SSI are obviously more than compensated by the increase in the downwelling terrestrial irradiance from increased water vapor concentrations. However, the study shows that while water vapor is widely appreciated as a greenhouse gas, water vapor impacts on the atmospheric transmission of solar radiation and the future of global dimming/brightening should not be overlooked.

Aerosol forcing in the Climate Model Intercomparison Project (CMIP5) simulations by HadGEM2-ES and the role of ammonium nitrate

The latest Hadley Centre climate model, HadGEM2-ES, includes Earth system components such as interactive chemistry and eight species of tropospheric aerosols. It has been run for the period 1860–2100 in support of the fifth phase of the Climate Model Intercomparison Project (CMIP5). Anthropogenic aerosol emissions peak between 1980 and 2020, resulting in a present-day all-sky top of the atmosphere aerosol forcing of −1.6 and −1.4 W m−2 with and without ammonium nitrate aerosols, respectively, for the sum of direct and first indirect aerosol forcings. Aerosol forcing becomes significantly weaker in the 21st century, being weaker than −0.5 W m−2 in 2100 without nitrate. However, nitrate aerosols become the dominant species in Europe and Asia and decelerate the decrease in global mean aerosol forcing. Considering nitrate aerosols makes aerosol radiative forcing 2–4 times stronger by 2100 depending on the representative concentration pathway, although this impact is lessened when changes in the oxidation properties of the atmosphere are accounted for. Anthropogenic aerosol residence times increase in the future in spite of increased precipitation, as cloud cover and aerosol-cloud interactions decrease in tropical and midlatitude regions. Deposition of fossil fuel black carbon onto snow and ice surfaces peaks during the 20th century in the Arctic and Europe but keeps increasing in the Himalayas until the middle of the 21st century. Results presented here confirm the importance of aerosols in influencing the Earth's climate, albeit with a reduced impact in the future, and suggest that nitrate aerosols will partially replace sulphate aerosols to become an important anthropogenic species in the remainder of the 21st century.

Updated estimate of aerosol direct radiative forcing from satellite observations and comparison against the Hadley Centre climate model

The fourth assessment report of the Intergovernmental Panel on Climate Change (IPCC) includes a comparison of observation-based and modeling-based estimates of the aerosol direct radiative forcing. In this comparison, satellite-based studies suggest a more negative aerosol direct radiative forcing than modeling studies. A previous satellite-based study, part of the IPCC comparison, uses aerosol optical depths and accumulation-mode fractions retrieved by the Moderate Resolution Imaging Spectroradiometer (MODIS) at collection 4. The latest version of MODIS products, named collection 5, improves aerosol retrievals. Using these products, the direct forcing in the shortwave spectrum defined with respect to present-day natural aerosols is now estimated at −1.30 and −0.65 Wm−2 on a global clear-sky and all-sky average, respectively, for 2002. These values are still significantly more negative than the numbers reported by modeling studies. By accounting for differences between present-day natural and preindustrial aerosol concentrations, sampling biases, and investigating the impact of differences in the zonal distribution of anthropogenic aerosols, good agreement is reached between the direct forcing derived from MODIS and the Hadley Centre climate model HadGEM2-A over clear-sky oceans. Results also suggest that satellite estimates of anthropogenic aerosol optical depth over land should be coupled with a robust validation strategy in order to refine the observation-based estimate of aerosol direct radiative forcing. In addition, the complex problem of deriving the aerosol direct radiative forcing when aerosols are located above cloud still needs to be addressed.

Estimates of aerosol radiative forcing from the MACC re-analysis

The European Centre for Medium-range Weather Forecast (ECMWF) provides an aerosol re-analysis starting from year 2003 for the Monitoring Atmospheric Composition and Climate (MACC) project. The re-analysis assimilates total aerosol optical depth retrieved by the Moderate Resolution Imaging Spectroradiometer (MODIS) to correct for model departures from observed aerosols. The reanalysis therefore combines satellite retrievals with the full spatial coverage of a numerical model. Re-analysed products are used here to estimate the shortwave direct and first indirect radiative forcing of anthropogenic aerosols over the period 2003–2010, using methods previously applied to satellite retrievals of aerosols and clouds. The best estimate of globally-averaged, all-sky direct radiative forcing is −0.7±0.3Wm−2. The standard deviation is obtained by a Monte-Carlo analysis of uncertainties, which accounts for uncertainties in the aerosol anthropogenic fraction, aerosol absorption, and cloudy-sky effects. Further accounting for differences between the present-day natural and pre-industrial aerosols provides a direct radiative forcing estimate of −0.4±0.3Wm−2. The best estimate of globally-averaged, all-sky first indirect radiative forcing is −0.6±0.4Wm−2. Its standard deviation accounts for uncertainties in the aerosol anthropogenic fraction, and in cloud albedo and cloud droplet number concentration susceptibilities to aerosol changes. The distribution of first indirect radiative forcing is asymmetric and is bounded by −0.1 and −2.0Wm−2. In order to decrease uncertainty ranges, better observational constraints on aerosol absorption and sensitivity of cloud droplet number concentrations to aerosol changes are required.

Radiative forcing of the direct aerosol effect from AeroCom Phase II simulations

We report on the AeroCom Phase II direct aerosol effect (DAE) experiment where 16 detailed global aerosol models have been used to simulate the changes in the aerosol distribution over the industrial era. All 16 models have estimated the radiative forcing (RF) of the anthropogenic DAE, and have taken into account anthropogenic sulphate, black carbon (BC) and organic aerosols (OA) from fossil fuel, biofuel, and biomass burning emissions. In addition several models have simulated the DAE of anthropogenic nitrate and anthropogenic influenced secondary organic aerosols (SOA). The model simulated all-sky RF of the DAE from total anthropogenic aerosols has a range from −0.58 to −0.02Wm−2, with a mean of −0.27Wm−2 for the 16 models. Several models did not include nitrate or SOA and modifying the estimate by accounting for this with information from the other AeroCom models reduces the range and slightly strengthens the mean. Modifying the model estimates for missing aerosol components and for the time period 1750 to 2010 results in a mean RF for the DAE of −0.35Wm−2. Compared to AeroCom Phase I (Schulz et al., 2006) we find very similar spreads in both total DAE and aerosol component RF. However, the RF of the total DAE is stronger negative and RF from BC from fossil fuel and biofuel emissions are stronger positive in the present study than in the previous AeroCom study.We find a tendency for models having a strong (positive) BC RF to also have strong (negative) sulphate or OA RF. This relationship leads to smaller uncertainty in the total RF of the DAE compared to the RF of the sum of the individual aerosol components. The spread in results for the individual aerosol components is substantial, and can be divided into diversities in burden, mass extinction coefficient (MEC), and normalized RF with respect to AOD. We find that these three factors give similar contributions to the spread in results.

Parameterization of contrails in the UK Met Office Climate Model

Persistent contrails are believed to currently have a relatively small but significant positive radiative forcing on climate. With air travel predicted to continue its rapid growth over the coming years, the contrail warming effect on climate is expected to increase. Nevertheless, there remains a high level of uncertainty in the current estimates of contrail radiative forcing. Contrail formation depends mostly on the aircraft flying in cold and moist enough air masses. Most studies to date have relied on simple parameterizations using averaged meteorological conditions. In this paper we take into account the short‐term variability in background cloudiness by developing an on‐line contrail parameterization for the UK Met Office climate model. With this parameterization, we estimate that for the air traffic of year 2002 the global mean annual linear contrail coverage was approximately 0.11%. Assuming a global mean contrail optical depth of 0.2 or smaller and assuming hexagonal ice crystals, the corresponding contrail radiative forcing was calculated to be less than 10 mW m−2 in all‐sky conditions. We find that the natural cloud masking effect on contrails may be significantly higher than previously believed. This new result is explained by the fact that contrails seem to preferentially form in cloudy conditions, which ameliorates their overall climate impact by approximately 40%.

Estimates of global multicomponent aerosol optical depth and direct radiative perturbation in the Laboratoire de Météorologie Dynamique general circulation model

The global cycle of multicomponent aerosols including sulfate, black carbon (BC),organic matter (OM), mineral dust, and sea salt is simulated in the Laboratoire de Me´te´orologie Dynamique general circulation model (LMDZT GCM). The seasonal open biomass burning emissions for simulation years 2000–2001 are scaled from climatological emissions in proportion to satellite detected fire counts. The emissions of dust and sea salt are parameterized online in the model. The comparison of model-predicted monthly mean aerosol optical depth (AOD) at 500 nm with Aerosol Robotic Network (AERONET) shows good agreement with a correlation coefficient of 0.57(N = 1324) and 76% of data points falling within a factor of 2 deviation. The correlation coefficient for daily mean values drops to 0.49 (N = 23,680). The absorption AOD (ta at 670 nm) estimated in the model is poorly correlated with measurements (r = 0.27, N = 349). It is biased low by 24% as compared to AERONET. The model reproduces the prominent features in the monthly mean AOD retrievals from Moderate Resolution Imaging Spectroradiometer (MODIS). The agreement between the model and MODIS is better over source and outflow regions (i.e., within a factor of 2).There is an underestimation of the model by up to a factor of 3 to 5 over some remote oceans. The largest contribution to global annual average AOD (0.12 at 550 nm) is from sulfate (0.043 or 35%), followed by sea salt (0.027 or 23%), dust (0.026 or 22%),OM (0.021 or 17%), and BC (0.004 or 3%). The atmospheric aerosol absorption is predominantly contributed by BC and is about 3% of the total AOD. The globally and annually averaged shortwave (SW) direct aerosol radiative perturbation (DARP) in clear-sky conditions is �2.17 Wm�2 and is about a factor of 2 larger than in all-sky conditions (�1.04 Wm�2). The net DARP (SW + LW) by all aerosols is �1.46 and �0.59 Wm�2 in clear- and all-sky conditions, respectively. Use of realistic, less absorbing in SW, optical properties for dust results in negative forcing over the dust-dominated regions.

Host model uncertainties in aerosol radiative forcing estimates: results from the AeroCom Prescribed intercomparison study

Simulated multi-model “diversity” in aerosol direct radiative forcing estimates is often perceived as a measure of aerosol uncertainty. However, current models used for aerosol radiative forcing calculations vary considerably in model components relevant for forcing calculations and the associated “host-model uncertainties” are generally convoluted with the actual aerosol uncertainty. In this AeroCom Prescribed intercomparison study we systematically isolate and quantify host model uncertainties on aerosol forcing experiments through prescription of identical aerosol radiative properties in twelve participating models. Even with prescribed aerosol radiative properties, simulated clear-sky and all-sky aerosol radiative forcings show significant diversity. For a purely scattering case with globally constant optical depth of 0.2, the global-mean all-sky top-of-atmosphere radiative forcing is −4.47Wm−2 and the inter-model standard deviation is 0.55Wm−2, corresponding to a relative standard deviation of 12 %. For a case with partially absorbing aerosol with an aerosol optical depth of 0.2 and single scattering albedo of 0.8, the forcing changes to 1.04Wm−2, and the standard deviation increases to 1.01W−2, corresponding to a significant relative standard deviation of 97 %. However, the top-of-atmosphere forcing variability owing to absorption (subtracting the scattering case from the case with scattering and absorption) is low, with absolute (relative) standard deviations of 0.45Wm−2 (8 %) clear-sky and 0.62Wm−2 (11 %) all-sky. Scaling the forcing standard deviation for a purely scattering case to match the sulfate radiative forcing in the Aero- Com Direct Effect experiment demonstrates that host model uncertainties could explain about 36% of the overall sulfate forcing diversity of 0.11Wm−2 in the AeroCom Direct Radiative Effect experiment.

A multimodel assessment of the influence of regional anthropogenic emission reductions on aerosol direct radiative forcing and the role of intercontinental transport

In this study, we assess changes of aerosol optical depth (AOD) and direct radiative forcing (DRF) in response to the reduction of anthropogenic emissions in four major pollution regions in the Northern Hemisphere by using results from nine global models in the framework of the Hemispheric Transport of Air Pollution (HTAP). DRF at top of atmosphere (TOA) and surface is estimated based on AOD results from the HTAP models and AOD-normalized DRF (NDRF) from a chemical transport model. The multimodel results show that, on average, a 20% reduction of anthropogenic emissions in North America, Europe, East Asia, and South Asia lowers the global mean AOD (all-sky TOA DRF) by 9.2% (9.0%), 3.5% (3.0%), and 9.4% (10.0%) for sulfate, particulate organic matter (POM), and black carbon (BC), respectively. Global annual average TOA all-sky forcing efficiency relative to particle or gaseous precursor emissions from the four regions (expressed as multimodel mean ± one standard deviation) is  ±3.5 ±0.8,  ±4.0 ±1.7, and 29.5 ±18.1mWm ±2 per Tg for sulfate (relative to SO2), POM, and BC, respectively. The impacts of the regional emission reductions on AOD and DRF extend well beyond the source regions because of intercontinental transport (ICT). On an annual basis, ICT accounts for 11 ±5% to 31 ±9% of AOD and DRF in a receptor region at continental or subcontinental scale, with domestic emissions accounting for the remainder, depending on regions and species. For sulfate AOD, the largest ICT contribution of 31 ±9% occurs in South Asia, which is dominated by the emissions from Europe. For BC AOD, the largest ICT contribution of 28 ±18% occurs in North America, which is dominated by the emissions from East Asia. The large spreads among models highlight the need to improve aerosol processes in models, and evaluate and constrain models with observations.

Global radiative and climate effect of the water vapour continuum at visible and near-infrared wavelengths

Recent laboratory measurements show that absorption by the water vapour continuum in near-infrared windows may be about an order of magnitude higher than assumed in many radiation codes. The radiative impact of the continuum at visible and near-infrared wavelengths is examined for the present day and for a possible future warmer climate (with a global-mean total column water increase of 33%). The calculations use a continuum model frequently used in climate models (‘CKD’) and a continuum model where absorption is enhanced at wavelengths greater than 1 µm based on recent measurements (‘CAVIAR’). The continuum predominantly changes the partitioning between solar radiation absorbed by the surface and the atmosphere; changes in top-of-atmosphere net irradiances are smaller. The global-mean clear-sky atmospheric absorption is enhanced by 1.5 W m−2 (about 2%) and 2.8 W m−2 (about 3.5%) for CKD and CAVIAR respectively, relative to a hypothetical no-continuum case, with all-sky enhancements about 80% of these values. The continuum is, in relative terms, more important for radiation budget changes between the present day and a possible future climate. Relative to the no-continuum case, the increase in global-mean clear-sky absorption is 8% higher using CKD and almost 20% higher using CAVIAR; all-sky enhancements are about half these values. The effect of the continuum is estimated for the solar component of the water vapour feedback, the reduction in downward surface irradiance and precipitation change in a warmer world. For CKD and CAVIAR respectively, and relative to the no-continuum case, the solar component of the water vapour feedback is enhanced by about 4 and 9%, the change in clear-sky downward surface irradiance is 7 and 18% more negative, and the global-mean precipitation response decreases by 1 and 4%. There is a continued need for improved continuum measurements, especially at atmospheric temperatures and at wavelengths below 2 µm.