ALBERT

Description: Tidally induced velocity variations of the Beardmore Glacier, Antarctica, and their representation in satellite measurements of ice velocity The Cryosphere, 7, 1375-1384, 2013 Author(s): O. J. Marsh, W. Rack, D. Floricioiu, N. R. Golledge, and W. Lawson Ocean tides close to the grounding line of outlet glaciers around Antarctica have been shown to directly influence ice velocity, both linearly and non-linearly. These fluctuations can be significant and have the potential to affect satellite measurements of ice discharge, which assume displacement between satellite passes to be consistent and representative of annual means. Satellite observations of horizontal velocity variation in the grounding zone are also contaminated by vertical tidal effects, the importance of which is highlighted here in speckle tracking measurements. Eight TerraSAR-X scenes from the grounding zone of the Beardmore Glacier are analysed in conjunction with GPS measurements to determine short-term and decadal trends in ice velocity. Diurnal tides produce horizontal velocity fluctuations of 〉50% on the ice shelf, recorded in the GPS data 4 km downstream of the grounding line. This variability decreases rapidly to 〈5% only 15 km upstream of the grounding line. Daily fluctuations are smoothed to 〈1% in the 11-day repeat pass TerraSAR-X imagery, but fortnightly variations over this period are still visible and show that satellite-velocity measurements can be affected by tides over longer periods. The measured tidal displacement observed in radar look direction over floating ice also allows the grounding line to be identified, using differential speckle tracking where phase information cannot be easily unwrapped.

Description: Snow on the Ross Ice Shelf: comparison of reanalyses and observations from automatic weather stations The Cryosphere, 7, 1399-1410, 2013 Author(s): L. Cohen and S. Dean Snow accumulation measurements from automatic weather stations (AWS) around the Ross Ice Shelf (RIS), Antarctica, are used to provide a new set of ground-based observations which are compared to precipitation from the ECMWF ERA-Interim and NCEP/NCAR Reanalysis-2 datasets. The high temporal resolution of the AWS snow accumulation measurements allow for an event-based comparison of reanalyses precipitation to the in situ observations. Snow accumulation records from nine AWS provide multiple years of accumulation data between 2008 and 2012 over a relatively large, homogeneous region of Antarctica, and also provide the basis for a statistical evaluation of accumulation and precipitation events. The complex effects of wind on snow accumulation (which can both limit and enhance accumulation) complicate the use of the accumulation measurements, but this analysis shows that they can provide a valuable source of ground-based observations for comparisons to modelled precipitation on synoptic timescales. The analysis shows that ERA-Interim reproduces more precipitation events than NCEP-2, and these events correspond to an average 8.2% more precipitation. Significant correlations between reanalyses and AWS event sizes are seen at several stations and show that ERA-Interim consistently produces larger precipitation events than NCEP-2.

Description: A general treatment of snow microstructure exemplified by an improved relation for thermal conductivity The Cryosphere, 7, 1473-1480, 2013 Author(s): H. Löwe, F. Riche, and M. Schneebeli Finding relevant microstructural parameters beyond density is a longstanding problem which hinders the formulation of accurate parameterizations of physical properties of snow. Towards a remedy, we address the effective thermal conductivity tensor of snow via anisotropic, second-order bounds. The bound provides an explicit expression for the thermal conductivity and predicts the relevance of a microstructural anisotropy parameter Q , which is given by an integral over the two-point correlation function and unambiguously defined for arbitrary snow structures. For validation we compiled a comprehensive data set of 167 snow samples. The set comprises individual samples of various snow types and entire time series of metamorphism experiments under isothermal and temperature gradient conditions. All samples were digitally reconstructed by micro-computed tomography to perform microstructure-based simulations of heat transport. The incorporation of anisotropy via Q considerably reduces the root mean square error over the usual density-based parameterization. The systematic quantification of anisotropy via the two-point correlation function suggests a generalizable route to incorporate microstructure into snowpack models. We indicate the inter-relation of the conductivity to other properties and outline a potential impact of Q on dielectric constant, permeability and adsorption rate of diffusing species in the pore space.

Description: Quantification of model uncertainty in aerosol optical thickness retrieval from Ozone Monitoring Instrument (OMI) measurements Atmospheric Measurement Techniques Discussions, 6, 8509-8541, 2013 Author(s): A. Määttä, M. Laine, J. Tamminen, and J. P. Veefkind We study uncertainty quantification in remote sensing of aerosols in the atmosphere with top of the atmosphere reflectance measurements from the nadir-viewing Ozone Monitoring Instrument (OMI). Focus is on the uncertainty in aerosol model selection of pre-calculated aerosol models and on the statistical modelling of the model inadequacies. The aim is to apply statistical methodologies that improve the uncertainty estimates of the aerosol optical thickness (AOT) retrieval by propagating model selection and model error related uncertainties more realistically. We utilise Bayesian model selection and model averaging methods for the model selection problem and use Gaussian processes to model the smooth systematic discrepancies from the modelled to observed reflectance. The systematic model error is learned from an ensemble of operational retrievals. The operational OMI multi-wavelength aerosol retrieval algorithm OMAERO is used for cloud free, over land pixels of the OMI instrument with the additional Bayesian model selection and model discrepancy techniques. The method is demonstrated with four examples with different aerosol properties: weakly absorbing aerosols, forest fires over Greece and Russia, and Sahara dessert dust. The presented statistical methodology is general; it is not restricted to this particular satellite retrieval application.

Description: Influence of high-order mechanics on simulation of glacier response to climate change: insights from Haig Glacier, Canadian Rocky Mountains The Cryosphere, 7, 1527-1541, 2013 Author(s): S. Adhikari and S. J. Marshall Evolution of glaciers in response to climate change has mostly been simulated using simplified dynamical models. Because these models do not account for the influence of high-order physics, corresponding results may exhibit some biases. For Haig Glacier in the Canadian Rocky Mountains, we test this hypothesis by comparing simulation results obtained from 3-D numerical models that deal with different assumptions concerning physics, ranging from simple shear deformation to comprehensive Stokes flow. In glacier retreat scenarios, we find a minimal role of high-order mechanics in glacier evolution, as geometric effects at our site (the presence of an overdeepened bed) result in limited horizontal movement of ice (flow speed on the order of a few meters per year). Consequently, high-order and reduced models all predict that Haig Glacier ceases to exist by ca. 2080 under ongoing climate warming. The influence of high-order mechanics is evident, however, in glacier advance scenarios, where ice speeds are greater and ice dynamical effects become more important. Although similar studies on other glaciers are essential to generalize such findings, we advise that high-order mechanics are important and therefore should be considered while modeling the evolution of active glaciers. Reduced model predictions may be adequate for other glaciologic and topographic settings, particularly where flow speeds are low and where mass balance changes dominate over ice dynamics in determining glacier geometry.

Description: Meteorological drivers of ablation processes on a cold glacier in the semi-arid Andes of Chile The Cryosphere, 7, 1513-1526, 2013 Author(s): S. MacDonell, C. Kinnard, T. Mölg, L. Nicholson, and J. Abermann Meteorological and surface change measurements collected during a 2.5 yr period are used to calculate surface mass and energy balances at 5324 m a.s.l. on Guanaco Glacier, a cold-based glacier in the semi-arid Andes of Chile. Meteorological conditions are marked by extremely low vapour pressures (annual mean of 1.1 hPa), strong winds (annual mean of 10 m s −1 ), shortwave radiation receipt persistently close to the theoretical site maximum during cloud-free days (mean annual 295 W m −2 ; summer hourly maximum 1354 W m −2 ) and low precipitation rates (mean annual 45 mm w.e.). Snowfall occurs sporadically throughout the year and is related to frontal events in the winter and convective storms during the summer months. Net shortwave radiation provides the greatest source of energy to the glacier surface, and net longwave radiation dominates energy losses. The turbulent latent heat flux is always negative, which means that the surface is always losing mass via sublimation, which is the main form of ablation at the site. Sublimation rates are most strongly correlated with net shortwave radiation, incoming shortwave radiation, albedo and vapour pressure. Low glacier surface temperatures restrict melting for much of the period, however episodic melting occurs during the austral summer, when warm, humid, calm and high pressure conditions restrict sublimation and make more energy available for melting. Low accumulation (131 mm w.e. over the period) and relatively high ablation (1435 mm w.e.) means that mass change over the period was negative (−1304 mm w.e.), which continued the negative trend recorded in the region over the last few decades.

Description: Antarctic ice-mass balance 2003 to 2012: regional reanalysis of GRACE satellite gravimetry measurements with improved estimate of glacial-isostatic adjustment based on GPS uplift rates The Cryosphere, 7, 1499-1512, 2013 Author(s): I. Sasgen, H. Konrad, E. R. Ivins, M. R. Van den Broeke, J. L. Bamber, Z. Martinec, and V. Klemann We present regional-scale mass balances for 25 drainage basins of the Antarctic Ice Sheet (AIS) from satellite observations of the Gravity and Climate Experiment (GRACE) for time period January 2003 to September 2012. Satellite gravimetry estimates of the AIS mass balance are strongly influenced by mass movement in the Earth interior caused by ice advance and retreat during the last glacial cycle. Here, we develop an improved glacial-isostatic adjustment (GIA) estimate for Antarctica using newly available GPS uplift rates, allowing us to more accurately separate GIA-induced trends in the GRACE gravity fields from those caused by current imbalances of the AIS. Our revised GIA estimate is considerably lower than previous predictions, yielding an estimate of apparent mass change of 53 ± 18 Gt yr −1 . Therefore, our AIS mass balance of −114 ± 23 Gt yr −1 is less negative than previous GRACE estimates. The northern Antarctic Peninsula and the Amundsen Sea sector exhibit the largest mass loss (−26 ± 3 Gt yr −1 and −127 ± 7 Gt yr −1 , respectively). In contrast, East Antarctica exhibits a slightly positive mass balance (26 ± 13 Gt yr −1 ), which is, however, mostly the consequence of compensating mass anomalies in Dronning Maud and Enderby Land (positive) and Wilkes and George V Land (negative) due to interannual accumulation variations. In total, 6% of the area constitutes about half the AIS imbalance, contributing 151 ± 7 Gt yr −1 (ca. 0.4 mm yr −1 ) to global mean sea-level change. Most of this imbalance is caused by ice-dynamic speed-up expected to prevail in the near future.

Description: The Greenhouse Gas Climate Change Initiative (GHG-CCI): comparative validation of GHG-CCI SCIAMACHY/ENVISAT and TANSO-FTS/GOSAT CO 2 and CH 4 retrieval algorithm products with measurements from the TCCON network Atmospheric Measurement Techniques Discussions, 6, 8679-8741, 2013 Author(s): B. Dils, M. Buchwitz, M. Reuter, O. Schneising, H. Boesch, R. Parker, S. Guerlet, I. Aben, T. Blumenstock, J. P. Burrows, A. Butz, N. M. Deutscher, C. Frankenberg, F. Hase, O. P. Hasekamp, J. Heymann, M. De Mazière, J. Notholt, R. Sussmann, T. Warneke, D. Griffith, V. Sherlock, and D. Wunch Column-averaged dry-air mole fractions of carbon dioxide and methane have been retrieved from spectra acquired by the TANSO-FTS and SCIAMACHY instruments on board GOSAT and ENVISAT using a range of European retrieval algorithms. These retrievals have been compared with data from ground-based high-resolution Fourier Transform Spectrometers (FTS) from the Total Carbon Column Observing Network (TCCON). The participating algorithms are the Weighting Function Modified Differential Optical Absorption Spectroscopy (DOAS) algorithm (WFMD, University of Bremen), the Bremen Optimal Estimation DOAS algorithm (BESD, University of Bremen), the Iterative Maximum A Posteriori DOAS (IMAP, Jet Propulsion Laboratory (JPL) and Netherlands Institute for Space Research algorithm (SRON)), the proxy and full-physics versions of SRON's RemoTeC algorithm (SRPR and SRFP respectively) and the proxy and full-physics versions of the University of Leicester's adaptation of the OCO (Orbiting Carbon Observatory) algorithm (OCPR and OCFP respectively). The goal of this algorithm inter-comparison was to identify strengths and weaknesses of the various so-called Round Robin data sets generated with the various algorithms so as to determine which of the competing algorithms would proceed to the next round of the European Space Agency's (ESA) Greenhouse Gas Climate Change Initiative (GHG-CCI) project, which is the generation of the so-called Climate Research Data Package (CRDP), which is the first version of the Essential Climate Variable (ECV) "Greenhouse Gases" (GHG). For CO 2 , all algorithms reach the precision requirements for inverse modelling ( 〈 8 ppb), with only WFMD having a lower precision (4.7 ppm) than the other algorithm products (2.4–2.5 ppm). When looking at the seasonal relative accuracy (SRA, variability of the bias in space and time), none of the algorithms have reached the demanding 〈 0.5 ppm threshold. For CH 4 , the precision for both SCIAMACHY products (50.2 ppb for IMAP and 76.4 ppb for WFMD) fail to meet the 〈 34 ppb threshold, but note that this work focusses on the period after the 2005 SCIAMACHY detector degradation. The GOSAT X CH 4 precision ranges between 18.1 and 14.0 ppb. Looking at the SRA, all GOSAT algorithm products reach the 〈 10 ppm threshold (values ranging between 5.4 and 6.2 ppb). For SCIAMACHY, IMAP and WFMD have a SRA of 17.2 ppb and 10.5 ppb respectively.

Description: The effect of phase partitioning of semivolatile compounds on the measured CCN activity of aerosol particles Atmospheric Measurement Techniques Discussions, 6, 8413-8433, 2013 Author(s): S. Romakkaniemi, A. Jaatinen, A. Laaksonen, A. Nenes, and T. Raatikainen The effect of inorganic semivolatile aerosol compounds on the CCN activity of aerosol particles was studied by using a computational model for a DMT-CCN counter, a cloud parcel model for condensation kinetics and experiments to quantify the modelled results. Concentrations of water vapour and semivolatiles as well as aerosol trajectories in the CCN column were calculated by a computational fluid dynamics model. These trajectories and vapour concentrations were then used as an input for the cloud parcel model to simulate mass transfer kinetics of water and semivolatiles between aerosol particles and the gas phase. Two different questions were studied: (1) how big fraction of semivolatiles is evaporated from particles before activation in the CCN counter? (2) How much the CCN activity can be increased due to condensation of semivolatiles prior to the maximum water supersaturation in the case of high semivolatile concentration in the gas phase? The results show that, to increase the CCN activity of aerosol particles, a very high gas phase concentration (as compared to typical ambient conditions) is needed. We used nitric acid as a test compound. A concentration of several ppb or higher is needed for measurable effect. In the case of particle evaporation, we used ammonium nitrate as a test compound and found that it partially evaporates before maximum supersaturation is reached in the CCN counter, thus causing an underestimation of CCN activity. The effect of evaporation is clearly visible in all supersaturations, leading to an underestimation of the critical dry diameter by 10 to 15 nanometres in the case of ammonium nitrate particles in different supersaturations. This result was also confirmed by measurements in supersaturations between 0.1 and 0.7%.

Description: Evidence of meltwater retention within the Greenland ice sheet The Cryosphere, 7, 1433-1445, 2013 Author(s): A. K. Rennermalm, L. C. Smith, V. W. Chu, J. E. Box, R. R. Forster, M. R. Van den Broeke, D. Van As, and S. E. Moustafa Greenland ice sheet mass losses have increased in recent decades with more than half of these attributed to surface meltwater runoff. However, the magnitudes of englacial storage, firn retention, internal refreezing and other hydrologic processes that delay or reduce true water export to the global ocean remain less understood, partly due to a scarcity of in situ measurements. Here, ice sheet surface meltwater runoff and proglacial river discharge between 2008 and 2010 near Kangerlussuaq, southwestern Greenland were used to establish sub- and englacial meltwater storage for a small ice sheet watershed (36–64 km 2 ). This watershed lacks significant potential meltwater storage in firn, surface lakes on the ice sheet and in the proglacial area, and receives limited proglacial precipitation. Thus, ice sheet surface runoff not accounted for by river discharge can reasonably be attributed to retention in sub- and englacial storage. Evidence for meltwater storage within the ice sheet includes (1) characteristic dampened daily river discharge amplitudes relative to ice sheet runoff; (2) three cold-season river discharge anomalies at times with limited ice sheet surface melt, demonstrating that meltwater may be retained up to 1–6 months; (3) annual ice sheet watershed runoff is not balanced by river discharge, and while near water budget closure is possible as much as 54% of melting season ice sheet runoff may not escape to downstream rivers; (4) even the large meltwater retention estimate (54%) is equivalent to less than 1% of the ice sheet volume, which suggests that storage in en- and subglacial cavities and till is plausible. While this study is the first to provide evidence for meltwater retention and delayed release within the Greenland ice sheet, more information is needed to establish how widespread this is along the Greenland ice sheet perimeter.

Description: High-resolution provenance of desert dust deposited on Mt. Elbrus, Caucasus in 2009–2012 using snow pit and firn core records The Cryosphere, 7, 1481-1498, 2013 Author(s): S. Kutuzov, M. Shahgedanova, V. Mikhalenko, P. Ginot, I. Lavrentiev, and S. Kemp The first record of dust deposition events on Mt. Elbrus, Caucasus Mountains derived from a snow pit and a shallow firn core is presented for the 2009–2012 period. A combination of isotopic analysis, SEVIRI red-green-blue composite imagery, MODIS atmospheric optical depth fields derived using the Deep Blue algorithm, air mass trajectories derived using the HYSPLIT model and analyses of meteorological data enabled identification of dust source regions with high temporal (hours) and spatial (ca. 20–100 km) resolution. Seventeen dust deposition events were detected; fourteen occurred in March–June, one in February and two in October. Four events originated in the Sahara, predominantly in northeastern Libya and eastern Algeria. Thirteen events originated in the Middle East, in the Syrian Desert and northern Mesopotamia, from a mixture of natural and anthropogenic sources. Dust transportation from Sahara was associated with vigorous Saharan depressions, strong surface winds in the source region and mid-tropospheric southwesterly flow with daily winds speeds of 20–30 m s −1 at 700 hPa level. Although these events were less frequent than those originating in the Middle East, they resulted in higher dust concentrations in snow. Dust transportation from the Middle East was associated with weaker depressions forming over the source region, high pressure centred over or extending towards the Caspian Sea and a weaker southerly or southeasterly flow towards the Caucasus Mountains with daily wind speeds of 12–18 m s −1 at 700 hPa level. Higher concentrations of nitrates and ammonium characterised dust from the Middle East deposited on Mt. Elbrus in 2009 indicating contribution of anthropogenic sources. The modal values of particle size distributions ranged between 1.98 μm and 4.16 μm. Most samples were characterised by modal values of 2.0–2.8 μm with an average of 2.6 μm and there was no significant difference between dust from the Sahara and the Middle East.

Description: Changes in glacier equilibrium-line altitude in the western Alps from 1984 to 2010: evaluation by remote sensing and modeling of the morpho-topographic and climate controls The Cryosphere, 7, 1455-1471, 2013 Author(s): A. Rabatel, A. Letréguilly, J.-P. Dedieu, and N. Eckert We present time series of equilibrium-line altitude (ELA) measured from the end-of-summer snow line altitude computed using satellite images, for 43 glaciers in the western Alps over the 1984–2010 period. More than 120 satellite images acquired by Landsat, SPOT and ASTER were used. In parallel, changes in climate variables, summer cumulative positive degree days (CPDD) and winter precipitation, were analyzed over the same time period using 22 weather stations located inside and around the study area. Assuming a continuous linear trend over the study period: (1) the average ELA of the 43 glaciers increased by about 170 m; (2) summer CPDD increased by about 150 PDD at 3000 m a.s.l.; and (3) winter precipitation remained rather stationary. Summer CPDD showed homogeneous spatial and temporal variability; winter precipitation showed homogeneous temporal variability, but some stations showed a slightly different spatial pattern. Regarding ELAs, temporal variability between the 43 glaciers was also homogeneous, but spatially, glaciers in the southern part of the study area differed from glaciers in the northern part, mainly due to a different precipitation pattern. A sensitivity analysis of the ELAs to climate and morpho-topographic variables (elevation, aspect, latitude) highlighted the following: (1) the average ELA over the study period of each glacier is strongly controlled by morpho-topographic variables; and (2) the interannual variability of the ELA is strongly controlled by climate variables, with the observed increasing trend mainly driven by increasing temperatures, even if significant nonlinear, low-frequency fluctuations appear to be driven by winter precipitation anomalies. Finally, we used an expansion of Lliboutry's approach to reconstruct fluctuations in the ELA of any glacier of the study area with respect to morpho-topographic and climate variables, by quantifying their respective weight and the related uncertainties in a consistent manner within a hierarchical Bayesian framework. This method was tested and validated using the ELA measured on the satellite images.

Description: A note on the water budget of temperate glaciers The Cryosphere, 7, 1557-1564, 2013 Author(s): J. Oerlemans In this note, the total dissipative melting in temperate glaciers is studied. The analysis is based on the notion that the dissipation is determined by the loss of potential energy due to the downward motion of mass (ice, snow, meltwater and rain). A mathematical formulation of the dissipation is developed and applied to a simple glacier geometry. In the next step, meltwater production resulting from enhanced ice motion during a glacier surge is calculated. The amount of melt energy available follows directly from the lowering of the centre of gravity of the glacier. To illustrate the concept, schematic calculations are presented for a number of glaciers with different geometric characteristics. Typical dissipative melt rates, expressed as water-layer depth averaged over the glacier, range from a few centimetres per year for smaller glaciers to half a metre per year for Franz Josef Glacier, one of the most active glaciers in the world (in terms of mass turnover). The total generation of meltwater during a surge is typically half a metre. For Variegated Glacier a value of 70 cm is found, for Kongsvegen 20 cm. These values refer to water layer depth averaged over the entire glacier. The melt \textit{rate} depends on the duration of the surge. It is generally an order of magnitude greater than water production by `normal' dissipation. On the other hand, the additional basal melt rate during a surge is comparable in magnitude with the water input from meltwater and precipitation. This suggests that enhanced melting during a surge does not grossly change the total water budget of a glacier. Basal water generated by enhanced sliding is an important ingredient in many theories of glacier surges. It provides a positive feedback mechanism that actually makes the surge happen. The results found here suggest that this can only work if water generated by enhanced sliding accumulates in a part of the glacier base where surface meltwater and rain have no or very limited access. This finding seems compatible with the fact that, on many glaciers, surges are initiated in the lower accumulation zone.

Description: Spectral Aerosol Extinction Monitoring System (SÆMS): setup, observational products, and comparisons Atmospheric Measurement Techniques Discussions, 6, 8647-8677, 2013 Author(s): A. Skupin, A. Ansmann, R. Engelmann, and H. Baars A Spectral Aerosol Extinction Monitoring System (SÆMS) is presented that allows us to continuously measure the spectral extinction coefficient of atmospheric aerosol particles along an about 2.7 km long optical path at 30–50 m height above ground at Leipzig (51.3° N, 12.4° E), Germany. The fully automated instrument measures the ambient aerosol extinction coefficients from 300–1000 nm. The main goal of SÆMS observations are long-term studies of the relationship between particle extinction and relative humidity from below 40 % to almost 100 %. The setup is presented and observations (a case study and statistical results for 2009) are discussed in terms of time series of 550 nm particle optical depth, Ångström exponent, and particle size distribution retrieved from the spectrally resolved extinction. The SÆMS measurements are compared with simultaneously performed EARLINET lidar, AERONET photometer, and in situ aerosol observations of particle size distribution and related extinction coefficients at the roof of our institute. Consistency between the different measurements is found which corroborates the quality of the SÆMS observations.

Description: Validation of SCIAMACHY O 2 A band cloud heights using Cloudnet radar/lidar measurements Atmospheric Measurement Techniques Discussions, 6, 8603-8645, 2013 Author(s): P. Wang and P. Stammes For the first time two SCIAMACHY O 2 A band cloud height products are validated using ground-based radar/lidar measurements between January 2003 and December 2011. The products are the ESA Level 2 (L2) version 5.02 cloud top height and the FRESCO (Fast Retrieval Scheme for Clouds from the Oxygen A band) version 6 cloud height. The radar/lidar profiles are obtained at the Cloudnet sites of Cabauw and Lindenberg, and are averaged for one hour centered at the SCIAMACHY overpass time to achieve an optimal temporal and spatial match. In total we have about 220 cases of single layer clouds and 200 cases of multi-layer clouds. The FRESCO cloud height and ESA L2 cloud top height are compared with the Cloudnet cloud top height and Cloudnet cloud middle height. We find that the ESA L2 cloud top height has a better agreement with the Cloudnet cloud top height than the Cloudnet cloud middle height. The ESA L2 cloud top height is on average 0.44 km higher than the Cloudnet cloud top height, with a standard deviation of 3.07 km. The FRESCO cloud height is closer to the Cloudnet cloud middle height than the Cloudnet cloud top height. The mean difference between the FRESCO cloud height and the Cloudnet cloud middle height is −0.14 km with a standard deviation of 1.88 km. The SCIAMACHY cloud height products are further compared to the Cloudnet cloud top height and the Cloudnet cloud middle height in 1 km bins. For single layer clouds, the difference between the ESA L2 cloud top height and the Cloudnet cloud top height is less than 1 km for each cloud bin at 3–7 km, which is 24 % percent of the data. The difference between the FRESCO cloud height and the Cloudnet cloud middle height is less than 1 km for each cloud bin at 0–6 km, which is 85 % percent of the data. The results are similar for multi-layer clouds, but the percentage of cases having a bias within 1 km is smaller than for single layer clouds. Since globally about 60 % of all clouds are low clouds and 42 % are single-layer low clouds, we expect that globally for a large percentage of cases the FRESCO cloud height would be close to the cloud middle height.

Description: Characteristics of cloud liquid water path from SEVIRI on the Meteosat Second Generation 2 satellite for several cloud types Atmospheric Measurement Techniques Discussions, 6, 8743-8782, 2013 Author(s): A. Kniffka, M. Stengel, M. Lockhoff, R. Bennartz, and R. Hollmann In this study the temporal and spatial characteristics of liquid water path (LWP) of low, middle level and high clouds are analysed using space-based observations of the Spinning Enhanced Visible and Infrared Imager (SEVIRI) onboard the Meteosat Second Generation 2 (MSG2) satellite. Both geophysical quantities are part of the dataset CLAAS (CLoud property dAtAset using SEVIRI) and are generated by EUMETSAT's Satellite Application Facility on Climate Monitoring (CM SAF). In this article we focus on the statistical properties of LWP retrieved at daylight associated with the individual cloud type. Our results reveal that each cloud type possesses a characteristic LWP distribution. These frequency distributions are constant with time in the entire SEVIRI field of view, but vary for smaller regions like Central Europe. The average LWP is higher over land than over sea, in case of low clouds 15–27% for 2009 and the variance of the frequency distributions is enhanced. Also, the average diurnal cycle of LWP is related to cloud type where most pronounced diurnal variations were detected for middle level clouds. With SEVIRI it is possible to distinguish between intrinsic LWP variability and variations driven by cloud amount. The relative amplitude of the intrinsic diurnal cycle can exceed the cloud amount driven amplitude.

Description: Global glacier changes: a revised assessment of committed mass losses and sampling uncertainties The Cryosphere, 7, 1565-1577, 2013 Author(s): S. H. Mernild, W. H. Lipscomb, D. B. Bahr, V. Radić, and M. Zemp Most glaciers and ice caps (GIC) are out of balance with the current climate. To return to equilibrium, GIC must thin and retreat, losing additional mass and raising sea level. Because glacier observations are sparse and geographically biased, there is an undersampling problem common to all global assessments. Here, we further develop an assessment approach based on accumulation-area ratios (AAR) to estimate committed mass losses and analyze the undersampling problem. We compiled all available AAR observations for 144 GIC from 1971 to 2010, and found that most glaciers and ice caps are farther from balance than previously believed. Accounting for regional and global undersampling errors, our model suggests that GIC are committed to additional losses of 32 ± 12% of their area and 38 ± 16% of their volume if the future climate resembles the climate of the past decade. These losses imply global mean sea-level rise of 163 ± 69 mm, assuming total glacier volume of 430 mm sea-level equivalent. To reduce the large uncertainties in these projections, more long-term glacier measurements are needed in poorly sampled regions.

Description: Remote sensing of atmospheric trace gas columns: an efficient approach for regularization and calculation of total column averaging kernels Atmospheric Measurement Techniques Discussions, 6, 4999-5031, 2013 Author(s): T. Borsdorff, O. P. Hasekamp, A. Wassmann, and J. Landgraf A concept is proposed to retrieve the vertical column densities of atmospheric trace gases from remote sensing measurements. It combines the numerical simplicity of a least-squares profile scaling retrieval with the numerically robust calculation of the total column averaging kernel using an analytic expression. The approach enables calculation of the total column averaging kernel on arbitrary vertical grids. Formally, the proposed method is equivalent to Tikhonov regularization of the first kind with an infinite regularization strength. Due to its efficiency it is particularly suited for implementation in operational data processing with high demands on processing time. To demonstrate the method, we apply it to CO column retrieval from simulated measurements in the 2.3 μm spectral region and to O 3 column retrieval from the UV, which represents ideal measurements of a series of space-borne spectrometers like SCIAMACHY, TROPOMI, GOME, and GOME-2. For both spectral ranges, we consider clear-sky and cloudy scenes where clouds are modelled as an elevated Lambertian surface. Here, the smoothing error for the clear-sky and cloudy atmosphere is significant and reaches several percent, depending on the reference profile which is used for scaling. This underlines the importance of the column averaging kernel for a proper interpretation of retrieved column densities. Furthermore, we show that the total column smoothing error is affected by a discretization error when total column averaging kernels are not represented on a fine enough vertical grid. For both retrievals this effect becomes negligible by using a vertical grid with 20–40 equally thick layers between 0 and 50 km.

Description: Mixing layer height retrievals by multichannel microwave radiometer observations Atmospheric Measurement Techniques Discussions, 6, 4971-4998, 2013 Author(s): D. Cimini, F. De Angelis, J.-C. Dupont, S. Pal, and M. Haeffelin The mixing layer height (MLH) is a key parameter for boundary layer studies, including meteorology, air quality, and climate. MLH estimates are inferred from in situ radiosonde measurements or remote sensing observations from instruments like lidar, wind profiling radar, or sodar. Methods used to estimate MLH from radiosonde profiles are also used with atmospheric temperature and humidity profiles retrieved by microwave radiometers (MWR). This paper proposes an alternative approach to estimate MLH from MWR data, based on direct observations (brightness temperatures, Tb) instead of retrieved profiles. To our knowledge, MLH estimates directly from Tb observations has never been attempted before. The method consists of a multivariate linear regression trained with an a priori set of collocated MWR Tb observations (multi-frequency and multi-angle) and MLH estimates from a state-of-the-art lidar system. Results show that the method is able to follow both the diurnal cycle and the day-to-day variability as suggested by the lidar measurements, and also it can detect low MLH values that are below the full overlap limit (~ 200 m) of the lidar system used. Statistics of the comparison between MWR- and reference lidar-based MLH retrievals show mean difference within 10 m, RMS within 340 m, and correlation coefficient higher than 0.77. Monthly mean analysis for day-time MLH from MWR, lidar, and radiosonde shows consistent seasonal variability, peaking at ~ 1200–1400 m in June and decreasing down to ~ 600 m in October. Conversely, night-time monthly mean MLH from all methods are within 300–500 m without any significant seasonal variability. The proposed method provides results that are more consistent with radiosonde estimates than MLH estimates from MWR retrieved profiles. MLH monthly mean values agree well within 1 std with bulk Richardson number method applied at radiosonde profiles at 11:00 and 23:00 UTC. The method described herewith operates continuously and it is expected to work with analogous performances for the entire diurnal cycle, except during considerable precipitation, demonstrating new potential for atmospheric observation by ground-based microwave radiometry.

Description: Stratospheric aerosol particle size information in Odin-OSIRIS limb scatter spectra Atmospheric Measurement Techniques Discussions, 6, 5065-5099, 2013 Author(s): L. A. Rieger, A. E. Bourassa, and D. A. Degenstein The Optical Spectrograph and InfraRed Imaging System (OSIRIS) on-board the Odin satellite has now taken over a decade of limb scatter measurements that have been used to retrieve the Version 5 stratospheric aerosol extinction product. This product is retrieved using a representative particle size distribution to calculate scattering cross sections and scattering phase functions for the forward model calculations. In this work the information content of OSIRIS measurements with respect to stratospheric aerosol is systematically examined for the purpose of retrieving particle size information along with the extinction coefficient. The benefit of using measurements at different wavelengths and scattering angles in the retrieval is studied and it is found that incorporation of the 1530 nm radiance measurement is key for a robust retrieval of particle size information. It is also found that using OSIRIS measurements at different solar geometries simultaneously provides little additional benefit. Based on these results, an improved aerosol retrieval algorithm is developed that couples the retrieval of aerosol extinction and mode radius of a log-normal particle size distribution. Comparison of these results with coincident measurements from SAGE III show agreement in retrieved extinction to within approximately 10% over the bulk of the aerosol layer, which is comparable to Version 5. The retrieved particle size, when converted to Ångström coefficient, shows good qualitative agreement with SAGE II measurements made at somewhat shorter wavelengths.

Description: SAGE version 7.0 algorithm: application to SAGE II Atmospheric Measurement Techniques Discussions, 6, 5101-5171, 2013 Author(s): R. P. Damadeo, J. M. Zawodny, L. W. Thomason, and N. Iyer This paper details the SAGE version 7.0 algorithm and how it is applied to SAGE II. Changes made between the previous (v6.2) and current (v7.0) versions are described and their impacts on the data products explained for both coincident event comparisons and time-series analysis. Users of the data will notice a general improvement in all of the SAGE II data products, which are now in better agreement with more modern data sets (e.g. SAGE III) and more robust for use with trend studies.

Description: Recent changes in spring snowmelt timing in the Yukon River basin detected by passive microwave satellite data The Cryosphere, 7, 905-916, 2013 Author(s): K. A. Semmens and J. M. Ramage Spring melt is a significant feature of high latitude snowmelt dominated drainage basins influencing hydrological and ecological processes such as snowmelt runoff and green-up. Melt duration, defined as the transition period from snowmelt onset until the end of the melt refreeze, is characterized by high diurnal amplitude variations (DAV) where the snowpack is melting during the day and refreezing at night, after which the snowpack melts constantly until depletion. Determining trends for this critical period is necessary for understanding how the Arctic is changing with rising temperatures and provides a baseline from which to assess future change. To study this dynamic period, brightness temperature ( T b ) data from the Special Sensor Microwave Imager (SSM/I) 37 V-GHz frequency from 1988 to 2010 were used to assess snowmelt timing trends for the Yukon River basin, Alaska/Canada. Annual T b and DAV for 1434 Equal-Area Scalable Earth (EASE)-Grid pixels (25 km resolution) were processed to determine melt onset and melt refreeze dates from T b and DAV thresholds previously established in the region. Temporal and spatial trends in the timing of melt onset and melt refreeze, and the duration of melt were analyzed for the 13 sub-basins of the Yukon River basin with three different time interval approaches. Results show a lengthening of the melt period for the majority of the sub-basins with a significant trend toward later end of melt refreeze after which the snowpack melts day and night leading to snow clearance, peak discharge, and green-up. Earlier melt onset trends were also found in the higher elevations and northernmost sub-basins (Porcupine, Chandalar, and Koyukuk rivers). Latitude and elevation displayed the dominant controls on melt timing variability and spring solar flux was highly correlated with melt timing in middle (∼600–1600 m) elevations.

Description: Seasonal controls on snow distribution and aerial ablation at the snow-patch and landscape scales, McMurdo Dry Valleys, Antarctica The Cryosphere, 7, 917-931, 2013 Author(s): J. W. Eveland, M. N. Gooseff, D. J. Lampkin, J. E. Barrett, and C. D. Takacs-Vesbach Accumulated snow in the McMurdo Dry Valleys, while limited, has great ecological significance to subnivian soil environments. Though sublimation dominates the ablation process in this region, measurable increases in soil moisture and insulation from temperature extremes provide more favorable conditions with respect to subnivian soil communities. While precipitation is not substantial, significant amounts of snow can accumulate, via wind transport, in topographic lees along the valley bottoms, forming thousands of discontinuous snow patches. These patches have the potential to act as significant sources of local meltwater, controlling biogeochemical cycling and the landscape distribution of microbial communities. Therefore, determining the spatial and temporal dynamics of snow at multiple scales is imperative to understanding the broader ecological role of snow in this region. High-resolution satellite imagery acquired during the 2009–2010 and 2010–2011 austral summers was used to quantify the distribution of snow across Taylor and Wright valleys. Extracted snow-covered area from the imagery was used as the basis for assessing inter-annual variability and seasonal controls on accumulation and ablation of snow at multiple scales. In addition to landscape analyses, fifteen 1 km 2 plots (3 in each of 5 study regions) were selected to assess the prevalence of snow cover at finer spatial scales, referred to herein as the snow-patch scale. Results confirm that snow patches tend to form in the same locations each year with some minor deviations observed. At the snow-patch scale, neighboring patches often exhibit considerable differences in aerial ablation rates, and particular snow patches do not reflect trends for snow-covered area observed at the landscape scale. These differences are presumably related to microtopographic influences acting on individual snow patches, such as wind sheltering and differences in snow depth due to the underlying topography. This highlights the importance of both the landscape and snow-patch scales in assessing the effects of snow cover on biogeochemical cycling and microbial communities.

Description: Technical Note: Aeolian dust proxies produce visible luminescence upon intense laser-illumination that results from incandescence of internally mixed carbon Atmospheric Measurement Techniques Discussions, 6, 5173-5194, 2013 Author(s): L. Ma, T. Cao, and J. E. Thompson Mineral dust mimics dispersed in air produced visible luminescence between 550–800 nm when illuminated with a high peak power (MW range) Nd:YAG laser beam at 532 or 1064 nm. The luminescence persists for a few microseconds after the laser pulse and the measured emission spectrum is roughly consistent with a blackbody emitter at ≈4300 K. Both observations are consistent with assigning laser-induced incandescence (LII) as the source of the luminescence. However, light emission intensity from the mineral dust proxies is 240–4600 less intense than incandescence from fresh kerosene soot on a per-mass basis at laser pulse energies 90% on average. Heating to 350 °C reduced emission by 45–72%. Since black carbon soot and char (BC) oxidizes at elevated temperatures and BC is known to be present in soils, we conclude emission of light from the mineral dust aerosol proxies is likely a result of black carbon or char internally mixed within the soil dust sample. The reduction in LII response for samples heated to temperatures of 250–350 °C may result from partial oxidation of BC, but alternatively, could implicate a role for carbon present within organic molecules. The study suggests laser-induced incandescence measurements may allow quantitation of black carbon in soils and that soil dust is not truly an interferent in BC analysis by LII, but rather, a BC containing material.

Description: Ice-shelf buttressing and the stability of marine ice sheets The Cryosphere, 7, 647-655, 2013 Author(s): G. H. Gudmundsson Ice-shelf buttressing and the stability of marine-type ice sheets are investigated numerically. Buttressing effects are analysed for a situation where a stable grounding line is located on a bed sloping upwards in the direction of flow. Such grounding-line positions are known to be unconditionally unstable in the absence of transverse flow variations. It is shown that ice-shelf buttressing can restore stability under these conditions. Ice flux at the grounding line is, in general, not a monotonically increasing function of ice thickness. This, possibly at first somewhat counterintuitive result, is found to be fully consistent with recent theoretical work. Grounding lines on retrograde slopes are conditionally stable, and the stability regime is a non-trivial function of bed and ice-shelf geometry. The stability of grounding lines cannot be assessed from considerations of local bed slope only.

Description: Evidence and analysis of 2012 Greenland records from spaceborne observations, a regional climate model and reanalysis data The Cryosphere, 7, 615-630, 2013 Author(s): M. Tedesco, X. Fettweis, T. Mote, J. Wahr, P. Alexander, J. E. Box, and B. Wouters A combined analysis of remote sensing observations, regional climate model (RCM) outputs and reanalysis data over the Greenland ice sheet provides evidence that multiple records were set during summer 2012. Melt extent was the largest in the satellite era (extending up to ∼97% of the ice sheet) and melting lasted up to ∼2 months longer than the 1979–2011 mean. Model results indicate that near surface temperature was ∼3 standard deviations (σ) above the 1958–2011 mean, while surface mass balance (SMB) was ∼3σ below the mean and runoff was 3.9σ above the mean over the same period. Albedo, exposure of bare ice and surface mass balance also set new records, as did the total mass balance with summer and annual mass changes of, respectively, −627 Gt and −574 Gt, 2σ below the 2003–2012 mean. We identify persistent anticyclonic conditions over Greenland associated with anomalies in the North Atlantic Oscillation (NAO), changes in surface conditions (e.g., albedo, surface temperature) and preconditioning of surface properties from recent extreme melting as major driving mechanisms for the 2012 records. Less positive if not increasingly negative SMB will likely occur should these characteristics persist.

Description: Retrieval of aerosol parameters from the oxygen A band in the presence of chlorophyll fluorescence Atmospheric Measurement Techniques Discussions, 6, 3181-3213, 2013 Author(s): A. F. J. Sanders and J. F. de Haan We have investigated precision of retrieved parameters for a generic aerosol retrieval algorithm over vegetated land using the O 2 A band. Chlorophyll fluorescence is taken into account in the forward model. Fluorescence emissions are modeled as isotropic contributions to the upwelling radiance field at the surface and they are retrieved along with aerosol parameters. Precision is calculated by propagating measurement noise using the forward model's derivatives. We assume that measurement noise is dominated by shot noise; thus, results apply to grating spectrometers in particular. In a number of retrieval simulations, we describe precision for various atmospheric states, observation geometries and spectral resolutions of the instrument. Our results show that aerosol optical thickness, aerosol pressure, fluorescence emission and surface albedo can be simultaneously retrieved from the O 2 A band. We also show that most of the fluorescence signal is provided by filling-in of the O 2 A band and to a lesser extent by filling-in of Fraunhofer lines.

Description: Balanced conditions or slight mass gain of glaciers in the Lahaul and Spiti region (northern India, Himalaya) during the nineties preceded recent mass loss The Cryosphere, 7, 569-582, 2013 Author(s): C. Vincent, Al. Ramanathan, P. Wagnon, D. P. Dobhal, A. Linda, E. Berthier, P. Sharma, Y. Arnaud, M. F. Azam, P. G. Jose, and J. Gardelle The volume change of the Chhota Shigri Glacier (India, 32° 20 N, 77° 30' E) between 1988 and 2010 has been determined using in situ geodetic measurements. This glacier has experienced only a slight mass loss between 1988 and 2010 (−3.8 ± 2.0 m w.e. (water equivalent) corresponding to −0.17 ± 0.09 m w.e. yr −1 ). Using satellite digital elevation models (DEM) differencing and field measurements, we measure a negative mass balance (MB) between 1999 and 2010 (−4.8 ± 1.8 m w.e. corresponding to −0.44 ± 0.16 m w.e. yr −1 ). Thus, we deduce a slightly positive or near-zero MB between 1988 and 1999 (+1.0 ± 2.7 m w.e. corresponding to +0.09 ± 0.24 m w.e. yr −1 ). Furthermore, satellite DEM differencing reveals that the MB of the Chhota Shigri Glacier (−0.39 ± 0.15 m w.e. yr −1 ) has been only slightly less negative than the MB of a 2110 km 2 glaciarized area in the Lahaul and Spiti region (−0.44 ± 0.09 m w.e. yr −1 ) during 1999−2011. Hence, we conclude that the ice wastage is probably moderate in this region over the last 22 yr, with near equilibrium conditions during the nineties, and an ice mass loss after. The turning point from balanced to negative mass budget is not known but lies probably in the late nineties and at the latest in 1999. This positive or near-zero MB for Chhota Shigri Glacier (and probably for the surrounding glaciers of the Lahaul and Spiti region) during at least part of the 1990s contrasts with a recent compilation of MB data in the Himalayan range that indicated ice wastage since 1975. However, in agreement with this compilation, we confirm more negative balances since the beginning of the 21st century.

Description: McClear: a new model estimating downwelling solar radiation at ground level in clear-sky conditions Atmospheric Measurement Techniques Discussions, 6, 3367-3405, 2013 Author(s): M. Lefèvre, A. Oumbe, P. Blanc, B. Espinar, B. Gschwind, Z. Qu, L. Wald, M. Schroedter-Homscheidt, C. Hoyer-Klick, A. Arola, A. Benedetti, J. W. Kaiser, and J.-J. Morcrette A new fast clear-sky model called McClear was developed to estimate the downwelling shortwave direct and global irradiances received at ground level under clear skies. McClear implements a fully physical modelling replacing empirical relations or simpler models used before. It exploits the recent results on aerosol properties, and total column content in water vapor and ozone produced by the MACC project (Monitoring Atmosphere Composition and Climate). It accurately reproduces the irradiance computed by the libRadtran reference radiative transfer model with a computational speed approximately 10 5 times greater by adopting the abaci, or look-up tables, approach combined with interpolation functions. It is therefore suited for geostationary satellite retrievals or numerical weather prediction schemes with many pixels or grid points, respectively. McClear irradiances were compared to 1 min measurements made in clear-sky conditions in several stations within the Baseline Surface Radiation Network in various climates. For global, respectively direct, irradiance, the correlation coefficient ranges between 0.95 and 0.99, resp. 0.86 and 0.99. The bias is comprised between −14 and 25 W m −2 , resp. −49 and +33 W m −2 . The RMSE ranges between 20 W m −2 (3% of the mean observed irradiance) and 36 W m −2 (5%), resp. 33 W m −2 (5%) and 64 W m −2 (10%). These results are much better than those from state-of-the-art models. This work demonstrates the quality of the McClear model combined with MACC products, and indirectly the quality of the aerosol properties modeled by the MACC reanalysis.

Description: Experimental quantification of contact freezing in an electrodynamic balance Atmospheric Measurement Techniques Discussions, 6, 3407-3437, 2013 Author(s): N. Hoffmann, A. Kiselev, D. Rzesanke, D. Duft, and T. Leisner Heterogeneous nucleation of ice in a supercooled water droplet induced by an external contact with a dry aerosol particle has long been known to be more effective than freezing induced by the same nucleus immersed in the droplet. However, the experimental quantification of contact freezing is challenging. Here we report an experimental method allowing to determine the temperature dependent ice nucleation probability of size selected aerosol particles. The method uses supercooled charged water droplets suspended in a laminar flow of air containing aerosol particles as contact freezing nuclei. The rate of droplet–particle collisions is calculated numerically with account for Coulomb attraction, drag force and induced dipole interaction between charged droplet and aerosol particles. The calculation is verified by direct counting of aerosol particles collected by a levitated droplet. By repeating the experiment on individual droplets for a sufficient number of times, we are able to reproduce the statistical freezing behavior of a large ensemble of supercooled droplets and measure the average rate of freezing events. The freezing rate is equal to the product of the droplet–particle collision rate and the probability of freezing on a single contact, the latter being a function of temperature, size and composition of the contact ice nuclei. Based on these observations, we show that for the types of particles investigated so far, contact freezing is the dominating freezing mechanism on the time scale of our experiment.

Description: The geomorphological effect of cornice fall avalanches in the Longyeardalen valley, Svalbard The Cryosphere, 7, 1361-1374, 2013 Author(s): M. Eckerstorfer, H. H. Christiansen, L. Rubensdotter, and S. Vogel The study of snow avalanches and their geomorphological effect in the periglacial parts of the cryosphere is important for enhanced geomorphological process understanding as well as hazard-related studies. Only a few field studies, and particularly few in the High Arctic, have quantified avalanche sedimentation. Snow avalanches are traditionally ranked behind rockfall in terms of their significance for mass-wasting processes of rockslopes. Cornice fall avalanches are at present the most dominant snow avalanche type at two slope systems, called Nybyen and Larsbreen, in the valley Longyeardalen in central Svalbard. Both slope systems are on northwest-facing lee slopes underneath a large summit plateau, with annual cornices forming on the top. High-frequency and magnitude cornice fall avalanching is observed by daily automatic time-lapse photography. In addition, rock debris sedimentation by cornice fall avalanches was measured directly in permanent sediment traps or by snow inventories. The results from a maximum of seven years of measurements in a total of 13 catchments show maximum mean rock debris sedimentation rates ranging from 8.2 to 38.7 kg m −2 at Nybyen, and from 0.8 to 55.4 kg m −2 at Larsbreen. Correspondingly, avalanche fan surfaces accreted from 2.6 to 8.8 mm yr −1 at Nybyen, and from 0.2 to 13.9 mm yr −1 at Larsbreen. This comparably efficient rockslope mass wasting is due to collapsing cornices producing cornice fall avalanches containing large amounts of rock debris throughout the entire winter. The rock debris of different origin stems from the plateau crests, the adjacent free rock face and the transport pathway, accumulating distinct avalanche fans at both slope systems. Cornice fall avalanche sedimentation also contributed to the development of a rock glacier at the Larsbreen site during the Holocene. We have recorded present maximum rockwall retreat rates of 0.9 mm yr −1 at Nybyen, but as much as 6.7 mm yr −1 at Larsbreen, while average Holocene rockwall retreat rates of 1.1 mm yr −1 at Nybyen have been determined earlier. As cornice fall avalanches are the dominant type of avalanche in central Svalbard, the related geomorphological effect is assumed to be of significance at periglacial landscape scale. A climate-induced shift in prevailing winter wind direction could change the rockslope sedimentation effectively by changing the snow avalanche activity.

Description: Separating mixtures of aerosol types in airborne High Spectral Resolution Lidar data Atmospheric Measurement Techniques Discussions, 6, 8269-8309, 2013 Author(s): S. P. Burton, M. A. Vaughan, R. A. Ferrare, and C. A. Hostetler Knowledge of aerosol type is important for source attribution and for determining the magnitude and assessing the consequences of aerosol radiative forcing. However, atmospheric aerosol is frequently not a single pure type, but instead occurs as a mixture of types, and this mixing affects the optical and radiative properties of the aerosol. This paper extends the work of earlier researchers by using the aerosol intensive parameters measured by the NASA Langley Research Center airborne High Spectral Resolution Lidar (HSRL-1) to develop a comprehensive and unified set of rules for characterizing the external mixing of several key aerosol intensive parameters: extinction-to-backscatter ratio (i.e. lidar ratio), backscatter color ratio, and depolarization ratio. We present the mixing rules in a particularly simple form that leads easily to mixing rules for the covariance matrices that describe aerosol distributions, rather than just scalar values of measured parameters. These rules can be applied to infer mixing ratios from the lidar-observed aerosol parameters, even for cases without significant depolarization. We demonstrate our technique with measurement curtains from three HSRL-1 flights which exhibit mixing between two aerosol types, urban pollution plus dust, marine plus dust, and smoke plus marine. For these cases, we infer a time-height cross-section of mixing ratio along the flight track, and partition aerosol extinction into portions attributed to the two pure types.

Description: Lidar-based remote sensing of atmospheric boundary layer height over land and ocean Atmospheric Measurement Techniques Discussions, 6, 8311-8338, 2013 Author(s): T. Luo, R. Yuan, and Z. Wang Atmospheric boundary layer (ABL) processes are important in climate, weather and air quality. A better understanding of the structure and the behavior of the ABL is required for understanding and modeling of the chemistry and dynamics of the atmosphere on all scales. Based on the systematic variations of ABL structures over different surfaces, different lidar-based methods were developed and evaluated to determine the boundary layer height and mixing layer height over land and ocean. With Atmospheric Radiation Measurement Program (ARM) Climate Research Facility (ACRF) micropulse lidar (MPL) and radiosonde measurements, diurnal and season cycles of atmospheric boundary layer depth and ABL vertical structure over ocean (TWP_C2 cite) and land (SGP_C1) are analyzed. The new methods are also applied to satellite lidar measurements. The derived global marine boundary layer structure database shows good agreement with marine ABL stratiform cloud top height.

Description: Fast reconstruction of hyperspectral radiative transfer simulations by using small spectral subsets: application to the oxygen A band Atmospheric Measurement Techniques Discussions, 6, 8339-8370, 2013 Author(s): A. Hollstein and R. Lindstrot Hyperspectral radiative transfer simulations are a versatile tool in remote sensing but can pose a major computational burden. We describe a simple method to construct hyperspectral simulation results by using only a small spectral subsample of the simulated wavelength range, thus leading to major speedups in such simulations. This is achieved by computing principal components for a small number of representative hyperspectral spectra and then deriving a reconstruction matrix for a specific spectral subset of channels to compute the hyperspectral data. The method is applied and discussed in detail using the example of top of atmosphere radiances in the oxygen A band, leading to speedups in the range of one to two orders of magnitude when compared to radiative transfer simulations at full spectral resolution.

Description: An assessment of cloud top thermodynamic phase products obtained from A-Train passive and active sensors Atmospheric Measurement Techniques Discussions, 6, 8371-8411, 2013 Author(s): S. Zeng, J. Riedi, F. Parol, C. Cornet, and F. Thieuleux The A-Train observations provide an unprecedented opportunity for the production of high quality dataset describing cloud properties. We illustrate in this study the use of one year of coincident POLDER (Polarization and Directionality of the Earth Reflectance), MODIS (MODerate Resolution Imaging Spectroradiometer) and CALIOP (Cloud-Aerosol Lidar with Orthogonal Polarization) observations to establish a reference dataset for the description of cloud top thermodynamic phase at global scale. We present the results of an extensive comparison between POLDER and MODIS cloud top phase products and discuss those in view of cloud vertical structure and optical properties derived simultaneously from collocated CALIOP active measurements. These results allow to identify and quantify potential biases present in the 3 considered dataset. Among those, we discuss the impacts of observation geometry, thin cirrus in multilayered and single layered cloud systems, supercooled liquid droplets, aerosols, fractional cloud cover and snow/ice or bright surfaces on global statistics of cloud phase derived from POLDER and MODIS passive measurements. Based on these analysis we define criteria for the selection of high confidence cloud phase retrievals which in turn can serve for the establishment of a reference cloud phase product. This high confidence joint product derived from POLDER/PARASOL and MODIS/Aqua can be used in the future as a benchmark for the evaluation of other cloud climatologies, for the assessment of cloud phase representation in models and the development of better cloud phase parametrization in the general circulation models (GCMs).

Description: Heterogeneity in glacier response in the upper Shyok valley, northeast Karakoram The Cryosphere, 7, 1385-1398, 2013 Author(s): R. Bhambri, T. Bolch, P. Kawishwar, D. P. Dobhal, D. Srivastava, and B. Pratap Glaciers in the Karakoram show long-term irregular behaviour with comparatively frequent and sudden advances. A glacier inventory of the upper Shyok valley situated in northeast Karakoram has been generated for the year 2002 using Landsat ETM+ and SRTM3 DEM as baseline data for the investigations and subsequent change analysis. The upper Shyok valley contained 2123 glaciers (larger than 0.02 km 2 in size) with an area of 2977.9 ± 95.3 km 2 in 2002. Out of these, 18 glaciers with an area of 1004.1 ± 32.1 km 2 showed surge-type behaviour. Change analysis based on Hexagon KH-9 (years 1973 and 1974) and Landsat TM/ETM+ (years 1989, 2002 and 2011) images had to be restricted to a subset of 136 glaciers (covering an area of 1609.7 ± 51.5 km 2 in 2002) due to adverse snow conditions. The area of the investigated glaciers, including the 18 surge-type glaciers identified, showed no significant changes during all studied periods. However, the analysis provides a hint that the overall glacier area slightly decreased until about 1989 (area 1973: 1613.6 ± 43.6 km 2 ; area 1989: 1602.0 ± 33.6 km 2 ) followed by an increase (area 2002: 1609.7 ± 51.5; area 2011: 1615.8 ± 35.5 km 2 ). Although the overall change in area is insignificant, advances in glacier tongues since the end of the 1980s are clearly visible. Detailed estimations of length changes for individual glaciers since the 1970s and for Central Rimo Glacier since the 1930s confirm the irregular retreat and advance.

Description: Odin-OSIRIS detection of the Chelyabinsk meteor Atmospheric Measurement Techniques Discussions, 6, 8435-8443, 2013 Author(s): L. A. Rieger, A. E. Bourassa, and D. A. Degenstein On 15 February 2013 an 11 000 ton meteor entered Earth's atmosphere south east of Chelyabinsk creating a large fireball at 23 km altitude. The resulting stratospheric aerosol loading was detected by the Ozone Mapping and Profiler Suite (OMPS) in a high altitude polar belt. This work confirms the presence and lifetime of the stratospheric debris using the Optical Spectrograph and InfraRed Imaging System (OSIRIS) onboard the Odin satellite. Although OSIRIS coverage begins in mid-March, the measurements show a belt of enhanced scattering near 35 km altitude between 50° N and 70° N. Initially, enhancements show increased scattering of up to 15% over the background conditions, decaying in intensity and dropping in altitude until they are indistinguishable from background conditions by mid-May.

Description: Weak layer fracture: facets and depth hoar The Cryosphere, 7, 1447-1453, 2013 Author(s): I. Reiweger and J. Schweizer Understanding failure initiation within weak snow layers is essential for modeling and predicting dry-snow slab avalanches. We therefore performed laboratory experiments with snow samples containing a weak layer consisting of either faceted crystals or depth hoar. During these experiments the samples were loaded with different loading rates and at various tilt angles until fracture. The strength of the samples decreased with increasing loading rate and increasing tilt angle. Additionally, we took pictures of the side of four samples with a high-speed video camera and calculated the displacement using a particle image velocimetry (PIV) algorithm. The fracture process within the weak layer could thus be observed in detail. Catastrophic failure started due to a shear fracture just above the interface between the depth hoar layer and the underlying crust.

Description: Scatter of mass changes estimates at basin scale for Greenland and Antarctica The Cryosphere, 7, 1411-1432, 2013 Author(s): V. R. Barletta, L. S. Sørensen, and R. Forsberg During the last decade, the GRACE mission has provided valuable data for determining the mass changes of the Greenland and Antarctic ice sheets. Yet, discrepancies still exist in the published mass balance results, and comprehensive analyses on the sources of errors and discrepancies are lacking. Here, we present monthly mass changes together with trends derived from GRACE data at basin scale for both the Greenland and Antarctic ice sheets, and we assess the variability and errors for each of the possible sources of discrepancies, and we do this in an unprecedented systematic way, taking into account mass inference methods, data sets and background models. We find a very good agreement between the monthly mass change results derived from two independent methods, which represents a cross validation. For the monthly solutions, we find that most of the scatter is caused by the use of the two different data sets rather than the two different methods applied. Besides the well-known GIA trend uncertainty, we find that the geocenter motion and the recent de-aliasing corrections significantly impact the trends, with contributions of +13.2 Gt yr −1 and −20 Gt yr −1 , respectively, for Antarctica, which is more affected by these than Greenland. We show differences between the use of release RL04 and the new RL05 and confirm a lower noise content in the new release. The overall scatter of the solutions well exceeds the uncertainties propagated from the data errors and the leakage (as done in the past); hence we calculate new sound total errors for the monthly solutions and the trends. We find that the scatter in the monthly solutions caused by applying different estimates of geocenter motion time series (degree-1 corrections) is significant – contributing with up to 40% of the total error. For the whole GRACE period (2003–2011) our trend estimate for Greenland is −234 ± 20 Gt yr −1 and −83 ± 36 Gt yr −1 for Antarctica (−111 ± 15 Gt yr −1 in the western part). We also find a clear (with respect to our errors) increase of mass loss in the last four years.

Description: Probabilistic approach to cloud and snow detection on AVHRR imagery Atmospheric Measurement Techniques Discussions, 6, 8445-8507, 2013 Author(s): J. P. Musial, F. Hüsler, M. Sütterlin, C. Neuhaus, and S. Wunderle The derivation of probability estimates complementary to geophysical data sets has gained special attention over the last years. The information about a confidence level of provided physical quantities is required to construct an error budget of higher level products and to correctly interpret final results of a particular analysis. Regarding the generation of products based on satellite data the common input consists of a cloud mask which allows discrimination between surface and cloud signals. Further the surface information is divided between snow and snow-free components. At any step of this discrimination process a misclassification in a cloud/snow mask propagates to higher level products and may alter their usability. Within this scope a novel Probabilistic Cloud Mask (PCM) algorithm suited for the 1×1 km Advanced Very High Resolution Radiometer (AVHRR) data is proposed which provides three types of probability estimates between: cloudy/clear-sky, cloudy/snow and clear-sky/snow conditions. As opposed to the majority of available techniques which are usually based on a decision-tree approach in the PCM algorithm all spectral, angular and ancillary information is used in a single step to retrieve the probability estimates from the pre-computed Look Up Tables (LUTs). Moreover, the issue of derivation of a single threshold value for a spectral test was overcome by the concept of multidimensional information space which is divided into small bins by an extensive set of thresholds. The discrimination between snow and ice clouds and detection of broken, thin clouds was enhanced by means of the Invariant Coordinate System (ICS) transformation. The study area covers a wide range of environmental conditions spanning from Iceland through central Europe to northern parts of Africa which exhibit diverse difficulties for cloud/snow masking algorithms. The retrieved PCM cloud classification was compared to the PPSv2012 and MOD35 collection 6 cloud masks, SYNOP weather reports, CALIPSO vertical feature mask version 3 and to MOD10A1 collection 5 snow mask. The outcomes of conducted analyses proved fine detection skills of the PCM method with comparable or better results than the reference PPS algorithm.

Description: Retrieval techniques for airborne imaging of methane concentrations using high spatial and moderate spectral resolution: application to AVIRIS Atmospheric Measurement Techniques Discussions, 6, 8543-8588, 2013 Author(s): A. K. Thorpe, C. Frankenberg, and D. A. Roberts Two quantitative retrieval techniques were evaluated to estimate methane (CH 4 ) enhancement in concentrated plumes using high spatial and moderate spectral resolution data from the Airborne Visible/Infrared Imaging Spectrometer (AVIRIS). An Iterative Maximum a Posteriori Differential Optical Absorption Spectroscopy (IMAP-DOAS) algorithm performed well for an ocean scene containing natural CH 4 emissions from the Coal Oil Point (COP) seep field near Santa Barbara, California. IMAP-DOAS retrieval precision errors are expected to equal between 0.31 to 0.61 ppm CH 4 over the lowest atmospheric layer (height up to 1.04 km), corresponding to about a 30 to 60 ppm error for a 10 m thick plume. However, IMAP-DOAS results for a terrestrial scene were adveresly influenced by the underlying landcover. A hybrid approach using Singular Value Decomposition (SVD) was particularly effective for terrestrial surfaces because it could better account for spectral variability in surface reflectance. Using this approach, a CH 4 plume was observed immediately downwind of two hydrocarbon storage tanks at the Inglewood Oil Field in Los Angeles, California, with a maximum near surface enhancement of 8.45 ppm above background. At COP, the distinct plume had a maximum enhancement of 2.85 ppm CH 4 above background and was consistent with known seep locations and local wind direction. A sensitivity analysis also indicates CH 4 sensitivity should be more than doubled for the next generation AVIRIS sensor (AVIRISng) due to improved spectral resolution and sampling. AVIRIS-like sensors offer the potential to better constrain emissions on local and regional scales, including sources of increasing concern like industrial point source emissions and fugitive CH 4 from the oil and gas industry.

Description: Effects of solar activity and geomagnetic field on noise in CALIOP profiles above the South Atlantic Anomaly Atmospheric Measurement Techniques Discussions, 6, 8589-8602, 2013 Author(s): V. Noel, H. Chepfer, C. Hoareau, M. Reverdy, and G. Cesana By documenting noise levels in 6.5 yr of nighttime measurements by the spaceborne lidar CALIOP above the South Atlantic Anomaly (SAA), we show they contain information about the evolution of upwelling high-energy radiation levels in the area. We find the amount of noisy profiles is influenced by the 11 yr cycle of solar activity, fluctuates by ±5% between 2006 and 2013, and is anticorrelated with solar activity with a 1 yr lag. The size of the SAA grows as solar activity decreases, and an overall westward shift of the SAA region is detectable. We predict SAA noise levels will increase anew after 2014, and will affect future spaceborne lidar missions most near 2020. In other areas, supposedly unaffected by incoming sunlight, nighttime noise levels are much weaker but follow the same 11 yr cycle, superimposed with a one-year cycle that affects both hemispheres similarly and could be attributed to geomagnetic activity.

Description: Are seasonal calving dynamics forced by buttressing from ice mélange or undercutting by melting? Outcomes from full-Stokes simulations of Store Glacier, West Greenland The Cryosphere, 8, 2353-2365, 2014 Author(s): J. Todd and P. Christoffersen We use a full-Stokes 2-D model (Elmer/Ice) to investigate the flow and calving dynamics of Store Glacier, a fast-flowing outlet glacier in West Greenland. Based on a new, subgrid-scale implementation of the crevasse depth calving criterion, we perform two sets of simulations: one to identify the primary forcing mechanisms and another to constrain future stability. We find that the mixture of icebergs and sea ice, known as ice mélange or sikussak, is principally responsible for the observed seasonal advance of the ice front. On the other hand, the effect of submarine melting on the calving rate of Store Glacier appears to be limited. Sensitivity analysis demonstrates that the glacier's calving dynamics are sensitive to seasonal perturbation, but are stable on interannual timescales due to the strong topographic control on the flow regime. Our results shed light on the dynamics of calving glaciers and may help explain why neighbouring glaciers do not necessarily respond synchronously to changes in atmospheric and oceanic forcing.

Description: Deglaciation of the Caucasus Mountains, Russia/Georgia, in the 21st century observed with ASTER satellite imagery and aerial photography The Cryosphere, 8, 2367-2379, 2014 Author(s): M. Shahgedanova, G. Nosenko, S. Kutuzov, O. Rototaeva, and T. Khromova Changes in the map area of 498 glaciers located on the Main Caucasus ridge (MCR) and on Mt. Elbrus in the Greater Caucasus Mountains (Russia and Georgia) were assessed using multispectral ASTER and panchromatic Landsat imagery with 15 m spatial resolution in 1999/2001 and 2010/2012. Changes in recession rates of glacier snouts between 1987–2001 and 2001–2010 were investigated using aerial photography and ASTER imagery for a sub-sample of 44 glaciers. In total, glacier area decreased by 4.7 ± 2.1% or 19.2 ± 8.7 km 2 from 407.3 ± 5.4 km 2 to 388.1 ± 5.2 km 2 . Glaciers located in the central and western MCR lost 13.4 ± 7.3 km 2 (4.7 ± 2.5%) in total or 8.5 km 2 (5.0 ± 2.4%) and 4.9 km 2 (4.1 ± 2.7%) respectively. Glaciers on Mt. Elbrus, although located at higher elevations, lost 5.8 ± 1.4 km 2 (4.9 ± 1.2%) of their total area. The recession rates of valley glacier termini increased between 1987–2000/01 and 2000/01–2010 (2000 for the western MCR and 2001 for the central MCR and Mt.~Elbrus) from 3.8 ± 0.8, 3.2 ± 0.9 and 8.3 ± 0.8 m yr −1 to 11.9 ± 1.1, 8.7 ± 1.1 and 14.1 ± 1.1 m yr −1 in the central and western MCR and on Mt. Elbrus respectively. The highest rate of increase in glacier termini retreat was registered on the southern slope of the central MCR where it has tripled. A positive trend in summer temperatures forced glacier recession, and strong positive temperature anomalies in 1998, 2006, and 2010 contributed to the enhanced loss of ice. An increase in accumulation season precipitation observed in the northern MCR since the mid-1980s has not compensated for the effects of summer warming while the negative precipitation anomalies, observed on the southern slope of the central MCR in the 1990s, resulted in stronger glacier wastage.

Description: Glacier-like forms on Mars The Cryosphere, 8, 2047-2061, 2014 Author(s): B. Hubbard, C. Souness, and S. Brough More than 1300 glacier-like forms (GLFs) are located in Mars' mid-latitudes. These GLFs are predominantly composed of ice–dust mixtures and are visually similar to terrestrial valley glaciers, showing signs of downhill viscous deformation and an expanded former extent. However, several fundamental aspects of their behavior are virtually unknown, including temporal and spatial variations in mass balance, ice motion, landscape erosion and deposition, and hydrology. Here, we investigate the physical glaciology of martian GLFs. We use satellite images of specific examples and case studies to build on existing knowledge relating to (i) GLF current and former extent, exemplified via a GLF located in Phlegra Montes; (ii) indicators of GLF motion, focusing on the presence of surface crevasses on several GLFs; (iii) processes of GLF debris transfer, focusing on mapping and interpreting boulder trains on one GLF located in Protonilus Mensae, the analysis of which suggests a best-estimate mean GLF flow speed of 7.5 mm a −1 ; and (iv) GLF hydrology, focusing on supra-GLF gulley networks. On the basis of this information, we summarize the current state of knowledge of the glaciology of martian GLFs and identify future research avenues.

Description: Modeling the elastic transmission of tidal stresses to great distances inland in channelized ice streams The Cryosphere, 8, 2007-2029, 2014 Author(s): J. Thompson, M. Simons, and V. C. Tsai Geodetic surveys suggest that ocean tides can modulate the motion of Antarctic ice streams, even at stations many tens of kilometers inland from the grounding line. These surveys suggest that ocean tidal stresses can perturb ice stream motion at distances about an order of magnitude farther inland than tidal flexure of the ice stream alone. Recent models exploring the role of tidal perturbations in basal shear stress are primarily one- or two-dimensional, with the impact of the ice stream margins either ignored or parameterized. Here, we use two- and three-dimensional finite-element modeling to investigate transmission of tidal stresses in ice streams and the impact of considering more realistic, three-dimensional ice stream geometries. Using Rutford Ice Stream as a real-world comparison, we demonstrate that the assumption that elastic tidal stresses in ice streams propagate large distances inland fails for channelized glaciers due to an intrinsic, exponential decay in the stress caused by resistance at the ice stream margins. This behavior is independent of basal conditions beneath the ice stream and cannot be fit to observations using either elastic or nonlinear viscoelastic rheologies without nearly complete decoupling of the ice stream from its lateral margins. Our results suggest that a mechanism external to the ice stream is necessary to explain the tidal modulation of stresses far upstream of the grounding line for narrow ice streams. We propose a hydrologic model based on time-dependent variability in till strength to explain transmission of tidal stresses inland of the grounding line. This conceptual model can reproduce observations from Rutford Ice Stream.

Description: Fluctuations of a Greenlandic tidewater glacier driven by changes in atmospheric forcing: observations and modelling of Kangiata Nunaata Sermia, 1859–present The Cryosphere, 8, 2031-2045, 2014 Author(s): J. M. Lea, D. W. F. Mair, F. M. Nick, B. R. Rea, D. van As, M. Morlighem, P. W. Nienow, and A. Weidick Many tidewater glaciers in Greenland are known to have undergone significant retreat during the last century following their Little Ice Age maxima. Where it is possible to reconstruct glacier change over this period, they provide excellent records for comparison to climate records, as well as calibration/validation for numerical models. These glacier change records therefore allow for tests of numerical models that seek to simulate tidewater glacier behaviour over multi-decadal to centennial timescales. Here we present a detailed record of behaviour from Kangiata Nunaata Sermia (KNS), SW Greenland, between 1859 and 2012, and compare it against available oceanographic and atmospheric temperature data between 1871 and 2012. We also use these records to evaluate the ability of a well-established one-dimensional flow-band model to replicate behaviour for the observation period. The record of terminus change demonstrates that KNS has advanced/retreated in phase with atmosphere and ocean climate anomalies averaged over multi-annual to decadal timescales. Results from an ensemble of model runs demonstrate that observed dynamics can be replicated. Model runs that provide a reasonable match to observations always require a significant atmospheric forcing component, but do not necessarily require an oceanic forcing component. Although the importance of oceanic forcing cannot be discounted, these results demonstrate that changes in atmospheric forcing are likely to be a primary driver of the terminus fluctuations of KNS from 1859 to 2012. We propose that the detail and length of the record presented makes KNS an ideal site for model validation exercises investigating links between climate, calving rates, and tidewater glacier dynamics.

Description: Potential of the TROPOspheric Monitoring Instrument (TROPOMI) onboard the Sentinel-5 Precursor for the monitoring of terrestrial chlorophyll fluorescence Atmospheric Measurement Techniques Discussions, 7, 12545-12588, 2014 Author(s): L. Guanter, I. Aben, P. Tol, J. M. Krijger, A. Hollstein, P. Köhler, A. Damm, J. Joiner, C. Frankenberg, and J. Landgraf Global monitoring of sun-induced chlorophyll fluorescence (SIF) can improve our knowledge about the photosynthetic functioning of terrestrial ecosystems. The feasibility of SIF retrievals from spaceborne atmospheric spectrometers has been demonstrated by a number of studies in the last years. In this work, we investigate the potential of the upcoming TROPOspheric Monitoring Instrument (TROPOMI) onboard the Sentinel-5 Precursor satellite mission for SIF retrieval. TROPOMI will sample the 675–775 nm spectral window with a spectral resolution of 0.5 nm and a pixel size of 7 km × 7 km. We use an extensive set of simulated TROPOMI data in order to assess the uncertainty of single SIF retrievals and subsequent spatio-temporal composites. Our results illustrate the enormous improvement in SIF monitoring achievable with TROPOMI with respect to comparable spectrometers currently in-flight, such as the Global Ozone Monitoring Experiment-2 (GOME-2) instrument. We find that TROPOMI can reduce global uncertainties in SIF mapping by more than a factor 2 with respect to GOME-2, which comes together with an about 5-fold improvement in spatial sampling. Finally, we discuss the potential of TROPOMI to accurately map other important vegetation parameters, such as leaf photosynthetic pigments and proxies for canopy structure, which will complement SIF retrievals for a self-contained description of vegetation condition and functioning.

Description: The Level 2 research product algorithms for the Superconducting Submillimeter-Wave Limb-Emission Sounder (SMILES) Atmospheric Measurement Techniques Discussions, 4, 3593-3645, 2011 Author(s): P. Baron, J. Urban, H. Sagawa, J. Möller, D. P. Murtagh, J. Mendrok, E. Dupuy, T. O. Sato, S. Ochiai, K. Suzuki, T. Manabe, T. Nishibori, K. Kikuchi, R. Sato, M. Takayanagi, Y. Murayama, M. Shiotani, and Y. Kasai This paper describes the algorithms of the level-2 research (L2r) processing chain developed for the Superconducting Submillimeter-Wave Limb-Emission Sounder (SMILES). The chain has been developed in parallel to the operational chain for conducting researches on calibration and retrieval algorithms. L2r chain products are available to the scientific community. The objective of version 2 is the retrieval of the vertical distribution of trace gases in the altitude range of 18–90 km. An theoretical error analysis is conducted to estimate the retrieval feasibility of key parameters of the processing: line-of-sight elevation tangent altitudes (or angles), temperature and O 3 profiles. The line-of-sight tangent altitudes are retrieved between 20 and 50 km from the strong ozone (O 3 ) line at 625.371 GHz, with low correlation with the O 3 volume-mixing ratio and temperature retrieved profiles. Neglecting the non-linearity of the radiometric gain in the calibration procedure is the main systematic error. It is large for the retrieved temperature (between 5–10 K). Therefore, atmospheric pressure can not be derived from the retrieved temperature, and, then, in the altitude range where the line-of-sight tangent altitudes are retrieved, the retrieved trace gases profiles are found to be better represented on pressure levels than on altitude levels. The error analysis for the retrieved HOCl profile demonstrates that best results for inverting weak lines can be obtained by using narrow spectral windows. Future versions of the L2r algorithms will improve the temperature/pressure retrievals and also provide information in the upper tropospheric/lower stratospheric region (e.g., water vapor, ice content, O 3 ) and on stratospheric and mesospheric line-of-sight winds.

Description: Near infrared nadir sounding of vertical column densities: methodology and application to SCIAMACHY Atmospheric Measurement Techniques Discussions, 4, 3685-3737, 2011 Author(s): S. Gimeno García, F. Schreier, G. Lichtenberg, and S. Slijkhuis Nadir observations with the shortwave infrared channels of SCIAMACHY onboard the ENVISAT satellite can be used to derive information on atmospheric gases such as CO, CH 4 , N 2 O, CO 2 , and H 2 O. For the operational level 1b–2 processing of SCIAMACHY data a new retrieval code BIRRA (Beer InfraRed Retrieval Algorithm) has been developed: BIRRA performs a nonlinear least squares fit of the measured radiance, where molecular concentration vertical profiles are scaled to fit the observed data. Here we present the forward modeling (radiative transfer) and inversion (least squares optimization) fundamentals of the code along with the further processing steps required to generate higher level products such as global distributions and time series. Moreover, various aspects of level 1 (observed spectra) and auxiliary input data relevant for successful retrievals are discussed. BIRRA is currently used for operational analysis of carbon monoxide vertical column densities from SCIAMACHY channel 8 observations, and is being prepared for methane retrievals using channel 6 spectra. A set of representative CO retrievals and first CH 4 results are presented to demonstrate BIRRA's capabilities.

Description: Volatilizable biogenic organic compounds (VBOCs) with two dimensional gas chromatography-time of flight mass spectrometry (GC × GC-TOFMS): sampling methods, VBOC complexity, and chromatographic retention data Atmospheric Measurement Techniques Discussions, 4, 3647-3684, 2011 Author(s): J. F. Pankow, W. Luo, A. N. Melnychenko, K. C. Barsanti, L. M. Isabelle, C. Chen, A. B. Guenther, and T. N. Rosenstiel Two dimensional gas chromatography (GC × GC) with detection by time-of-flight mass spectrometry (TOFMS) was applied in the rapid analysis of air samples containing highly complex mixtures of volatilizable biogenic organic compounds (VBOCs). VBOC analytical methodologies are briefly reviewed, and optimal conditions are discussed for sampling with both adsorption/thermal desorption (ATD) cartridges and solid-phase microextraction (SPME) fibers. Air samples containing VBOC emissions from leaves of two tree species ( Cedrus atlantica and Calycolpus moritzianus ) were obtained by both ATD and SPME. The optimized gas chromatographic conditions utilized a 45 m, 0.25 mm I.D. low-polarity primary column (DB-VRX, 1.4 μm film) and a 1.5 m, 0.25 mm I.D. polar secondary column (Stabilwax® 0.25 μm film). Excellent separation was achieved in a 36 min temperature programmed GC × GC chromatogram. Thousands of VBOC peaks were present in the sample chromatograms; hundreds of tentative identifications by NIST mass spectral matching are provided. Very few of the tentatively identified compounds are currently available as authentic standards. Method detection limit values for a 5 l ATD sample were 3.5 pptv (10 ng m −3 ) for isoprene, methyl vinyl ketone, and methacrolein, and ~1.5 pptv (~10 ng m −3 ) for monoterpenes and sesquiterpenes. Kovats-type chromatographic retention index values on the primary column and relative retention time values on the secondary column are provided for 21 standard compounds and for 417 tentatively identified VBOCs. 19 of the 21 authentic standard compounds were found in one of the Cedrus atlantica SPME samples. In addition, easily quantifiable levels of at least 13 sesquiterpenes were found in an ATD sample obtained from a branch enclosure of Calycolpus moritzianus . Overall, the results obtained via GC × GC-TOFMS highlight an extreme, and largely uncharacterized diversity of VBOCs, consistent with the hypothesis that sesquiterpenes and other compounds beyond the current list of typically determined VBOC analytes may well be important contributors to global atmospheric levels of organic particulate matter.

Description: Correcting orbital drift signal in the time series of AVHRR derived convective cloud fraction using rotated empirical orthogonal function Atmospheric Measurement Techniques Discussions, 4, 3877-3890, 2011 Author(s): A. Devasthale, K. Karlsson, J. Quaas, and H. Grassl The AVHRRs instruments onboard the series of NOAA satellites offer the longest available meteorological data records from space. These satellites have drifted in orbit resulting in shifts in the local time sampling during the life span of sensors onboard. Depending on the amplitude of a diurnal cycle of the geophysical parameters derived, orbital drift may cause spurious trends in their time series. We investigate tropical deep convective clouds, which show pronounced diurnal cycle amplitude, to bracket an upper bound of the impact of orbital drift on their time series. We carry out a rotated empirical orthogonal function analysis and show that the REOFs are useful in delineating orbital drift signal and, more importantly, in correcting this signal in the time series of convective cloud amount. These results will help facilitate the derivation of homogenized data series of cloud amount from NOAA satellite sensors and ultimately analyzing trends from them. However, we suggest detailed comparison of various methods and their rigorous testing before applying final orbital drift corrections.

Description: Sensitivity studies for a space-based methane lidar mission Atmospheric Measurement Techniques Discussions, 4, 3545-3592, 2011 Author(s): C. Kiemle, M. Quatrevalet, G. Ehret, A. Amediek, A. Fix, and M. Wirth Methane is the third most important greenhouse gas in the atmosphere after water vapour and carbon dioxide. A major handicap to quantify the emissions at the Earth's surface in order to better understand biosphere-atmosphere exchange processes and potential climate feedbacks is the lack of accurate and global observations of methane. Space-based integrated path differential absorption (IPDA) lidar has potential to fill this gap, and a Methane Remote Lidar Mission (MERLIN) on a small satellite in Polar orbit was proposed by DLR and CNES in the frame of a German-French climate monitoring initiative. System simulations are used to identify key performance parameters and to find an advantageous instrument configuration, given the environmental, technological, and budget constraints. The sensitivity studies use representative averages of the atmospheric and surface state to estimate the measurement precision, i.e. the random uncertainty due to instrument noise. Key performance parameters for MERLIN are average laser power, telescope size, orbit height, surface reflectance, and detector noise. A modest-size lidar instrument with 0.45 W average laser power and 0.55 m telescope diameter on a 506 km orbit could provide 50-km averaged methane column measurement along the sub-satellite track with a precision of about 1 % over vegetation. The use of a methane absorption trough at 1.65 μm improves the near-surface measurement sensitivity and vastly relaxes the wavelength stability requirement that was identified as one of the major technological risks in the pre-phase A studies for A-SCOPE, a space-based IPDA lidar for carbon dioxide at the European Space Agency. Minimal humidity and temperature sensitivity at this wavelength position will enable accurate measurements in tropical wetlands, key regions with largely uncertain methane emissions. In contrast to actual passive remote sensors, measurements in Polar Regions will be possible and biases due to aerosol layers and thin ice clouds will be minimised.

Description: Information operator approach applied to the retrieval of the vertical distribution of atmospheric constituents from ground-based high-resolution FTIR measurements Atmospheric Measurement Techniques Discussions, 4, 3739-3785, 2011 Author(s): C. Senten, M. De Mazière, G. Vanhaelewyn, and C. Vigouroux The analysis of high spectral resolution Fourier Transform infrared (FTIR) solar absorption spectra is an important issue in remote sensing. If this is done carefully, one can obtain information, not only about the total column abundances, but also about the vertical distribution of various constituents in the atmosphere. This work introduces the application of the information operator approach for extracting vertical profile information from ground-based FTIR measurements. The algorithm is implemented and tested within the well-known retrieval code SFIT2, adapting the optimal estimation method such as to take into account only the significant contributions to the solution. In particular, we demonstrate the feasibility of the method in an application to ground-based FTIR spectra taken in the frame of the Network for the Detection of Atmospheric Composition Change (NDACC) at Ile de La Réunion (21° S, 55° E). A thorough comparison is made between the original optimal estimation method and this alternative retrieval algorithm, regarding information content, retrieval robustness and corresponding full error budget evaluation for the target species ozone (O 3 ), nitrous oxide (N 2 O), methane (CH 4 ), and carbon monoxide (CO). For O 3 and CH 4 , a comparison with the Tikhonov regularization method has also been made. It is shown that the information operator approach performs well and in most cases yields both a better accuracy and stability than the optimal estimation method. Additionally, the information operator approach has the advantage of being less sensitive to the choice of a priori information. The Tikhonov regularization results seem to be situated between both methods' results, as to profile retrievals, error budgets and column stability.

Description: Opportunistic validation of sulfur dioxide in the Sarychev Peak volcanic eruption cloud Atmospheric Measurement Techniques Discussions, 4, 3861-3875, 2011 Author(s): S. A. Carn and T. M. Lopez We report attempted validation of Ozone Monitoring Instrument (OMI) sulfur dioxide (SO 2 ) retrievals in the stratospheric volcanic cloud from Sarychev Peak (Kurile Islands) in June 2009, through opportunistic deployment of a ground-based ultraviolet (UV) spectrometer (FLYSPEC) as the volcanic cloud drifted over Central Alaska. The volcanic cloud altitude (~12–14 km) was constrained using coincident CALIPSO lidar observations. By invoking some assumptions about the spatial distribution of SO 2 , we derive averages of FLYSPEC vertical SO 2 columns for comparison with OMI SO 2 measurements. Despite limited data, we find minimum OMI-FLYSPEC differences of ~5–6 % which support the validity of the operational OMI SO 2 algorithm. These measurements represent the first attempt to validate SO 2 in a stratospheric volcanic cloud using a mobile ground-based instrument, and demonstrate the need for a network of rapidly deployable instruments for validation of space-based volcanic SO 2 measurements.

Description: Inversion of tropospheric profiles of aerosol extinction and HCHO and NO 2 mixing ratios from MAX-DOAS observations in Milano during the summer of 2003 and comparison with independent data sets Atmospheric Measurement Techniques Discussions, 4, 3891-3964, 2011 Author(s): T. Wagner, S. Beirle, T. Brauers, T. Deutschmann, U. Frieß, C. Hak, J. D. Halla, K. P. Heue, W. Junkermann, X. Li, U. Platt, and I. Pundt-Gruber We present aerosol and trace gas profiles derived from MAX-DOAS observations. Our inversion scheme is based on simple profile parameterisations used as input for an atmospheric radiative transfer model (forward model). From a least squares fit of the forward model to the MAX-DOAS measurements, two profile parameters are retrieved including integrated quantities (aerosol optical depth or trace gas vertical column density), and parameters describing the height and shape of the respective profiles. From these results, the aerosol extinction and trace gas mixing ratios can also be calculated. We apply the profile inversion to MAX-DOAS observations during a measurement campaign in Milano, Italy, September 2003, which allowed simultaneous observations from three telescopes (directed to north, west, south). Profile inversions for aerosols and trace gases were possible on 23 days. Especially in the middle of the campaign (17–20 September 2003), enhanced values of aerosol optical depth and NO 2 and HCHO mixing ratios were found. The retrieved layer heights were typically similar for HCHO and aerosols. For NO 2 , lower layer heights were found, which increased during the day. The MAX-DOAS inversion results are compared to independent measurements: (1) aerosol optical depth measured at an AERONET station at Ispra; (2) near-surface NO 2 and HCHO (formaldehyde) mixing ratios measured by long path DOAS and Hantzsch instruments at Bresso; (3) vertical profiles of HCHO and aerosols measured by an ultra light aircraft. Depending on the viewing direction, the aerosol optical depths from MAX-DOAS are either smaller or larger than those from AERONET observations. Similar comparison results are found for the MAX-DOAS NO 2 mixing ratios versus long path DOAS measurements. In contrast, the MAX-DOAS HCHO mixing ratios are generally higher than those from long path DOAS or Hantzsch instruments. The comparison of the HCHO and aerosol profiles from the aircraft showed reasonable agreement with the respective MAX-DOAS layer heights. From the comparison of the results for the different telescopes, it was possible to investigate the internal consistency of the MAX-DOAS observations. As part of our study, a cloud classification algorithm was developed (based on the MAX-DOAS zenith viewing directions), and the effects of clouds on the profile inversion were investigated. Different effects of clouds on aerosols and trace gas retrievals were found: while the aerosol optical depth is systematically underestimated and the HCHO mixing ratio is systematically overestimated under cloudy conditions, the NO 2 mixing ratios are only slightly affected. These findings are in basic agreement with radiative transfer simulations.

Description: Modeling the ascent of sounding balloons: derivation of the vertical air motion Atmospheric Measurement Techniques Discussions, 4, 3965-4012, 2011 Author(s): A. Gallice, F. G. Wienhold, C. R. Hoyle, F. Immler, and T. Peter A new model to describe the ascent of sounding balloons in the troposphere and lower stratosphere (up to ~30–35 km altitude) is presented. Contrary to previous models, detailed account is taken of both the variation of the drag coefficient with altitude and the heat imbalance between the balloon and the atmosphere. To compensate for the lack of data on the drag coefficient of sounding balloons, a reference curve for the relationship between drag coefficient and Reynolds number is derived from a dataset of flights launched during the Lindenberg Upper Air Methods Intercomparisons (LUAMI) campaign. The transfer of heat from the surrounding air into the balloon is accounted for by solving the radial heat diffusion equation inside the balloon. The potential applications of the model include the forecast of the trajectory of sounding balloons, which can be used to increase the accuracy of the match technique, and the derivation of the air vertical velocity. The latter is obtained by subtracting the ascent rate of the balloon in still air calculated by the model from the actual ascent rate. This technique is shown to provide an approximation for the vertical air motion with an uncertainty error of 0.5 m s −1 in the troposphere and 0.2 m s −1 in the stratosphere. An example of extraction of the air vertical velocity is provided in this paper. We show that the air vertical velocities derived from the balloon soundings in this paper are in general agreement with small-scale atmospheric velocity fluctuations related to gravity waves, mechanical turbulence, or other small-scale air motions measured during the SUCCESS campaign (Subsonic Aircraft: Contrail and Cloud Effects Special Study) in the orographically unperturbed mid-latitude middle troposphere.

Description: The surface energy balance of a polygonal tundra site in northern Siberia – Part 2: Winter The Cryosphere, 5, 509-524, 2011 Author(s): M. Langer, S. Westermann, S. Muster, K. Piel, and J. Boike In this study, we present the winter time surface energy balance at a polygonal tundra site in northern Siberia based on independent measurements of the net radiation, the sensible heat flux and the ground heat flux from two winter seasons. The latent heat flux is inferred from measurements of the atmospheric turbulence characteristics and a model approach. The long-wave radiation is found to be the dominant factor in the surface energy balance. The radiative losses are balanced to about 60 % by the ground heat flux and almost 40 % by the sensible heat fluxes, whereas the contribution of the latent heat flux is small. The main controlling factors of the surface energy budget are the snow cover, the cloudiness and the soil temperature gradient. Large spatial differences in the surface energy balance are observed between tundra soils and a small pond. The ground heat flux released at a freezing pond is by a factor of two higher compared to the freezing soil, whereas large differences in net radiation between the pond and soil are only observed at the end of the winter period. Differences in the surface energy balance between the two winter seasons are found to be related to differences in snow depth and cloud cover which strongly affect the temperature evolution and the freeze-up at the investigated pond.

Description: Potential and limitations of the MAX-DOAS method to retrieve the vertical distribution of tropospheric nitrogen dioxide Atmospheric Measurement Techniques Discussions, 4, 4013-4072, 2011 Author(s): T. Vlemmix, A. J. M. Piters, A. J. C. Berkhout, L. F. L. Gast, P. Wang, and P. F. Levelt Muliple Axis Differential Optical Absorption Spectroscopy (MAX-DOAS) instruments can measure from the ground the absorption by nitrogen dioxide (NO 2 ) of scattered sunlight seen in multiple viewing directions. This paper studies the potential of this technique to derive the vertical distribution of NO 2 in the troposphere. Such profile information is essential in validation studies in which MAX-DOAS retrievals play a role. The retrieval algorithm used is based on a pre-calculated look-up table and assumes homogeneous mixing of aerosols and NO 2 in layers extending from the surface to a variable height. Two retrieval models are compared: one including and one excluding an elevated NO 2 layer at a fixed altitude in the free troposphere. An ensemble technique is applied to derive retrieved model uncertainties. Sensitivity studies demonstrate that MAX-DOAS based retrievals can make a distinction between an NO 2 layer that extends from the surface to a certain height (having a constant mixing ratio, or a mixing ratio that decreases with altitude) and an elevated NO 2 layer. The height of the elevated NO 2 layer can only be retrieved accurately when the aerosol extinction profile is known and the measurement noise is low. The uncertainty in this elevated NO 2 layer height provides the main source of uncertainty in the retrieval of the free tropospheric contribution to the tropospheric NO 2 column. A comparison was performed with independent data, based on observations done at the CINDI campaign, held in the Netherlands in 2009. Comparison with lidar partial tropospheric NO 2 columns showed a correlation of 0.78, and an average difference of 0.1× 10 15 molec cm −2 . The diurnal evolution of the NO 2 volume mixing ratio measured by in-situ monitors at the surface and averaged over five days with cloud-free mornings, compares quite well to the MAX-DOAS retrieval: a correlation was found of 0.8, and an average difference of 0.2 ppb.

Description: Water isotopic ratios from a continuously melted ice core sample Atmospheric Measurement Techniques Discussions, 4, 4073-4104, 2011 Author(s): V. Gkinis, T. J. Popp, T. Blunier, M. Bigler, S. Schüpbach, and S. J. Johnsen A new technique for on-line high resolution isotopic analysis of liquid water, tailored for ice core studies is presented. We build an interface between an Infra Red Cavity Ring Down Spectrometer (IR-CRDS) and a Continuous Flow Analysis (CFA) system. The system offers the possibility to perform simultaneuous water isotopic analysis of δ 18 O and δD on a continuous stream of liquid water as generated from a continuously melted ice rod. Injection of sub μl amounts of liquid water is achieved by pumping sample through a fused silica capillary and instantaneously vaporizing it with 100 % efficiency in a home made oven at a temperature of 170 °C. A calibration procedure allows for proper reporting of the data on the VSMOW scale. We apply the necessary corrections based on the assessed performance of the system regarding instrumental drifts and dependance on humidity levels. The melt rates are monitored in order to assign a depth scale to the measured isotopic profiles. Application of spectral methods yields the combined uncertainty of the system at below 0.1 ‰ and 0.5 ‰ for δ 18 O and δD, respectively. This performance is comparable to that achieved with mass spectrometry. Dispersion of the sample in the transfer lines limits the resolution of the technique. In this work we investigate and assess these dispersion effects. By using an optimal filtering method we show how the measured profiles can be corrected for the smoothing effects resulting from the sample dispersion. Considering the significant advantages the technique offers, i.e. simultaneuous measurement of δ 18 O and δD, potentially in combination with chemical components that are traditionally measured on CFA systems, notable reduction on analysis time and power consumption, we consider it as an alternative to traditional isotope ratio mass spectrometry with the possibility to be deployed for field ice core studies. We present data acquired in the framework of the NEEM deep ice core drilling project in Greenland, during the 2010 field season.

Description: A full Stokes ice flow model for the vicinity of Dome Fuji, Antarctica, with induced anisotropy and fabric evolution The Cryosphere, 5, 495-508, 2011 Author(s): H. Seddik, R. Greve, T. Zwinger, and L. Placidi A three-dimensional, thermo-mechanically coupled ice flow model with induced anisotropy has been applied to a ~200 × 200 km domain around the Dome Fuji drill site, Antarctica. The model ("Elmer/Ice") is based on the open-source multi-physics package Elmer ( http://www.csc.fi/elmer/ ) and solves the full Stokes equations. Flow-induced anisotropy in ice is accounted for by an implementation of the C ontinuum-mechanical, A nisotropic F low model, based on an anisotropic F low E nhancement factor ("CAFFE model"). Steady-state simulations for present-day climate conditions are conducted. The main findings are: (i) the flow regime at Dome Fuji is a complex superposition of vertical compression, horizontal extension and bed-parallel shear; (ii) for an assumed geothermal heat flux of 60 mW m −2 the basal temperature at Dome Fuji reaches the pressure melting point and the basal melting rate is ~0.35 mm a −1 ; (iii) in agreement with observational data, the fabric shows a strong single maximum at Dome Fuji, and the age of the ice is decreased compared to an isotropic scenario; (iv) as a consequence of spatially variable basal melting conditions, the basal age tends to be smaller where the ice is thicker and larger where the ice is thinner. The latter result is of great relevance for the consideration of a future drill site in the area.

Description: A 3-D tomographic trajectory retrieval for the air-borne limb-imager GLORIA Atmospheric Measurement Techniques Discussions, 4, 3805-3859, 2011 Author(s): J. Ungermann, J. Blank, J. Lotz, K. Leppkes, T. Guggenmoser, M. Kaufmann, P. Preusse, U. Naumann, and M. Riese Infrared limb sounding from aircraft can provide 2-D curtains of multiple trace gas species. However, conventional limb sounders view perpendicular to the aircraft axis and are unable to resolve the observed airmass along their line-of-sight. GLORIA (Gimballed Limb Observer for Radiance Imaging of the Atmosphere) is a new remote sensing instrument able to adjust its horizontal view angle with respect to the aircraft flight direction from 45° to 135°. This will allow for tomographic measurements of mesoscale structures for a wide variety of atmospheric constituents. Many flights of the GLORIA instrument will not follow closed curves that allow measuring an airmass from all directions. Consequently, it is examined by means of simulations, what results can be expected from tomographic evaluation of measurements made during a straight flight. It is demonstrated that the achievable resolution and stability is enhanced compared to conventional retrievals. In a second step, it is shown that the incorporation of channels exhibiting different optical depth can greatly enhance the 3-D retrieval quality enabling the exploitation of previously unused spectral samples. A second problem for tomographic retrievals is that advection, which can be neglected for conventional retrievals, plays an important role for the time-scales involved in a tomographic measurement flight. This paper presents a method to diagnose the effect of a time-varying atmosphere on a 3-D retrieval and demonstrates an effective way to compensate for effects of advection by incorporating wind-fields from meteorological datasets as a priori information.

Description: High-resolution air quality monitoring from space: a fast retrieval scheme for CO from hyperspectral infrared measurements Atmospheric Measurement Techniques Discussions, 4, 3787-3803, 2011 Author(s): N. Smith, H.-L. Huang, E. Weisz, H. J. Annegarn, and R. B. Pierce The first results of the Fast Linear Inversion Trace gas System (FLITS) retrieval scheme are presented here for CO from IASI (Infrared Atmospheric Sounding Interferometer) measurements using RAQMS (Real time Air Quality Modelling System) as atmospheric background. FLITS is a simple linear inversion scheme with a stable performance that retrieves total column CO concentrations (molec cm −2 ) at single field-of-view (FOV) irrespective of cloud cover. A case study is presented here for a biomass burning plume over the Pacific on 29 March 2010. For each FOV a single tropospheric CO density, vertically integrated over 200–800 hPa, is retrieved with 12 channels in the spectral range 2050–2225 cm −1 . Despite variations in cloud cover and temperature, the degrees of freedom for signal (DFS) of the solution ranges between 0.8 and 0.95. In addition, the retrieval error is at least half the background error of 10 %, with dominant contribution from uncertainty in the measurement and temperature. With its stability and processing speed, FLITS meet two of the key requirements for operational processing. We conclude that the linear combination of space-borne measurements with a chemical transport model in the FLITS retrieval scheme holds potential for real-time air quality monitoring and evaluation of pollutant transport at high spatial resolution.

Description: Early in-flight detection of SO 2 via Differential Optical Absorption Spectroscopy: a feasible aviation safety measure to prevent potential encounters with volcanic plumes Atmospheric Measurement Techniques Discussions, 4, 2827-2881, 2011 Author(s): L. Vogel, B. Galle, C. Kern, H. Delgado Granados, V. Conde, P. Norman, S. Arellano, O. Landgren, P. Lübcke, J. M. Alvarez Nieves, L. Cárdenas Gonzáles, and U. Platt Volcanic ash constitutes a risk to aviation, mainly due to its ability to cause jet engines to fail. Other risks include the possibility of abrasion of windshields and potentially serious damage to avionic systems. These hazards have been widely recognized since the early 1980s, when volcanic ash provoked several incidents of engine failure in commercial aircraft. In addition to volcanic ash, volcanic gases also pose a threat. Prolonged and/or cumulative exposure to sulphur dioxide (SO 2 ) or sulphuric acid (H 2 SO 4 ) aerosols potentially affects e.g. windows, air frame and may cause permanent damage to engines. SO 2 receives most attention among the gas species commonly found in volcanic plumes because its presence above the lower troposphere is a clear proxy for a volcanic cloud and indicates that fine ash could also be present. Up to now, remote sensing of SO 2 via Differential Optical Absorption Spectroscopy (DOAS) in the ultraviolet spectral region has been used to measure volcanic clouds from ground based, airborne and satellite platforms. Attention has been given to volcanic emission strength, chemistry inside volcanic clouds and measurement procedures were adapted accordingly. Here we present a set of experimental and model results, highlighting the feasibility of DOAS to be used as an airborne early detection system of SO 2 in two spatial dimensions. In order to prove our new concept, simultaneous airborne and ground-based measurements of the plume of Popocatépetl volcano, Mexico, were conducted in April 2010. The plume extended at an altitude around 5250 m above sea level and was approached and traversed at the same altitude with several forward looking DOAS systems aboard an airplane. These DOAS systems measured SO 2 in the flight direction and at ± 40 mrad (2.3°) angles relative to it in both, horizontal and vertical directions. The approaches started at up to 25 km distance to the plume and SO 2 was measured at all times well above the detection limit. In combination with radiative transfer studies, this study indicates that an extended volcanic cloud with a concentration of 10 12 molecules cm −3 at typical flight levels of 10 km can be detected unambiguously at distances of up to 80 km away. This range provides enough time (approx. 5 min) for pilots to take action to avoid entering a volcanic cloud in the flight path, suggesting that this technique can be used as an effective aid to prevent dangerous aircraft encounters with potentially ash rich volcanic clouds.

Description: Catalytic oxidation of H 2 on platinum: a method for in situ calibration of hygrometers Atmospheric Measurement Techniques Discussions, 4, 3083-3095, 2011 Author(s): A. W. Rollins, T. D. Thornberry, R.-S. Gao, B. D. Hall, and D. W. Fahey Standard reference samples of water vapor suitable for in situ calibration of atmospheric hygrometers are not currently widespread, leading to difficulties in unifying the calibrations of these hygrometers and potentially contributing to measurement discrepancies. We describe and evaluate a system for reliably and quantitatively converting mixtures of H 2 in air to H 2 O on a heated platinum (Pt) surface, providing a compact, portable, adjustable source of water vapor. The technique is shown to be accurate and can be used to easily and predictably produce a wide range of water vapor concentrations (≈1 ppm−2 %) on demand. The result is a H 2 O standard that is suitable for in situ calibration of hygrometers, with an accuracy nearly that of the available H 2 standards (≈±2 %).

Description: A versatile and reproducible automatic injection system for liquid standard introduction: application to in-situ calibration Atmospheric Measurement Techniques Discussions, 4, 3233-3249, 2011 Author(s): G. Isaacman, N. M. Kreisberg, D. R. Worton, S. V. Hering, and A. H. Goldstein The quantitation of trace organic compounds in ambient organic aerosol is difficult due to the chemical complexity of these mixtures, but is needed to provide insight into their sources and formation processes. Compound-level characterization of organic aerosols is typically performed through sample collection followed by gas or liquid chromatography. With these methods, introduction of liquid standards has long been used as an effective means of quantifying trace compounds, but automating this technique for use with in-situ instrumentation has not previously been achieved. Here we develop an automatic injection system (AutoInject) for the introduction of liquids into a custom collection and analysis cell for improved quantitation in chromatographic measurements. The system consists of chilled reservoirs containing liquid standards from which a sample loop is loaded and then injected into the cell. The AutoInject is shown to have reproducibility over 106 injections with a relative standard deviation of 1.5 %, and have negligible injection-to-injection carryover. A 6-port selector allows injection of different liquid standards separately or simultaneously. Additionally, automatic injection of multiple sample loops is shown to generate a linear multi-point calibration curve. Tests conducted in this work focus on use with the T hermal desorption A erosol G as chromatograph (TAG), but the flexibility of the system allows it to be used for a variety of applications.

Description: Scale-dependent measurement and analysis of ground surface temperature variability in alpine terrain The Cryosphere, 5, 431-443, 2011 Author(s): S. Gubler, J. Fiddes, M. Keller, and S. Gruber Measurements of environmental variables are often used to validate and calibrate physically-based models. Depending on their application, the models are used at different scales, ranging from few meters to tens of kilometers. Environmental variables can vary strongly within the grid cells of these models. Validating a model with a single measurement is therefore delicate and susceptible to induce bias in further model applications. To address the question of uncertainty associated with scale in permafrost models, we present data of 390 spatially-distributed ground surface temperature measurements recorded in terrain of high topographic variability in the Swiss Alps. We illustrate a way to program, deploy and refind a large number of measurement devices efficiently, and present a strategy to reduce data loss reported in earlier studies. Data after the first year of deployment is presented. The measurements represent the variability of ground surface temperatures at two different scales ranging from few meters to some kilometers. On the coarser scale, the dependence of mean annual ground surface temperature on elevation, slope, aspect and ground cover type is modelled with a multiple linear regression model. Sampled mean annual ground surface temperatures vary from −4 °C to 5 °C within an area of approximately 16 km 2 subject to elevational differences of approximately 1000 m. The measurements also indicate that mean annual ground surface temperatures vary up to 6 °C (i.e., from −2 °C to 4 °C) even within an elevational band of 300 m. Furthermore, fine-scale variations can be high (up to 2.5 °C) at distances of less than 14 m in homogeneous terrain. The effect of this high variability of an environmental variable on model validation and applications in alpine regions is discussed.

Description: Spatial resolution of tropical terrestrial CO 2 fluxes inferred using space-borne column CO 2 sampled in different earth orbits: the role of spatial error correlations Atmospheric Measurement Techniques Discussions, 4, 3251-3276, 2011 Author(s): P. I. Palmer, L. Feng, and H. Bösch We use realistic numerical experiments to assess the sensitivity of 8-day CO 2 flux estimates, inferred from space-borne short-wave infrared measurements of column-averaged CO 2 dry air mixing ratio X CO 2 , to the choice of Earth observing orbit. We focus on three orbits: (1) a low-inclination circular orbit used by the NASA Tropical Rainfall Measuring Mission (TRMM); (2) a sun-synchronous orbit used by the Japanese Greenhouse Gases Observing SATellite (GOSAT) and proposed for the NASA Orbiting Carbon Observatory (OCO-2) instrument; and (3) a precessing orbit used by the International Space Station (ISS). For each orbit, we assume an instrument based on the specification of the OCO-2; for GOSAT we use the relevant instrument specification. Sun-synchronous orbits offer near global coverage within a few days but have implications for the density of clear-sky measurements. The TRMM and ISS orbits intensively sample tropical latitudes, with sun-lit clear-sky measurements evenly distributed between a.m./p.m. For a specified spatial resolution for inferred fluxes, we find there is a critical number of measurements beyond which there is a disproportionately small decrease in flux uncertainty. We also find that including spatial correlations for measurements and model errors (of length 300 km) reduces the effectiveness of high measurement density for flux estimation and so should be considered when deciding sampling strategies. We show that cloud-free data from the TRMM orbit generally can improve the spatial resolution of CO 2 fluxes achieved by OCO-2 over tropical South America, for example, from 950 km to 630 km, and that combining data from these low-inclination and sun-synchronous orbits have the potential to reduce this spatial length further. Decreasing the length of the error correlations to 50 km, reflecting anticipated future improvements to transport models, results in CO 2 flux estimates on spatial scales that approach those observed by regional aircraft.

Description: Validation of first chemistry mode retrieval results from new limb- imaging FTS GLORIA with correlative MIPAS-STR observations Atmospheric Measurement Techniques Discussions, 7, 12691-12717, 2014 Author(s): W. Woiwode, O. Suminska-Ebersoldt, H. Oelhaf, M. Höpfner, G. V. Belyaev, A. Ebersoldt, F. Friedl-Vallon, J.-U. Grooß, T. Gulde, M. Kaufmann, A. Kleinert, M. Krämer, E. Kretschmer, T. Kulessa, G. Maucher, T. Neubert, C. Piesch, P. Preusse, M. Riese, H. Rongen, C. Sartorius, G. Schardt, A. Schönfeld, D. Schuettemeyer, M. K. Sha, F. Stroh, J. Ungermann, C. M. Volk, and J. Orphal We report first chemistry mode retrieval results from the new airborne limb- imaging infrared FTS (Fourier transform spectrometer) GLORIA and comparisons with observations by the conventional airborne limb- scanning infrared FTS MIPAS-STR. For GLORIA, the flights aboard the high-altitude research aircraft M55 Geophysica during the ESSenCe campaign (ESa Sounder Campaign 2011) were the very first in field deployment after several years of development. The simultaneous observations of GLORIA and MIPAS-STR during the flight on 16 December 2011 inside the polar vortex and under the conditions of optically partially transparent polar stratospheric clouds (PSCs) provided us the unique opportunity to compare the observations by two different infrared FTS generations directly. The retrieval results of temperature, HNO 3 , O 3 , H 2 O, CFC-11 and CFC-12 show reasonable agreement of GLORIA with MIPAS-STR and collocated in-situ observations. For the horizontally binned hyperspectral limb-images, the GLORIA sampling outnumbered the horizontal cross-track sampling of MIPAS-STR by up to one order of magnitude. Depending on the target parameter, typical vertical resolutions of 0.5 to 2.0 km were obtained for GLORIA and are typically by factors of 2 to 4 better compared to MIPAS-STR. While the improvement of the performance, characterisation and data processing of GLORIA are subject of ongoing work, the presented first results already demonstrate the considerable gain in sampling and vertical resolution achieved with GLORIA.

Description: A new approach to mapping permafrost and change incorporating uncertainties in ground conditions and climate projections The Cryosphere, 8, 2177-2194, 2014 Author(s): Y. Zhang, I. Olthof, R. Fraser, and S. A. Wolfe Spatially detailed information on permafrost distribution and change with climate is important for land use planning, infrastructure development, and environmental assessments. However, the required soil and surficial geology maps in the North are coarse, and projected climate scenarios vary widely. Considering these uncertainties, we propose a new approach to mapping permafrost distribution and change by integrating remote sensing data, field measurements, and a process-based model. Land cover types from satellite imagery are used to capture the general land conditions and to improve the resolution of existing permafrost maps. For each land cover type, field observations are used to estimate the probabilities of different ground conditions. A process-based model is used to quantify the evolution of permafrost for each ground condition under three representative climate scenarios (low, medium, and high warming). From the model results, the probability of permafrost occurrence and the most likely permafrost conditions are determined. We apply this approach at 20 m resolution to a large area in Northwest Territories, Canada. Mapped permafrost conditions are in agreement with field observations and other studies. The data requirements, model robustness, and computation time are reasonable, and this approach may serve as a practical means to mapping permafrost and changes at high resolution in other regions.

Description: Comparison of nitric oxide measurements in the mesosphere and lower thermosphere from ACE-FTS, MIPAS, SCIAMACHY, and SMR Atmospheric Measurement Techniques Discussions, 7, 12735-12794, 2014 Author(s): S. Bender, M. Sinnhuber, T. von Clarmann, G. Stiller, B. Funke, M. López-Puertas, J. Urban, K. Pérot, K. A. Walker, and J. P. Burrows We compare the nitric oxide measurements in the mesosphere and lower thermosphere (60 to 150 km) from four instruments: ACE-FTS, MIPAS, SCIAMACHY, and SMR. We use the daily zonal mean data in that altitude range for the years 2004–2010 (ACE-FTS), 2005–2012 (MIPAS), 2008–2012 (SCIAMACHY), and 2003–2012 (SMR). We first compare the data qualitatively with respect to the morphology, focussing on the major features, and then compare the time series directly and quantitatively. In three geographical regions, we compare the vertical density profiles on coincident measurement days. Since none of the instruments delivers continuous daily measurements in this altitude region, we carried out a multi-linear regression analysis. This regression analysis considers annual and semi-annual variability in form of harmonic terms and inter-annual variability by responding linearly to the solar Lyman-α radiation index and the geomagnetic Kp index. This analysis helps to find similarities and differences in the individual data sets with respect to the inter-annual variations caused by geomagnetic and solar variability. We find that the data sets are consistent and that they only disagree on minor aspects. SMR and ACE-FTS deliver the longest time series in the mesosphere and they both agree remarkably well. The shorter time series from MIPAS and SCIAMACHY also agree with them where they overlap. The data agree within ten to twenty percent when the number densities are large, but they can differ by 50 to 100% in some cases.

Description: GPS radio occultation with TerraSAR-X and TanDEM-X: sensitivity of lower troposphere sounding to the Open-Loop Doppler model Atmospheric Measurement Techniques Discussions, 7, 12719-12733, 2014 Author(s): F. Zus, G. Beyerle, S. Heise, T. Schmidt, and J. Wickert The Global Positioning System (GPS) radio occultation (RO) technique provides valuable input for numerical weather prediction and is considered as a data source for climate related research. Numerous studies outline the high precision and accuracy of RO atmospheric soundings in the upper troposphere and lower stratosphere. In this altitude region (8–25 km) RO atmospheric soundings are considered to be free of any systematic error. In the tropical (30° S–30° N) Lower (

Description: Seasonal cycle and long-term trend of solar energy fluxes through Arctic sea ice The Cryosphere, 8, 2219-2233, 2014 Author(s): S. Arndt and M. Nicolaus Arctic sea ice has not only decreased in volume during the last decades, but has also changed in its physical properties towards a thinner and more seasonal ice cover. These changes strongly impact the energy budget, and might affect the ice-associated ecosystems. In this study, we quantify solar shortwave fluxes through sea ice for the entire Arctic during all seasons. To focus on sea-ice-related processes, we exclude fluxes through open water, scaling linearly with sea ice concentration. We present a new parameterization of light transmittance through sea ice for all seasons as a function of variable sea ice properties. The maximum monthly mean solar heat flux under the ice of 30 × 10 5 Jm −2 occurs in June, enough heat to melt 0.3 m of sea ice. Furthermore, our results suggest that 96% of the annual solar heat input through sea ice occurs during only a 4-month period from May to August. Applying the new parameterization to remote sensing and reanalysis data from 1979 to 2011, we find an increase in transmitted light of 1.5% yr −1 for all regions. This corresponds to an increase in potential sea ice bottom melt of 63% over the 33-year study period. Sensitivity studies reveal that the results depend strongly on the timing of melt onset and the correct classification of ice types. Assuming 2 weeks earlier melt onset, the annual transmitted solar radiation to the upper ocean increases by 20%. Continuing the observed transition from a mixed multi-year/first-year sea ice cover to a seasonal ice cover results in an increase in light transmittance by an additional 18%.

Description: Post-LIA glacier changes along a latitudinal transect in the Central Italian Alps The Cryosphere, 8, 2235-2252, 2014 Author(s): R. Scotti, F. Brardinoni, and G. B. Crosta The variability of glacier response to atmospheric temperature rise in different topo-climatic settings is still a matter of debate. To address this question in the Central Italian Alps, we compile a post-LIA (Little Ice Age) multitemporal glacier inventory (1860–1954–1990–2003–2007) along a latitudinal transect that originates north of the continental divide in the Livigno Mountains and extends south through the Disgrazia and Orobie ranges, encompassing continental-to-maritime climatic settings. In these sub-regions, we examine the area change of 111 glaciers. Overall, the total glacierized area has declined from 34.1 to 10.1 km 2 , with a substantial increase in the number of small glaciers due to fragmentation. The average annual decrease (AAD) in glacier area has risen by about 1 order of magnitude from 1860–1990 (Livigno: 0.45; Orobie: 0.42; and Disgrazia: 0.39 % a −1 ) to 1990–2007 (Livigno: 3.08; Orobie: 2.44; and Disgrazia: 2.27 % a −1 ). This ranking changes when considering glaciers smaller than 0.5 km 2 only (i.e., we remove the confounding caused by large glaciers in Disgrazia), so that post-1990 AAD follows the latitudinal gradient and Orobie glaciers stand out (Livigno: 4.07; Disgrazia: 3.57; and Orobie: 2.47 % a −1 ). More recent (2007–2013) field-based mass balances in three selected small glaciers confirm post-1990 trends showing the consistently highest retreat in continental Livigno and minimal area loss in maritime Orobie, with Disgrazia displaying transitional behavior. We argue that the recent resilience of glaciers in Orobie is a consequence of their decoupling from synoptic atmospheric temperature trends, a decoupling that arises from the combination of local topographic configuration (i.e., deep, north-facing cirques) and high winter precipitation, which ensures high snow-avalanche supply, as well as high summer shading and sheltering. Our hypothesis is further supported by the lack of correlations between glacier change and glacier attributes in Orobie, as well as by the higher variability in ELA,sub〉0 positioning, post-LIA glacier change, and interannual mass balances, as we move southward along the transect.

Description: A method to improve the determination of wave perturbations close to the tropopause by using a digital filter Atmospheric Measurement Techniques Discussions, 4, 1181-1197, 2011 Author(s): P. Alexander, A. de la Torre, P. Llamedo, R. Hierro, T. Schmidt, A. Haser, and J. Wickert Calculations of gravity wave activity all over the globe derived from GPS radio occultation temperature profiles led some years ago to the following question: are the wave amplitude enhancements systematically observed around tropopause levels due to physical processes or are they a simple artifact generated by any digital filter used to isolate the wave components? The latter alternative has been found to be the correct one. This has been concluded after almost a decade of work on global wave climatologies obtained from GPS radio occultation satellite data, which allowed to analyze, for the first time, a large amount of atmospheric profiles including both the troposphere and the stratosphere. We present a new filtering method which can be equally applied to temperature or refractivity profiles. The suggested technique significantly reduces artificial enhancements around the tropopause, which represents an improvement in comparison to previous applications of standard filters.

Description: An algorithm for retrieving black carbon optical parameters from thermal-optical (OC/EC) instruments Atmospheric Measurement Techniques Discussions, 4, 1233-1254, 2011 Author(s): A. Andersson, R. J. Sheesley, E. N. Kirillova, and Ö. Gustafsson Through absorption of sun light atmospheric black carbon (BC) is expected to influence regional/global climate by warming the atmosphere and dimming the surface. To evaluate the impact of these effects it is of interest to examine both the radiative properties of BC and the concentrations in air. Building on recent developments we present a novel application for combining these two aspects using the common thermal-optical (OC/EC) instrument. By correlating the OC/EC laser transmission with the FID-carbon detection non-carbon contributions to the light attenuation are detected. Such analysis allows the calculation of mass absorption cross-sections (MACs) for BC, corrected for certain in-organic components. This approach has been applied to data from two S S Asian and two SN S European sites, including a time series analysis for one of the SN S European sites. Taken together this study demonstrates broad applicability for this method while providing new insights into the optical properties of BC.

Description: Melting trends over the Greenland ice sheet (1958–2009) from spaceborne microwave data and regional climate models The Cryosphere, 5, 359-375, 2011 Author(s): X. Fettweis, M. Tedesco, M. van den Broeke, and J. Ettema To study near-surface melt changes over the Greenland ice sheet (GrIS) since 1979, melt extent estimates from two regional climate models were compared with those obtained from spaceborne microwave brightness temperatures using two different remote sensing algorithms. The results from the two models were consistent with those obtained with the remote sensing algorithms at both daily and yearly time scales, encouraging the use of the models for analyzing melting trends before the satellite era (1958–1979), when forcing data is available. Differences between satellite-derived and model-simulated results still occur and are used here to identify (i) biases in the snow models (notably in the albedo parametrization, in the thickness of a snow layer, in the maximum liquid water content within the snowpack and in the snowfall impacting the bare ice appearance in summer) and (ii) limitations in the use of passive microwave data for snowmelt detection at the edge of the ice sheet due to mixed pixel effect (e.g., tundra or rock nearby the ice sheet). The results from models and spaceborne microwave sensors confirm a significant (p-value = 0.01) increase in GrIS surface melting since 1979. The melt extent recorded over the last years (1998, 2003, 2005 and 2007) is unprecedented in the last 50 yr with the cumulated melt area in the 2000's being, on the average, twice that of the 1980's.

Description: Determination of field scale ammonia emissions for common slurry spreading practice with two independent methods Atmospheric Measurement Techniques Discussions, 4, 2635-2687, 2011 Author(s): J. Sintermann, C. Ammann, U. Kuhn, C. Spirig, R. Hirschberger, A. Gärtner, and A. Neftel At a cropland and a grassland site field scale ammonia (NH 3 ) emissions from slurry application were determined simultaneously by two approaches based on (i) eddy covariance (EC) flux measurements using high temperature Chemical Ionisation Mass Spectrometry (HT-CIMS) and on (ii) backward Lagrangian Stochastic (bLS) dispersion modelling using concentration measurements by three optical open path Fourier Transform Infrared (FTIR) systems. Slurry was spread on the fields in sequential tracks over a period of one to two hours. In order to calculate field emissions, measured EC/HT-CIMS fluxes were combined with flux footprint analysis of individual slurry spreading tracks to parameterise the NH 3 volatilisation with a bi-exponential time dependence. Accordingly, track-resolved concentration footprints for the FTIR measurements were calculated using bLS. Comparison of concentrations calculated from the parameterised fluxes with concentrations measured by impingers showed that the EC/HT-CIMS emissions on the two fertilisations corresponded to the impinger concentrations within 10 % while the bLS/FTIR results showed larger deviations. For both events, the period during fertilisation and the subsequent two hours contributed by more than 80 % to the total field emissions. Averaged over the two measurement methods, the cumulated emissions of the first day amounted to 17 ± 3 % loss of applied total ammoniacal nitrogen over the cropland and 16 ± 3 % over the grassland field.

Description: Empirical analysis and modeling of errors of atmospheric profiles from GPS radio occultation Atmospheric Measurement Techniques Discussions, 4, 2599-2633, 2011 Author(s): B. Scherllin-Pirscher, A. K. Steiner, G. Kirchengast, Y.-H. Kuo, and U. Foelsche The utilization of radio occultation (RO) data in atmospheric studies requires precise knowledge of error characteristics. We present results of an empirical error analysis of GPS radio occultation (RO) bending angle, refractivity, dry pressure, dry geopotential height, and dry temperature. We find very good agreement between data characteristics of different missions (CHAMP, GRACE-A, and Formosat-3/COSMIC (F3C)). In the global mean, observational errors (standard deviation from "true" profiles at mean tangent point location) agree within 0.3 % in bending angle, 0.1 % in refractivity, and 0.2 K in dry temperature at all altitude levels between 4 km and 35 km. Above ≈20 km, the observational errors show a strong seasonal dependence at high latitudes. Larger errors occur in hemispheric wintertime and are associated mainly with background data used in the retrieval process. The comparison between UCAR and WEGC results (both data centers have independent inversion processing chains) reveals different magnitudes of observational errors in atmospheric parameters, which are attributable to different background fields used. Based on the empirical error estimates, we provide a simple analytical error model for GPS RO atmospheric parameters and account for vertical, latitudinal, and seasonal variations. In the model, which spans the altitude range from 4 km to 35 km, a constant error is adopted around the tropopause region amounting to 0.8 % for bending angle, 0.35 % for refractivity, 0.15 % for dry pressure, 10 m for dry geopotential height, and 0.7 K for dry temperature. Below this region the observational error increases following an inverse height power-law and above it increases exponentially. The observational error model is the same for UCAR and WEGC data but due to somewhat different error characteristics below about 10 km and above about 20 km some parameters have to be adjusted. Overall, the observational error model is easily applicable and adjustable to individual error characteristics.

Description: Isotope hydrological studies of the perennial ice deposit of Saarhalle, Mammuthöhle, Dachstein Mts, Austria The Cryosphere, 5, 291-298, 2011 Author(s): Z. Kern, I. Fórizs, R. Pavuza, M. Molnár, and B. Nagy A 5.28 m-long ice core was extracted from a major cave ice body in the Mammuthöhle cave system. The upper ~1.2 m of ice most likely originate from precipitation fallen before the 1960s (based on

Description: Co-registration and bias corrections of satellite elevation data sets for quantifying glacier thickness change The Cryosphere, 5, 271-290, 2011 Author(s): C. Nuth and A. Kääb There are an increasing number of digital elevation models (DEMs) available worldwide for deriving elevation differences over time, including vertical changes on glaciers. Most of these DEMs are heavily post-processed or merged, so that physical error modelling becomes difficult and statistical error modelling is required instead. We propose a three-step methodological framework for assessing and correcting DEMs to quantify glacier elevation changes: (i) remove DEM shifts, (ii) check for elevation-dependent biases, and (iii) check for higher-order, sensor-specific biases. A simple, analytic and robust method to co-register elevation data is presented in regions where stable terrain is either plentiful (case study New Zealand) or limited (case study Svalbard). The method is demonstrated using the three global elevation data sets available to date, SRTM, ICESat and the ASTER GDEM, and with automatically generated DEMs from satellite stereo instruments of ASTER and SPOT5-HRS. After 3-D co-registration, significant biases related to elevation were found in some of the stereoscopic DEMs. Biases related to the satellite acquisition geometry (along/cross track) were detected at two frequencies in the automatically generated ASTER DEMs. The higher frequency bias seems to be related to satellite jitter , most apparent in the back-looking pass of the satellite. The origins of the more significant lower frequency bias is uncertain. ICESat-derived elevations are found to be the most consistent globally available elevation data set available so far. Before performing regional-scale glacier elevation change studies or mosaicking DEMs from multiple individual tiles (e.g. ASTER GDEM), we recommend to co-register all elevation data to ICESat as a global vertical reference system.

Description: Present dynamics and future prognosis of a slowly surging glacier The Cryosphere, 5, 299-313, 2011 Author(s): G. E. Flowers, N. Roux, S. Pimentel, and C. G. Schoof Glacier surges are a well-known example of an internal dynamic oscillation whose occurrence is not a direct response to the external climate forcing, but whose character (i.e. period, amplitude, mechanism) may depend on the glacier's environmental or climate setting. We examine the dynamics of a small (∼5 km 2 ) valley glacier in Yukon, Canada, where two previous surges have been photographically documented and an unusually slow surge is currently underway. To characterize the dynamics of the present surge, and to speculate on the future of this glacier, we employ a higher-order flowband model of ice dynamics with a regularized Coulomb-friction sliding law in both diagnostic and prognostic simulations. Diagnostic (force balance) calculations capture the measured ice-surface velocity profile only when non-zero basal water pressures are prescribed over the central region of the glacier, coincident with where evidence of the surge has been identified. This leads to sliding accounting for 50–100% of the total surface motion in this region. Prognostic simulations, where the glacier geometry evolves in response to a prescribed surface mass balance, reveal a significant role played by a bedrock ridge beneath the current equilibrium line of the glacier. Ice thickening occurs above the ridge in our simulations, until the net mass balance reaches sufficiently negative values. We suggest that the bedrock ridge may contribute to the propensity for surges in this glacier by promoting the development of the reservoir area during quiescence, and may permit surges to occur under more negative balance conditions than would otherwise be possible. Collectively, these results corroborate our interpretation of the current glacier flow regime as indicative of a slow surge that has been ongoing for some time, and support a relationship between surge incidence or character and the net mass balance. Our results also highlight the importance of glacier bed topography in controlling ice dynamics, as observed in many other glacier systems.

Description: Remote sensing of aerosols by using polarized, directional and spectral measurements within the A-Train: the PARASOL mission Atmospheric Measurement Techniques Discussions, 4, 2037-2069, 2011 Author(s): D. Tanré, F. M. Bréon, J. L. Deuzé, O. Dubovik, F. Ducos, P. François, P. Goloub, M. Herman, A. Lifermann, and F. Waquet The aerosol remote sensing from space has started in the 1980's using observations provided by geostationary satellites or by polar orbiting platforms not specifically designed for observing aerosols. As a result, the number of retrieved parameters was limited and retrievals in the visible restricted over ocean. Over land, because of the important surface contribution, the aerosol detection was performed in the UV (or in the dark blue) where most of the earth surfaces are dark enough but with overlap of multiple aerosols parameters, content, altitude and absorption. Instruments dedicated to aerosol monitoring are recently available and the POLDER instrument on board the PARASOL mission is one of them. By measuring the wavelength, angular and polarization properties of the radiance at the top of the atmosphere, in coordination with the other A-Train instruments, PARASOL can better quantify aerosol optical depths (AOD) and improve the derivation of the radiative and physical properties. The instrument, the inversion schemes and the list of aerosol parameters are described. Examples of retrieved aerosol parameters are provided as well as innovative approaches and further inversion techniques.

Description: Atmospheric influences on infrared-laser signals used for occultation measurements between Low Earth Orbit satellites Atmospheric Measurement Techniques Discussions, 4, 2689-2747, 2011 Author(s): S. Schweitzer, G. Kirchengast, and V. Proschek LEO-LEO infrared-laser occultation (LIO) is a new occultation technique between Low Earth Orbit (LEO) satellites, which applies signals in the short wave infrared spectral range (SWIR) within 2 μm to 2.5 μm. It is part of the LEO-LEO microwave and infrared-laser occultation (LMIO) method, recently introduced by Kirchengast and Schweitzer (2011), that enables to retrieve thermodynamic profiles (pressure, temperature, humidity) and accurate altitude levels from microwave signals and profiles of greenhouse gases and further variables such as line-of-sight wind speed from simultaneously measured LIO signals. For enabling trace species retrieval based on differential transmission, the LIO signals are spectrally located as pairs, one in the centre of a suitable absorption line of a target species (absorption signal) and one close by but outside of any absorption lines (reference signal). Due to the novelty of the LMIO method, detailed knowledge of atmospheric influences on LIO signals and of their suitability for accurate trace species retrieval did not yet exist. Here we discuss the atmospheric influences on the transmission and differential transmission of LIO signals. Refraction effects, trace species absorption (by target species, and cross-sensitivity to foreign species), aerosol extinction and Rayleigh scattering are studied in detail. The influences of clouds, turbulence, wind, scattered solar radiation and terrestrial thermal radiation are discussed as well. We show that the influence of defocusing, foreign species absorption, aerosols and turbulence is observable, but can be rendered small to negligible by use of the differential transmission principle and by a design with close frequency spacing of absorption and reference signals within 0.5 %. The influences of Rayleigh scattering and thermal radiation on the received signal intensities are found negligible. Cloud-scattered solar radiation can be observable under bright-day conditions but this influence can be as well made negligible by a design with a close time spacing (within 5 ms) of interleaved laser-pulse and background signals. Cloud extinction loss generally blocks SWIR signals, except very thin or sub-visible cirrus clouds, which can be addressed by a design allowing retrieval of a cloud layering profile from reference signals and its use in trace species retrieval when scanning through intermittent upper tropospheric cloudiness. Wind can have a small influence via Doppler shift resulting in a slightly modified trace species absorption in comparison to calm air, which can be made negligible by using a simultaneously retrieved wind speed profile or a moderately accurate (to about 10 m s −1 ) background wind profile. Considering all these influences, we conclude that the set of SWIR channels proposed for implementing the LMIO method (Kirchengast et al., 2010; Kirchengast and Schweitzer, 2011) provides adequate sensitivity to accurately retrieve eight greenhouse gas/isotope trace species of key importance to climate and atmospheric chemistry (H 2 O, 12 CO 2 , 13 CO 2 , C 18 OO, CH 4 , N 2 O, O 3 , CO) in the upper troposphere/lower stratosphere region outside clouds under all atmospheric conditions. Two further isotope species (HDO, H 2 18 O) can be retrieved in the upper troposphere.

Description: Effective density of Aquadag and fullerene soot black carbon reference materials used for SP2 calibration Atmospheric Measurement Techniques Discussions, 4, 4937-4955, 2011 Author(s): M. Gysel, M. Laborde, J. S. Olfert, R. Subramanian, and A. J. Gröhn The mass and effective density of black carbon (BC) particles generated from aqueous suspensions of Aquadag and fullerene soot was measured and parametrized as a function of their mobility diameter. The measurements were made by two independent research groups by operating a differential mobility analyser (DMA) in series with an aerosol particle mass analyser (APM) or a Couette centrifugal particle mass analyser (CPMA). Consistent and reproducible results were found in this study for different production lots of Aquadag, indicating that the effective density of these particles is a stable quantity and largely unaffected by differences in aerosol generation procedures and suspension treatments. The effective density of fullerene soot particles from one production lot was also found to be stable and independent of suspension treatments. Some difference to previous literature data was observed for both Aquadag and fullerene soot at larger particle diameters. Knowledge of the exact relationship between mobility diameter and particle mass is of great importance, as DMAs are commonly used to size-select particles from BC reference materials for calibration of single particle soot photometers (SP2), which quantitatively detect the BC mass in single particles.

Description: Two instruments based on differential optical absorption spectroscopy (DOAS) to measure accurate ammonia concentrations in the atmosphere Atmospheric Measurement Techniques Discussions, 4, 5037-5078, 2011 Author(s): H. Volten, J. B. Bergwerff, M. Haaima, D. E. Lolkema, A. J. C. Berkhout, G. R. van der Hoff, C. J. M. Potma, R. J. Wichink Kruit, W. A. J. van Pul, and D. P. J. Swart We present two Differential Optical Absorption Spectroscopy (DOAS) instruments built at RIVM, the RIVM DOAS and the miniDOAS. Both instruments provide virtually interference free measurements of NH 3 concentrations in the atmosphere, since they measure over an open path, without suffering from inlet problems or interference problems by ammonium aerosols dissociating on tubes or filters. They measure concentrations up to at least 200 μg m −3 , have a fast response, low maintenance demands, and a high up-time. The RIVM DOAS has a high accuracy of typically 0.15 μg m −3 for ammonia over 5-min averages and over a total light path of 100 m. The miniDOAS has been developed for application in measurement networks such as the Dutch National Air Quality Monitoring Network (LML). Compared to the RIVM DOAS it has a similar accuracy, but is significantly reduced in size, costs, and handling complexity. The RIVM DOAS and miniDOAS results showed excellent agreement ( R 2 = 0.996) during a field measurement campaign in Vredepeel, the Netherlands. This measurement site is located in an agricultural area and is characterized by highly variable, but on average high ammonia concentrations in the air. The RIVM-DOAS and miniDOAS results were compared to the results of the AMOR instrument, a continuous-flow wet denuder system, which is currently used in the LML. Averaged over longer time spans of typically a day the (mini)DOAS and AMOR results agree reasonably well, although an offset of the AMOR values compared to the (mini)DOAS results exists. On short time scales the (mini)DOAS shows a faster response and does not show the memory effects due to inlet tubing and transport of absorption fluids encountered by the AMOR. Due to its high accuracy, high uptime, low maintenance and its open path, the (mini)DOAS shows a good potential for flux measurements by using two (or more) systems in a gradient set-up and applying the aerodynamic gradient technique.

Description: Towards space based verification of CO 2 emissions from strong localized sources: fossil fuel power plant emissions as seen by a CarbonSat constellation Atmospheric Measurement Techniques Discussions, 4, 5147-5182, 2011 Author(s): V. A. Velazco, M. Buchwitz, H. Bovensmann, M. Reuter, O. Schneising, J. Heymann, T. Krings, K. Gerilowski, and J. P. Burrows Carbon dioxide (CO 2 ) is the most important man-made greenhouse gas (GHG) that cause global warming. With electricity generation through fossil-fuel power plants now as the economic sector with the largest source of CO 2 , power plant emissions monitoring has become more important than ever in the fight against global warming. In a previous study done by Bovensmann et al. (2010), random and systematic errors of power plant CO 2 emissions have been quantified using a single overpass from a proposed CarbonSat instrument. In this study, we quantify errors of power plant annual emission estimates from a hypothetical CarbonSat and constellations of several CarbonSats while taking into account that power plant CO 2 emissions are time-dependent. Our focus is on estimating systematic errors arising from the sparse temporal sampling as well as random errors that are primarily dependent on wind speeds. We used hourly emissions data from the US Environmental Protection Agency (EPA) combined with assimilated and re-analyzed meteorological fields from the National Centers of Environmental Prediction (NCEP). CarbonSat orbits were simulated as a sun-synchronous low-earth orbiting satellite (LEO) with an 828-km orbit height, local time ascending node (LTAN) of 13:30 (01:30 p.m.) and achieves global coverage after 5 days. We show, that despite the variability of the power plant emissions and the limited satellite overpasses, one CarbonSat can verify reported US annual CO 2 emissions from large power plants (≥5 Mt CO 2 yr −1 ) with a systematic error of less than ~4.9 % for 50 % of all the power plants. For 90 % of all the power plants, the systematic error was less than ~12.4 %. We additionally investigated two different satellite configurations using a combination of 5 CarbonSats. One achieves global coverage everyday but only samples the targets at fixed local times. The other configuration samples the targets five times at two-hour intervals approximately every 6th day but only achieves global coverage after 5 days. From the statistical analyses, we found, as expected, that the random errors improve by approximately a factor of two if 5 satellites are used. On the other hand, more satellites do not result in a large reduction of the systematic error. The systematic error is somewhat smaller for the CarbonSat constellation configuration achieving global coverage everyday. Finally, we recommend the CarbonSat constellation configuration that achieves daily global coverage.

Description: In-situ multispectral and bathymetric measurements over a supraglacial lake in western Greenland using a remotely controlled watercraft The Cryosphere, 5, 445-452, 2011 Author(s): M. Tedesco and N. Steiner Supraglacial lakes form from meltwater on the Greenland ice sheet in topographic depressions on the surface, affecting both surface and sub-glacial processes. As the reflectance in the visible and near-infrared regions of a column of water is modulated by its height, retrieval techniques using spaceborne remote sensing data (e.g. Landsat, MODIS) have been proposed in the literature for the detection of lakes and estimation of their volume. These techniques require basic assumptions on the spectral properties of the water as well as the bottom of the lake, among other things. In this study, we report results obtained from the analysis of concurrent in-situ multi-spectral and depth measurements collected over a supraglacial lake during early July 2010 in West Greenland (Lake Olivia, 69°36'35" N, 49°29'40" W) and aim to assess some of the underlying hypotheses in remote sensing based bathymetric approaches. In particular, we focus our attention on the analysis of the lake bottom albedo and of the water attenuation coefficient. The analysis of in-situ data (collected by means of a remotely controlled boat equipped with a GPS, a sonar and a spectrometer) highlights the exponential trend of the water-leaving reflectance with lake depth. The values of the attenuation factor obtained from in-situ data are compared with those computed using approaches proposed in the literature. Also, the values of the lake bottom albedo from in-situ measurements are compared with those obtained from the analysis of reflectance of shallow waters. Finally, we quantify the error between in-situ measured and satellite-estimated lake depth values for the lake under study.

Description: Tropospheric and total ozone columns over Paris (France) measured using medium-resolution ground-based solar-absorption Fourier-transform infrared spectroscopy Atmospheric Measurement Techniques Discussions, 4, 3337-3358, 2011 Author(s): C. Viatte, B. Gaubert, M. Eremenko, F. Hase, M. Schneider, T. Blumenstock, M. Ray, P. Chelin, J.-M. Flaud, and J. Orphal Ground-based Fourier-transform infrared (FTIR) solar absorption spectroscopy is a powerful remote sensing technique providing information on the vertical distribution of various atmospheric constituents. This work presents the first evaluation of a mid-resolution ground-based FTIR to measure tropospheric ozone, independently of stratospheric ozone. This is demonstrated using a new atmospheric observatory (named OASIS for "Observations of the Atmosphere by Solar absorption Infrared Spectroscopy"), installed in Créteil (France). Indeed, the information content of OASIS ozone retrievals is clearly sufficient to monitor separately tropospheric (from the surface up to 8 km) and stratospheric ozone. Daily mean tropospheric ozone columns derived from the Infrared Atmospheric Sounding Interferometer (IASI) and from OASIS measurements have been compared for summer 2009 and a good agreement of −5.6 (±16.1) % is observed. Also, a qualitative comparison between in-situ surface ozone measurements and OASIS data clearly shows OASIS's capacity to monitor seasonal tropospheric ozone variations, as well as ozone pollution episodes in summer 2009 around Paris. Two extreme pollution events were identified (on the 1 July and 6 August 2009) for which ozone partial columns from OASIS and predictions from a regional air-quality model (CHIMERE) were compared by respecting temporal and spatial coincidence criteria. Quantitatively, an average bias of 0.2 %, a mean square error deviation of 7.6 %, and a correlation coefficient of 0.91 was found between CHIMERE and OASIS. This demonstrates that a mid-resolution FTIR instrument in ground-based solar absorption geometry is a promising technique for monitoring tropospheric ozone.

Description: Measurements of Humidity in the Atmosphere and Validation Experiments (MOHAVE)-2009: overview of campaign operations and results Atmospheric Measurement Techniques Discussions, 4, 3277-3336, 2011 Author(s): T. Leblanc, T. D. Walsh, I. S. McDermid, G. C. Toon, J.-F. Blavier, B. Haines, W. G. Read, B. Herman, E. Fetzer, S. Sander, T. Pongetti, D. N. Whiteman, T. G. McGee, L. Twigg, G. Sumnicht, D. Venable, M. Calhoun, A. Dirisu, D. Hurst, A. Jordan, E. Hall, L. Miloshevich, H. Vömel, C. Straub, N. Kampfer, G. E. Nedoluha, R. M. Gomez, K. Holub, S. Gutman, J. Braun, T. Vanhove, G. Stiller, and A. Hauchecorne The Measurements of Humidity in the Atmosphere and Validation Experiment (MOHAVE) 2009 campaign took place on 11–27 October 2009 at the JPL Table Mountain Facility in California (TMF). The main objectives of the campaign were to (1) validate the water vapor measurements of several instruments, including, three Raman lidars, two microwave radiometers, two Fourier-Transform spectrometers, and two GPS receivers (column water), (2) cover water vapor measurements from the ground to the mesopause without gaps, and (3) study upper tropospheric humidity variability at timescales varying from a few minutes to several days. A total of 58 radiosondes and 20 Frost-Point hygrometer sondes were launched. Two types of radiosondes were used during the campaign. Non negligible differences in the readings between the two radiosonde types used (Vaisala RS92 and InterMet iMet-1) made a small, but measurable impact on the derivation of water vapor mixing ratio by the Frost-Point hygrometers. As observed in previous campaigns, the RS92 humidity measurements remained within 5 % of the Frost-point in the lower and mid-troposphere, but were too dry in the upper troposphere. Over 270 h of water vapor measurements from three Raman lidars (JPL and GSFC) were compared to RS92, CFH, and NOAA-FPH. The JPL lidar profiles reached 20 km when integrated all night, and 15 km when integrated for 1 h. Excellent agreement between this lidar and the frost-point hygrometers was found throughout the measurement range, with only a 3 % (0.3 ppmv) mean wet bias for the lidar in the upper troposphere and lower stratosphere (UTLS). The other two lidars provided satisfactory results in the lower and mid-troposphere (2–5 % wet bias over the range 3–10 km), but suffered from contamination by fluorescence (wet bias ranging from 5 to 50 % between 10 km and 15 km), preventing their use as an independent measurement in the UTLS. The comparison between all available stratospheric sounders allowed to identify only the largest biases, in particular a 10 % dry bias of the Water Vapor Millimeter-wave Spectrometer compared to the Aura-Microwave Limb Sounder. No other large, or at least statistically significant, biases could be observed. Total Precipitable Water (TPW) measurements from six different co-located instruments were available. Several retrieval groups provided their own TPW retrievals, resulting in the comparison of 10 different datasets. Agreement within 7 % (0.7 mm) was found between all datasets. Such good agreement illustrates the maturity of these measurements and raises confidence levels for their use as an alternate or complementary source of calibration for the Raman lidars. Tropospheric and stratospheric ozone and temperature measurements were also available during the campaign. The water vapor and ozone lidar measurements, together with the advected potential vorticity results from the high-resolution transport model MIMOSA, allowed the identification and study of a deep stratospheric intrusion over TMF. These observations demonstrated the lidar strong potential for future long-term monitoring of water vapor in the UTLS.

Description: A method to resolve the phase state of aerosol particles Atmospheric Measurement Techniques Discussions, 4, 6229-6248, 2011 Author(s): E. Saukko, H. Kuuluvainen, and A. Virtanen The phase state of atmospheric aerosols has impact on their chemical aging and their deliquescence and thus their ability to act as cloud condensation nuclei (CCN). The phase change of particles can be induced by the deliquescence or efflorescence of water or by chemical aging. Existing methods, such as tandem differential mobility analysis rely on the size change of particles related to the water uptake or release related to deliquescence and efflorescence. To address the need to study the phase change induced by mass-preserving and nearly mass-preserving processes a new method has been developed. The method relies on the physical impaction of particles on a smooth substrate and subsequent counting of bounced particles by condensation particle counter (CPC). The connection between the bounce probability and physical properties of particles is so far qualitative. To evaluate the performance of this method, the phase state of ammonium sulfate and levoglucosan, crystalline and amorphous solid, in the presence of water vapor was studied. The results show a marked difference in particle bouncing properties between substances – not only at the critical relative humidity level, but also on the slope of the bouncing probability with respect to humidity. This suggests that the method can be used to differentiate between amorphous and crystalline substances as well as to differentiate between liquid and solid phases.

Description: Fast simulators for satellite cloud optical centroid pressure retrievals, 1. evaluation of OMI cloud retrievals Atmospheric Measurement Techniques Discussions, 4, 6185-6228, 2011 Author(s): J. Joiner, A. P. Vasilkov, P. Gupta, P. K. Bhartia, P. Veefkind, M. Sneep, J. de Haan, I. Polonsky, and R. Spurr The cloud Optical Centroid Pressure (OCP) is a satellite-derived parameter that is commonly used in trace-gas retrievals to account for the effects of clouds on near-infrared through ultraviolet radiance measurements. A fast simulator is desirable in order to further expand the use of cloud OCP retrievals into the operational and climate communities for applications such as data assimilation and evaluation of cloud vertical structure in general circulation models. In this paper, we develop and validate a fast simulator that provides estimates of the cloud OCP given a vertical profile of optical extinction. We use a pressure-weighting scheme where the weights depend upon optical parameters of clouds and/or aerosol. A cloud weighting function is easily extracted using this formulation. We then use the fast simulator to compare two different satellite cloud OCP retrievals from the Ozone Monitoring Instrument (OMI) with estimates based on collocated cloud extinction profiles from a combination of CloudSat radar and MODIS visible radiance data. These comparisons are made over a wide range of conditions in order to provide a comprehensive validation of the OMI cloud OCP retrievals. We find generally good agreement between OMI cloud OCPs and those predicted by CloudSat. However, the OMI cloud OCPs from the two independent algorithms agree better with each other than either does with the estimates from CloudSat/MODIS. Differences between OMI cloud OCPs and those based on CloudSat/MODIS may result from undetected snow/ice at the surface, cloud 3-D effects, and the fact that CloudSat only observes a relatively small fraction of OMI pixel.

Description: Quantification of gas-phase glyoxal and methylglyoxal via the Laser-Induced Phosphorescence of (methyl)GLyOxal Spectrometry (LIPGLOS) method Atmospheric Measurement Techniques Discussions, 4, 6159-6183, 2011 Author(s): S. B. Henry, A. Kammrath, and F. N. Keutsch Glyoxal and methylglyoxal are key products of oxidative photochemistry in the lower troposphere. Reliable measurements of such compounds are critical for testing our understanding of volatile organic compound (VOC) processing in this region. We present a new method for obtaining sensitive, high time resolution, in situ measurements of these compounds via laser-induced phosphorescent decays. By exploiting the unique phosphorescent lifetimes for each molecule, this method achieves speciation and high-sensitivity quantification of both molecules (3σ limits of detection of 11 ppt v in 5 min for glyoxal and 243 ppt v in 5 min for methylglyoxal). Additionally, this method enables the simultaneous measurement of both glyoxal and methylglyoxal using a single, non-wavelength-tunable light source, which will allow for the development of inexpensive and turnkey instrumentation. The simplicity and affordability of this new instrumentation would enable the construction of a long-term, spatially distributed database of these two key species. This chemical map can be used to constrain or drive regional or global models as well as provide verification of satellite observations.

Description: Temporal co-registration for TROPOMI cloud clearing Atmospheric Measurement Techniques Discussions, 4, 6249-6272, 2011 Author(s): I. Genkova, J. Robaidek, R. Roebeling, M. Sneep, and P. Veefkind The TROPOspheric Monitoring Instrument (TROPOMI) is planed for launch in 2014 on board of the Sentinel 5 Precursor (S5P) and is anticipated to provide high-quality and timely information on the global atmospheric composition for climate and air quality applications. TROPOMI will observe key atmospheric constituents such as ozone, nitrogen dioxide, sulfur dioxide, carbon monoxide, methane, formaldehyde and aerosol properties. The retrieval algorithms for the anticipated products require cloud information on a pixel basis. Most of them will use the cloud properties derived from TROPOMI's own measurements, such as the O 2 A-band measurements. However, the methane and the aerosol retrievals require very precise cloud clearing, which is difficult to achieve at the TROPOMI spatial resolution (7 × 7 km 2 ) and without thermal IR measurements. The current payload of the Sentinel 5 Precursor (S-5P) does not include a cloud imager, thus it is planned to fly the S5P mission in a constellation with another instrument yielding an accurate cloud mask. The cloud imagery data will be provided by the US NPOESS Preparatory Project (NPP) mission which will have the Visible Infrared Imager Radiometer Suite (VIIRS) on board (Scalione, 2004). VIIRS will have 22 bands in the VIS and IR spectral ranges, and will deliver data with two spatial resolutions: imagery resolution bands with a nominal pixel size of 370 m at nadir, and moderate resolution bands with nominal pixel size 740 m at nadir. The instrument is combining fine spatial resolution with high-accuracy calibration similar or superior to AVHRR. This paper presents results from investigating the temporal co-registration requirements for suitable time differences between the VIIRS measurements of clouds and the TROPOMI methane and aerosol measurements, so that the former could be used for cloud clearing. The temporal co-registration is studied using Meteosat Second Generation (MSG) Spinning Enhanced Visible and Infrared Imager (SEVIRI) data with 15 min temporal resolution (Veefkind, 2008a), and with data from the Geostationary Operational Environmental Satellite-10 (GOES-10) having 1 min temporal resolution. The aim is to understand and assess the relation between the amount of allowed cloud contamination and the required time difference between the two satellites' overflights. Quantitative analysis shows that a time difference of approximately 5 min is sufficient (in most conditions) to use the cloud information from the first instrument for cloud clearing in the retrievals using data from the second instrument. In recent years the A-train constellation demonstrated the benefit of flying satellites in formation. Therefore this study's findings will be useful and applicable for designing future Low Earth Orbit (LEO) satellite constellations.

Description: 1-D-Var retrieval of daytime total columnar water vapour from MERIS measurements Atmospheric Measurement Techniques Discussions, 4, 6811-6844, 2011 Author(s): R. Lindstrot, R. Preusker, H. Diedrich, L. Doppler, R. Bennartz, and J. Fischer A new scheme for the retrieval of total columnar water vapour from measurements of MERIS (Medium Resolution Imaging Spectrometer) on ENVISAT (ENVIronmental SATellite) is presented. The algorithm is based on a fast forward model of the water vapour transmittance around 900 nm, including a correction for atmospheric scattering. It provides the water vapour column amount for cloud-free scenes above land and ocean at a spatial resolution of 0.25 km × 0.3 km and 1 km × 1.2 km, depending on whether applied to the "full resolution" or the operational "reduced resolution" mode of MERIS. Uncertainties are provided on a pixel-by-pixel basis, taking into account all relevant sources of error. An extensive validation against various sources of ground-based reference data reveals a high accuracy of MERIS water vapour above land (root mean square deviations between 1 mm and 2.7 mm), apart from a wet bias of MERIS between 5 and 10% that is found in all comparison studies. This wet bias might be caused by spectroscopic uncertainties, such as the description of the water vapour continuum. Above ocean the accuracy is reduced, due to the uncertainty introduced by the unknown atmospheric scattering. Consequently, an increased root mean square deviation of ≥5 mm was found by comparing MERIS total columnar water vapour above ocean against SSM/I and ENVISAT MWR data. An increased wet bias of 2–3 mm is found over ocean, potentially du to a not properly working atmospheric correction scheme.

Description: Initial investigation of the wavelength dependence of optical properties measured with a new multi-pass aerosol extinction differential optical absorption spectrometer (AE-DOAS) Atmospheric Measurement Techniques Discussions, 4, 6315-6349, 2011 Author(s): R. T. Chartier and M. E. Greenslade Atmospheric aerosols directly affect climate by scattering and absorbing radiation. The magnitude of the impact is dependent upon the wavelength of light, but is often estimated near 550 nm. When light scattering and absorption by aerosols is approximated, the wavelength dependence of the refractive index for specific components is lost. As a result, climate models would have inherent uncertainties for aerosol contributions to radiative forcing when considering the entire solar spectrum. An aerosol extinction differential optical absorption spectrometer has been developed to directly measure aerosol extinction at mid-ultraviolet to near infrared wavelengths. The instrument consists of a spectrometer coupled to a closed White-type multi-pass gas cell with an adjustable path length of up to approximately 20 m. Laboratory measurements of various gases are compared with known absorption cross sections. Additionally, the extinction of monodisperse samples of polystyrene latex spheres are measured and compared to Mie theory generated with refractive index values from the literature to validate the new instrument. The polystyrene experiments also emphasize the ability of the new instrument to retrieve the wavelength dependent refractive index, especially in the ultraviolet wavelength regions where variability is expected. The spectrometer will be a significant advancement for determining wavelength dependent complex refractive indices in future laboratory studies as well as provide the ability to monitor ambient aerosol light extinction.

Description: Temperature variability and offset in steep alpine rock and ice faces The Cryosphere, 5, 977-988, 2011 Author(s): A. Hasler, S. Gruber, and W. Haeberli The thermal condition of high-alpine mountain flanks can be an important determinant of climate change impact on slope stability and correspondingly down-slope hazard regimes. In this study we analyze time-series from 17 shallow temperature-depth profiles at two field sites in steep bedrock and ice. Extending earlier studies that revealed the topographic variations in temperatures, we demonstrate considerable differences of annual mean temperatures for variable surface characteristics and depths within the measured profiles. This implies that measurements and model related to compact and near-vertical bedrock temperatures may deviate considerably from conditions in the majority of bedrock slopes in mountain ranges that are usually non-vertical and fractured. For radiation-exposed faces mean annual temperatures at depth are up to 3 °C lower and permafrost is likely to exist at lower elevations than reflected by estimates based on near-vertical homogeneous cases. Retention of a thin snow cover and ventilation effects in open clefts are most likely responsible for this cooling. The measurements presented or similar data could be used in the future to support the development and testing of models related to the thermal effect of snow-cover and fractures in steep bedrock.

Description: On the influence of model physics on simulations of Arctic and Antarctic sea ice The Cryosphere, 5, 687-699, 2011 Author(s): F. Massonnet, T. Fichefet, H. Goosse, M. Vancoppenolle, P. Mathiot, and C. König Beatty Two hindcast (1983–2007) simulations are performed with the global, ocean-sea ice models NEMO-LIM2 and NEMO-LIM3 driven by atmospheric reanalyses and climatologies. The two simulations differ only in their sea ice component, while all other elements of experimental design (resolution, initial conditions, atmospheric forcing) are kept identical. The main differences in the sea ice models lie in the formulation of the subgrid-scale ice thickness distribution, of the thermodynamic processes, of the sea ice salinity and of the sea ice rheology. To assess the differences in model skill over the period of investigation, we develop a set of metrics for both hemispheres, comparing the main sea ice variables (concentration, thickness and drift) to available observations and focusing on both mean state and seasonal to interannual variability. Based upon these metrics, we discuss the physical processes potentially responsible for the differences in model skill. In particular, we suggest that (i) a detailed representation of the ice thickness distribution increases the seasonal to interannual variability of ice extent, with spectacular improvement for the simulation of the recent observed summer Arctic sea ice retreats, (ii) the elastic-viscous-plastic rheology enhances the response of ice to wind stress, compared to the classical viscous-plastic approach, (iii) the grid formulation and the air-sea ice drag coefficient affect the simulated ice export through Fram Strait and the ice accumulation along the Canadian Archipelago, and (iv) both models show less skill in the Southern Ocean, probably due to the low quality of the reanalyses in this region and to the absence of important small-scale oceanic processes at the models' resolution (~1°).

Description: Retrieval of sulphur dioxide from the infrared atmospheric sounding interferometer (IASI) Atmospheric Measurement Techniques Discussions, 4, 7241-7275, 2011 Author(s): L. Clarisse, D. Hurtmans, C. Clerbaux, J. Hadji-Lazaro, Y. Ngadi, and P.-F. Coheur Thermal infrared sounding of sulphur dioxide (SO 2 ) from space has gained appreciation and popularity as a valuable complement to ultraviolet sounding. There are several strong absorption bands of SO 2 in the infrared, and atmospheric sounders, primarily designed for weather forecasting, have therefore often the ability to globally monitor SO 2 abundances. Most of the observed SO 2 is found in volcanic plumes. In this paper we outline a novel algorithm for the sounding of SO 2 above ~500 hPa altitude using high resolution infrared sounders and apply it to measurements of the infrared atmospheric sounding interferometer (IASI). The main features of the algorithm are a wide applicable total column range (over 4 orders of magnitude, from 0.5 to 5000 dobson units), a low theoretical uncertainty (3–5%) and near real time applicability. We make an error analysis and demonstrate the algorithm on the recent eruptions of Sarychev, Kasatochi, Grimsvötn, Puyehue-Cordón Caulle and Nabro.

Description: Glacier contribution to streamflow in two headwaters of the Huasco River, Dry Andes of Chile The Cryosphere, 5, 1099-1113, 2011 Author(s): S. Gascoin, C. Kinnard, R. Ponce, S. Lhermitte, S. MacDonell, and A. Rabatel Quantitative assessment of glacier contribution to present-day streamflow is a prerequisite to the anticipation of climate change impact on water resources in the Dry Andes. In this paper we focus on two glaciated headwater catchments of the Huasco Basin (Chile, 29° S). The combination of glacier monitoring data for five glaciers (Toro 1, Toro 2, Esperanza, Guanaco, Estrecho and Ortigas) with five automatic streamflow records at sites with glacier coverage of 0.4 to 11 % allows the estimation of the mean annual glacier contribution to discharge between 2003/2004 and 2007/2008 hydrological years. In addition, direct manual measurements of glacier runoff were conducted in summer at the snouts of four glaciers, which provide the instantaneous contribution of glacier meltwater to stream runoff during summer. The results show that the mean annual glacier contribution to streamflow ranges between 3.3 and 23 %, which is greater than the glaciated fraction of the catchments. We argue that glacier contribution is partly enhanced by the effect of snowdrift from the non-glacier area to the glacier surface. Glacier mass loss is evident over the study period, with a mean of −0.84 m w.e. yr −1 for the period 2003/2004–2007/2008, and also contributes to increase glacier runoff. An El Niño episode in 2002 resulted in high snow accumulation, modifying the hydrological regime and probably reducing the glacier contribution in favor of seasonal snowmelt during the subsequent 2002/2003 hydrological year. At the hourly timescale, summertime glacier contributions are highly variable in space and time, revealing large differences in effective melting rates between glaciers and glacierets (from 1 mm w.e. h −1 to 6 mm w.e. h −1 ).