Impact of a Reduced Arctic Sea Ice Cover on Ocean and Atmospheric Properties

The Arctic sea ice cover declined over the last few decades and reached a record minimum in 2007, with a slight recovery thereafter. Inspired by this the authors investigate the response of atmospheric and oceanic properties to a 1-yr period of reduced sea ice cover. Two ensembles of equilibrium and transient simulations are produced with the Community Climate System Model. A sea ice change is induced through an albedo change of 1 yr. The sea ice area and thickness recover in both ensembles after 3 and 5 yr, respectively. The sea ice anomaly leads to changes in ocean temperature and salinity to a depth of about 200 m in the Arctic Basin. Further, the salinity and temperature changes in the surface layer trigger a “Great Salinity Anomaly” in the North Atlantic that takes roughly 8 yr to travel across the North Atlantic back to high latitudes. In the atmosphere the changes induced by the sea ice anomaly do not last as long as in the ocean. The response in the transient and equilibrium simulations, while similar overall, differs in specific regional and temporal details. The surface air temperature increases over the Arctic Basin and the anomaly extends through the whole atmospheric column, changing the geopotential height fields and thus the storm tracks. The patterns of warming and thus the position of the geopotential height changes vary in the two ensembles. While the equilibrium simulation shifts the storm tracks to the south over the eastern North Atlantic and Europe, the transient simulation shifts the storm tracks south over the western North Atlantic and North America. The authors propose that the overall reduction in sea ice cover is important for producing ocean anomalies; however, for atmospheric anomalies the regional location of the sea ice anomalies is more important. While observed trends in Arctic sea ice are large and exceed those simulated by comprehensive climate models, there is little evidence based on this particular model that the seasonal loss of sea ice (e.g., as occurred in 2007) would constitute a threshold after which the Arctic would exhibit nonlinear, irreversible, or strongly accelerated sea ice loss. Caution should be exerted when extrapolating short-term trends to future sea ice behavior.

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.

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.

Atmospheric water parameters measured by a ground-based microwave radiometer and compared with the WRF model

The microwave radiometer TROWARA measures integrated water vapour (IWV) and integrated cloud liquid water (ILW) at Bern since 1994 with a time resolution of 7 s. In this study, we compare TROWARA measurements with a simulation of summer 2012 in Switzerland performed with the Weather Research and Forecasting (WRF) model. It is found that the WRF model agrees very well with TROWARA’s IWV variations with a mean bias of only 0.7 mm. The ILW distribution of the WRF model, although similar in shape to TROWARA’s distribution, overestimates the fraction of clear sky periods (83% compared to 60%).

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.

Accumulation rates and predominant atmospheric sources of natural and anthropogenic Hg and Pb on the Faroe Islands

A monolith representing 5420 14C yr of peat accumulation was collected from a blanket bog at Myrarnar, Faroe Islands. The maximum Hg concentration (498 ng/g at a depth of 4.5 cm) coincides with the maximum concentration of anthropogenic Pb (111 μg/g). Age dating of recent peat accumulation using 210Pb (CRS model) shows that the maxima in Hg and Pb concentrations occur at AD 1954 ± 2. These results, combined with the isotopic composition of Pb in that sample (206Pb/207Pb = 1.1720 ± 0.0017), suggest that coal burning was the dominant source of both elements. From the onset of peat accumulation (ca. 4286 BC) until AD 1385, the ratios Hg/Br and Hg/Se were constant (2.2 ± 0.5 × 10-4 and 8.5 ± 1.8 × 10-3, respectively). Since then, Hg/Br and Hg/Se values have increased, also reaching their maxima in AD 1954. The age date of the maximum concentrations of anthropogenic Hg and Pb in the Faroe Islands is consistent with a previous study of peat cores from Greenland and Denmark (dated using the atmospheric bomb pulse curve of 14C), which showed maximum concentrations in AD 1953. The average rate of atmospheric Hg accumulation from 1520 BC to AD 1385 was 1.27 ± 0.38 μg/m2/yr. The Br and Se concentrations and the background Hg/Br and Hg/Se ratios were used to calculate the average rate of natural Hg accumulation for the same period, 1.32 ± 0.36 μg/m2/yr and 1.34 ± 0.29 μg/m2/yr, respectively. These fluxes are similar to the preanthropogenic rates obtained using peat cores from Switzerland, southern Greenland, southern Ontario, Canada, and the northeastern United States. Episodic volcanic emissions and the continual supply of marine aerosols to the Faroe Islands, therefore, have not contributed significantly to the Hg inventory or the Hg accumulation rates, relative to these other areas. The maximum rate of Hg accumulation was 34 μg/m2/yr. The greatest fluxes of anthropogenic Hg accumulation calculated using Br and Se, respectively, were 26 and 31 μg/m2/yr. The rate of atmospheric Hg accumulation in 1998 (16 μg/m2/yr) is comparable to the values recently obtained by atmospheric transport modeling for Denmark, the Faroe Islands, and Greenland.

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.