ALBERT

Description: Influence of microphysical schemes on atmospheric water in the Weather Research and Forecasting model Geoscientific Model Development Discussions, 6, 4563-4601, 2013 Author(s): F. Cossu and K. Hocke This study examines how different microphysical parameterization schemes influence orographically-induced precipitation and the distributions of hydrometeors and water vapour for mid-latitude summer conditions in the Weather Research and Forecasting (WRF) model. A high-resolution two-dimensional idealized simulation is used to assess the differences between the schemes in which a moist air flow is interacting with a bell-shaped 2 km high mountain. Periodic lateral boundary conditions are chosen to recirculate atmospheric water in the domain. It is found that the 13 selected microphysical schemes conserve the water in the model domain. The gain or loss of water is less than 0.81% over a simulation time interval of 61 days. The differences of the microphysical schemes in terms of the distributions of water vapour, hydrometeors and accumulated precipitation are presented and discussed. The Kessler scheme, the only scheme without ice-phase processes, shows final values of cloud liquid water 14 times greater than the other schemes. The differences among the other schemes are not as extreme, but still they differ up to 79% in water vapour, up to 10 times in hydrometeors and up to 64% in accumulated precipitation at the end of the simulation. The microphysical schemes also differ in the surface evaporation rate. The WRF single-moment 3-class scheme has the highest surface evaporation rate compensated by the highest precipitation rate. The different distributions of hydrometeors and water vapour of the microphysical schemes induce differences up to 49 W m −2 in the downwelling shortwave radiation and up to 33 W m −2 in the downwelling longwave radiation.

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: 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: AnaWEGE: a weather generator based on analogues of atmospheric circulation Geoscientific Model Development Discussions, 6, 4745-4774, 2013 Author(s): P. Yiou This paper presents a stochastic weather generator based on analogues of circulation (AnaWEGE). Analogues of circulation have been a promising paradigm to analyse climate variability and its extremes. The weather generator uses precomputed analogues of sea-level pressure over the North Atlantic. The stochastic rules of the generator constrain the continuity in time of the simulations. The generator then simulates spatially coherent time series of a climate variable, drawn from meteorological observations. The weather generator is tested for European temperatures, and for winter and summer seasons. The biases in temperature quantiles and autocorrelation are rather small compared to observed variability. The ability of simulating extremely hot summers and cold winters is also assessed.

Description: Methodological aspects of a pattern-scaling approach to produce global fields of monthly means of daily maximum and minimum temperature Geoscientific Model Development Discussions, 6, 4833-4882, 2013 Author(s): S. Kremser, G. E. Bodeker, and J. Lewis A Climate Pattern-Scaling Model (CPSM) that simulates global patterns of climate change, for a prescribed emissions scenario, is described. A CPSM works by quantitatively establishing the statistical relationship between a climate variable at a specific location (e.g. daily maximum surface temperature, T max ) and one or more predictor time series (e.g. global mean surface temperature, T global ) – referred to as the "training" of the CPSM. This training uses a regression model to derive fit-coefficients that describe the statistical relationship between the predictor time series and the target climate variable time series. Once that relationship has been determined, and given the predictor time series for any greenhouse gas (GHG) emissions scenario, the change in the climate variable of interest can be reconstructed – referred to as the "application" of the CPSM. The advantage of using a CPSM rather than a typical atmosphere-ocean global climate model (AOGCM) is that the predictor time series required by the CPSM can usually be generated quickly using a simple climate model (SCM) for any prescribed GHG emissions scenario and then applied to generate global fields of the climate variable of interest. The training can be performed either on historical measurements or on output from an AOGCM. Using model output from 21st century simulations has the advantage that the climate change signal is more pronounced than in historical data and therefore a more robust statistical relationship is obtained. The disadvantage of using AOGCM output is that the CPSM training might be compromised by any AOGCM inadequacies. For the purposes of exploring the various methodological aspects of the CPSM approach, AOGCM output was used in this study to train the CPSM. These investigations of the CPSM methodology focus on monthly mean fields of daily temperature extremes ( T max and T min ). Key conclusions are: (1) overall, the CPSM trained on simulations based on the Representative Concentration Pathway (RCP) 8.5 emissions scenario is able to reproduce AOGCM simulations of T max and T min based on predictor time series from an RCP 4.5 emissions scenario; (2) access to hemisphere average land and ocean temperatures as predictors improves the variance that can be explained, particularly over the oceans; (3) regression model fit-coefficients derived from individual simulations based on the RCP 2.6, 4.5 and 8.5 emissions scenarios agree well over most regions of the globe (the Arctic is the exception); (4) training the CPSM on concatenated time series from an ensemble of simulations does not result in fit-coefficients that explain significantly more of the variance than an approach that weights results based on single simulation fits; and (5) the inclusion of a linear time dependence in the regression model fit-coefficients improves the variance explained, primarily over the oceans.

Description: Modeling different freeze/thaw processes in heterogeneous landscapes of the Arctic polygonal tundra using an ecosystem model Geoscientific Model Development Discussions, 6, 4883-4932, 2013 Author(s): given_name prefix surname suffix, S. Yi, K. Wischnewski, M. Langer, S. Muster, and J. Boike Freeze/thaw (F/T) processes can be quite different under the various land surface types found in the heterogeneous polygonal tundra of the Arctic. Proper simulation of these different processes is essential for accurate prediction of the release of greenhouse gases under a warming climate scenario. In this study we have modified the dynamic organic soil version of the Terrestrial Ecosystem Model (DOS-TEM) to simulate F/T processes beneath the polygon rims, polygon centers (with and without water), and lakes that are common features in Arctic lowland regions. We first verified the F/T algorithm in the DOS-TEM against analytical solutions, and then compared the results with in situ measurements from Samoylov Island, Siberia. In the final stage, we examined the different responses of the F/T processes for different water levels at the various land surface types. The simulations revealed that (1) the DOS-TEM was very efficient and its results compared very well with analytical solutions for idealized cases, (2) the simulations compared reasonably well with in situ measurements although there were a number of model limitations and uncertainties, (3) the DOS-TEM was able to successfully simulate the differences in F/T dynamics under different land surface types, and (4) permafrost beneath water bodies was found to respond highly sensitive to changes in water depths between 1 and 2 m. Our results indicate that water is very important in the thermal processes simulated by the DOS-TEM; the heterogeneous nature of the landscape and different water depths therefore need to be taken into account when simulating methane emission responses to a warming climate.

Description: A standard test case suite for two-dimensional linear transport on the sphere: results from a collection of state-of-the-art schemes Geoscientific Model Development Discussions, 6, 4983-5076, 2013 Author(s): P. H. Lauritzen, P. A. Ullrich, C. Jablonowski, P. A. Bosler, D. Calhoun, A. J. Conley, T. Enomoto, L. Dong, S. Dubey, O. Guba, A. B. Hansen, E. Kaas, J. Kent, J.-F. Lamarque, M. J. Prather, D. Reinert, V. V. Shashkin, W. C. Skamarock, B. Sørensen, M. A. Taylor, and M. A. Tolstykh Recently, a standard test case suite for 2-D linear transport on the sphere was proposed to assess important aspects of accuracy in geophysical fluid dynamics with a "minimal" set of idealized model configurations/runs/diagnostics. Here we present results from 19 state-of-the-art transport scheme formulations based on finite-difference/finite-volume methods as well as emerging (in the context of atmospheric/oceanographic sciences) Galerkin methods. Discretization grids range from traditional regular latitude-longitude grids to more isotropic domain discretizations such as icosahedral and cubed-sphere tessellations of the sphere. The schemes are evaluated using a wide range of diagnostics in idealized flow environments. Accuracy is assessed in single- and two-tracer configurations using conventional error norms as well as novel diagnostics designed for climate and climate-chemistry applications. In addition, algorithmic considerations that may be important for computational efficiency are reported on. The latter is inevitably computing platform dependent, The ensemble of results from a wide variety of schemes presented here helps shed light on the ability of the test case suite diagnostics and flow settings to discriminate between algorithms and provide insights into accuracy in the context of global atmospheric/ocean modeling. A library of benchmark results is provided to facilitate scheme intercomparison and model development. Simple software and data-sets are made available to facilitate the process of model evaluation and scheme intercomparison.

Description: On the analytic approximation of bulk collision rates of non-spherical hydrometeors Geoscientific Model Development Discussions, 6, 5077-5116, 2013 Author(s): A. Seifert, U. Blahak, and R. Buhr Analytic approximations of the binary collision rates of hydrometeors are derived for use in bulk microphysical parameterizations. Special attention is given to non-spherical hydrometeors like raindrops and snowflakes. The terminal fall velocity of these particles cannot be sufficiently well approximated by power law relations which are used in most microphysical parameterizations and therefore an improved formulation is needed. The analytic approximations of the bulk collision rates given in this paper are an alternative to look-up tables and can replace the Wisner approximation which is used in many atmospheric models.

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: Ensemble initialization of the oceanic component of a coupled model through bred vectors at seasonal-to-interannual time scales Geoscientific Model Development Discussions, 6, 5189-5214, 2013 Author(s): J. Baehr and R. Piontek We evaluate the ensemble spread at seasonal-to-interannual time scales for two perturbation techniques implemented into the ocean component of a coupled model: (1) lagged initial conditions as commonly used for decadal predictions, (2) bred vectors as commonly used for weather and seasonal forecasting. We show that relative to an uninitialized reference simulation the implementation for bred vectors can improve the ensemble spread compared to lagged initialization at time scales from one months up to three years. As bred vectors have so far mostly been used at short time scales, we initially focus on the implementation of the bred vectors into the ocean component. We introduce a depth-dependent vertical rescaling norm, accounting for the vertical dependence of the variability, and extending the commonly used upper-ocean rescaling norm to the full water column. We further show that it is sufficient for the (sub-surface) ocean to breed temperature and salinity (i.e., scalar quantities), and rely on the governing physics to carry the temperature and salinity perturbations to the flow field. Using these bred vectors with a rescaling interval of 12 months, we initialize hindcast simulations and compare them to hindcast simulations initialized with lagged initial conditions. We quantify the ensemble spread by analyzing Talagrand diagrams and spread-error ratios. For both temperature and salinity, the lagged initialized ensemble is particularly under-dispersive for the first few months of predictable lead time. The ensemble initialized with bred vectors improves the spread for temperature and salinity for the 0–700 m and 1000–3500 m means, compared to the lagged ensemble at lead times of several months to one year. As the lead time increases to years, the differences between the two ensemble initialization techniques becomes more difficult to discern. While the results need to be confirmed in an initialized framework, the present analysis represents a first step towards an improved ensemble generation at the transition from seasonal-to-interannual time scales, in particular at lead times up to one year.

Description: Understanding the performance of the FLake model over the African Great Lakes Geoscientific Model Development Discussions, 6, 5141-5187, 2013 Author(s): W. Thiery, A. Martynov, F. Darchambeau, J.-P. Descy, P.-D. Plisnier, L. Sushama, and N. P. M. van Lipzig The ability of the one-dimensional lake model FLake to represent the mixolimnion temperatures for tropical conditions was tested for three locations in East Africa: Lake Kivu, Lake Tanganyika's northern and southern basins. Meteorological observations from surrounding Automatic Weather Stations were corrected and used to drive FLake, whereas a comprehensive set of water temperature profiles served to evaluate the model at each site. Careful forcing data correction and model configuration allowed to reproduce the observed mixed layer seasonality at Lake Kivu and Lake Tanganyika (northern and southern basins), with correct representation of both the mixed layer depth and temperature structure. At Lake Kivu, mixolimnion temperatures predicted by FLake were found sensitive both to minimal variations in the external parameters (lake depth and water transparency) as to small changes in the meteorological driving data, in particular wind velocity. In each case, small modifications may already lead to a regime switch from the correctly represented seasonal mixed layer deepening to either completely mixed (down to the model lake bottom) or permanently stratified (from ~10 m downwards) conditions. In contrast, model temperatures are found robust close to the surface, with acceptable predictions of near-surface water temperatures even when the seasonal mixing regime is not reproduced. FLake can thus be a suitable tool to parameterize tropical lake water surface temperatures within atmospheric prediction models, but may be less appropriate, in its current form, to study complex limnological processes within tropical lakes. Furthermore, a study of different initial conditions showed that for tropical lakes lacking reliable initial data, a fully mixed, artificially warm initialisation is to be preferred, but only if the model is allowed to spin up until convergence is reached. Finally, FLake was used to attribute the seasonal mixing cycle at Lake Kivu to variations in the near-surface meteorological conditions. It was found that the annual mixing down to 60 m during the main dry season is primarily due to enhanced lake evaporation and secondarily due to the decreased incoming long wave radiation, both causing a significant heat loss from the lake surface and associated mixolimnion cooling.

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: 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: 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: 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: The mid-Pliocene climate simulated by FGOALS-g2 Geoscientific Model Development Discussions, 6, 2403-2428, 2013 Author(s): W. Zheng, Z. Zhang, L. Chen, and Y. Yu Within the framework of Pliocene Model Intercomparison Project (PlioMIP), the mid-Pliocene (3.264–3.025 Ma) climate simulated by the Flexible Global Ocean-Atmosphere-Land System model grid-point version 2 (FGOALS-g2) are analyzed in this study. Results show that the model reproduces the large-scale features of the global warming over the land and ocean. The simulated mid-Pliocene global annual mean surface air temperature (TAS) and sea surface temperature (SST) are 4.17 and 2.62°C warmer than the pre-Industrial simulation, respectively. In particular, the feature of larger warming over mid-high latitudes is well captured. In the simulated warm mid-Pliocene climate, the Atlantic Meridional Overturning Circulation (AMOC) and El Niño-Southern Oscillation (ENSO) become weaker.

Description: The potential of an observational data set for calibration of a computationally expensive computer model Geoscientific Model Development Discussions, 6, 2369-2401, 2013 Author(s): D. J. McNeall, P. G. Challenor, J. R. Gattiker, and E. J. Stone We measure the potential of an observational data set to constrain a set of inputs to a complex and computationally expensive computer model. We use each member in turn of an ensemble of output from a computationally expensive model, corresponding to some observable part of a modelled system, as a proxy for an observational data set. We argue that our ability to constrain uncertain parameter inputs to a model using its own output as data, provides a maximum bound for our ability to constrain the model inputs using observations of the real system. The ensemble provides a set of known parameter input and model output pairs, which we use to build a computationally efficient statistical proxy for the full computer model, termed an emulator. We use the emulator to find and rule out ''implausible" values for the inputs of held-out ensemble members, given the computer model output. As we know the true values of the inputs for the ensemble, we can compare our constraint of the model inputs with the true value of the input for any ensemble member. Measures of the quality of constraint have the potential to inform strategy for data collection campaigns, before any real-world data is collected, as well as acting as an effective sensitivity analysis. We use an ensemble of the ice sheet model Glimmer to demonstrate our measures of quality of constraint. The ensemble has 250 model runs with 5 uncertain input parameters, and an output variable representing the pattern of the thickness of ice over Greenland. We have an observation of historical ice sheet thickness that directly matches the output variable, and offers an opportunity to constrain the model. We show that different ways of summarising our output variable (ice volume, ice surface area and maximum ice thickness) offer different potential constraints on individual input parameters. We show that combining the observational data gives increased power to constrain the model. We investigate the impact of uncertainty in observations or in model biases on our measures, showing that even a modest uncertainty can seriously degrade the potential of the observational data to constrain the model.

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: 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: On searching for optimized set of physical parameterization schemes in a multi-physics land surface process model Geoscientific Model Development Discussions, 6, 4511-4530, 2013 Author(s): S. Hong, X. Yu, S. K. Park, Y.-S. Choi, and B. Myoung Optimization of land surface models has been very challenging due to the increasing complexity of such models. Typical parameter calibration techniques often limit the solution of the spatiotemporal discrepancy in the modeling performance levels especially for regional applications. Thus, in this study, an attempt was made to perform scheme-based model optimization by designing a framework for coupling a micro-genetic algorithm (micro-GA) with the Noah land surface model that has multiple physics options (Noah-MP). Micro-GA controls the scheme selections in 10 different land surface parameterization fields in Noah-MP in order to extract the optimal scheme combination for a certain region. This coupling framework was successfully applied to the optimization of the surface water partitioning in the Korean Peninsula, promising not only the effectiveness of the scheme-based optimization but also model diagnosis capability by exploring the scheme sensitivity during the micro-GA evolution process. Then, the method was applied to four different regions in East Asia that have different climatic characteristics. The results indicate that (1) the optimal scheme combinations vary with the regions, (2) schemes related to the surface water partitioning are important for the modeling accuracy, and (3) specialized post-parameter optimization for each region may be required.

Description: Atmospheric inverse modeling with known physical bounds: an example from trace gas emissions Geoscientific Model Development Discussions, 6, 4531-4562, 2013 Author(s): S. M. Miller, A. M. Michalak, and P. J. Levi Many inverse problems in the atmospheric sciences involve parameters with known physical constraints. Examples include non-negativity (e.g., emissions of some urban air pollutants) or upward limits implied by reaction or solubility constants. However, probabilistic inverse modeling approaches based on Gaussian assumptions cannot incorporate such bounds and thus often produce unrealistic results. The atmospheric literature lacks consensus on the best means to overcome this problem, and existing atmospheric studies rely on a limited number of the possible methods with little examination of the relative merits of each. This paper investigates the applicability of several approaches to bounded inverse problems and is also the first application of Markov chain Monte Carlo (MCMC) to estimation of atmospheric trace gas fluxes. The approaches discussed here are broadly applicable. A common method of data transformations is found to unrealistically skew estimates for the examined example application. The method of Lagrange multipliers and two MCMC methods yield more realistic and accurate results. In general, the examined MCMC approaches produce the most realistic result but can require substantial computational time. Lagrange multipliers offer an appealing alternative for large, computationally intensive problems when exact uncertainty bounds are less central to the analysis. A synthetic data inversion of US anthropogenic methane emissions illustrates the strengths and weaknesses of each approach.

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: A database and tool for boundary conditions for regional air quality modeling: description and evaluation Geoscientific Model Development Discussions, 6, 4665-4704, 2013 Author(s): B. H. Henderson, F. Akhtar, H. O. T. Pye, S. L. Napelenok, and W. T. Hutzell Transported air pollutants receive increasing attention as regulations tighten and global concentrations increase. The need to represent international transport in regional air quality assessments requires improved representation of boundary concentrations. Currently available observations are too sparse vertically to provide boundary information, particularly for ozone precursors, but global simulations can be used to generate spatially and temporally varying Lateral Boundary Conditions (LBC). This study presents a public database of global simulations designed and evaluated for use as LBC for air quality models (AQMs). The database covers the contiguous United States (CONUS) for the years 2000–2010 and contains hourly varying concentrations of ozone, aerosols, and their precursors. The database is complimented by a tool for configuring the global results as inputs to regional scale models (e.g., Community Multiscale Air Quality or Comprehensive Air quality Model with extensions). This study also presents an example application based on the CONUS domain, which is evaluated against satellite retrieved ozone vertical profiles. The results show performance is largely within uncertainty estimates for the Tropospheric Emission Spectrometer (TES) with some exceptions. The major difference shows a high bias in the upper troposphere along the southern boundary in January. This publication documents the global simulation database, the tool for conversion to LBC, and the fidelity of concentrations on the boundaries. This documentation is intended to support applications that require representation of long-range transport of air pollutants.

Description: Effects of vegetation structure on biomass accumulation in a Balanced Optimality Structure Vegetation Model (BOSVM v1.0) Geoscientific Model Development Discussions, 6, 4603-4663, 2013 Author(s): Z. Yin, S. C. Dekker, B. J. J. M. van den Hurk, and H. A. Dijkstra A myriad of interactions exist between vegetation and local climate for arid and semi-arid regions. Vegetation function, structure and individual behavior have large impacts on carbon-water-energy balances, which consequently influence local climate variability that, in turn, feeds back to the vegetation. In this study, a conceptual vegetation structure scheme is formulated and tested in a new carbon-water-energy coupled model to explore the importance of vegetation structure and vegetation adaptation to water stress on equilibrium biomass states. Surface energy, water and carbon fluxes are simulated for a range of vegetation structures across a precipitation gradient in West Africa and optimal vegetation structures that maximizes biomass for each precipitation regime are determined. Two different strategies of vegetation adaptation to water stress are included. Under dry conditions vegetation tries to maximize the Water Use Efficiency and Leaf Area Index as it tries to maximize carbon gain. However, an important negative feedback mechanism is found as the vegetation also tries to minimize its cover to optimize the surrounding bare ground area from which water can be extracted, thereby forming patches of vertical vegetation. Under larger precipitation, a positive feedback mechanism is found in which vegetation tries to maximize its cover as it then can reduce water loss from bare soil while having maximum carbon gain due to a large Leaf Area Index. The competition between vegetation and bare soil determines a transition between a "survival" state to a "growing" state.

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: A distributed computing approach to improve the performance of the Parallel Ocean Program (v2.1) Geoscientific Model Development Discussions, 6, 4705-4744, 2013 Author(s): B. van Werkhoven, J. Maassen, M. Kliphuis, H. A. Dijkstra, S. E. Brunnabend, M. van Meersbergen, F. J. Seinstra, and H. E. Bal The Parallel Ocean Program (POP) is used in many strongly eddying ocean circulation simulations. Ideally one would like to do thousand-year long simulations, but the current performance of POP prohibits this type of simulations. In this work, using a new distributed computing approach, two innovations to improve the performance of POP are presented. The first is a new block partitioning scheme for the optimization of the load balancing of POP such that it can be run efficiently in a multi-platform setting. The second is an implementation of part of the POP model code on Graphics Processing Units. We show that the combination of both innovations leads to a substantial performance increase also when running POP simultaneously over multiple computational platforms.

Description: Inherently mass-conservative version of the semi-Lagrangian Absolute Vorticity (SL-AV) atmospheric model dynamical core Geoscientific Model Development Discussions, 6, 4809-4832, 2013 Author(s): V. V. Shashkin and M. A. Tolstykh The semi-Lagrangian Absolute Vorticity (SL-AV) atmospheric model is the global semi-Lagrangian hydrostatic model used for operational medium-range and seasonal forecasts at Hydrometeorological centre of Russia. The distinct feature of SL-AV dynamical core is the semi-implicit semi-Lagrangian vorticity-divergence formulation on the unstaggered grid. Semi-implicit semi-Lagrangian approach allows for long time steps while violates the global and local mass-conservation. In particular, the total mass in simulations with semi-Lagrangian models can drift significantly if no aposteriori mass-fixing algorithms are applied. However, the global mass-fixing algorithms degrade the local mass conservation. The inherently mass-conservative version of SL-AV model dynamical core presented in the article ensures global and local mass conservation without mass-fixing algorithms. The mass conservation is achieved with the introduction of the finite-volume semi-Lagrangian discretization for continuity equation based on the 3-D extension of the conservative cascade semi-Lagrangian transport scheme (CCS). The numerical experiments show that the presented new version of SL-AV dynamical core combines the accuracy and stability of the standard SL-AV dynamical core with the mass-conservation properties. The results of the mountain induced Rossby wave test and baroclinic instability test for mass-conservative dynamical core are found to be in agreement with the results available in literature.

Description: A fast input/output library for high resolution climate models Geoscientific Model Development Discussions, 6, 4775-4807, 2013 Author(s): X. Huang, W. Wang, H. Fu, G. Yang, B. Wang, and C. Zhang We describe the design and implementation of Climate Fast Input/Output (CFIO), a fast input/output (I/O) library for high resolution climate models. CFIO provides a simple method for modelers to overlap the I/O phase with the computing phase automatically, so as to shorten the running time of numerical simulations. To minimize the code modifications required for porting, CFIO provides similar interfaces and features to Parallel network Common Data Form (PnetCDF), which is one of the most widely used I/O libraries in climate models. We deployed CFIO in three high resolution climate models, including two ocean models (POP and LICOM) and one sea ice model (CICE). The experimental results show that CFIO improves the performance of climate models significantly versus the original serial I/O approach. When running with CFIO at 0.1° resolution with about 1000 CPU cores, we managed to reduce the running time by factors of 7.9, 4.6 and 2.0 for POP, CICE, and LICOM respectively. We also compared the performance of CFIO against PnetCDF in different scenarios. For scenarios with both data output operations and computations, CFIO decreases the I/O overhead by a factor of 5.1 compared to PnetCDF.

Description: Application and evaluation of McICA scheme with new radiation code in BCC_AGCM2.0.1 Geoscientific Model Development Discussions, 6, 4933-4982, 2013 Author(s): given_name prefix surname suffix, H. Zhang, X. Jing, and J. Li This research incorporates the Monte Carlo Independent Column Approximation (McICA) scheme with the correlated k-distribution BCC-RAD radiation model into the climate model BCC_AGCM2.0.1 and examines the impacts on modeled climate through several simulations with variations in cloud structures. Results from experiments with consistent sub-grid cloud structures show that both clear-sky radiation fluxes and cloud radiative forcings (CRFs) calculated by the new scheme are mostly improved relative to those calculated from the original one. The modeled atmospheric temperature and specific humidity are also improved due to changes in the radiative heating rates. The vertical overlap of fractional clouds and horizontal distribution of cloud condensation are important for computing CRFs. The maximum changes in seasonal CRF using the general overlap assumption (GenO) with different decorrelation depths ( L cf ) are larger than 10 and 20 Wm 2 for longwave (LW) CRF and shortwave (SW) CRF, respectively, mostly located in the Tropics and mid-latitude storm tracks. Larger (smaller) L cf in the Tropics (mid-latitude storm tracks) yield better cloud fraction and CRF compared with observations. The inclusion of an observation-based horizontal inhomogeneity of cloud condensation has a distinct impact on LW CRF and SW CRF, with global means of ∼1.2 Wm −2 and ∼3.7 Wm −2 at the top of atmosphere, respectively, making these much closer to observations. These results prove the reliability of the new model configuration to be used in BCC_AGCM2.0.1 for climate simulations, and also indicate that more detailed real-world information on cloud structures should be obtained to constrain cloud settings in McICA in the future.

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: 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: A regional climate modelling projection ensemble experiment – NARCliM Geoscientific Model Development Discussions, 6, 5117-5139, 2013 Author(s): J. P. Evans, F. Ji, C. Lee, P. Smith, D. Argüeso, and L. Fita Including the impacts of climate change in decision making and planning processes is a challenge facing many regional governments including the New South Wales (NSW) and Australian Capital Territory (ACT) governments in Australia. NARCliM (NSW/ACT Regional Climate Modelling project) is a regional climate modelling project that aims to provide a comprehensive and consistent set of climate projections that can be used by all relevant government departments when considering climate change. To maximise end user engagement and ensure outputs are relevant to the planning process, a series of stakeholder workshops were run to define key aspects of the model experiment including spatial resolution, time slices, and output variables. As with all such experiments, practical considerations limit the number of ensembles members that can be simulated such that choices must be made concerning which Global Climate Models (GCMs) to downscale from, and which Regional Climate Models (RCMs) to downscale with. Here a methodology for making these choices is proposed that aims to sample the uncertainty in both GCMs and RCMs, as well as spanning the range of future climate projections present in the full GCM ensemble. The created ensemble provides a more robust view of future regional climate changes.

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: Efficient performance of the Met Office Unified Model v8.2 on Intel Xeon partially used nodes Geoscientific Model Development Discussions, 7, 7395-7425, 2014 Author(s): I. Bermous The atmospheric Unified Model (UM) developed at the UK Met Office is used for weather and climate prediction by forecast teams at a number of international meteorological centres and research institutes on a wide variety of hardware and software environments. Over its 25 year history the UM sources have been optimised for a better application performance on a number of HPC systems including NEC SX vector architecture systems and recently the IBM Power6/Power7 platforms. Understanding the influence of the compiler flags, MPI libraries and run configurations is crucial to achieving the shortest elapsed times for a UM application on any particular HPC system. These aspects are very important for applications that must run within operational time frames. Driving the current study is the HPC industry trend since 1980 for processor arithmetic performance to increase at a faster rate than memory bandwidth. This gap has been growing especially fast for multicore processors in the past 10 years and it can have significant implication for the performance and performance scaling of memory bandwidth intensive applications, such as the UM. Analysis of partially used nodes on Intel Xeon clusters is provided in this paper for short and medium range weather forecasting systems using global and limited-area configurations. It is shown that on the Intel Xeon based clusters the fastest elapsed times and the most efficient system usage can be achieved using partially committed nodes.

Description: Simulation of groundwater and surface water over the continental US using a hyperresolution, integrated hydrologic model Geoscientific Model Development Discussions, 7, 7317-7349, 2014 Author(s): R. M. Maxwell, L. E. Condon, and S. J. Kollet Interactions between surface and groundwater systems are well-established theoretically and observationally. While numerical models that solve both surface and subsurface flow equations in a single framework (matrix) are increasingly being applied, computational limitations have restricted their use to local and regional studies. Regional or watershed, scale simulations have been effective tools in understanding hydrologic processes, however there are still many questions, such as the adaptation of water resources to anthropogenic stressors and climate variability, that need to be answered across large spatial extents at high resolution. In response to this "grand challenge" in hydrology, we present the results of a parallel, integrated hydrologic model simulating surface and subsurface flow at high spatial resolution (1 km) over much of continental North America (~ 6 300 000 or 6.3 million km 2 ). These simulations provide predictions of hydrologic states and fluxes, namely water table depth and streamflow, at unprecedented scale and resolution. The physically-based modeling approach used here requires limited parameterizations and relies only on more fundamental inputs, such as topography, hydrogeologic properties and climate forcing. Results are compared to observations and provide mechanistic insight into hydrologic process interaction. This study demonstrates both the feasibility of continental scale integrated models and their utility for improving our understanding of large-scale hydrologic systems; the combination of high resolution and large spatial extent facilitates novel analysis of scaling relationships using model outputs.

Description: EDDA: integrated simulation of debris flow erosion, deposition and property changes Geoscientific Model Development Discussions, 7, 7267-7316, 2014 Author(s): H. X. Chen and L. M. Zhang Debris flow material properties change during the initiation, transportation and deposition processes, which influences the runout characteristics of the debris flow. A quasi-three-dimensional depth-integrated numerical model, EDDA, is presented in this paper to simulate debris flow erosion, deposition and induced material property changes. The model considers changes in debris flow density, yield stress and dynamic viscosity during the flow process. The yield stress of debris flow mixture is determined at limit equilibrium using the Mohr–Coulomb equation, which is applicable to clear water flow, hyper-concentrated flow and fully developed debris flow. To assure numerical stability and computational efficiency at the same time, a variable time stepping algorithm is developed to solve the governing differential equations. Four numerical tests are conducted to validate the model. The first two tests involve a one-dimensional dam-break water flow and a one-dimensional debris flow with constant properties. The last two tests involve erosion and deposition, and the movement of multi-directional debris flows. The changes in debris flow mass and properties due to either erosion or deposition are shown to affect the runout characteristics significantly. The model is also applied to simulate a large-scale debris flow in Xiaojiagou Ravine to test the performance of the model in catchment-scale simulations. The results suggest that the model estimates well the volume, inundated area, and runout distance of the debris flow. The model is intended for use as a module in a real-time debris flow warning system.

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: Northern Hemisphere storminess in the Norwegian Earth System Model (NorESM1-M) Geoscientific Model Development Discussions, 7, 8975-9015, 2014 Author(s): E. M. Knudsen and J. E. Walsh Metrics of storm activity in Northern Hemisphere high- and midlatitudes are evaluated from historical output and future projections by the Norwegian Earth System Model (NorESM1-M) coupled global climate model. The European Re-Analysis Interim (ERA-Interim) and the Community Climate System Model (CCSM4), a global climate model of the same vintage as NorESM1-M, provide benchmarks for comparison. The focus is on the autumn and early winter (September through December), the period when the ongoing and projected Arctic sea ice retreat is greatest. Storm tracks derived from a vorticity-based algorithm for storm identification are reproduced well by NorESM1-M, although the tracks are somewhat better resolved in the higher-resolution ERA-Interim and CCSM4. The tracks are projected to shift polewards in the future as climate changes under the Representative Concentration Pathway (RCP) forcing scenarios. Cyclones are projected to become generally more intense in the high-latitudes, especially over the Alaskan region, although in some other areas the intensity is projected to decrease. While projected changes in track density are less coherent, there is a general tendency towards less frequent storms in midlatitudes and more frequent storms in high-latitudes, especially the Baffin Bay/Davis Strait region. Autumn precipitation is projected to increase significantly across the entire high-latitudes. Together with the projected increases in storm intensity and sea level and the loss of sea ice, this increase in precipitation implies a greater vulnerability to coastal flooding and erosion, especially in the Alaskan region. The projected changes in storm intensity and precipitation (as well as sea ice and sea level pressure) scale generally linearly with the RCP value of the forcing and with time through the 21st century.

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: The 1-way on-line coupled atmospheric chemistry model system MECO(n) – Part 1: The limited-area atmospheric chemistry model COSMO/MESSy Geoscientific Model Development Discussions, 4, 1305-1358, 2011 Author(s): A. Kerkweg and P. Jöckel The numerical weather prediction model of the Consortium for Small Scale Modelling (COSMO), maintained by the German weather service (DWD), is connected with the Modular Earth Submodel System (MESSy). This effort is undertaken in preparation of a~new, limited-area atmospheric chemistry model. This model is as consistent as possible, with respect to atmospheric chemistry and related processes, with a previously developed global atmospheric chemistry general circulation model: the ECHAM/MESSy Atmospheric Chemistry (EMAC) model. The combined system constitutes a new research tool, bridging the global to the meso-γ scale for atmospheric chemistry research. MESSy provides the infrastructure and includes, among others, the process and diagnostic submodels for atmospheric chemistry simulations. Furthermore, MESSy is highly flexible allowing model setups with tailor made complexity, depending on the scientific question. Here, the connection of the MESSy infrastructure to the COSMO model is documented. Previously published prototype submodels for simplified tracer studies are generalised to be plugged-in and used in the global and the limited-area model. They are used to evaluate the tracer transport characteristics of the new COSMO/MESSy model system, an important prerequisite for future atmospheric chemistry applications. A supplementary document with further details on the technical implementation of the MESSy interface into COSMO with a complete list of modifications to the COSMO code is provided.

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: Carbon monoxide as a tracer for tropical troposphere to stratosphere transport in the Chemical Lagrangian Model of the Stratosphere (CLaMS) Geoscientific Model Development Discussions, 4, 1185-1211, 2011 Author(s): R. Pommrich, R. Müller, J.-U. Grooß, P. Konopka, G. Günther, H.-C. Pumphrey, S. Viciani, F. D'Amato, and M. Riese Variations in the mixing ratio of trace gases of tropospheric origin entering the stratosphere in the tropics are of interest for assessing both troposphere to stratosphere transport fluxes in the tropics and the impact on the composition of the tropical lower stratosphere of quasi-horizontal in-mixing into the tropical tropopause layer from the mid-latitude stratosphere. Here, we present a simplified chemistry scheme for the Chemical Lagrangian Model of the Stratosphere (CLaMS) for the simulation, at comparatively low numerical cost, of CO, ozone, and long-lived trace substances (CH 4 , N 2 O, CCl 3 F, and CO 2 ) in the lower tropical stratosphere. The boundary conditions at the ground are represented for the long-lived trace substances CH 4 , N 2 O, CCl 3 F, and CO 2 based on ground-based measurements. The boundary condition for CO in the free troposphere is deduced from MOPITT measurements. We find that the zonally averaged tropical CO anomaly patterns simulated by this model version of CLaMS are in good agreement with observations. The introduction of a new scheme in the ECMWF integrated forecast system (Tompkins et al., 2007) for the ice supersaturation after September 2006, results in a somewhat less good agreement between observed and simulated CO patterns in the tropical lower stratosphere after this date.

Description: The Met Office Unified Model Global Atmosphere 3.0/3.1 and JULES Global Land 3.0/3.1 configurations Geoscientific Model Development Discussions, 4, 1213-1271, 2011 Author(s): D. N. Walters, M. J. Best, A. C. Bushell, D. Copsey, J. M. Edwards, P. D. Falloon, C. M. Harris, A. P. Lock, J. C. Manners, C. J. Morcrette, M. J. Roberts, R. A. Stratton, S. Webster, J. M. Wilkinson, M. R. Willett, I. A. Boutle, P. D. Earnshaw, P. G. Hill, C. MacLachlan, G. M. Martin, W. Moufouma-Okia, M. D. Palmer, J. C. Petch, G. G. Rooney, A. A. Scaife, and K. D. Williams We describe Global Atmosphere 3.0 (GA3.0): a configuration of the Met Office Unified Model (MetUM) developed for use across climate research and weather prediction activities. GA3.0 has been formulated by converging the development paths of the Met Office's weather and climate global atmospheric model components such that wherever possible, atmospheric processes are modelled or parametrized seamlessly across spatial resolutions and timescales. This unified development process will provide the Met Office and its collaborators with regular releases of a configuration that has been evaluated, and can hence be applied, over a variety of modelling régimes. We also describe Global Land 3.0 (GL3.0): a configuration of the JULES community land surface model developed for use with GA3.0. This paper provides a comprehensive technical and scientific description of the GA3.0 and GL3.0 (and related GA3.1 and GL3.1) configurations and presents the results of some initial evaluations of their performance in various applications. It is to be the first in a series of papers describing each subsequent Global Atmosphere release; this will provide a single source of reference for established users and developers as well as researchers requiring access to a current, but trusted, global MetUM setup.

Description: FAMOUS, faster: using parallel computing techniques to accelerate the FAMOUS/HadCM3 climate model with a focus on the radiative transfer algorithm Geoscientific Model Development Discussions, 4, 1273-1303, 2011 Author(s): P. Hanappe, A. Beurivé, F. Laguzet, L. Steels, N. Bellouin, O. Boucher, Y. H. Yamazaki, T. Aina, and M. Allen We have optimised the atmospheric radiation algorithm of the FAMOUS climate model on several hardware platforms. The optimisation involved translating the Fortran code to C and restructuring the algorithm around the computation of a single air column. A task queue and a thread pool are used to distribute the computation to several processors. Finally, four air columns are packed together in a single data structure and computed simultaneously using Single Instruction Multiple Data operations. The modified algorithm runs more than 50 times faster on the CELL's Synergistic Processing Elements than on its main PowerPC processing element. On Intel-compatible processors, the new radiation code runs 4 times faster and on graphics processors, using OpenCL, more than 2.5 times faster, as compared to the original code. Because the radiation code takes more than 60 % of the total CPU time, FAMOUS executes more than twice as fast. Our version of the algorithm returns bit-wise identical results, which demonstrates the robustness of our approach.

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: Analyzing numerics of bulk microphysics schemes in Community models: warm rain processes Geoscientific Model Development Discussions, 4, 1403-1434, 2011 Author(s): I. Sednev and S. Menon In the last decade there has been only one study that discussed time integration scheme (TIS) applied to advance governing differential equations in bulk microphysics (BLK) schemes. Recently, Morrison and Gettelman (2008) examine numerical aspects of double-moment BLK scheme with diagnostic treatment of precipitating hydrometeors implemented into Community Atmosphere Model, version 3 (CAM) to find an acceptable level of accuracy and numerical stability. However, stability condition for their explicit non-positive definite TIS was not defined. It is conventionally thought that the Weather Research and Forecasting (WRF) model can be applied for a broad range of spatial scales from large eddy up to global scale simulations if time steps used for model integration satisfy to a certain limit imposed mainly by dynamics. However, numerics used in WRF BLK schemes has never been analyzed in detail. To improve creditability of BLK schemes we derive a general analytical stability and positive definiteness criteria for explicit Eulerian time integration scheme used to advanced finite-difference equations that govern warm rain formation processes in microphysics packages in Community models (CAM and WRF) and define well-behaved, conditionally well-behaved, and non-well-behaved Explicit Eulerian Bulk Microphysics Code (EEBMPC) classes. We highlight that source codes of BLK schemes, originally developed for use in cloud-resolving models, implemented in Community models belong to conditionally well-behaved EEBMPC class and exhibit better performance for finer spatial resolutions when time steps do not exceed seconds or tenths of seconds. For coarser spatial resolutions used in regional and global scale simulations time steps are usually increased from hundredths up to thousands of seconds. This might lead to a degradation of conditionally well-behaved EEBMPCs ability to calculate the amount of precipitation as well as its spatial and temporal distribution since both stability and positive definiteness conditions are not met in the TIS. The correction through the so called “mass conservation” technique commonly used in many models with bulk microphysics is a main characteristic of non-well-behaved EEBMPC, whose utilization leads to erroneous conclusions regarding relative importance of different microphysical processes. Moreover, surface boundary conditions for ocean, land, lake, and sea ice models are dependent on the precipitation and its spatial and temporal distribution. Uncertainties in calculations of temporal and spatial patterns of accumulated precipitation influence the global water cycle. In fact, numerics in non-well-behaved EEBMPCs, which are used in Community Earth System Model, act as a hidden climate forcing agent, if relatively long time steps are used for the host model integration. By analyzing numerics of warm rain processes in EEBMPCs implemented in Community models we provide general guidelines regarding appropriate choice of integration time steps for use in these models.

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 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: The 1-way on-line coupled atmospheric chemistry model system MECO(n) – Part 2: On-line coupling Geoscientific Model Development Discussions, 4, 1359-1402, 2011 Author(s): A. Kerkweg and P. Jöckel A new, highly flexible model system for the seamless dynamical down-scaling of meteorological and chemical processes from the global to the meso-γ scale is presented. A global model and a cascade of an arbitrary number of limited-area model instances run concurrently in the same parallel environment, in which the coarser grained instances provide the boundary data for the finer grained instances. Thus, disk-space intensive and time consuming intermediate and pre-processing steps are entirely avoided and the time interpolation errors of common off-line nesting approaches are minimised. More specifically, the regional model COSMO of the German Weather Service (DWD) is nested on-line into the atmospheric general circulation model ECHAM5 within the Modular Earth Submodel System (MESSy) framework. ECHAM5 and COSMO have previously been equipped with the MESSy infrastructure, implying that the same process formulations (MESSy submodels) are available for both models. This guarantees the highest degree of achievable consistency, between both, the meteorological and chemical conditions at the domain boundaries of the nested limited-area model, and between the process formulations on all scales. The on-line nesting of the different models is established by a client-server approach with the newly developed Multi-Model-Driver (MMD), an additional component of the MESSy infrastructure. With MMD an arbitrary number of model instances can be run concurrently within the same message passing interface (MPI) environment, the respective coarser model (either global or regional) is the server for the nested finer (regional) client model, i.e., it provides the data required to calculate the initial and boundary fields to the client model. On-line nesting means that the coupled (client-server) models exchange their data via the computer memory, in contrast to the data exchange via files on disk in common off-line nesting approaches. MMD consists of a library (Fortran95 and some parts in C) which is based on the MPI standard and two new MESSy submodels, MMDSERV and MMDCLNT (both Fortran95) for the server and client models, respectively. MMDCLNT contains a further sub-submodel, INT2COSMO, for the interpolation of the coarse grid data provided by the server models (either ECHAM5/MESSy or COSMO/MESSy) to the grid of the respective client model (COSMO/MESSy). INT2COSMO is based on the off-line pre-processing tool INT2LM provided by the DWD. The new achievements allow the setup of model cascades for zooming (down-scaling) from the global scale to the lower edge of the meso-γ scale (≈1 km) with a very high degree of consistency between the different models and between the chemical and meteorological boundary conditions.

Description: Development and evaluation of an Earth-system model – HadGEM2 Geoscientific Model Development Discussions, 4, 997-1062, 2011 Author(s): W. J. Collins, N. Bellouin, M. Doutriaux-Boucher, N. Gedney, P. Halloran, T. Hinton, J. Hughes, C. D. Jones, M. Joshi, S. Liddicoat, G. Martin, F. O'Connor, J. Rae, C. Senior, S. Sitch, I. Totterdell, A. Wiltshire, and S. Woodward We describe here the development and evaluation of an Earth system model suitable for centennial-scale climate prediction. The principal new components added to the physical climate model are the terrestrial and ocean ecosystems and gas-phase tropospheric chemistry, along with their coupled interactions. The individual Earth system components are described briefly and the relevant interactions between the components are explained. Because the multiple interactions could lead to unstable feedbacks, we go through a careful process of model spin up to ensure that all components are stable and the interactions balanced. This spun-up configuration is evaluated against observed data for the Earth system components and is generally found to perform very satisfactorily. The reason for the evaluation phase is that the model is to be used for the core climate simulations carried out by the Met Office Hadley Centre for the Coupled Model Intercomparison Project (CMIP5), so it is essential that addition of the extra complexity does not detract substantially from its climate performance. Localised changes in some specific meteorological variables can be identified, but the impacts on the overall simulation of present day climate are slight. This model is proving valuable both for climate predictions, and for investigating the strengths of biogeochemical feedbacks.

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: 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 multi-layer land surface energy budget model for implicit coupling with global atmospheric simulations Geoscientific Model Development Discussions, 7, 8649-8701, 2014 Author(s): J. Ryder, J. Polcher, P. Peylin, C. Ottlé, Y. Chen, E. van Gorsel, V. Haverd, M. J. McGrath, K. Naudts, J. Otto, A. Valade, and S. Luyssaert In Earth system modelling, a description of the energy budget of the vegetated surface layer is fundamental as it determines the meteorological conditions in the planetary boundary layer and as such contributes to the atmospheric conditions and its circulation. The energy budget in most Earth system models has long been based on a "big-leaf approach", with averaging schemes that represent in-canopy processes. Such models have difficulties in reproducing consistently the energy balance in field observations. We here outline a newly developed numerical model for energy budget simulation, as a component of the land surface model ORCHIDEE-CAN (Organising Carbon and Hydrology In Dynamic Ecosystems – CANopy). This new model implements techniques from single-site canopy models in a practical way. It includes representation of in-canopy transport, a multilayer longwave radiation budget, height-specific calculation of aerodynamic and stomatal conductance, and interaction with the bare soil flux within the canopy space. Significantly, it avoids iterations over the height of tha canopy and so maintains implicit coupling to the atmospheric model LMDz. As a first test, the model is evaluated against data from both an intensive measurement campaign and longer term eddy covariance measurements for the intensively studied Eucalyptus stand at Tumbarumba, Australia. The model performs well in replicating both diurnal and annual cycles of fluxes, as well as the gradients of sensible heat fluxes. However, the model overestimates sensible heat flux against an underestimate of the radiation budget. Improved performance is expected through the implementation of a more detailed calculation of stand albedo and a more up-to-date stomatal conductance calculation.

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: 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: iGen: the automated generation of a parameterisation of entrainment in marine stratocumulus Geoscientific Model Development Discussions, 4, 971-995, 2011 Author(s): D. F. Tang and S. Dobbie In a previous paper we described a new technique for automatically generating parameterisations using a program called iGen. iGen generates parameterisations by analysing the source code of a high resolution model that resolves the physics to be parameterised. In order to demonstrate that this technique scales up to deal with models of realistic complexity we have used iGen to generate a parameterisation of entrainment in marine stratocumulus. We present details of our technique in which iGen was used to analyse the source code of a cloud resolving model and generate a parameterisation of the mean and standard deviation of entrainment velocity in marine stratocumulus in terms of the large-scale state of the boundary layer. The parameterisation was tested against results from the DYCOMS-II intercomparison of cloud resolving models and iGen's parameterisation of mean entrainment velocity was found to be 5.27 × 10 −3 ± 0.62 × 10 −3 m s −1 compared to 5.2 × 10 −3 ± 0.8 × 10 −3 m s −1 for the DYCOMS-II ensemble of cloud resolving models.

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: The HadGEM2-ES implementation of CMIP5 centennial simulations Geoscientific Model Development Discussions, 4, 689-763, 2011 Author(s): C. D. Jones, J. K. Hughes, N. Bellouin, S. C. Hardiman, G. S. Jones, J. Knight, S. Liddicoat, F. M. O'Connor, R. J. Andres, C. Bell, K.-O. Boo, A. Bozzo, N. Butchart, P. Cadule, K. D. Corbin, M. Doutriaux-Boucher, P. Friedlingstein, J. Gornall, L. Gray, P. R. Halloran, G. Hurtt, W. Ingram, J.-F. Lamarque, R. M. Law, M. Meinshausen, S. Osprey, E. J. Palin, L. Parsons Chini, T. Raddatz, M. Sanderson, A. A. Sellar, A. Schurer, P. Valdes, N. Wood, S. Woodward, M. Yoshioka, and M. Zerroukat The scientific understanding of the Earth's climate system, including the central question of how the climate system is likely to respond to human-induced perturbations, is comprehensively captured in GCMs and Earth System Models(ESM). Diagnosing the simulated climate response, and comparing responses across different models, is crucially dependent on transparent assumptions of how the GCM/ESM has been driven – especially because the implementation can involve subjective decisions and may differ between modelling groups performing the same experiment. This paper outlines the climate forcings and setup of the Met Office Hadley Centre ESM, HadGEM2-ES for the CMIP5 set of centennial experiments. We document the prescribed greenhouse gas concentrations, aerosol precursors, stratospheric and tropospheric ozone assumptions, as well as implementation of land-use change and natural forcings for the HadGEM2-ES historical and future experiments following the Representative Concentration Pathways. In addition, we provide details of how HadGEM2-ES ensemble members were initialised from the control run and how the palaeoclimate and AMIP experiments, as well as the "emission-driven" RCP experiments were performed.

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: The HadGEM2 family of Met Office Unified Model Climate configurations Geoscientific Model Development Discussions, 4, 765-841, 2011 Author(s): The HadGEM2 Development Team: Martin, G. M., N. Bellouin, W. J. Collins, I. D. Culverwell, P. R. Halloran, S. C. Hardiman, T. J. Hinton, C. D. Jones, R. E. McDonald, A. J. McLaren, F. M. O'Connor, M. J. Roberts, J. M. Rodriguez, S. Woodward, M. J. Best, M. E. Brooks, A. R. Brown, N. Butchart, C. Dearden, S. H. Derbyshire, I. Dharssi, M. Doutriaux-Boucher, J. M. Edwards, P. D. Falloon, N. Gedney, L. J. Gray, H. T. Hewitt, M. Hobson, M. R. Huddleston, J. Hughes, S. Ineson, W. J. Ingram, P. M. James, T. C. Johns, C. E. Johnson, A. Jones, C. P. Jones, M. M. Joshi, A. B. Keen, S. Liddicoat, A. P. Lock, A. V. Maidens, J. C. Manners, S. F. Milton, J. G. L. Rae, J. K. Ridley, A. Sellar, C. A. Senior, I. J. Totterdell, A. Verhoef, P. L. Vidale, and A. Wiltshire We describe the HadGEM2 family of climate configurations of the Met Office Unified Model, MetUM. The concept of a model "family" comprises a range of specific model configurations incorporating different levels of complexity but with a common physical framework. The HadGEM2 family of configurations includes atmosphere and ocean components, with and without a vertical extension to include a well-resolved stratosphere, and an Earth-System (ES) component which includes dynamic vegetation, ocean biology and atmospheric chemistry. The HadGEM2 physical model includes improvements designed to address specific systematic errors encountered in the previous climate configuration, HadGEM1, namely Northern Hemisphere continental temperature biases and tropical sea surface temperature biases and poor variability. Targeting these biases was crucial in order that the ES configuration could represent important biogeochemical climate feedbacks. Detailed descriptions and evaluations of particular HadGEM2 family members are included in a number of other publications, and the discussion here is limited to a summary of the overall performance using a set of model metrics which compare the way in which the various configurations simulate present-day climate and its variability.

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: Towards an online-coupled chemistry-climate model: evaluation of COSMO-ART Geoscientific Model Development Discussions, 4, 1809-1874, 2011 Author(s): C. Knote, D. Brunner, H. Vogel, J. Allan, A. Asmi, M. Äijälä, S. Carbone, H. D. van der Gon, J. L. Jimenez, A. Kiendler-Scharr, C. Mohr, L. Poulain, A. S. H. Prévôt, E. Swietlicki, and B. Vogel The online-coupled, regional chemistry transport model COSMO-ART is evaluated for periods in all seasons against several measurement datasets to assess its ability to represent gaseous pollutants and ambient aerosol characteristics over the European domain. Measurements used in the comparison include long-term station observations, satellite and ground-based remote sensing products, and complex datasets of aerosol chemical composition and number size distribution from recent field campaigns. This is the first time these comprehensive measurements of aerosol characteristics in Europe are used to evaluate a regional chemistry transport model. We show a detailed analysis of the simulated size-resolved chemical composition under different meteorological conditions. The model is able to represent trace gas concentrations with good accuracy and reproduces bulk aerosol properties rather well though with a clear tendency to underestimate both total mass (PM 10 and PM 2.5 ) and aerosol optical depth. We find indications of an overestimation of shipping emissions. Time evolution of aerosol chemical composition is captured, although some biases are found in relative composition. Nitrate aerosol components are on average overestimated, and sulfates underestimated. The accuracy of simulated organics depends strongly on season and location. While strongly underestimated during summer, organic mass is comparable in spring and autumn. We see indications for an overestimated fractional contribution of primary organic matter in urban areas and an underestimation of SOA at many locations. Aerosol number concentrations can be simulated well, size distributions are comparable. Our work sets the basis for subsequent studies of aerosol characteristics and climate impacts with COSMO-ART, and highlights areas where improvements are necessary for current regional modeling systems in general.

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: Mapping the climate: guidance on appropriate techniques to map climate variables and their uncertainty Geoscientific Model Development Discussions, 4, 1875-1906, 2011 Author(s): N. R. Kaye, A. Hartley, and D. Hemming Maps are a crucial asset in communicating climate science to a diverse audience, and there is a wealth of software available to analyse and visualise climate information. However, this availability makes it easy to create poor maps as users often lack an underlying cartographic knowledge. Furthermore, communicating and visualising uncertainties in climate data and climate change projections, using for example ensemble based approaches, presents additional challenges for mapping that require careful consideration. This paper assesses a range of techniques for mapping uncertainties, comparing "intrinsic" approaches that use colour in much the same way as conventional thematic maps, and "extrinsic" approaches that incorporate additional geometry such as points or features. We proposes that, unlike traditional cartography, where many known standards allow maps to be interpreted easily, there is no standard mapping approach used to represent uncertainty (in climate or other information). Consequently, a wide range of techniques have been applied for this purpose, and users may spend unnecessary time trying to understand the mapping approach rather than interpreting the information presented. We use cartographic knowledge and lessons learned from mapping other information to propose a suitable mapping technique that represents both magnitude and uncertainty in climate data. This technique adjusts the hue of a small palette of colours to show the mean or median of a climate variable, and the saturation of the colour to illustrate a measure of uncertainty. It is designed to be easy to replicate, visible to colour blind people and intuitive to understand. This technique may be utilised to map a wide range of climate data, and it is proposed that it would provide a consistent approach suitable for mapping information for the Fifth Assessment Report of the Intergovernmental Panel on Climate Change (IPCC AR5).

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: 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 new version of the CNRM Chemistry-Climate Model, CNRM-CCM: description and improvements from the CCMVal-2 simulations Geoscientific Model Development Discussions, 4, 1129-1183, 2011 Author(s): M. Michou, D. Saint-Martin, H. Teyssèdre, A. Alias, F. Karcher, D. Olivié, A. Voldoire, B. Josse, V.-H. Peuch, H. Clark, J. N. Lee, and F. Chéroux This paper presents a new version of the Météo-France CNRM Chemistry-Climate Model, so-called CNRM-CCM. It includes some fundamental changes from the previous version (CNRM-ACM) which was extensively evaluated in the context of the CCMVal-2 validation activity. The most notable changes concern the radiative code of the GCM, and the inclusion of the detailed stratospheric chemistry of our Chemistry-Transport model MOCAGE on-line within the GCM. A 47-yr transient simulation (1960–2006) is the basis of our analysis. CNRM-CCM generates satisfactory dynamical and chemical fields in the stratosphere. Several shortcomings of CNRM-ACM simulations for CCMVal-2 that resulted from an erroneous representation of the impact of volcanic aerosols as well as from transport deficiencies have been eliminated. Remaining problems concern the upper stratosphere (5 to 1 hPa) where temperatures are too high, and where there are biases in the NO 2 , N 2 O 5 and O 3 mixing ratios. In contrast, temperatures at the tropical tropopause are too cold. These issues are addressed through the implementation of a more accurate radiation scheme at short wavelengths. Despite these problems we show that this new CNRM CCM is a useful tool to study chemistry-climate applications.

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: Pliocene Model Intercomparison Project: implementation strategy and mid-Pliocene Global climatology using GENESIS v3.0 GCM Geoscientific Model Development Discussions, 4, 2577-2603, 2011 Author(s): S. J. Koenig, R. M. DeConto, and D. Pollard The mid-Pliocene Warm Period (3.29 to 2.97 Ma BP) has been identified as an analogue for the future, with the potential to help understand climate processes in a warmer than modern world. Sets of climate proxies, combined to provide boundary conditions for Global Climate Model (GCM) simulations of the mid-Pliocene, form the basis for the international, data-driven Pliocene Model Intercomparison Project (PlioMIP). Here, we outline the strategy for implementing pre-industrial (modern) and mid-Pliocene forcings and boundary conditions into the GENESIS version 3 GCM, as part of PlioMIP. We describe the prescription of greenhouse gas concentrations and orbital parameters and the implementation of geographic boundary conditions such as land-ice-sea distribution, topography, sea surface temperatures, sea ice extent, vegetation, soils, and ice sheets. We further describe model-specific details including spin-up and integration times. In addition, the global climatology of the mid-Pliocene as simulated by the GENESIS v3 GCM is analyzed and compared to the pre-industrial control simulation. The simulated climate of the mid-Pliocene warm interval is found to differ considerably from pre-industrial. We identify model sensitivity to imposed forcings, and internal feedbacks that collectively affect both local and far-field responses. Our analysis points out the need to assess both the direct impacts of external forcings and the combined effects of indirect, internal feedbacks. This paper provides the basis for assessing model biases within the PlioMIP framework, and will be useful for comparisons with other studies of mid-Pliocene climates.

Description: Development of the high-order decoupled direct method in three dimensions for particulate matter: enabling advanced sensitivity analysis in air quality models Geoscientific Model Development Discussions, 4, 2605-2633, 2011 Author(s): W. Zhang, S. L. Capps, Y. Hu, A. Nenes, S. L. Napelenok, and A. G. Russell The high-order decoupled direct method in three dimensions for particular matter (HDDM-3D/PM) has been implemented in the Community Multiscale Air Quality (CMAQ) model to enable advanced sensitivity analysis. The major effort of this work is to develop high-order DDM sensitivity analysis of ISORROPIA, the inorganic aerosol module of CMAQ. A case-specific approach has been applied, and the sensitivities of activity coefficients and water content are explicitly computed. Stand-alone tests are performed for ISORROPIA by comparing the sensitivities (first- and second-order) computed by HDDM and the brute force (BF) approximations. Similar comparison has also been carried out for CMAQ results simulated using a week-long winter episode for a continental US domain. Second-order sensitivities of aerosol species (e.g., sulfate, nitrate, and ammonium) with respect to domain-wide SO 2 , NO x , and NH 3 emissions show agreement with BF results, yet exhibit less noise in locations where BF results are demonstrably inaccurate. Second-order sensitivity analysis elucidates nonlinear responses of secondary inorganic aerosols to their precursors and competing species that have not yet been well-understood with other approaches. Including second-order sensitivity coefficients in the Taylor series projection of the nitrate concentrations with a 50% reduction in domain-wide NO x emission shows a statistically significant improvement compared to the first-order Taylor series projection.

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: Pliocene Ice Sheet Modelling Intercomparison Project (PLISMIP) – experimental design Geoscientific Model Development Discussions, 4, 2661-2686, 2011 Author(s): A. M. Dolan, S. J. Koenig, D. J. Hill, A. M. Haywood, and R. M. DeConto During the mid-Pliocene Warm Period (3.264 to 3.025 million yr ago), global mean temperature was similar to that predicted for the next century and atmospheric carbon dioxide concentrations were slightly higher. Sea level was also higher than today, implying a reduction in the extent of the ice sheets. Thus, the mid-Pliocene Warm Period provides a unique testing ground to investigate the stability of the Earth's ice sheets and their contribution to sea level in a warmer-than-modern world. Climate models and ice sheet models can be used to enhance our understanding of ice sheet stability, however, uncertainties associated with different ice-sheet modelling frameworks/approaches mean that a rigorous comparison of numerical ice sheet model simulations for the Pliocene is essential. As an extension to the Pliocene Model Intercomparison Project (PlioMIP; Haywood et al., 2010, 2011a), the Pliocene Ice Sheet Modelling Intercomparison Project (PLISMIP) will address these uncertainties. Here we outline the PLISMIP experimental design and initialisation conditions that have been adopted to simulate the Greenland and Antarctic ice sheets under present day and warm mid-Pliocene conditions. Not only will this project provide a new benchmark in the simulation of ice sheets in a past warm period, but the analysis of model sensitivity to various uncertainties could directly inform future predictions of ice sheet and sea level change.