Attributing trends in extremely hot days to changes in atmospheric dynamics

This paper presents a method for attributing regional trends in the frequency of extremely hot days (EHDs) to changes in the frequency of the atmospheric patterns that characterize such extraordinary events. The study is applied to mainland Spain and the Balearic Islands for the extended summers of the period 1958–2008, where significant and positive trends in maximum temperature (Tx) have been reported during the second half of the past century. First, the study area was split into eight regions attending to their different temporal variability of the daily Tx series obtained from the Spain02 gridded data set using a clustering procedure. Second, the large-scale atmospheric situations causing EHDs are defined by circulation types (CTs). The obtainment of the CTs differs from the majority of CT classifications proposed in the literature. It is based on regional series and on a previous characterization of the main atmospheric situations obtained using only some days classified as extremes in the different regions. Three different atmospheric fields (SLP, T850, and Z500) from ECMWF reanalysis and analysis data and combinations of them (SLP–T850, SLP–Z500, and T850–Z500) are used to produce six different CT classifications. Subsequently, links between EHD occurrence in the different regions and CT for all days have been established. Finally, a simple model to relate the trends in EHDsfor each region to the changes in the CT frequency appearance has been formulated. Most regions present positive and significant trends in the occurrence of EHDs. The CT classifications using two variables perform better. In particular, SLP–T850 is the best for characterizing the atmospheric situations leading to EHD occurrences for most of the regions. Only a small number of CTs have significant trends in their frequency and are associated with high efficiency causing EHD occurrences in most regions simultaneously, especially in the northern and central regions. Attribution results show that changes in circulation can only explain some part of the regional EHD trends. The percentage of the trend attributable to changes in atmospheric dynamics varies from 15 to 50 %, depends on the region and is sensitive to the selected large-scale variables.

Natural and anthropogenic methane fluxes in Eurasia: a mesoscale quantification by generalized atmospheric inversion

Eight surface observation sites providing quasi-continuous measurements of atmospheric methane mixingratios have been operated since the mid-2000’s in Siberia. For the first time in a single work, we assimilate 1 year of these in situ observations in an atmospheric inversion. Our objective is to quantify methane surface fluxes from anthropogenic and wetland sources at the mesoscale in the Siberian lowlands for the year 2010. To do so, we first inquire about the way the inversion uses the observations and the way the fluxes are constrained by the observation sites. As atmospheric inver- sions at the mesoscale suffer from mis-quantified sources of uncertainties, we follow recent innovations in inversion techniques and use a new inversion approach which quantifies the uncertainties more objectively than the previous inversion systems. We find that, due to errors in the representation of the atmospheric transport and redundant pieces of information, only one observation every few days is found valuable by the inversion. The remaining high-resolution quasi-continuous signal is representative of very local emission patterns difficult to analyse with a mesoscale system. An analysis of the use of information by the inversion also reveals that the observation sites constrain methane emissions within a radius of 500 km. More observation sites than the ones currently in operation are then necessary to constrain the whole Siberian lowlands. Still, the fluxes within the constrained areas are quantified with objectified uncertainties. Finally, the tolerance intervals for posterior methane fluxes are of roughly 20 % (resp. 50 %) of the fluxes for anthropogenic (resp. wetland) sources. About 50–70 % of Siberian lowlands emissions are constrained by the inversion on average on an annual basis. Extrapolating the figures on the constrained areas to the whole Siberian lowlands, we find a regional methane budget of 5–28 TgCH4 for the year 2010, i.e. 1–5 % of the global methane emissions. As very few in situ observations are available in the region of interest, observations of methane total columns from the Greenhouse Gas Observing SATellite (GOSAT) are tentatively used for the evaluation of the inversion results, but they exhibit only a marginal signal from the fluxes within the region of interest.

Evaluating atmospheric methane inversion model results for Pallas, northern Finland

A state-of-the-art inverse model, CarbonTracker Data Assimilation Shell (CTDAS), was used to optimize estimates of methane (CH4) surface fluxes using atmospheric observations of CH4 as a constraint. The model consists of the latest version of the TM5 atmospheric chemistry-transport model and an ensemble Kalman filter based data assimilation system. The model was constrained by atmospheric methane surface concentrations, obtained from the World Data Centre for Greenhouse Gases (WDCGG). Prior methane emissions were specified for five sources: biosphere, anthropogenic, fire, termites and ocean, of which bio-sphere and anthropogenic emissions were optimized. Atmospheric CH 4 mole fractions for 2007 from northern Finland calculated from prior and optimized emissions were compared with observations. It was found that the root mean squared errors of the posterior esti - mates were more than halved. Furthermore, inclusion of NOAA observations of CH 4 from weekly discrete air samples collected at Pallas improved agreement between posterior CH 4 mole fraction estimates and continuous observations, and resulted in reducing optimized biosphere emissions and their uncertainties in northern Finland.

Impact of oceanic circulation changes on atmospheric δ¹³ CO₂

δ¹³ CO₂ measured in Antarctic ice cores provides constraints on oceanic and terrestrial carbon cycle processes linked with millennial-scale changes in atmospheric CO₂. However, the interpretation of δ¹³ CO₂ is not straight-forward. Using carbon isotope-enabled versions of the LOVECLIM and Bern3D models, we perform a set of sensitivity experiments in which the formation rates of North Atlantic Deep Water (NADW), North Pacific Deep Water (NPDW), Antarctic Bottom Water (AABW), and Antarctic Intermediate Water (AAIW) are varied. We study the impact of these circulation changes on atmospheric δ¹³ CO₂ as well as on the oceanic δ¹³ CO₂ distribution. In general, we find that the formation rates of AABW, NADW, NPDW, and AAIW are negatively correlated with changes in δ¹³ CO₂: namely, strong oceanic ventilation decreases atmospheric δ¹³ CO₂. However, since large-scale oceanic circulation reorganizations also impact nutrient utilization and the Earth’s climate, the relationship between atmospheric δ¹³ CO₂ levels and ocean ventilation rate is not unequivocal. In both models atmospheric δ¹³ CO₂ is very sensitive to changes in AABW formation rates: increased AABW formation enhances the transport of low δ¹³ CO₂ waters to the surface and decreases atmospheric δ¹³ CO₂. By contrast, the impact of NADW changes on atmospheric δ¹³ CO₂ is less robust and might be model dependent. This results from complex interplay between global climate, carbon cycle, and the formation rate of NADW, a water body characterized by relatively high δ¹³ CO₂.

Trajectory mapping of middle atmospheric water vapor by a mini network of NDACC instruments

The important task to observe the global coverage of middle atmospheric trace gases like water vapor or ozone usually is accomplished by satellites. Climate and atmospheric studies rely upon the knowledge of trace gas distributions throughout the stratosphere and mesosphere. Many of these gases are currently measured from satellites, but it is not clear whether this capability will be maintained in the future. This could lead to a significant knowledge gap of the state of the atmosphere. We explore the possibilities of mapping middle atmospheric water vapor in the Northern Hemisphere by using Lagrangian trajectory calculations and water vapor profile data from a small network of five ground-based microwave radiometers. Four of them are operated within the frame of NDACC (Network for the Detection of Atmospheric Composition Change). Keeping in mind that the instruments are based on different hardware and calibration setups, a height-dependent bias of the retrieved water vapor profiles has to be expected among the microwave radiometers. In order to correct and harmonize the different data sets, the Microwave Limb Sounder (MLS) on the Aura satellite is used to serve as a kind of traveling standard. A domain-averaging TM (trajectory mapping) method is applied which simplifies the subsequent validation of the quality of the trajectory-mapped water vapor distribution towards direct satellite observations. Trajectories are calculated forwards and backwards in time for up to 10 days using 6 hourly meteorological wind analysis fields. Overall, a total of four case studies of trajectory mapping in different meteorological regimes are discussed. One of the case studies takes place during a major sudden stratospheric warming (SSW) accompanied by the polar vortex breakdown; a second takes place after the reformation of stable circulation system. TM cases close to the fall equinox and June solstice event from the year 2012 complete the study, showing the high potential of a network of ground-based remote sensing instruments to synthesize hemispheric maps of water vapor.

Zeeman effect in atmospheric O₂ measured by ground-based microwave radiometry

In this work we study the Zeeman effect on stratospheric O₂ using ground-based microwave radiometer measurements. The interaction of the Earth magnetic field with the oxygen dipole leads to a splitting of O₂ energy states, which polarizes the emission spectra. A special campaign was carried out in order to measure this effect in the oxygen emission line centered at 53.07 GHz. Both a fixed and a rotating mirror were incorporated into the TEMPERA (TEMPERature RAdiometer) in order to be able to measure under different observational angles. This new configuration allowed us to change the angle between the observational path and the Earth magnetic field direction. Moreover, a high-resolution spectrometer (1 kHz) was used in order to measure for the first time the polarization state of the radiation due to the Zeeman effect in the main isotopologue of oxygen from ground-based microwave measurements. The measured spectra showed a clear polarized signature when the observational angles were changed, evidencing the Zeeman effect in the oxygen molecule. In addition, simulations carried out with the Atmospheric Radiative Transfer Simulator (ARTS) allowed us to verify the microwave measurements showing a very good agreement between model and measurements. The results suggest some interesting new aspects for research of the upper atmosphere.

Hygroscopic growth of atmospheric aerosol particles based on active remote sensing and radiosounding measurements: selected cases in southeastern Spain

A new methodology based on combining active and passive remote sensing and simultaneous and collocated radiosounding data to study the aerosol hygroscopic growth effects on the particle optical and microphysical properties is presented. The identification of hygroscopic growth situations combines the analysis of multispectral aerosol particle backscatter coefficient and particle linear depolarization ratio with thermodynamic profiling of the atmospheric column. We analyzed the hygroscopic growth effects on aerosol properties, namely the aerosol particle backscatter coefficient and the volume concentration profiles, using data gathered at Granada EARLINET station. Two study cases, corresponding to different aerosol loads and different aerosol types, are used for illustrating the potential of this methodology. Values of the aerosol particle backscatter coefficient enhancement factors range from 2.1 ± 0.8 to 3.9 ± 1.5, in the ranges of relative humidity 60–90 and 40–83%, being similar to those previously reported in the literature. Differences in the enhancement factor are directly linked to the composition of the atmospheric aerosol. The largest value of the aerosol particle backscatter coefficient enhancement factor corresponds to the presence of sulphate and marine particles that are more affected by hygroscopic growth. On the contrary, the lowest value of the enhancement factor corresponds to an aerosol mixture containing sulphates and slight traces of mineral dust. The Hänel parameterization is applied to these case studies, obtaining results within the range of values reported in previous studies, with values of the γ exponent of 0.56 ± 0.01 (for anthropogenic particles slightly influenced by mineral dust) and 1.07 ± 0.01 (for the situation dominated by anthropogenic particles), showing the convenience of this remote sensing approach for the study of hygroscopic effects of the atmospheric aerosol under ambient unperturbed conditions. For the first time, the retrieval of the volume concentration profiles for these cases using the Lidar Radiometer Inversion Code (LIRIC) allows us to analyze the aerosol hygroscopic growth effects on aerosol volume concentration, observing a stronger increase of the fine mode volume concentration with increasing relative humidity.

Comparison of co-located independent ground-based middle atmospheric wind and temperature measurements with numerical weather prediction models

High-resolution, ground-based and independent observations including co-located wind radiometer, lidar stations, and infrasound instruments are used to evaluate the accuracy of general circulation models and data-constrained assimilation systems in the middle atmosphere at northern hemisphere midlatitudes. Systematic comparisons between observations, the European Centre for Medium-Range Weather Forecasts (ECMWF) operational analyses including the recent Integrated Forecast System cycles 38r1 and 38r2, the NASA’s Modern-Era Retrospective Analysis for Research and Applications (MERRA) reanalyses, and the free-running climate Max Planck Institute–Earth System Model–Low Resolution (MPI-ESM-LR) are carried out in both temporal and spectral dom ains. We find that ECMWF and MERRA are broadly consistent with lidar and wind radiometer measurements up to ~40 km. For both temperature and horizontal wind components, deviations increase with altitude as the assimilated observations become sparser. Between 40 and 60 km altitude, the standard deviation of the mean difference exceeds 5 K for the temperature and 20 m/s for the zonal wind. The largest deviations are observed in winter when the variability from large-scale planetary waves dominates. Between lidar data and MPI-ESM-LR, there is an overall agreement in spectral amplitude down to 15–20 days. At shorter time scales, the variability is lacking in the model by ~10 dB. Infrasound observations indicate a general good agreement with ECWMF wind and temperature products. As such, this study demonstrates the potential of the infrastructure of the Atmospheric Dynamics Research Infrastructure in Europe project that integrates various measurements and provides a quantitative understanding of stratosphere-troposphere dynamical coupling for numerical weather prediction applications.

Fossil and non-fossil source contributions to atmospheric carbonaceous aerosols during extreme spring grassland fires in Eastern Europe

In early spring the Baltic region is frequently affected by high-pollution events due to biomass burning in that area. Here we present a comprehensive study to investigate the impact of biomass/grass burning (BB) on the evolution and composition of aerosol in Preila, Lithuania, during springtime open fires. Non-refractory submicron particulate matter (NR-PM1) was measured by an Aerodyne aerosol chemical speciation monitor (ACSM) and a source apportionment with the multilinear engine (ME-2) running the positive matrix factorization (PMF) model was applied to the organic aerosol fraction to investigate the impact of biomass/grass burning. Satellite observations over regions of biomass burning activity supported the results and identification of air mass transport to the area of investigation. Sharp increases in biomass burning tracers, such as levoglucosan up to 683 ngm-3 and black carbon (BC) up to 17 μgm-3 were observed during this period. A further separation between fossil and non-fossil primary and secondary contributions was obtained by coupling ACSM PMF results and radiocarbon (14C) measurements of the elemental (EC) and organic (OC) carbon fractions. Non-fossil organic carbon (OCnf/ was the dominant fraction of PM1, with the primary (POCnf/ and secondary (SOCnf/ fractions contributing 26–44% and 13–23% to the total carbon (TC), respectively. 5–8% of the TC had a primary fossil origin (POCf/, whereas the contribution of fossil secondary organic carbon (SOCf/ was 4–13 %. Nonfossil EC (ECnf/ and fossil EC (ECf/ ranged from 13–24 and 7–13 %, respectively. Isotope ratios of stable carbon and nitrogen isotopes were used to distinguish aerosol particles associated with solid and liquid fossil fuel burning.

Evolution of the stable carbon isotope composition of atmospheric CO₂ over the last glacial cycle

We present new δ¹³C measurements of atmospheric CO₂ covering the last glacial/interglacial cycle, complementing previous records covering Terminations I and II. Most prominent in the new record is a significant depletion in δ¹³C(atm) of 0.5‰ occurring during marine isotope stage (MIS) 4, followed by an enrichment of the same magnitude at the beginning of MIS 3. Such a significant excursion in the record is otherwise only observed at glacial terminations, suggesting that similar processes were at play, such as changing sea surface temperatures, changes in marine biological export in the Southern Ocean (SO) due to variations in aeolian iron fluxes, changes in the Atlantic meridional overturning circulation, upwelling of deep water in the SO, and long-term trends in terrestrial carbon storage. Based on previous modeling studies, we propose constraints on some of these processes during specific time intervals. The decrease in δ¹³C(atm) at the end of MIS 4 starting approximately 64 kyr B.P. was accompanied by increasing [CO₂]. This period is also marked by a decrease in aeolian iron flux to the SO, followed by an increase in SO upwelling during Heinrich event 6, indicating that it is likely that a large amount of δ¹³C-depleted carbon was transferred to the deep oceans previously, i.e., at the onset of MIS 4. Apart from the upwelling event at the end of MIS 4 (and potentially smaller events during Heinrich events in MIS 3), upwelling of deep water in the SO remained reduced until the last glacial termination, whereupon a second pulse of isotopically light carbon was released into the atmosphere.

Atmospheric histories and global emissions of halons H-1211 (CBrClF 2 ), H-1301 (CBrF 3 ), and H-2402 (CBrF 2 CBrF 2 )

We report ground-based atmospheric measurements and emission estimates for the halons H-1211 (CBrClF₂), H-1301 (CBrF₃), and H-2402 (CBrF₂CBrF₂) from the AGAGE (Advanced Global Atmospheric Gases Experiment) and the National Oceanic and Atmospheric Administration global networks. We also include results from archived air samples in canisters and from polar firn in both hemispheres, thereby deriving an atmospheric record of nearly nine decades (1930s to present). All three halons were absent from the atmosphere until ~1970, when their atmospheric burdens started to increase rapidly. In recent years H-1211 and H-2402 mole fractions have been declining, but H-1301 has continued to grow. High-frequency observations show continuing emissions of H-1211 and H-1301 near most AGAGE sites. For H-2402 the only emissions detected were derived from the region surrounding the Sea of Japan/East Sea. Based on our observations, we derive global emissions using two different inversion approaches. Emissions for H-1211 declined from a peak of 11 kt yr⁻¹ (late 1990s) to 3.9 kt yr⁻¹ at the end of our record (mean of 2013–2015), for H-1301 from 5.4 kt yr⁻¹ (late 1980s) to 1.6 kt yr⁻¹, and for H-2402 from 1.8 kt yr⁻¹ (late 1980s) to 0.38 kt yr⁻¹. Yearly summed halon emissions have decreased substantially; nevertheless, since 2000 they have accounted for ~30% of the emissions of all major anthropogenic ozone depletion substances, when weighted by ozone depletion potentials.

Validation of the Swiss methane emission inventory by atmospheric observations and inverse modelling

Atmospheric inverse modelling has the potential to provide observation-based estimates of greenhouse gas emissions at the country scale, thereby allowing for an independent validation of national emission inventories. Here, we present a regional-scale inverse modelling study to quantify the emissions of methane (CH₄) from Switzerland, making use of the newly established CarboCount-CH measurement network and a high-resolution Lagrangian transport model. In our reference inversion, prior emissions were taken from the "bottom-up" Swiss Greenhouse Gas Inventory (SGHGI) as published by the Swiss Federal Office for the Environment in 2014 for the year 2012. Overall we estimate national CH₄ emissions to be 196 ± 18 Gg yr⁻¹ for the year 2013 (1σ uncertainty). This result is in close agreement with the recently revised SGHGI estimate of 206 ± 33 Gg yr⁻¹ as reported in 2015 for the year 2012. Results from sensitivity inversions using alternative prior emissions, uncertainty covariance settings, large-scale background mole fractions, two different inverse algorithms (Bayesian and extended Kalman filter), and two different transport models confirm the robustness and independent character of our estimate. According to the latest SGHGI estimate the main CH₄ source categories in Switzerland are agriculture (78 %), waste handling (15 %) and natural gas distribution and combustion (6 %). The spatial distribution and seasonal variability of our posterior emissions suggest an overestimation of agricultural CH₄ emissions by 10 to 20 % in the most recent SGHGI, which is likely due to an overestimation of emissions from manure handling. Urban areas do not appear as emission hotspots in our posterior results, suggesting that leakages from natural gas distribution are only a minor source of CH₄ in Switzerland. This is consistent with rather low emissions of 8.4 Gg yr⁻¹ reported by the SGHGI but inconsistent with the much higher value of 32 Gg yr⁻¹ implied by the EDGARv4.2 inventory for this sector. Increased CH₄ emissions (up to 30 % compared to the prior) were deduced for the north-eastern parts of Switzerland. This feature was common to most sensitivity inversions, which is a strong indicator that it is a real feature and not an artefact of the transport model and the inversion system. However, it was not possible to assign an unambiguous source process to the region. The observations of the CarboCount-CH network provided invaluable and independent information for the validation of the national bottom-up inventory. Similar systems need to be sustained to provide independent monitoring of future climate agreements.

M-AARE - Coupling atmospheric, hydrological, hydrodynamic and damage models in the Aare river basin, Switzerland

The triggering mechanism and the temporal evolution of large flood events, especially of worst-case scenarios, are not yet fully understood. Consequently, the cumulative losses of extreme floods are unknown. To study the link between weather conditions, discharges and flood losses it is necessary to couple atmospheric, hydrological, hydrodynamic and damage models. The objective of the M-AARE project is to test the potentials and opportunities of a model chain that relates atmospheric conditions to flood losses or risks. The M-AARE model chain is a set of coupled models consisting of four main components: the precipitation module, the hydrology module, the hydrodynamic module, and the damage module. The models are coupled in a cascading framework with harmonized time-steps. First exploratory applications show that the one way coupling of the WRF-PREVAH-BASEMENT models has been achieved and provides promising new insights for a better understanding of key aspects in flood risk analysis.