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: Calibration of the passive cavity aerosol spectrometer probe for airborne determination of the size distribution Atmospheric Measurement Techniques, 6, 2349-2358, 2013 Author(s): Y. Cai, J. R. Snider, and P. Wechsler This work describes calibration methods for the particle sizing and particle concentration systems of the passive cavity aerosol spectrometer probe (PCASP). Laboratory calibrations conducted over six years, in support of the deployment of a PCASP on a cloud physics research aircraft, are analyzed. Instead of using the many calibration sizes recommended by the PCASP manufacturer, a relationship between particle diameter and scattered light intensity is established using three sizes of mobility-selected polystyrene latex particles, one for each amplifier gain stage. In addition, studies of two factors influencing the PCASP's determination of the particle size distribution – amplifier baseline and particle shape – are conducted. It is shown that the PCASP-derived size distribution is sensitive to adjustments of the sizing system's baseline voltage, and that for aggregates of spheres, a PCASP-derived particle size and a sphere-equivalent particle size agree within uncertainty dictated by the PCASP's sizing resolution. Robust determinations of aerosol concentration, and size distribution, also require calibration of the PCASP's aerosol flowrate sensor. Sensor calibrations, calibration drift, and the sensor's non-linear response are documented.

Description: Experimental quantification of contact freezing in an electrodynamic balance Atmospheric Measurement Techniques, 6, 2373-2382, 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 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 to determine the temperature-dependent ice nucleation probability of size-selected aerosol particles. The method is based on the suspension of supercooled charged water droplets 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 timescale of our experiment.

Description: Validation and empirical correction of MODIS AOT and AE over ocean Atmospheric Measurement Techniques, 6, 2455-2475, 2013 Author(s): N. A. J. Schutgens, M. Nakata, and T. Nakajima We present a validation study of Collection 5 MODIS level 2 Aqua and Terra AOT (aerosol optical thickness) and AE (Ångström exponent) over ocean by comparison to coastal and island AERONET (AErosol RObotic NETwork) sites for the years 2003–2009. We show that MODIS (MODerate-resolution Imaging Spectroradiometer) AOT exhibits significant biases due to wind speed and cloudiness of the observed scene, while MODIS AE, although overall unbiased, exhibits less spatial contrast on global scales than the AERONET observations. The same behaviour can be seen when MODIS AOT is compared against Maritime Aerosol Network (MAN) data, suggesting that the spatial coverage of our datasets does not preclude global conclusions. Thus, we develop empirical correction formulae for MODIS AOT and AE that significantly improve agreement of MODIS and AERONET observations. We show these correction formulae to be robust. Finally, we study random errors in the corrected MODIS AOT and AE and show that they mainly depend on AOT itself, although small contributions are present due to wind speed and cloud fraction in AOT random errors and due to AE and cloud fraction in AE random errors. Our analysis yields significantly higher random AOT errors than the official MODIS error estimate (0.03 + 0.05 τ), while random AE errors are smaller than might be expected. This new dataset of bias-corrected MODIS AOT and AE over ocean is intended for aerosol model validation and assimilation studies, but also has consequences as a stand-alone observational product. For instance, the corrected dataset suggests that much less fine mode aerosol is transported across the Pacific and Atlantic oceans.

Description: A neural network algorithm for cloud fraction estimation using NASA-Aura OMI VIS radiance measurements Atmospheric Measurement Techniques, 6, 2301-2309, 2013 Author(s): G. Saponaro, P. Kolmonen, J. Karhunen, J. Tamminen, and G. de Leeuw The discrimination of cloudy from cloud-free pixels is required in almost any estimate of a parameter retrieved from satellite data in the ultraviolet (UV), visible (VIS) or infrared (IR) parts of the electromagnetic spectrum. In this paper we report on the development of a neural network (NN) algorithm to estimate cloud fractions using radiances measured at the top of the atmosphere with the NASA-Aura Ozone Monitoring Instrument (OMI). We present and discuss the results obtained from the application of two different types of neural networks, i.e., extreme learning machine (ELM) and back propagation (BP). The NNs were trained with an OMI data sets existing of six orbits, tested with three other orbits and validated with another two orbits. The results were evaluated by comparison with cloud fractions available from the MODerate Resolution Imaging Spectrometer (MODIS) flying on Aqua in the same constellation as OMI, i.e., with minimal time difference between the OMI and MODIS observations. The results from the ELM and BP NNs are compared. They both deliver cloud fraction estimates in a fast and automated way, and they both performs generally well in the validation. However, over highly reflective surfaces, such as desert, or in the presence of dust layers in the atmosphere, the cloud fractions are not well predicted by the neural network. Over ocean the two NNs work equally well, but over land ELM performs better.

Description: Validation of stratospheric and mesospheric ozone observed by SMILES from International Space Station Atmospheric Measurement Techniques, 6, 2311-2338, 2013 Author(s): Y. Kasai, H. Sagawa, D. Kreyling, E. Dupuy, P. Baron, J. Mendrok, K. Suzuki, T. O. Sato, T. Nishibori, S. Mizobuchi, K. Kikuchi, T. Manabe, H. Ozeki, T. Sugita, M. Fujiwara, Y. Irimajiri, K. A. Walker, P. F. Bernath, C. Boone, G. Stiller, T. von Clarmann, J. Orphal, J. Urban, D. Murtagh, E. J. Llewellyn, D. Degenstein, A. E. Bourassa, N. D. Lloyd, L. Froidevaux, M. Birk, G. Wagner, F. Schreier, J. Xu, P. Vogt, T. Trautmann, and M. Yasui We observed ozone (O 3 ) in the vertical region between 250 and 0.0005 hPa (~ 12–96 km) using the Superconducting Submillimeter-Wave Limb-Emission Sounder (SMILES) on the Japanese Experiment Module (JEM) of the International Space Station (ISS) between 12 October 2009 and 21 April 2010. The new 4 K superconducting heterodyne receiver technology of SMILES allowed us to obtain a one order of magnitude better signal-to-noise ratio for the O 3 line observation compared to past spaceborne microwave instruments. The non-sun-synchronous orbit of the ISS allowed us to observe O 3 at various local times. We assessed the quality of the vertical profiles of O 3 in the 100–0.001 hPa (~ 16–90 km) region for the SMILES NICT Level 2 product version 2.1.5. The evaluation is based on four components: error analysis; internal comparisons of observations targeting three different instrumental setups for the same O 3 625.371 GHz transition; internal comparisons of two different retrieval algorithms; and external comparisons for various local times with ozonesonde, satellite and balloon observations (ENVISAT/MIPAS, SCISAT/ACE-FTS, Odin/OSIRIS, Odin/SMR, Aura/MLS, TELIS). SMILES O 3 data have an estimated absolute accuracy of better than 0.3 ppmv (3%) with a vertical resolution of 3–4 km over the 60 to 8 hPa range. The random error for a single measurement is better than the estimated systematic error, being less than 1, 2, and 7%, in the 40–1, 80–0.1, and 100–0.004 hPa pressure regions, respectively. SMILES O 3 abundance was 10–20% lower than all other satellite measurements at 8–0.1 hPa due to an error arising from uncertainties of the tangent point information and the gain calibration for the intensity of the spectrum. SMILES O 3 from observation frequency Band-B had better accuracy than that from Band-A. A two month period is required to accumulate measurements covering 24 h in local time of O 3 profile. However such a dataset can also contain variation due to dynamical, seasonal, and latitudinal effects.

Description: The Atmospheric radiation measurement (ARM) program network of microwave radiometers: instrumentation, data, and retrievals Atmospheric Measurement Techniques, 6, 2359-2372, 2013 Author(s): M. P. Cadeddu, J. C. Liljegren, and D. D. Turner The Climate Research Facility of the US Department of Energy's Atmospheric Radiation Measurement (ARM) Program operates a network of ground-based microwave radiometers. Data and retrievals from these instruments have been available to the scientific community for almost 20 yr. In the past five years the network has expanded to include a total of 22 microwave radiometers deployed in various locations around the world. The new instruments cover a frequency range between 22 and 197 GHz and are consistently and automatically calibrated. The latest addition to the network is a new generation of three-channel radiometers, currently in the early stage of deployment at all ARM sites. The network has been specifically designed to achieve increased accuracy in the retrieval of precipitable water vapor (PWV) and cloud liquid water path (LWP) with the long-term goal of providing the scientific community with reliable, calibrated radiometric data and retrievals of important geophysical quantities with well-characterized uncertainties. The radiometers provide high-quality, continuous datasets that can be utilized in a wealth of applications and scientific studies. This paper presents an overview of the microwave instrumentation, calibration procedures, data, and retrievals that are available for download from the ARM data archive.

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: Airborne sun photometer PLASMA: concept, measurements, comparison of aerosol extinction vertical profile with lidar Atmospheric Measurement Techniques, 6, 2383-2389, 2013 Author(s): Y. Karol, D. Tanré, P. Goloub, C. Vervaerde, J. Y. Balois, L. Blarel, T. Podvin, A. Mortier, and A. Chaikovsky A 15-channel airborne sun-tracking photometer has been developed. The instrument provides aerosol optical depths over a wide spectral range (0.34–2.25 μm) with an accuracy (ΔAOD) of approximately 0.01. Taking measurements at different altitudes allow us to derive the aerosol extinction vertical profile. Thanks to the wide spectral range of the instrument, information on the aerosol size distribution along the vertical is also available.

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: McClear: a new model estimating downwelling solar radiation at ground level in clear-sky conditions Atmospheric Measurement Techniques, 6, 2403-2418, 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. It is a fully physical model replacing empirical relations or simpler models used before. It exploits the recent results on aerosol properties, and total column content in water vapour 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 table, 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 at several stations within the Baseline Surface Radiation Network in various climates. The bias for global irradiance comprises between −6 and 25 W m −2 . The RMSE ranges from 20 W m −2 (3% of the mean observed irradiance) to 36 W m −2 (5%) and the correlation coefficient ranges between 0.95 and 0.99. The bias for the direct irradiance comprises between −48 and +33 W m −2 . The root mean square error (RMSE) ranges from 33 W m −2 (5%) to 64 W m −2 (10%). The correlation coefficient ranges between 0.84 and 0.98. This work demonstrates the quality of the McClear model combined with MACC products, and indirectly the quality of the aerosol properties modelled by the MACC reanalysis.

Description: Ten years of MIPAS measurements with ESA Level 2 processor V6 – Part 1: Retrieval algorithm and diagnostics of the products Atmospheric Measurement Techniques, 6, 2419-2439, 2013 Author(s): P. Raspollini, B. Carli, M. Carlotti, S. Ceccherini, A. Dehn, B. M. Dinelli, A. Dudhia, J.-M. Flaud, M. López-Puertas, F. Niro, J. J. Remedios, M. Ridolfi, H. Sembhi, L. Sgheri, and T. von Clarmann The MIPAS (Michelson Interferometer for Passive Atmospheric Sounding) instrument on the Envisat (Environmental satellite) satellite has provided vertical profiles of the atmospheric composition on a global scale for almost ten years. The MIPAS mission is divided in two phases: the full resolution phase, from 2002 to 2004, and the optimized resolution phase, from 2005 to 2012, which is characterized by a finer vertical and horizontal sampling attained through a reduction of the spectral resolution. While the description and characterization of the products of the ESA processor for the full resolution phase has been already described in previous papers, in this paper we focus on the performances of the latest version of the ESA (European Space Agency) processor, named ML2PP V6 (MIPAS Level 2 Prototype Processor), which has been used for reprocessing the entire mission. The ESA processor had to perform the operational near real time analysis of the observations and its products needed to be available for data assimilation. Therefore, it has been designed for fast, continuous and automated analysis of observations made in quite different atmospheric conditions and for a minimum use of external constraints in order to avoid biases in the products. The dense vertical sampling of the measurements adopted in the second phase of the MIPAS mission resulted in sampling intervals finer than the instantaneous field of view of the instrument. Together with the choice of a retrieval grid aligned with the vertical sampling of the measurements, this made ill-conditioned the retrieval problem of the MIPAS operational processor. This problem has been handled with minimal changes to the original retrieval approach but with significant improvements nonetheless. The Levenberg–Marquardt method, already present in the retrieval scheme for its capability to provide fast convergence for nonlinear problems, is now also exploited for the reduction of the ill-conditioning of the inversion. An expression specifically designed for the regularizing Levenberg–Marquardt method has been implemented for the computation of the covariance matrices and averaging kernels of the retrieved products. The regularization of the Levenberg–Marquardt method is controlled by the convergence criteria and is deliberately kept weak. The resulting oscillations of the retrieved profile are a posteriori damped by an innovative self-adapting Tikhonov regularization. The convergence criteria and the weakness of the self-adapting regularization ensure that minimum constraints are used and the best vertical resolution obtainable from the measurements is achieved in all atmospheric conditions. Random and systematic errors, as well as vertical and horizontal resolution are compared in the two phases of the mission for all products, namely: temperature, H 2O , O 3 , HNO 3 , CH 4 , N 2O , NO 2 , CFC-11, CFC-12, N 2 O 5 and ClONO 2 . The use in the two phases of the mission of different optimized sets of spectral intervals ensures that, despite the different spectral resolutions, comparable performances are obtained in the whole MIPAS mission in terms of random and systematic errors, while the vertical resolution and the horizontal resolution are significantly better in the case of the optimized resolution measurements.

Description: Ground-based stratospheric O 3 and HNO 3 measurements at Thule, Greenland: an intercomparison with Aura MLS observations Atmospheric Measurement Techniques, 6, 2441-2453, 2013 Author(s): I. Fiorucci, G. Muscari, L. Froidevaux, and M. L. Santee In response to the need for improving our understanding of the evolution and the interannual variability of the winter Arctic stratosphere, in January 2009 a Ground-Based Millimeter-wave Spectrometer (GBMS) was installed at the Network for the Detection of Atmospheric Composition Change (NDACC) site in Thule (76.5° N, 68.8° W), Greenland. In this work, stratospheric GBMS O 3 and HNO 3 vertical profiles obtained from Thule during the winters 2010 (HNO 3 only), 2011 and 2012 are characterized and intercompared with co-located measurements of the Aura Microwave Limb Sounder (MLS) experiment. Using a recently developed algorithm based on Optimal Estimation, we find that the GBMS O 3 retrievals show good sensitivity (〉 80%) to atmospheric variations between ~ 17 and ~ 50 km, where their 1σ uncertainty is estimated to be the larger of ~ 11% or 0.2 ppmv. Similarly, HNO 3 profiles can be considered for scientific use between ~ 17 and ~ 45 km altitude, with a 1σ uncertainty that amounts to the larger of 15% or 0.2 ppbv. Comparisons with Aura MLS version 3.3 observations show that, on average, GBMS O 3 mixing ratios are biased negatively with respect to MLS throughout the stratosphere, with differences ranging between ~ 0.3 ppmv (8%) and 0.9 ppmv (18%) in the 17–50 km vertical range. GBMS HNO 3 values display instead a positive bias with respect to MLS up to 26 km, reaching a maximum of ~ 1 ppbv (10%) near the mixing ratio profile peak. O 3 and HNO 3 values from the two datasets prove to be well correlated at all altitudes, although their correlations worsen at the lower end of the altitude ranges considered. Column contents of GBMS and MLS O 3 (from 20 km upwards) and HNO 3 (from 17 km upwards) correlate very well and indicate that GBMS measurements can provide valuable estimates of column interannual and seasonal variations for these compounds.

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: Inter-calibration of polar imager solar channels using SEVIRI Atmospheric Measurement Techniques, 6, 2495-2508, 2013 Author(s): J. F. Meirink, R. A. Roebeling, and P. Stammes Accurate calibration of satellite imagers is a prerequisite for using their measurements in climate applications. Here we present a method for the inter-calibration of geostationary and polar-orbiting imager solar channels based on regressions of collocated near-nadir reflectances. Specific attention is paid to correcting for differences in spectral response between instruments. The method is used to calibrate the solar channels of the Spinning Enhanced Visible and Infrared Imager (SEVIRI) on the geostationary Meteosat satellite with corresponding channels of the Moderate Resolution Imaging Spectroradiometer (MODIS) on the polar-orbiting Aqua satellite. The SEVIRI operational calibration is found to be stable during the years 2004 to 2009, but offset by −8, −6, and +3.5 % for channels 1 (0.6 μm), 2 (0.8 μm), and 3 (1.6 μm), respectively. These results are robust for a range of choices that can be made regarding data collocation and selection, as long as the viewing and illumination geometries of the two instruments are matched. Uncertainties in the inter-calibration method are estimated to be 1 % for channel 1 and 1.5 % for channels 2 and 3. A specific application of our method is the inter-calibration of polar imagers using SEVIRI as a transfer instrument. This offers an alternative to direct inter-calibration, which in general has to rely on high-latitude collocations. Using this method we have tied MODIS-Terra and Advanced Very High Resolution Radiometer (AVHRR) instruments on National Oceanic and Atmospheric Administration (NOAA) satellites 17 and 18 to MODIS-Aqua for the years 2007 to 2009. While reflectances of the two MODIS instruments differ less than 2 % for all channels considered, deviations of an existing AVHRR calibration from MODIS-Aqua reach −3.5 and +2.5 % for the 0.8 and 1.6 μm channels, respectively.

Description: Using ocean-glint scattered sunlight as a diagnostic tool for satellite remote sensing of greenhouse gases Atmospheric Measurement Techniques, 6, 2509-2520, 2013 Author(s): A. Butz, S. Guerlet, O. P. Hasekamp, A. Kuze, and H. Suto Spectroscopic measurements of sunlight backscattered by the Earth's surface is a technique widely used for remote sensing of atmospheric constituent concentrations from space. Thereby, remote sensing of greenhouse gases poses particularly challenging accuracy requirements for instrumentation and retrieval algorithms which, in general, suffer from various error sources. Here, we investigate a method that helps disentangle sources of error for observations of sunlight backscattered from the glint spot on the ocean surface. The method exploits the backscattering characteristics of the ocean surface, which is bright for glint geometry but dark for off-glint angles. This property allows for identifying a set of clean scenes where light scattering due to particles in the atmosphere is negligible such that uncertain knowledge of the lightpath can be excluded as a source of error. We apply the method to more than 3 yr of ocean-glint measurements by the Thermal And Near infrared Sensor for carbon Observation (TANSO) Fourier Transform Spectrometer (FTS) onboard the Greenhouse Gases Observing Satellite (GOSAT), which aims at measuring carbon dioxide (CO 2 ) and methane (CH 4 ) concentrations. The proposed method is able to clearly monitor recent improvements in the instrument calibration of the oxygen (O 2 ) A-band channel and suggests some residual uncertainty in our knowledge about the instrument. We further assess the consistency of CO 2 retrievals from several absorption bands between 6400 cm −1 (1565 nm) and 4800 cm −1 (2100 nm) and find that the absorption bands commonly used for monitoring of CO 2 dry air mole fractions from GOSAT allow for consistency better than 1.5 ppm. Usage of other bands reveals significant inconsistency among retrieved CO 2 concentrations pointing at inconsistency of spectroscopic parameters.

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: Parameterizing radiative transfer to convert MAX-DOAS dSCDs into near-surface box-averaged mixing ratios Atmospheric Measurement Techniques, 6, 1521-1532, 2013 Author(s): R. Sinreich, A. Merten, L. Molina, and R. Volkamer We present a novel parameterization method to convert multi-axis differential optical absorption spectroscopy (MAX-DOAS) differential slant column densities (dSCDs) into near-surface box-averaged volume mixing ratios. The approach is applicable inside the planetary boundary layer under conditions with significant aerosol load, and builds on the increased sensitivity of MAX-DOAS near the instrument altitude. It parameterizes radiative transfer model calculations and significantly reduces the computational effort, while retrieving ~ 1 degree of freedom. The biggest benefit of this method is that the retrieval of an aerosol profile, which usually is necessary for deriving a trace gas concentration from MAX-DOAS dSCDs, is not needed. The method is applied to NO 2 MAX-DOAS dSCDs recorded during the Mexico City Metropolitan Area 2006 (MCMA-2006) measurement campaign. The retrieved volume mixing ratios of two elevation angles (1° and 3°) are compared to volume mixing ratios measured by two long-path (LP)-DOAS instruments located at the same site. Measurements are found to agree well during times when vertical mixing is expected to be strong. However, inhomogeneities in the air mass above Mexico City can be detected by exploiting the different horizontal and vertical dimensions probed by the MAX-DOAS and LP-DOAS instruments. In particular, a vertical gradient in NO 2 close to the ground can be observed in the afternoon, and is attributed to reduced mixing coupled with near-surface emission inside street canyons. The existence of a vertical gradient in the lower 250 m during parts of the day shows the general challenge of sampling the boundary layer in a representative way, and emphasizes the need of vertically resolved measurements.

Description: Polarization data from SCIAMACHY limb backscatter observations compared to vector radiative transfer model simulations Atmospheric Measurement Techniques, 6, 1503-1520, 2013 Author(s): P. Liebing, K. Bramstedt, S. Noël, V. Rozanov, H. Bovensmann, and J. P. Burrows SCIAMACHY is a passive imaging spectrometer onboard ENVISAT designed to obtain trace gas abundances from measured radiances and irradiances in the UV to SWIR range in nadir-, limb- and occultation-viewing modes. Its grating spectrometer introduces a substantial sensitivity to the polarization of the incoming light with nonnegligible effects on the radiometric calibration. To be able to correct for the polarization sensitivity, SCIAMACHY utilizes broadband Polarization Measurement Devices (PMDs). While for the nadir-viewing mode the measured atmospheric polarization has been validated against POLDER data (Tilstra and Stammes, 2007, 2010), a similar validation study regarding the limb-viewing mode has not yet been performed. This paper aims at an assessment of the quality of the SCIAMACHY limb polarization data. Since limb polarization measurements by other air/spaceborne instruments in the spectral range of SCIAMACHY are not available, a comparison with radiative transfer simulations by SCIATRAN V3.1 (Rozanov et al., 2013) using a wide range of atmospheric parameters is performed. SCIATRAN is a vector radiative transfer model (VRTM) capable of performing calculations of the multiply scattered radiance in a spherically symmetric atmosphere. The study shows that the limb polarization data exhibit a large time-dependent bias that decreases with wavelength. Possible reasons for this bias are a still unknown combination of insufficient accuracy or inconsistencies of the on-ground calibration data, scan mirror degradation and stress induced changes of the polarization response of components inside the optical bench of the instrument. It is shown that it should in principle be feasible to recalibrate the effective polarization sensitivity of the instrument using the in-flight data and VRTM simulations.

Description: Evaluating calibration strategies for isotope ratio infrared spectroscopy for atmospheric 13 CO 2 / 12 CO 2 measurement Atmospheric Measurement Techniques, 6, 1491-1501, 2013 Author(s): X.-F. Wen, Y. Meng, X.-Y. Zhang, X.-M. Sun, and X. Lee Isotope ratio infrared spectroscopy (IRIS) provides an in situ technique for measuring δ 13 C in atmospheric CO 2 . A number of methods have been proposed for calibrating the IRIS measurements, but few studies have systematically evaluated their accuracy for atmospheric applications. In this study, we carried out laboratory and ambient measurements with two commercial IRIS analyzers and compared the accuracy of four calibration strategies. We found that calibration based on the 12 C and 13 C mixing ratios (Bowling et al., 2003) and on linear interpolation of the measured delta using the mixing ratio of the major isotopologue (Lee et al., 2005) yielded accuracy better than 0.06‰. Over a 7-day atmospheric measurement in Beijing, the two analyzers agreed to within −0.02 ± 0.18‰ after proper calibration. However, even after calibration the difference between the two analyzers showed a slight correlation with concentration, and this concentration dependence propagated through the Keeling analysis, resulting in a much larger difference of 2.44‰ for the Keeling intercept. The high sensitivity of the Keeling analysis to the concentration dependence underscores the challenge of IRIS for atmospheric research.

Description: Cavity ring-down spectroscopy sensor development for high-time-resolution measurements of gaseous elemental mercury in ambient air Atmospheric Measurement Techniques, 6, 1477-1489, 2013 Author(s): A. Pierce, D. Obrist, H. Moosmüller, X. Faïn, and C. Moore We describe further development of a previous laboratory prototype pulsed cavity ring-down spectroscopy (CRDS) sensor into a field-deployable system for high-time-resolution, continuous, and automated measurement of gaseous elemental mercury (GEM) concentrations in ambient air. We employed an external, isotopically enriched Hg cell for automated locking and stabilization of the laser wavelength on the GEM peak absorption during measurements. Further, we describe implementation of differential absorption measurements via a piezoelectric tuning element for pulse-by-pulse tuning of the laser wavelength onto and off of the GEM absorption line. This allowed us to continuously correct (at 25 Hz) for system baseline extinction losses unrelated to GEM absorption. Extensive measurement and calibration data obtained with the system were based on spike addition in both GEM-free air and ambient air. Challenges and interferences that occurred during measurements (particularly in ambient air) are discussed including temperature and ozone (O 3 ) concentration fluctuations, and steps taken to reduce these. CRDS data were highly linear ( r 2 ≥ 0.98) with data from a commercial Tekran 2537 Hg analyzer across a wide range of GEM concentrations (0 to 127 ng m −3 ) in Hg-free and ambient air. Measurements during periods of stable background GEM concentrations provided a conservative instrument sensitivity estimate of 0.35 ng m −3 for the CRDS system when time averaged for 5 min. This sensitivity, along with concentration patterns observed in ambient air (with the CRDS system and verified with the Tekran analyzer), showed that the sensor was capable of characterizing GEM fluctuations in ambient air. The value of fast-response GEM measurements was shown by a series of GEM spike additions – highlighting that high-temporal-resolution measurement allowed for detailed characterization of fast concentration fluctuations not possible with traditional analyzers.

Description: Broadband measurements of aerosol extinction in the ultraviolet spectral region Atmospheric Measurement Techniques, 6, 861-877, 2013 Author(s): R. A. Washenfelder, J. M. Flores, C. A. Brock, S. S. Brown, and Y. Rudich Aerosols influence the Earth's radiative budget by scattering and absorbing incoming solar radiation. The optical properties of aerosols vary as a function of wavelength, but few measurements have reported the wavelength dependence of aerosol extinction cross sections and complex refractive indices. We describe a new laboratory instrument to measure aerosol optical extinction as a function of wavelength, using cavity enhanced spectroscopy with a broadband light source. The instrument consists of two broadband channels which span the 360–390 and 385–420 nm spectral regions using two light emitting diodes (LED) and a grating spectrometer with charge-coupled device (CCD) detector. We determined aerosol extinction cross sections and directly observed Mie scattering resonances for aerosols that are purely scattering (polystyrene latex spheres and ammonium sulfate), slightly absorbing (Suwannee River fulvic acid), and strongly absorbing (nigrosin dye). We describe an approach for retrieving refractive indices as a function of wavelength from the measured extinction cross sections over the 360–420 nm wavelength region. The retrieved refractive indices for PSL and ammonium sulfate agree within uncertainty with the literature values for this spectral region. The refractive index determined for nigrosin is 1.78 (± 0.03) + 0.19 (± 0.08) i at 360 nm and 1.63 (± 0.03) + 0.21 (± 0.05) i at 420 nm. The refractive index determined for Suwannee River fulvic acid is 1.71 (± 0.02) + 0.07 (± 0.06) i at 360 nm and 1.66 (± 0.02) + 0.06 (± 0.04) i at 420 nm. These laboratory results support the potential for a field instrument capable of determining ambient aerosol optical extinction, average aerosol extinction cross section, and complex refractive index as a function of wavelength.

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: Climatologies from satellite measurements: the impact of orbital sampling on the standard error of the mean Atmospheric Measurement Techniques, 6, 937-948, 2013 Author(s): M. Toohey and T. von Clarmann Climatologies of atmospheric observations are often produced by binning measurements according to latitude and calculating zonal means. The uncertainty in these climatological means is characterised by the standard error of the mean (SEM). However, the usual estimator of the SEM, i.e., the sample standard deviation divided by the square root of the sample size, holds only for uncorrelated randomly sampled measurements. Measurements of the atmospheric state along a satellite orbit cannot always be considered as independent because (a) the time-space interval between two nearest observations is often smaller than the typical scale of variations in the atmospheric state, and (b) the regular time-space sampling pattern of a satellite instrument strongly deviates from random sampling. We have developed a numerical experiment where global chemical fields from a chemistry climate model are sampled according to real sampling patterns of satellite-borne instruments. As case studies, the model fields are sampled using sampling patterns of the Michelson Interferometer for Passive Atmospheric Sounding (MIPAS) and Atmospheric Chemistry Experiment Fourier-Transform Spectrometer (ACE-FTS) satellite instruments. Through an iterative subsampling technique, and by incorporating information on the random errors of the MIPAS and ACE-FTS measurements, we produce empirical estimates of the standard error of monthly mean zonal mean model O 3 in 5° latitude bins. We find that generally the classic SEM estimator is a conservative estimate of the SEM, i.e., the empirical SEM is often less than or approximately equal to the classic estimate. Exceptions occur only when natural variability is larger than the random measurement error, and specifically in instances where the zonal sampling distribution shows non-uniformity with a similar zonal structure as variations in the sampled field, leading to maximum sensitivity to arbitrary phase shifts between the sample distribution and sampled field. The occurrence of such instances is thus very sensitive to slight changes in the sampling distribution, and to the variations in the measured field. This study highlights the need for caution in the interpretation of the oft-used classically computed SEM, and outlines a relatively simple methodology that can be used to assess one component of the uncertainty in monthly mean zonal mean climatologies produced from measurements from satellite-borne instruments.

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: 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: A chemical analyzer for charged ultrafine particles Atmospheric Measurement Techniques, 6, 2339-2348, 2013 Author(s): S. G. Gonser and A. Held New particle formation is a frequent phenomenon in the atmosphere and of major significance for the Earth's climate and human health. To date the mechanisms leading to the nucleation of particles as well as to aerosol growth are not completely understood. A lack of appropriate measurement equipment for online analysis of the chemical composition of freshly nucleated particles is one major limitation. We have developed a Chemical Analyzer for Charged Ultrafine Particles (CAChUP) capable of analyzing particles with diameters below 30 nm. A bulk of size-separated particles is collected electrostatically on a metal filament, resistively desorbed and subsequently analyzed for its molecular composition in a time of flight mass spectrometer. We report on technical details as well as characterization experiments performed with the CAChUP. Our instrument was tested in the laboratory for its detection performance as well as for its collection and desorption capabilities. The manual application of defined masses of camphene (C 10 H 16 ) to the desorption filament resulted in a detection limit between 0.5 and 5 ng, and showed a linear response of the mass spectrometer. Flow tube experiments of 25 nm diameter secondary organic aerosol from ozonolysis of alpha-pinene also showed a linear relation between collection time and the mass spectrometer's signal intensity. The resulting mass spectra from the collection experiments are in good agreement with published work on particles generated by the ozonolysis of alpha-pinene. A sensitivity study shows that the current setup of CAChUP is ready for laboratory measurements and for the observation of new particle formation events in the field.

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: Comparison of AOD between CALIPSO and MODIS: significant differences over major dust and biomass burning regions Atmospheric Measurement Techniques, 6, 2391-2401, 2013 Author(s): X. Ma, K. Bartlett, K. Harmon, and F. Yu Cloud–Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO) provide global vertical profiles of aerosol optical properties for the first time. In this study, we employed about 6 yr (2006–2011) of CALIPSO level 3 monthly mean gridded aerosol optical depth (AOD) products (daytime and nighttime) for cloud-free conditions, to compare with the Moderate Resolution Imaging Spectroradiometer (MODIS) Terra/Aqua level 3 monthly mean AOD dataset for the same time period. While the spatial distribution and seasonal variability of CALIPSO AOD is generally consistent with that of MODIS, CALIPSO is overall lower than MODIS as MODIS has higher frequency than CALIPSO for most bins of AOD. The correlation between MODIS and CALIPSO is better over ocean than over land. We focused on four regions that have large systematic differences: two over dust regions (the Sahara and Northwest China) and two over biomass burning regions (South Africa and South America). It is found that CALIPSO AOD is significantly lower than MODIS AOD over dust regions during the whole time period, with a maximum difference of 0.3 over the Saharan region and 0.25 over Northwest China. For biomass burning regions, CALIPSO AOD is significantly higher than MODIS AOD over South Africa, with a maximum difference of 0.25. Additionally CALIPSO AOD is slightly higher than MODIS AOD over South America for most of the time period, with a few exceptions in 2006, 2007, and 2010, when biomass burning is significantly stronger than during other years. We analyzed the impact of the satellite spatial and temporal sampling issue by using level 2 CALIPSO and MODIS products, and these systematic differences can still be found. The results of this study indicate that systematic differences of CALIPSO relative to MODIS are closely associated with aerosol types, which vary by location and season. Large differences over dust and biomass burning regions may suggest that assumptions made in satellite retrievals, such as the assumed lidar ratios for CALIPSO retrievals over dust and biomass burning regions or the surface reflectance information and/or the aerosol model utilized by the MODIS algorithm, are not appropriate.

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: Extreme extension across Seram and Ambon, eastern Indonesia: evidence for Banda slab rollback Solid Earth, 4, 277-314, 2013 Author(s): J. M. Pownall, R. Hall, and I. M. Watkinson The island of Seram, which lies in the northern part of the 180°-curved Banda Arc, has previously been interpreted as a fold-and-thrust belt formed during arc-continent collision, which incorporates ophiolites intruded by granites thought to have been produced by anatexis within a metamorphic sole. However, new geological mapping and a re-examination of the field relations cause us to question this model. We instead propose that there is evidence for recent and rapid N–S extension that has caused the high-temperature exhumation of lherzolites beneath low-angle lithospheric detachment faults that induced high-temperature metamorphism and melting in overlying crustal rocks. These "Kobipoto Complex" migmatites include highly residual Al–Mg-rich garnet + cordierite + sillimanite + spinel + corundum granulites (exposed in the Kobipoto Mountains) which contain coexisting spinel + quartz, indicating that peak metamorphic temperatures likely approached 900 °C. Associated with these residual granulites are voluminous Mio-Pliocene granitic diatexites, or "cordierite granites", which crop out on Ambon, western Seram, and in the Kobipoto Mountains and incorporate abundant schlieren of spinel- and sillimanite-bearing residuum. Quaternary "ambonites" (cordierite + garnet dacites) emplaced on Ambon were also evidently sourced from the Kobipoto Complex migmatites as demonstrated by granulite-inherited xenoliths. Exhumation of the hot peridotites and granulite-facies Kobipoto Complex migmatites to shallower structural levels caused greenschist- to lower-amphibolite facies metapelites and amphibolites of the Tehoru Formation to be overprinted by sillimanite-grade metamorphism, migmatisation, and limited localised anatexis to form the Taunusa Complex. The extreme extension required to have driven Kobipoto Complex exhumation evidently occurred throughout Seram and along much of the northern Banda Arc. The lherzolites must have been juxtaposed against the crust at typical lithospheric mantle temperatures in order to account for such high-temperature metamorphism and therefore could not have been part of a cooled ophiolite. In central Seram, lenses of peridotites are incorporated with a major left-lateral strike-slip shear zone (the "Kawa Shear Zone"), demonstrating that strike-slip motions likely initiated shortly after the mantle had been partly exhumed by detachment faulting and that the main strike-slip faults may themselves be reactivated and steepened low-angle detachments. The geodynamic driver for mantle exhumation along the detachment faults and strike-slip faulting in central Seram is very likely the same; we interpret the extreme extension to be the result of eastward slab rollback into the Banda Embayment as outlined by the latest plate reconstructions for Banda Arc evolution.

Description: Microwave radiometer to retrieve temperature profiles from the surface to the stratopause Atmospheric Measurement Techniques, 6, 2477-2494, 2013 Author(s): O. Stähli, A. Murk, N. Kämpfer, C. Mätzler, and P. Eriksson TEMPERA (TEMPERature RAdiometer) is a new ground-based radiometer which measures in a frequency range from 51–57 GHz radiation emitted by the atmosphere. With this instrument it is possible to measure temperature profiles from ground to about 50 km. This is the first ground-based instrument with the capability to retrieve temperature profiles simultaneously for the troposphere and stratosphere. The measurement is done with a filterbank in combination with a digital fast Fourier transform spectrometer. A hot load and a noise diode are used as stable calibration sources. The optics consist of an off-axis parabolic mirror to collect the sky radiation. Due to the Zeeman effect on the emission lines used, the maximum height for the temperature retrieval is about 50 km. The effect is apparent in the measured spectra. The performance of TEMPERA is validated by comparison with nearby radiosonde and satellite data from the Microwave Limb Sounder on the Aura satellite. In this paper we present the design and measurement method of the instrument followed by a description of the retrieval method, together with a validation of TEMPERA data over its first year, 2012.

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 of nitric oxide in the mesosphere and lower thermosphere from SCIAMACHY limb spectra Atmospheric Measurement Techniques, 6, 2521-2531, 2013 Author(s): S. Bender, M. Sinnhuber, J. P. Burrows, M. Langowski, B. Funke, and M. López-Puertas We use the ultra-violet (UV) spectra in the range 230–300 nm from the SCanning Imaging Absorption spectroMeter for Atmospheric CHartographY (SCIAMACHY) to retrieve the nitric oxide (NO) number densities from atmospheric emissions in the gamma-bands in the mesosphere and lower thermosphere. Using 3-D ray tracing, a 2-D retrieval grid, and regularisation with respect to altitude and latitude, we retrieve a whole semi-orbit simultaneously for the altitude range from 60 to 160 km. We present details of the retrieval algorithm, first results, and initial comparisons to data from the Michelson Interferometer for Passive Atmospheric Sounding (MIPAS). Our results agree on average well with MIPAS data and are in line with previously published measurements from other instruments. For the time of available measurements in 2008–2011, we achieve a vertical resolution of 5–10 km in the altitude range 70–140 km and a horizontal resolution of about 9° from 60° S–60° N. With this we have independent measurements of the NO densities in the mesosphere and lower thermosphere with approximately global coverage. This data can be further used to validate climate models or as input for them.

Description: Development and field testing of a rapid and ultra-stable atmospheric carbon dioxide spectrometer Atmospheric Measurement Techniques, 7, 4445-4453, 2014 Author(s): B. Xiang, D. D. Nelson, J. B. McManus, M. S. Zahniser, R. A. Wehr, and S. C. Wofsy We present field test results for a new spectroscopic instrument to measure atmospheric carbon dioxide (CO 2 ) with high precision (0.02 μmol mol −1 , or ppm at 1 Hz) and demonstrate high stability (within 0.1 ppm over more than 8 months), without the need for hourly, daily, or even monthly calibration against high-pressure gas cylinders. The technical novelty of this instrument (ABsolute Carbon dioxide, ABC) is the spectral null method using an internal quartz reference cell with known CO 2 column density. Compared to a previously described prototype, the field instrument has better stability and benefits from more precise thermal control of the optics and more accurate pressure measurements in the sample cell (at the mTorr level). The instrument has been deployed at a long-term ecological research site (the Harvard Forest, USA), where it has measured for 8 months without on-site calibration and with minimal maintenance, showing drift bounds of less than 0.1 ppm. Field measurements agree well with those of a commercially available cavity ring-down CO 2 instrument (Picarro G2301) run with a standard calibration protocol. This field test demonstrates that ABC is capable of performing high-accuracy, unattended, continuous field measurements with minimal use of reference gas cylinders.

Description: Investigating uptake of N 2 O in agricultural soils using a high-precision dynamic chamber method Atmospheric Measurement Techniques, 7, 4455-4462, 2014 Author(s): N. J. Cowan, D. Famulari, P. E. Levy, M. Anderson, D. S. Reay, and U. M. Skiba Uptake (or negative flux) of nitrous oxide (N 2 O) in agricultural soils is a controversial issue which has proved difficult to investigate in the past due to constraints such as instrumental precision and methodological uncertainties. Using a recently developed high-precision quantum cascade laser gas analyser combined with a closed dynamic chamber, a well-defined detection limit of 4 μg N 2 O-N m −2 h −1 could be achieved for individual soil flux measurements. 1220 measurements of N 2 O flux were made from a variety of UK soils using this method, of which 115 indicated uptake by the soil (i.e. a negative flux in the micrometeorological sign convention). Only four of these apparently negative fluxes were greater than the detection limit of the method, which suggests that the vast majority of reported negative fluxes from such measurements are actually due to instrument noise. As such, we suggest that the bulk of negative N 2 O fluxes reported for agricultural fields are most likely due to limits in detection of a particular flux measurement methodology and not a result of microbiological activity consuming atmospheric N 2 O.

Description: Rapid, optical measurement of the atmospheric pressure on a fast research aircraft using open-path TDLAS Atmospheric Measurement Techniques, 7, 3653-3666, 2014 Author(s): B. Buchholz, A. Afchine, and V. Ebert Because of the high travel speed, the complex flow dynamics around an aircraft, and the complex dependency of the fluid dynamics on numerous airborne parameters, it is quite difficult to obtain accurate pressure values at a specific instrument location of an aircraft's fuselage. Complex simulations using computational fluid dynamics (CFD) models can in theory computationally "transfer" pressure values from one location to another. However, for long flight patterns, this process is inconvenient and cumbersome. Furthermore, these CFD transfer models require a local experimental validation, which is rarely available. In this paper, we describe an integrated approach for a spectroscopic, calibration-free, in-flight pressure determination in an open-path White cell on an aircraft fuselage using ambient, atmospheric water vapour as the "sensor species". The presented measurements are realised with the HAI (Hygrometer for Atmospheric Investigations) instrument, built for multiphase water detection via calibration-free TDLAS (tunable diode laser absorption spectroscopy). The pressure determination is based on raw data used for H 2 O concentration measurement, but with a different post-flight evaluation method, and can therefore be conducted at deferred time intervals on any desired flight track. The spectroscopic pressure is compared in-flight with the static ambient pressure of the aircraft avionic system and a micro-mechanical pressure sensor, located next to the open-path cell, over a pressure range from 150 to 800 hPa, and a water vapour concentration range of more than 3 orders of magnitude. The correlation between the micro-mechanical pressure sensor measurements and the spectroscopic pressure measurements shows an average deviation from linearity of only 0.14% and a small offset of 9.5 hPa. For the spectroscopic pressure evaluation we derive measurement uncertainties under laboratory conditions of 3.2 and 5.1% during in-flight operation on the HALO airplane. Under certain flight conditions we quantified, for the first time, stalling-induced, dynamic pressure deviations of up to 30% (at 200 hPa) between the avionic sensor and the optical and mechanical pressure sensors integrated in HAI. Such severe local pressure deviations from the typically used avionic pressure are important to take into account for other airborne sensors employed on such fast flying platforms as the HALO aircraft.

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: Changes in soil organic carbon and nitrogen capacities of Salix cheilophila Schneid along a revegetation chronosequence in semi-arid degraded sandy land of the Gonghe Basin, Tibet Plateau Solid Earth, 5, 1045-1054, 2014 Author(s): Y. Yu and Z. Q. Jia The Gonghe Basin is a sandified and desertified region of China, but the distribution of soil organic carbon (SOC) and total nitrogen (TN) along the cultivation chronosequence across this ecologically fragile region is not well understood. This study was carried out to understand the effects of restoration with Salix cheilophila for different periods of time (6, 11, 16, 21 years) to test whether it enhanced C and N storage. Soil samples, in four replications from seven depth increments (0–10, 10–20, 20–30, 30–50, 50–100, 100–150 and 150–200 cm), were collected in each stand. Soil bulk density, SOC, TN, aboveground biomass and root biomass were measured. Results indicated that changes occurred in both the upper and deeper soil layers with an increase in revegetation time. The 0–200 cm soil showed that the 6-year stand gained 3.89 Mg C ha −1 and 1.00 Mg N ha −1 , which accounted for 40.82% of the original SOC and 11.06% of the TN of the 0-year stand. The 11-year stand gained 7.82 Mg C ha −1 and 1.98 Mg N ha −1 in the 0–200 cm soil layers, accounting for 58.06% of the SOC and 19.80% of the TN of the 0-year stand. The 16-year stand gained 11.32 Mg C ha −1 and 3.30 Mg N ha −1 in the 0–200 cm soil layers, accounting for 66.71% of the SOC and 21.98% of the TN of the 0-year stand. The 21-year stand gained 13.05 Mg C ha −1 and 5.45 Mg N ha −1 from the same soil depth, accounting for 69.79% of the SOC and 40.47% of the TN compared with the 0-year stand. The extent of these changes depended on soil depth and plantation age. The results demonstrated that, as stand age increased, the storage of SOC and TN increased. These results further indicated that restoration with S. cheilophila has positive impacts on the Gonghe Basin and has increased the capacity of SOC sequestration and N storage. The shrub's role as carbon sink is compatible with system management and persistence. The findings are significant for assessing C and N sequestration accurately in semi-arid degraded high, cold sandy regions in the future.

Description: Mixing-layer height retrieval with ceilometer and Doppler lidar: from case studies to long-term assessment Atmospheric Measurement Techniques, 7, 3685-3704, 2014 Author(s): J. H. Schween, A. Hirsikko, U. Löhnert, and S. Crewell Aerosol signatures observed by ceilometers are frequently used to derive mixing-layer height (MLH) which is an essential variable for air quality modelling. However, Doppler wind lidar measurements of vertical velocity can provide a more direct estimation of MLH via simple thresholding. A case study reveals difficulties in the aerosol-based MLH retrieval during transition times when the mixing layer builds up in the morning and when turbulence decays in the afternoon. The difficulties can be explained by the fact that the aerosol distribution is related to the history of the mixing process and aerosol characteristics are modified by humidification. The results of the case study are generalized by evaluating one year of joint measurements by a Vaisala CT25K and a HALO Photonics Streamline wind lidar. On average the aerosol-based retrieval gives higher MLH than the wind lidar with an overestimation of MLH by about 300 m (600 m) in the morning (late afternoon). Also, the daily aerosol-based maximum MLH is larger and occurs later during the day and the average morning growth rates are smaller than those derived from the vertical wind. In fair weather conditions classified by less than 4 octa cloud cover the mean diurnal cycle of cloud base height corresponds well to the mixing-layer height showing potential for a simplified MLH estimation.

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: 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: Preface: Environmental benefits of biochar Solid Earth, 5, 1301-1303, 2014 Author(s): J. Paz-Ferreiro, A. Méndez, A. M. Tarquis, A. Cerdà, and G. Gascó

Description: Radio occultation bending angle anomalies during tropical cyclones Atmospheric Measurement Techniques, 4, 1053-1060, 2011 Author(s): R. Biondi, T. Neubert, S. Syndergaard, and J. K. Nielsen The tropical deep convection affects the radiation balance of the atmosphere changing the water vapor mixing ratio and the temperature of the upper troposphere lower stratosphere. The aim of this work is to better understand these processes and to investigate if severe storms leave a significant signature in radio occultation profiles in the tropical tropopause layer. Using tropical cyclone best track database and data from different GPS radio occultation missions (COSMIC, GRACE, CHAMP, SACC and GPSMET), we selected 1194 profiles in a time window of 3 h and a space window of 300 km from the eye of the cyclone. We show that the bending angle anomaly of a GPS radio occultation signal is typically larger than the climatology in the upper troposphere and lower stratosphere and that a double tropopause during deep convection can easily be detected using this technique. Comparisons with co-located radiosondes, climatology of tropopause altitudes and GOES analyses are also shown to support the hypothesis that the bending angle anomaly can be used as an indicator of convective towers. The results are discussed in connection to the GPS radio occultation receiver which will be part of the Atomic Clock Ensemble in Space (ACES) payload on the International Space Station.

Description: Assimilation of GPS radio occultation data at DWD Atmospheric Measurement Techniques, 4, 1105-1113, 2011 Author(s): H. Anlauf, D. Pingel, and A. Rhodin We describe the status of the assimilation of bending angles from GPS radio occultations in the 3D-Var for DWD's operational global forecast model GME ("Global Model for Europe"). Experiments show that the assimilation of GPSRO data leads to a significant reduction of biases in the analyses of temperature, humidity and wind in the upper troposphere and the stratosphere, as well as a better r. m. s. fit in the comparison to radiosondes. The impact on forecasts is most prominent in the data sparse Southern Hemisphere, but is also quite notable in the Northern Hemisphere extra-tropics. The positive results found in the impact experiments lead to the implementation of the assimilation of GPS radio occultations from GRACE-A, FORMOSAT-3/COSMIC and GRAS/MetOp-A into the operational suite on 3 August 2010. We also show some initial results from assimilation experiments using radio occultation data from the German research satellite TerraSAR-X.

Description: Exploring the potentials and limitations of the time-reversal imaging of finite seismic sources Solid Earth, 2, 95-105, 2011 Author(s): S. Kremers, A. Fichtner, G. B. Brietzke, H. Igel, C. Larmat, L. Huang, and M. Käser The characterisation of seismic sources with time-reversed wave fields is developing into a standard technique that has already been successful in numerous applications. While the time-reversal imaging of effective point sources is now well-understood, little work has been done to extend this technique to the study of finite rupture processes. This is despite the pronounced non-uniqueness in classic finite source inversions. The need to better constrain the details of finite rupture processes motivates the series of synthetic and real-data time reversal experiments described in this paper. We address questions concerning the quality of focussing in the source area, the localisation of the fault plane, the estimation of the slip distribution and the source complexity up to which time-reversal imaging can be applied successfully. The frequency band for the synthetic experiments is chosen such that it is comparable to the band usually employed for finite source inversion. Contrary to our expectations, we find that time-reversal imaging is useful only for effective point sources, where it yields good estimates of both the source location and the origin time. In the case of finite sources, however, the time-reversed field does not provide meaningful characterisations of the fault location and the rupture process. This result cannot be improved sufficiently with the help of different imaging fields, realistic modifications of the receiver geometry or weights applied to the time-reversed sources. The reasons for this failure are manifold. They include the choice of the frequency band, the incomplete recording of wave field information at the surface, the excitation of large-amplitude surface waves that deteriorate the depth resolution, the absence of a sink that should absorb energy radiated during the later stages of the rupture process, the invisibility of small slip and the neglect of prior information concerning the fault geometry and the inherent smoothness of seismologically inferred Earth models that prevents the beneficial occurrence of strong multiple-scattering. The condensed conclusion of our study is that the limitations of time-reversal imaging – at least in the frequency band considered here – start where the seismic source stops being effectively point-localised.

Description: CIAO: the CNR-IMAA advanced observatory for atmospheric research Atmospheric Measurement Techniques, 4, 1191-1208, 2011 Author(s): F. Madonna, A. Amodeo, A. Boselli, C. Cornacchia, V. Cuomo, G. D'Amico, A. Giunta, L. Mona, and G. Pappalardo Long-term observations of aerosol and clouds are of crucial importance to understand the weather climate system. At the Istituto di Metodologie per l'Analisi Ambientale of the Italian National Research Council (CNR-IMAA) an advanced atmospheric observatory, named CIAO, is operative. CIAO (CNR-IMAA Atmospheric Observatory) main scientific objective is the long term measurement for the climatology of aerosol and cloud properties. Its equipment addresses the state-of-the-art for the ground-based remote sensing of aerosol, water vapour and clouds including active and passive sensors, like lidars, ceilometers, radiometers, and a radar. This paper describes the CIAO infrastructure, its scientific activities as well as the observation strategy. The observation strategy is mainly organized in order to provide quality assured measurements for satellite validation and model evaluation and to fully exploit the synergy and integration of the active and passive sensors for the improvement of atmospheric profiling. Data quality is ensured both by the application of protocols and dedicated quality assurance programs mainly related to the projects and networks in which the infrastructure is involved. The paper also introduces examples of observations performed at CIAO and of the synergies and integration algorithms (using Raman lidar and microwave profiler data) developed and implemented at the observatory for the optimization and improvement of water vapour profiling. CIAO database represents an optimal basis to study the synergy between different sensors and to investigate aerosol-clouds interactions, and can give a significant contribution to the validation programs of the incoming new generation satellite missions.

Description: Airborne DOAS limb measurements of tropospheric trace gas profiles: case studies on the profile retrieval of O 4 and BrO Atmospheric Measurement Techniques, 4, 1241-1260, 2011 Author(s): C. Prados-Roman, A. Butz, T. Deutschmann, M. Dorf, L. Kritten, A. Minikin, U. Platt, H. Schlager, H. Sihler, N. Theys, M. Van Roozendael, T. Wagner, and K. Pfeilsticker A novel limb scanning mini-DOAS spectrometer for the detection of UV/vis absorbing radicals (e.g., O 3 , BrO, IO, HONO) was deployed on the DLR-Falcon (Deutsches Zentrum für Luft- und Raumfahrt) aircraft and tested during the ASTAR 2007 campaign (Arctic Study of Tropospheric Aerosol, Clouds and Radiation) that took place at Svalbard (78° N) in spring 2007. Our main objectives during this campaign were to test the instrument, and to perform spectral and profile retrievals of tropospheric trace gases, with particular interest on investigating the distribution of halogen compounds (e.g., BrO) during the so-called ozone depletion events (ODEs). In the present work, a new method for the retrieval of vertical profiles of tropospheric trace gases from tropospheric DOAS limb observations is presented. Major challenges arise from modeling the radiative transfer in an aerosol and cloud particle loaded atmosphere, and from overcoming the lack of a priori knowledge of the targeted trace gas vertical distribution (e.g., unknown tropospheric BrO vertical distribution). Here, those challenges are tackled by a mathematical inversion of tropospheric trace gas profiles using a regularization approach constrained by a retrieved vertical profile of the aerosols extinction coefficient E M . The validity and limitations of the algorithm are tested with in situ measured E M , and with an absorber of known vertical profile (O 4 ). The method is then used for retrieving vertical profiles of tropospheric BrO. Results indicate that, for aircraft ascent/descent observations, the limit for the BrO detection is roughly 1.5 pptv (pmol mol −1 ), and the BrO profiles inferred from the boundary layer up to the upper troposphere and lower stratosphere have around 10 degrees of freedom. For the ASTAR 2007 deployments during ODEs, the retrieved BrO vertical profiles consistently indicate high BrO mixing ratios (∼15 pptv) within the boundary layer, low BrO mixing ratios (≤1.5 pptv) in the free troposphere, occasionally enhanced BrO mixing ratios (∼1.5 pptv) in the upper troposphere, and increasing BrO mixing ratios with altitude in the lowermost stratosphere. These findings agree reasonably well with satellite and balloon-borne soundings of total and partial BrO atmospheric column densities.

Description: Processing of GRAS/METOP radio occultation data recorded in closed-loop and raw-sampling modes Atmospheric Measurement Techniques, 4, 1021-1026, 2011 Author(s): M. E. Gorbunov, K. B. Lauritsen, H.-H. Benzon, G. B. Larsen, S. Syndergaard, and M. B. Sørensen Instrument GRAS (Global Navigation Satellite System Receiver for Atmospheric Sounding) on-board of the Metop-A satellite was activated on 27 October 2006. Currently, Metop-A is a fully operational satellite with GRAS providing from 650–700 occultations per day. We describe our processing of GRAS data based on the modification of our OCC software, which was modified to become capable of reading and processing GRAS data. We perform a statistical comparison of bending angles and refractivities derived from GRAS data with those derived from ECMWF analyses. We conclude that GRAS data have error characteristics close to those of COSMIC data. In the height range 10–30 km, the systematic refractivity difference GRAS–ECMWF is of the order of 0.1–0.2 %, and the standard deviation is 0.3–0.6 %. In the lower troposphere GRAS refractivity and bending angle indicate a negative bias, which reaches its maximum value in the tropics. In particular the retrieved refractivity is biased by up to 2.5 %. The negative bias pattern is similar to that found in the statistical validation of COSMIC data. This makes it probable that the bias should not be attributed to the instrument design or hardware.

Description: Eight-component retrievals from ground-based MAX-DOAS observations Atmospheric Measurement Techniques, 4, 1027-1044, 2011 Author(s): H. Irie, H. Takashima, Y. Kanaya, K. F. Boersma, L. Gast, F. Wittrock, D. Brunner, Y. Zhou, and M. Van Roozendael We attempt for the first time to retrieve lower-tropospheric vertical profile information for 8 quantities from ground-based Multi-Axis Differential Optical Absorption Spectroscopy (MAX-DOAS) observations. The components retrieved are the aerosol extinction coefficients at two wavelengths, 357 and 476 nm, and NO 2 , HCHO, CHOCHO, H 2 O, SO 2 , and O 3 volume mixing ratios. A Japanese MAX-DOAS profile retrieval algorithm, version 1 (JM1), is applied to observations performed at Cabauw, the Netherlands (51.97° N, 4.93° E), in June–July 2009 during the Cabauw Intercomparison campaign of Nitrogen Dioxide measuring Instruments (CINDI). Of the retrieved profiles, we focus here on the lowest-layer data (mean values at altitudes 0–1 km), where the sensitivity is usually highest owing to the longest light path. In support of the capability of the multi-component retrievals, we find reasonable overall agreement with independent data sets, including a regional chemical transport model (CHIMERE) and in situ observations performed near the surface (2–3 m) and at the 200-m height level of the tall tower in Cabauw. Plumes of enhanced HCHO and SO 2 were likely affected by biogenic and ship emissions, respectively, and an improvement in their emission strengths is suggested for better agreement between CHIMERE simulations and MAX-DOAS observations. Analysis of air mass factors indicates that the horizontal spatial representativeness of MAX-DOAS observations is about 3–15 km (depending mainly on aerosol extinction), comparable to or better than the spatial resolution of current UV-visible satellite observations and model calculations. These demonstrate that MAX-DOAS provides multi-component data useful for the evaluation of satellite observations and model calculations and can play an important role in bridging different data sets having different spatial resolutions.

Description: Native American lithic procurement along the international border in the boot heel region of southwestern New Mexico Solid Earth, 2, 75-93, 2011 Author(s): K. E. Zeigler, P. Hogan, C. Hughes, and A. Kurota Multidisciplinary field projects can be very useful to a more fundamental understanding of the world around us, though these projects are not as common as they should be. In particular, the combination of archeology and geology combines our understanding of human behavior and human use of the landscape with an intimate knowledge of geologic processes and the materials available for human use in order to gain a broader understanding of human-Earth interaction. Here we present data from a cross-disciplinary project that uses a common dataset, archeological artifacts, to explore the anthropological and geologic implications of useage patterns. Archeological excavations and surveys conducted by the Office of Contract Archeology in 2007 along the route of the proposed international border fence reveal patterns of use of geologic materials by Archaic, Formative and Protohistoric Native Americans in the Boot Heel of southwestern New Mexico. Thousands of artifacts were recorded in multiple sites from Guadalupe Pass in the southern Peloncillo Mountains to the Carrizalillo Hills west of Columbus. We identified the lithologies of artifacts, ranging from projectile points to groundstones, and then constructed material movement maps based on either known procurement sites ("quarries") or outcrops identified as the closest source to a given site for each lithology. Not unexpectedly, the majority of the rock types utilized by native peoples are local siliceous volcanic materials. However, several artifacts constructed from obsidian were transported into the region from northern Mexico and eastern Arizona, indicating long-distance travel and/or trade routes. We also examine useage pattern difference between Archaic, Formative and Protohistoric sites. Additionally, a dramatic change in distribution of sources for geologic materials occurs between one pre-Spanish site and one post-Spanish site that are adjacent to one another.

Description: Preliminary validation of column-averaged volume mixing ratios of carbon dioxide and methane retrieved from GOSAT short-wavelength infrared spectra Atmospheric Measurement Techniques, 4, 1061-1076, 2011 Author(s): I. Morino, O. Uchino, M. Inoue, Y. Yoshida, T. Yokota, P. O. Wennberg, G. C. Toon, D. Wunch, C. M. Roehl, J. Notholt, T. Warneke, J. Messerschmidt, D. W. T. Griffith, N. M. Deutscher, V. Sherlock, B. Connor, J. Robinson, R. Sussmann, and M. Rettinger Column-averaged volume mixing ratios of carbon dioxide and methane retrieved from the Greenhouse gases Observing SATellite (GOSAT) Short-Wavelength InfraRed observation (GOSAT SWIR X CO 2 and X CH 4 ) were compared with the reference calibrated data obtained by ground-based high-resolution Fourier Transform Spectrometers (g-b FTSs) participating in the Total Carbon Column Observing Network (TCCON). Preliminary results are as follows: the GOSAT SWIR X CO 2 and X CH 4 (Version 01.xx) are biased low by 8.85 ± 4.75 ppm (2.3 ± 1.2 %) and 20.4 ± 18.9 ppb (1.2 ± 1.1 %), respectively. The standard deviation of the GOSAT SWIR X CO 2 and X CH 4 is about 1 % (1 σ) after correcting the negative biases of X CO 2 and X CH 4 by 8.85 ppm and 20.4 ppb, respectively. The latitudinal distributions of zonal means of the GOSAT SWIR X CO 2 and X CH 4 show similar features to those of the g-b FTS data except for the negative biases in the GOSAT data.

Description: Remote sensing of aerosols over snow using infrared AATSR observations Atmospheric Measurement Techniques, 4, 1133-1145, 2011 Author(s): L. G. Istomina, W. von Hoyningen-Huene, A. A. Kokhanovsky, E. Schultz, and J. P. Burrows Infrared (IR) retrievals of aerosol optical thickness (AOT) are challenging because of the low reflectance of aerosol layer at longer wavelengths. In this paper we present a closer analysis of this problem, performed with radiative transfer (RT) simulations for coarse and accumulation mode of four main aerosol components. It shows the strong angular dependence of aerosol IR reflectance at low solar elevations resulting from the significant asymmetry of aerosol phase function at these wavelengths. This results in detectable values of aerosol IR reflectance at certain non-nadir observation angles providing the advantage of multiangle remote sensing instruments for a retrieval of AOT at longer wavelengths. Such retrievals can be of importance e.g. in case of a very strong effect of the surface on the top of atmosphere (TOA) reflectance in the visible spectral range. In the current work, a new method to retrieve AOT of the coarse and accumulation mode particles over snow has been developed using the measurements of Advanced Along Track Scanning Radiometer (AATSR) on board the ENVISAT satellite. The algorithm uses AATSR channel at 3.7 μm and utilizes its dual-viewing observation technique, implying the forward view with an observation zenith angle of around 55 degrees and the nadir view. It includes cloud/snow discrimination, extraction of the atmospheric reflectance out of measured brightness temperature (BT) at 3.7 μm, and interpolation of look-up tables (LUTs) for a given aerosol reflectance. The algorithm uses LUTs, separately simulated with RT forward calculations. The resulting AOT at 500 nm is estimated from the value at 3.7 μm using a fixed Angström parameter. The presented method has been validated against ground-based Aerosol Robotic Network (AERONET) data for 4 high Arctic stations and shows good agreement. A case study has been performed at W-Greenland on 5 July 2008. The day before was characterized by a noticeable dust event. The retrieved AOT maps of the region show a clear increase of AOT in the Kangerlussuaq area. The area of increased AOT was detected on 5 July on the ice sheet east of Kangelussuaq, opposite to the observed north easterly wind at ground level. This position can be explained by a small scale atmospheric circulation transporting the mobilized mineral dust upslope, after its intrusion into the upper branch of the circulation. The performed study of atmospheric reflectance at 3.7 μm also shows possibilities for the detection and retrievals of cloud properties over snow surfaces.

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: Analytical system for stable carbon isotope measurements of low molecular weight (C 2 -C 6 ) hydrocarbons Atmospheric Measurement Techniques, 4, 1161-1175, 2011 Author(s): A. Zuiderweg, R. Holzinger, and T. Röckmann We present setup, testing and initial results from a new automated system for stable carbon isotope ratio measurements on C 2 to C 6 atmospheric hydrocarbons. The inlet system allows analysis of trace gases from air samples ranging from a few liters for urban samples and samples with high mixing ratios, to many tens of liters for samples from remote unpolluted regions with very low mixing ratios. The centerpiece of the sample preparation is the separation trap, which is used to separate CO 2 and methane from the compounds of interest. The main features of the system are (i) the capability to sample up to 300 l of air, (ii) long term (since May 2009) operational δ 13 C accuracy levels in the range 0.3–0.8 ‰ (1-σ), and (iii) detection limits of order 1.5–2.5 ngC (collected amount of substance) for all reported compounds. The first application of this system was the analysis of 21 ambient air samples taken during 48 h in August 2009 in Utrecht, the Netherlands. Results obtained are generally in good agreement with those from similar urban ambient air studies. Short sample intervals allowed by the design of the instrument help to illustrate the complex diurnal behavior of hydrocarbons in an urban environment, where diverse sources, dynamical processes, and chemical reactions are present.

Description: Source brightness fluctuation correction of solar absorption fourier transform mid infrared spectra Atmospheric Measurement Techniques, 4, 1045-1051, 2011 Author(s): T. Ridder, T. Warneke, and J. Notholt The precision and accuracy of trace gas observations using solar absorption Fourier Transform infrared spectrometry depend on the stability of the light source. Fluctuations in the source brightness, however, cannot always be avoided. Current correction schemes, which calculate a corrected interferogram as the ratio of the raw DC interferogram and a smoothed DC interferogram, are applicable only to near infrared measurements. Spectra in the mid infrared spectral region below 2000 cm −1 are generally considered uncorrectable, if they are measured with a MCT detector. Such measurements introduce an unknown offset to MCT interferograms, which prevents the established source brightness fluctuation correction. This problem can be overcome by a determination of the offset using the modulation efficiency of the instrument. With known modulation efficiency the offset can be calculated, and the source brightness correction can be performed on the basis of offset-corrected interferograms. We present a source brightness fluctuation correction method which performs the smoothing of the raw DC interferogram in the interferogram domain by an application of a running mean instead of high-pass filtering the corresponding spectrum after Fourier transformation of the raw DC interferogram. This smoothing can be performed with the onboard software of commercial instruments. The improvement of MCT spectra and subsequent ozone profile and total column retrievals is demonstrated. Application to InSb interferograms in the near infrared spectral region proves the equivalence with the established correction scheme.

Description: A re-evaluation of the Italian historical geomagnetic catalogue: implications for paleomagnetic dating at active Italian volcanoes Solid Earth, 2, 65-74, 2011 Author(s): F. D'Ajello Caracciolo, A. Pignatelli, F. Speranza, and A. Meloni Paleomagnetism is proving to represent one of the most powerful dating tools of volcanics emplaced in Italy during the last few centuries/millennia. This method requires that valuable proxies of the local geomagnetic field (paleo)secular variation ((P)SV) are available. To this end, we re-evaluate the whole Italian geomagnetic directional dataset, consisting of 833 and 696 declination and inclination measurements, respectively, carried out since 1640 AD at several localities. All directions were relocated via the virtual geomagnetic pole method to Stromboli (38.8° N, 15.2° E), the rough centre of the active Italian volcanoes. For declination-only measurements, missing inclinations were derived (always by pole method) by French data (for period 1670–1789), and by nearby Italian sites/years (for periods 1640–1657 and 1790–1962). Using post-1825 declination values, we obtain a 0.46 ± 0.19° yr −1 westward drift of the geomagnetic field for Italy. The original observation years were modified, considering such drift value, to derive at a drift-corrected relocated dataset. Both datasets were found to be in substantial agreement with directions derived from the field models by Jackson et al. (2000) and Pavon-Carrasco et al. (2009). However, the drift-corrected dataset minimizes the differences between the Italian data and both field models, and eliminates a persistent 1.6° shift of 1933–1962 declination values from Castellaccio with respect to other nearly coeval Italian data. The relocated datasets were used to calculate two post-1640 Italian SV curves, with mean directions calculated every 30 and 10 years before and after 1790, respectively. The curve comparison suggests that both available field models yield the best available SV curve to perform paleomagnetic dating of 1600–1800 AD Italian volcanics, while the Italian drift-corrected curve is probably preferable for the 19th century. For the 20th century, the global model by Jackson et al. (2000) yields more accurate inclination values, while the declinations from our drift-corrected curve seem to better represent the local field evolution, at least for the first half of the century.

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: An improved NO 2 retrieval for the GOME-2 satellite instrument Atmospheric Measurement Techniques, 4, 1147-1159, 2011 Author(s): A. Richter, M. Begoin, A. Hilboll, and J. P. Burrows Satellite observations of nitrogen dioxide (NO 2 ) provide valuable information on both stratospheric and tropospheric composition. Nadir measurements from GOME, SCIAMACHY, OMI, and GOME-2 have been used in many studies on tropospheric NO 2 burdens, the importance of different NO x emissions sources and their change over time. The observations made by the three GOME-2 instruments will extend the existing data set by more than a decade, and a high quality of the data as well as their good consistency with existing time series is of particular importance. In this paper, an improved GOME-2 NO 2 retrieval is described which reduces the scatter of the individual NO 2 columns globally but in particular in the region of the Southern Atlantic Anomaly. This is achieved by using a larger fitting window including more spectral points, and by applying a two step spike removal algorithm in the fit. The new GOME-2 data set is shown to have good consistency with SCIAMACHY NO 2 columns. Remaining small differences are shown to be linked to changes in the daily solar irradiance measurements used in both GOME-2 and SCIAMACHY retrievals. In the large retrieval window, a not previously identified spectral signature was found which is linked to deserts and other regions with bare soil. Inclusion of this empirically derived pseudo cross-section significantly improves the retrievals and potentially provides information on surface properties and desert aerosols. Using the new GOME-2 NO 2 data set, a long-term average of tropospheric columns was computed and high-pass filtered. The resulting map shows evidence for pollution from several additional shipping lanes, not previously identified in satellite observations. This illustrates the excellent signal to noise ratio achievable with the improved GOME-2 retrievals.

Description: Matching radiative transfer models and radiosonde data from the EPS/Metop Sodankylä campaign to IASI measurements Atmospheric Measurement Techniques, 4, 1177-1189, 2011 Author(s): X. Calbet, R. Kivi, S. Tjemkes, F. Montagner, and R. Stuhlmann Radiances observed from IASI are compared to calculated ones. Calculated radiances are obtained using several radiative transfer models (OSS, LBLRTM v11.3 and v11.6) on best estimates of the atmospheric state vectors. The atmospheric state vectors are derived from cryogenic frost point hygrometer and humidity dry bias corrected RS92 measurements flown on sondes launched 1 h and 5 min before IASI overpass time. The temperature and humidity best estimate profiles are obtained by interpolating or extrapolating these measurements to IASI overpass time. The IASI observed and calculated radiances match to within one sigma IASI instrument noise in the spectral region where water vapour is a strong absorber (wavenumber, ν, in the range of 1500 ≤ ν ≤ 1570 and 1615 ≤ ν ≤ 1800 cm −1 ).

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: Eddy covariance measurements with high-resolution time-of-flight aerosol mass spectrometry: a new approach to chemically resolved aerosol fluxes Atmospheric Measurement Techniques, 4, 1275-1289, 2011 Author(s): D. K. Farmer, J. R. Kimmel, G. Phillips, K. S. Docherty, D. R. Worsnop, D. Sueper, E. Nemitz, and J. L. Jimenez Although laboratory studies show that biogenic volatile organic compounds (VOCs) yield substantial secondary organic aerosol (SOA), production of biogenic SOA as indicated by upward fluxes has not been conclusively observed over forests. Further, while aerosols are known to deposit to surfaces, few techniques exist to provide chemically-resolved particle deposition fluxes. To better constrain aerosol sources and sinks, we have developed a new technique to directly measure fluxes of chemically-resolved submicron aerosols using the high-resolution time-of-flight aerosol mass spectrometer (HR-AMS) in a new, fast eddy covariance mode. This approach takes advantage of the instrument's ability to quantitatively identify both organic and inorganic components, including ammonium, sulphate and nitrate, at a temporal resolution of several Hz. The new approach has been successfully deployed over a temperate ponderosa pine plantation in California during the BEARPEX-2007 campaign, providing both total and chemically resolved non-refractory (NR) PM 1 fluxes. Average deposition velocities for total NR-PM 1 aerosol at noon were 2.05 ± 0.04 mm s −1 . Using a high resolution measurement of the NH 2 + and NH 3 + fragments, we demonstrate the first eddy covariance flux measurements of particulate ammonium, which show a noon-time deposition velocity of 1.9 ± 0.7 mm s −1 and are dominated by deposition of ammonium sulphate.

Description: Space and ground segment performance and lessons learned of the FORMOSAT-3/COSMIC mission: four years in orbit Atmospheric Measurement Techniques, 4, 1115-1132, 2011 Author(s): C.-J. Fong, D. Whiteley, E. Yang, K. Cook, V. Chu, B. Schreiner, D. Ector, P. Wilczynski, T.-Y. Liu, and N. Yen The FORMOSAT-3/COSMIC (Constellation Observing System for Meteorology, Ionosphere, and Climate) Mission consisting of six Low-Earth-Orbit (LEO) satellites is the world's first demonstration constellation using radio occultation signals from Global Positioning System (GPS) satellites. The atmospheric profiles derived by processing radio occultation signals are retrieved in near real-time for global weather/climate monitoring, numerical weather prediction, and space weather research. The mission has processed, on average, 1400 to 1800 high-quality atmospheric sounding profiles per day. The atmospheric radio occultation data are assimilated into operational numerical weather prediction models for global weather prediction, including typhoon/hurricane/cyclone forecasts. The radio occultation data has shown a positive impact on weather predictions at many national weather forecast centers. A follow-on mission was proposed that transitions the current experimental research mission into a significantly improved real-time operational mission, which will reliably provide 8000 radio occultation soundings per day. The follow-on mission, as planned, will consist of 12 LEO satellites (compared to 6 satellites for the current mission) with data latency requirement of 45 min (compared to 3 h for the current mission), which will provide greatly enhanced opportunities for operational forecasts and scientific research. This paper will address the FORMOSAT-3/COSMIC system and mission overview, the spacecraft and ground system performance after four years in orbit, the lessons learned from the encountered technical challenges and observations, and the expected design improvements for the spacecraft and ground system for FORMOSAT-7/COSMIC-2.

Description: Exploring Earth's atmosphere with radio occultation: contributions to weather, climate and space weather Atmospheric Measurement Techniques, 4, 1077-1103, 2011 Author(s): R. A. Anthes The launch of the proof-of-concept mission GPS/MET (Global Positioning System/Meteorology) in 1995 began a revolution in profiling Earth's atmosphere through radio occultation (RO). GPS/MET; subsequent single-satellite missions CHAMP (CHAllenging Minisatellite Payload), SAC-C (Satellite de Aplicaciones Cientificas-C), GRACE (Gravity Recovery and Climate Experiment), METOP-A, and TerraSAR-X (Beyerle et al., 2010); and the six-satellite constellation, FORMOSAT-3/COSMIC (Formosa Satellite mission {#}3/Constellation Observing System for Meteorology, Ionosphere, and Climate) have proven the theoretical capabilities of RO to provide accurate and precise profiles of electron density in the ionosphere and refractivity, containing information on temperature and water vapor, in the stratosphere and troposphere. This paper summarizes results from these RO missions and the applications of RO observations to atmospheric research and operational weather analysis and prediction.

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: Detection of HO 2 by laser-induced fluorescence: calibration and interferences from RO 2 radicals Atmospheric Measurement Techniques, 4, 1209-1225, 2011 Author(s): H. Fuchs, B. Bohn, A. Hofzumahaus, F. Holland, K. D. Lu, S. Nehr, F. Rohrer, and A. Wahner HO 2 concentration measurements are widely accomplished by chemical conversion of HO 2 to OH including reaction with NO and subsequent detection of OH by laser-induced fluorescence. RO 2 radicals can be converted to OH via a similar radical reaction sequence including reaction with NO, so that they are potential interferences for HO 2 measurements. Here, the conversion efficiency of various RO 2 species to HO 2 is investigated. Experiments were conducted with a radical source that produces OH and HO 2 by water photolysis at 185 nm, which is frequently used for calibration of LIF instruments. The ratio of HO 2 and the sum of OH and HO 2 concentrations provided by the radical source was investigated and was found to be 0.50 ± 0.02. RO 2 radicals are produced by the reaction of various organic compounds with OH in the radical source. Interferences via chemical conversion from RO 2 radicals produced by the reaction of OH with methane and ethane (H-atom abstraction) are negligible consistent with measurements in the past. However, RO 2 radicals from OH plus alkene- and aromatic-precursors including isoprene (mainly OH-addition) are detected with a relative sensitivity larger than 80 % with respect to that for HO 2 for the configuration of the instrument with which it was operated during field campaigns. Also RO 2 from OH plus methyl vinyl ketone and methacrolein exhibit a relative detection sensitivity of 60 %. Thus, previous measurements of HO 2 radical concentrations with this instrument were biased in the presence of high RO 2 radical concentrations from isoprene, alkenes or aromatics, but were not affected by interferences in remote clean environment with no significant emissions of biogenic VOCs, when the OH reactivity was dominated by small alkanes. By reducing the NO concentration and/or the transport time between NO addition and OH detection, interference from these RO 2 species are suppressed to values below 20 % relative to the HO 2 detection sensitivity. The HO 2 conversion efficiency is also smaller by a factor of four, but this is still sufficient for atmospheric HO 2 concentration measurements for a wide range of conditions.

Description: Diode laser-based cavity ring-down instrument for NO 3 , N 2 O 5 , NO, NO 2 and O 3 from aircraft Atmospheric Measurement Techniques, 4, 1227-1240, 2011 Author(s): N. L. Wagner, W. P. Dubé, R. A. Washenfelder, C. J. Young, I. B. Pollack, T. B. Ryerson, and S. S. Brown This article presents a diode laser-based, cavity ring-down spectrometer for simultaneous in situ measurements of four nitrogen oxide species, NO 3 , N 2 O 5 , NO, NO 2 , as well as O 3 , designed for deployment on aircraft. The instrument measures NO 3 and NO 2 by optical extinction at 662 nm and 405 nm, respectively; N 2 O 5 is measured by thermal conversion to NO 3 , while NO and O 3 are measured by chemical conversion to NO 2 . The instrument has several advantages over previous instruments developed by our group for measurement of NO 2 , NO 3 and N 2 O 5 alone, based on a pulsed Nd:YAG and dye laser. First, the use of continuous wave diode lasers reduces the requirements for power and weight and eliminates hazardous materials. Second, detection of NO 2 at 405 nm is more sensitive than our previously reported 532 nm instrument, and does not have a measurable interference from O 3 . Third, the instrument includes chemical conversion of NO and O 3 to NO 2 to provide measurements of total NO x (= NO + NO 2 ) and O x (= NO 2 + O 3 ) on two separate channels; mixing ratios of NO and O 3 are determined by subtraction of NO 2 . Finally, all five species are calibrated against a single standard based on 254 nm O 3 absorption to provide high accuracy. Disadvantages include an increased sensitivity to water vapor on the 662 nm NO 3 and N 2 O 5 channels and a modest reduction in sensitivity for these species compared to the pulsed laser instrument. The in-flight detection limit for both NO 3 and N 2 O 5 is 3 pptv (2 σ, 1 s) and for NO, NO 2 and O 3 is 140, 90, and 120 pptv (2 σ, 1 s) respectively. Demonstrated performance of the instrument in a laboratory/ground based environment is better by approximately a factor of 2–3. The NO and NO 2 measurements are less precise than research-grade chemiluminescence instruments. However, the combination of these five species in a single instrument, calibrated to a single analytical standard, provides a complete and accurate picture of nighttime nitrogen oxide chemistry. The instrument performance is demonstrated using data acquired during a recent field campaign in California.

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: Retrieval of water vapor vertical distributions in the upper troposphere and the lower stratosphere from SCIAMACHY limb measurements Atmospheric Measurement Techniques, 4, 933-954, 2011 Author(s): A. Rozanov, K. Weigel, H. Bovensmann, S. Dhomse, K.-U. Eichmann, R. Kivi, V. Rozanov, H. Vömel, M. Weber, and J. P. Burrows This study describes the retrieval of water vapor vertical distributions in the upper troposphere and lower stratosphere (UTLS) altitude range from space-borne observations of the scattered solar light made in limb viewing geometry. First results using measurements from SCIAMACHY (Scanning Imaging Absorption spectroMeter for Atmospheric CHartographY) aboard ENVISAT (Environmental Satellite) are presented here. In previous publications, the retrieval of water vapor vertical distributions has been achieved exploiting either the emitted radiance leaving the atmosphere or the transmitted solar radiation. In this study, the scattered solar radiation is used as a new source of information on the water vapor content in the UTLS region. A recently developed retrieval algorithm utilizes the differential absorption structure of the water vapor in 1353–1410 nm spectral range and yields the water vapor content in the 11–25 km altitude range. In this study, the retrieval algorithm is successfully applied to SCIAMACHY limb measurements and the resulting water vapor profiles are compared to in situ balloon-borne observations. The results from both satellite and balloon-borne instruments are found to agree typically within 10 %.

Description: Open Plot Project: an open-source toolkit for 3-D structural data analysis Solid Earth, 2, 53-63, 2011 Author(s): S. Tavani, P. Arbues, M. Snidero, N. Carrera, and J. A. Muñoz In this work we present the Open Plot Project, an open-source software for structural data analysis, including a 3-D environment. The software includes many classical functionalities of structural data analysis tools, like stereoplot, contouring, tensorial regression, scatterplots, histograms and transect analysis. In addition, efficient filtering tools are present allowing the selection of data according to their attributes, including spatial distribution and orientation. This first alpha release represents a stand-alone toolkit for structural data analysis. The presence of a 3-D environment with digitalising tools allows the integration of structural data with information extracted from georeferenced images to produce structurally validated dip domains. This, coupled with many import/export facilities, allows easy incorporation of structural analyses in workflows for 3-D geological modelling. Accordingly, Open Plot Project also candidates as a structural add-on for 3-D geological modelling software. The software (for both Windows and Linux O.S.), the User Manual, a set of example movies (complementary to the User Manual), and the source code are provided as Supplement. We intend the publication of the source code to set the foundation for free, public software that, hopefully, the structural geologists' community will use, modify, and implement. The creation of additional public controls/tools is strongly encouraged.

Description: Fast NO 2 retrievals from Odin-OSIRIS limb scatter measurements Atmospheric Measurement Techniques, 4, 965-972, 2011 Author(s): A. E. Bourassa, C. A. McLinden, C. E. Sioris, S. Brohede, A. F. Bathgate, E. J. Llewellyn, and D. A. Degenstein The feasibility of retrieving vertical profiles of NO 2 from space-based measurements of limb scattered sunlight has been demonstrated using several different data sets since the 1980's. The NO 2 data product routinely retrieved from measurements made by the Optical Spectrograph and InfraRed Imaging System (OSIRIS) instrument onboard the Odin satellite uses a spectral fitting technique over the 437 to 451 nm range, over which there are 36 individual wavelength measurements. In this work we present a proof of concept technique for the retrieval of NO 2 using only 4 of the 36 OSIRIS measurements in this wavelength range, which reduces the computational cost by almost an order of magnitude. The method is an adaptation of a triplet analysis technique that is currently used for the OSIRIS retrievals of ozone at Chappuis band wavelengths. The results obtained are shown to be in very good agreement with the spectral fit method, and provide an important alternative for applications where the computational burden is very high. Additionally this provides a baseline for future instrument design in terms of cost effectiveness and reducing spectral range requirements.

Description: Intercomparison of stratospheric gravity wave observations with AIRS and IASI Atmospheric Measurement Techniques, 7, 4517-4537, 2014 Author(s): L. Hoffmann, M. J. Alexander, C. Clerbaux, A. W. Grimsdell, C. I. Meyer, T. Rößler, and B. Tournier Gravity waves are an important driver for the atmospheric circulation and have substantial impact on weather and climate. Satellite instruments offer excellent opportunities to study gravity waves on a global scale. This study focuses on observations from the Atmospheric Infrared Sounder (AIRS) onboard the National Aeronautics and Space Administration Aqua satellite and the Infrared Atmospheric Sounding Interferometer (IASI) onboard the European MetOp satellites. The main aim of this study is an intercomparison of stratospheric gravity wave observations of both instruments. In particular, we analyzed AIRS and IASI 4.3 μm brightness temperature measurements, which directly relate to stratospheric temperature. Three case studies showed that AIRS and IASI provide a clear and consistent picture of the temporal development of individual gravity wave events. Statistical comparisons based on a 5-year period of measurements (2008–2012) showed similar spatial and temporal patterns of gravity wave activity. However, the statistical comparisons also revealed systematic differences of variances between AIRS and IASI that we attribute to the different spatial measurement characteristics of both instruments. We also found differences between day- and nighttime data that are partly due to the local time variations of the gravity wave sources. While AIRS has been used successfully in many previous gravity wave studies, IASI data are applied here for the first time for that purpose. Our study shows that gravity wave observations from different hyperspectral infrared sounders such as AIRS and IASI can be directly related to each other, if instrument-specific characteristics such as different noise levels and spatial resolution and sampling are carefully considered. The ability to combine observations from different satellites provides an opportunity to create a long-term record, which is an exciting prospect for future climatological studies of stratospheric gravity wave activity.

Description: Collection efficiency of the soot-particle aerosol mass spectrometer (SP-AMS) for internally mixed particulate black carbon Atmospheric Measurement Techniques, 7, 4507-4516, 2014 Author(s): M. D. Willis, A. K. Y. Lee, T. B. Onasch, E. C. Fortner, L. R. Williams, A. T. Lambe, D. R. Worsnop, and J. P. D. Abbatt The soot-particle aerosol mass spectrometer (SP-AMS) uses an intra-cavity infrared laser to vaporize refractory black carbon (rBC) containing particles, making the particle beam–laser beam overlap critical in determining the collection efficiency (CE) for rBC and associated non-refractory particulate matter (NR-PM). This work evaluates the ability of the SP-AMS to quantify rBC and NR-PM mass in internally mixed particles with different thicknesses of organic coating. Using apparent relative ionization efficiencies for uncoated and thickly coated rBC particles, we report measurements of SP-AMS sensitivity to NR-PM and rBC, for Regal Black, the recommended particulate calibration material. Beam width probe (BWP) measurements are used to illustrate an increase in sensitivity for highly coated particles due to narrowing of the particle beam, which enhances the CE of the SP-AMS by increasing the laser beam–particle beam overlap. Assuming complete overlap for thick coatings, we estimate CE for bare Regal Black particles of 0.6 ± 0.1, which suggests that previously measured SP-AMS sensitivities to Regal Black were underestimated by up to a factor of 2. The efficacy of the BWP measurements is highlighted by studies at a busy road in downtown Toronto and at a non-roadside location, which show particle beam widths similar to, but greater than that of bare Regal Black and coated Regal Black, respectively. Further BWP measurements at field locations will help to constrain the range of CE for fresh and aged rBC-containing particles. The ability of the SP-AMS to quantitatively assess the composition of internally mixed particles is validated through measurements of laboratory-generated organic coated particles, which demonstrate that the SP-AMS can quantify rBC and NR-PM over a wide range of particle compositions and rBC core sizes.

Description: Reference quality upper-air measurements: GRUAN data processing for the Vaisala RS92 radiosonde Atmospheric Measurement Techniques, 7, 4463-4490, 2014 Author(s): R. J. Dirksen, M. Sommer, F. J. Immler, D. F. Hurst, R. Kivi, and H. Vömel The GCOS (Global Climate Observing System) Reference Upper-Air Network (GRUAN) data processing for the Vaisala RS92 radiosonde was developed to meet the criteria for reference measurements. These criteria stipulate the collection of metadata, the use of well-documented correction algorithms, and estimates of the measurement uncertainty. An important and novel aspect of the GRUAN processing is that the uncertainty estimates are vertically resolved. This paper describes the algorithms that are applied in version 2 of the GRUAN processing to correct for systematic errors in radiosonde measurements of pressure, temperature, humidity, and wind, as well as how the uncertainties related to these error sources are derived. Currently, the RS92 is launched on a regular basis at 13 out of 15 GRUAN sites. An additional GRUAN requirement for performing reference measurements with the RS92 is that the manufacturer-prescribed procedure for the radiosonde's preparation, i.e. heated reconditioning of the sensors and recalibration during ground check, is followed. In the GRUAN processing however, the recalibration of the humidity sensors that is applied during ground check is removed. For the dominant error source, solar radiation, laboratory experiments were performed to investigate and model its effect on the RS92's temperature and humidity measurements. GRUAN uncertainty estimates are 0.15 K for night-time temperature measurements and approximately 0.6 K at 25 km during daytime. The other uncertainty estimates are up to 6% relative humidity for humidity, 10–50 m for geopotential height, 0.6 hPa for pressure, 0.4–1 m s −1 for wind speed, and 1° for wind direction. Daytime temperature profiles for GRUAN and Vaisala processing are comparable and consistent within the estimated uncertainty. GRUAN daytime humidity profiles are up to 15% moister than Vaisala processed profiles, of which two-thirds is due to the radiation dry bias correction and one-third is due to an additional calibration correction. Redundant measurements with frost point hygrometers (CFH and NOAA FPH) show that GRUAN-processed RS92 humidity profiles and frost point data agree within 15% in the troposphere. No systematic biases occur, apart from a 5% dry bias for GRUAN data around −40 °C at night.

Description: Deployment of a sequential two-photon laser-induced fluorescence sensor for the detection of gaseous elemental mercury at ambient levels: fast, specific, ultrasensitive detection with parts-per-quadrillion sensitivity Atmospheric Measurement Techniques, 7, 4251-4265, 2014 Author(s): D. Bauer, S. Everhart, J. Remeika, C. Tatum Ernest, and A. J. Hynes The operation of a laser-based sensor for gas-phase elemental mercury, Hg(0), is described. It utilizes sequential two-photon laser excitation with detection of blue-shifted laser-induced fluorescence (LIF) to provide a highly specific detection scheme that precludes detection of anything other than atomic mercury. It has high sensitivity, fast temporal resolution, and can be deployed for in situ measurements in the open atmosphere with essentially no perturbation of the environment. An ambient sample can also be pulled through a fluorescence cell, allowing for standard addition calibrations of the concentration. No type of preconcentration is required and there appears to be no significant interferences from other atmospheric constituents, including gas-phase oxidized mercury species. As a consequence, it is not necessary to remove oxidized mercury, commonly referred to as reactive gaseous mercury (RGM), from the air sample. The instrument has been deployed as part of an instrument intercomparison and compares well with conventional instrumentation that utilizes preconcentration on gold followed by analysis using cold-vapor atomic fluorescence spectroscopy (CVAFS). Currently, the achievable detection sensitivity is ~ 15 pg m −3 (~ 5 × 10 4 atoms cm −3 , ~ 2 ppq) at a sampling rate of 0.1 Hz, i.e., averaging 100 shots with a 10 Hz laser system. Preliminary results are described for a 50 Hz instrument that utilizes a modified excitation sequence and has monitored ambient elemental mercury with an effective sampling rate of 10 Hz. Additional work is required to produce the precision necessary to perform eddy correlation measurements. Addition of a pyrolysis channel should allow for the measurement of total gaseous mercury (TGM) and hence RGM (by difference) with good sensitivity and time resolution.

Description: Prerequisites for application of hyperbolic relaxed eddy accumulation on managed grasslands and alternative net ecosystem exchange flux partitioning Atmospheric Measurement Techniques, 7, 4237-4250, 2014 Author(s): M. Riederer, J. Hübner, J. Ruppert, W. A. Brand, and T. Foken Relaxed eddy accumulation is still applied in ecosystem sciences for measuring trace gas fluxes. On managed grasslands, the length of time between management events and the application of relaxed eddy accumulation has an essential influence on the determination of the proportionality factor b and thus on the resulting flux. In this study this effect is discussed for the first time. Also, scalar similarity between proxy scalars and scalars of interest is affected until the ecosystem has completely recovered. Against this background, CO 2 fluxes were continuously measured and 13 CO 2 isofluxes were determined with a high measurement precision on two representative days in summer 2010. Moreover, a common method for the partitioning of the net ecosystem exchange into assimilation and respiration based on temperature and light response was compared with an isotopic approach directly based on the isotope discrimination of the biosphere. This approach worked well on the grassland site and could enhance flux partitioning results by better reproducing the environmental conditions.

Description: Extending the satellite data record of tropospheric ozone profiles from Aura-TES to MetOp-IASI: characterisation of optimal estimation retrievals Atmospheric Measurement Techniques, 7, 4223-4236, 2014 Author(s): H. Oetjen, V. H. Payne, S. S. Kulawik, A. Eldering, J. Worden, D. P. Edwards, G. L. Francis, H. M. Worden, C. Clerbaux, J. Hadji-Lazaro, and D. Hurtmans We apply the Tropospheric Emission Spectrometer (TES) ozone retrieval algorithm to Infrared Atmospheric Sounding Instrument (IASI) radiances and characterise the uncertainties and information content of the retrieved ozone profiles. This study focuses on mid-latitudes for the year 2008. We validate our results by comparing the IASI ozone profiles to ozone sondes. In the sonde comparisons, we find a negative bias (1–10%) in the IASI profiles in the lower to mid-troposphere and a positive bias (up to 14%) in the upper troposphere/lower stratosphere (UTLS) region. For the described cases, the degrees of freedom for signal are on average 3.2, 0.3, 0.8, and 0.9 for the columns 0 km – top of atmosphere, (0–6), (0–11), and (8–16) km, respectively. We find that our biases with respect to sondes and our degrees of freedom for signal for ozone are comparable to previously published results from other IASI ozone algorithms. In addition to evaluating biases, we validate the retrieval errors by comparing predicted errors to the sample covariance matrix of the IASI observations themselves. For the predicted versus empirical error comparison, we find that these errors are consistent and that the measurement noise and the interference of temperature and water vapour on the retrieval together mostly explain the empirically derived random errors. In general, the precision of the IASI ozone profiles is better than 20%.

Description: Solar irradiances measured using SPN1 radiometers: uncertainties and clues for development Atmospheric Measurement Techniques, 7, 4267-4283, 2014 Author(s): J. Badosa, J. Wood, P. Blanc, C. N. Long, L. Vuilleumier, D. Demengel, and M. Haeffelin The fast development of solar radiation and energy applications, such as photovoltaic and solar thermodynamic systems, has increased the need for solar radiation measurement and monitoring, for not only the global but also the diffuse and direct components. End users look for the best compromise between getting close to state-of-the-art measurements and keeping low capital, maintenance and operating costs. Among the existing commercial options, SPN1 is a relatively low cost solar radiometer that estimates global and diffuse solar irradiances from seven thermopile sensors under a shading mask and without moving parts. This work presents a comprehensive study of SPN1 accuracy and sources of uncertainty, drawing on laboratory experiments, numerical modelling and comparison studies between measurements from this sensor and state-of-the art instruments for six diverse sites. Several clues are provided for improving the SPN1 accuracy and agreement with state-of-the art measurements.

Description: Differences in aerosol absorption Ångström exponents between correction algorithms for a particle soot absorption photometer measured on the South African Highveld Atmospheric Measurement Techniques, 7, 4285-4298, 2014 Author(s): J. Backman, A. Virkkula, V. Vakkari, J. P. Beukes, P. G. Van Zyl, M. Josipovic, S. Piketh, P. Tiitta, K. Chiloane, T. Petäjä, M. Kulmala, and L. Laakso Absorption Ångström exponents (AAEs) calculated from filter-based absorption measurements are often used to give information on the origin of the ambient aerosol, for example, to distinguish between urban pollution and biomass burning aerosol. Filter-based absorption measurements are widely used and are common at aerosol monitoring stations globally. Several correction algorithms are used to account for artefacts associated with filter-based absorption techniques. These algorithms are of profound importance when determining the absolute amount of absorption by the aerosol. However, this study shows that there are substantial differences between the AAEs calculated from these corrections. Depending on the used correction, AAEs can change by as much as 46%. The study also highlights that the difference between AAEs calculated using different corrections can lead to conflicting conclusions on the type of aerosol when using the same data set. The AAE ranged between 1.17 for non-corrected data to 1.96 for the correction that gave the greatest values. Furthermore, the study implies that the AAEs reported for a site depend on at which filter transmittance the filter is changed. In this work, the AAEs were calculated from data measured with a three-wavelength particle soot absorption photometer (PSAP) at Elandsfontein on the South African Highveld for 23 months. The sample air of the PSAP was diluted to prolong filter change intervals, by a factor of 15. The correlation coefficient between the dilution-corrected PSAP and a non-diluted Multi-Angle Absorption Photometer (MAAP) was 0.9. Thus, the study also shows that the applicability of the PSAP can be extended to remote sites that are not often visited or suffer from high levels of pollution.

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: Looking through the haze: evaluating the CALIPSO level 2 aerosol optical depth using airborne high spectral resolution lidar data Atmospheric Measurement Techniques, 7, 4317-4340, 2014 Author(s): R. R. Rogers, M. A. Vaughan, C. A. Hostetler, S. P. Burton, R. A. Ferrare, S. A. Young, J. W. Hair, M. D. Obland, D. B. Harper, A. L. Cook, and D. M. Winker The Cloud–Aerosol Lidar with Orthogonal Polarization (CALIOP) instrument onboard the Cloud–Aerosol Lidar and Pathfinder Satellite Observations (CALIPSO) spacecraft has provided over 8 yr of nearly continuous vertical profiling of Earth's atmosphere. In this paper we investigate the V3.01 and V3.02 CALIOP 532 nm aerosol layer optical depth (AOD) product (i.e the AOD of individual layers) and the column AOD product (i.e., the sum AOD of the complete column) using an extensive database of coincident measurements. The CALIOP AOD measurements and AOD uncertainty estimates are compared with collocated AOD measurements collected with the NASA High Spectral Resolution Lidar (HSRL) in the North American and Caribbean regions. In addition, the CALIOP aerosol lidar ratios are investigated using the HSRL measurements. In general, compared with the HSRL values, the CALIOP layer AOD are biased high by less than 50% for AOD 〈 0.3 with higher errors for higher AOD. Less than 60% of the HSRL AOD measurements are encompassed within the CALIOP layer 1 SD uncertainty range (around the CALIOP layer AOD), so an error estimate is created to encompass 68% of the HSRL data. Using this new metric, the CALIOP layer AOD error is estimated using the HSRL layer AOD as ±0.035 ± 0.05 · (HSRL layer AOD) at night and ±0.05 ± 0.05 · (HSRL layer AOD) during the daytime. Furthermore, the CALIOP layer AOD error is found to correlate with aerosol loading as well as aerosol subtype, with the AODs in marine and dust layers agreeing most closely with the HSRL values. The lidar ratios used by CALIOP for polluted dust, polluted continental, and biomass burning layers are larger than the values measured by the HSRL in the CALIOP layers, and therefore the AODs for these types retrieved by CALIOP were generally too large. We estimated the CALIOP column AOD error can be expressed as ±0.05 ± 0.07 · (HSRL column AOD) at night and ±0.08 ± 0.1 · (HSRL column AOD) during the daytime. Multiple sources of error contribute to both positive and negative errors in the CALIOP column AOD, including multiple layers in the column of different aerosol types, lidar ratio errors, cloud misclassification, and undetected aerosol layers. The undetected layers were further investigated and we found that the layer detection algorithm works well at night, although undetected aerosols in the free troposphere introduce a mean underestimate of 0.02 in the column AOD in the data set examined. The decreased signal-to-noise ratio (SNR) during the daytime led to poorer performance of the layer detection. This caused the daytime CALIOP column AOD to be less accurate than during the nighttime, because CALIOP frequently does not detect optically thin aerosol layers with AOD 〈 0.1. Given that the median vertical extent of aerosol detected within any column was 1.6 km during the nighttime and 1.5 km during the daytime, we can estimate the minimum extinction detection threshold to be 0.012 km −1 at night and 0.067 km −1 during the daytime in a layer median sense. This extensive validation of level 2 CALIOP AOD products extends previous validation studies to nighttime lighting conditions and provides independent measurements of the lidar ratio; thus, allowing the assessment of the effect on the CALIOP AOD of using inappropriate lidar ratio values in the extinction retrieval.

Description: The use of NO 2 absorption cross section temperature sensitivity to derive NO 2 profile temperature and stratospheric–tropospheric column partitioning from visible direct-sun DOAS measurements Atmospheric Measurement Techniques, 7, 4299-4316, 2014 Author(s): E. Spinei, A. Cede, W. H. Swartz, J. Herman, and G. H. Mount This paper presents a temperature sensitivity method (TESEM) to accurately calculate total vertical NO 2 column, atmospheric slant NO 2 profile-weighted temperature ( T ), and to separate stratospheric and tropospheric columns from direct-sun (DS), ground-based measurements using the retrieved T . TESEM is based on differential optical absorption spectroscopy (DOAS) fitting of the linear temperature-dependent NO 2 absorption cross section, σ ( T ), regression model (Vandaele et al., 2003). Separation between stratospheric and tropospheric columns is based on the primarily bimodal vertical distribution of NO 2 and an assumption that stratospheric effective temperature can be represented by temperature at 27 km ± 3 K, and tropospheric effective temperature is equal to surface temperature within 3–5 K. These assumptions were derived from the Global Modeling Initiative (GMI) chemistry-transport model (CTM) simulations over two northern midlatitude sites in 2011. TESEM was applied to the Washington State University Multi-Function DOAS instrument (MFDOAS) measurements at four midlatitude locations with low and moderate NO 2 anthropogenic emissions: (1) the Jet Propulsion Laboratory's Table Mountain Facility (JPL-TMF), CA, USA (34.38° N/117.68° W); (2) Pullman, WA, USA (46.73° N/117.17° W); (3) Greenbelt, MD, USA (38.99° N/76.84° W); and (4) Cabauw, the Netherlands (51.97° N/4.93° E) during July 2007, June–July 2009, July–August and October 2011, November 2012–May 2013, respectively. NO 2 T and total, stratospheric, and tropospheric NO 2 vertical columns were determined over each site.

Description: A Fabry–Perot interferometer-based camera for two-dimensional mapping of SO 2 distributions Atmospheric Measurement Techniques, 7, 3705-3715, 2014 Author(s): J. Kuhn, N. Bobrowski, P. Lübcke, L. Vogel, and U. Platt We examine a new imaging method for the remote sensing of volcanic gases, which relies on the regularly spaced narrow-band absorption structures in the UV–VIS of many molecules. A Fabry–Perot interferometer (FPI) is used to compare the scattered sunlight radiance at wavelengths corresponding to absorption bands with the radiance at wavelengths in between the bands, thereby identifying and quantifying the gas. In this first theoretical study, we present sample calculations for the detection of sulfur dioxide (SO 2 ). Optimum values for the FPI setup parameters are proposed. Furthermore, the performance of the FPI method is compared to SO 2 cameras. We show that camera systems using an FPI are far less influenced by changes in atmospheric radiative transfer (e.g., due to aerosol) and have a great potential as a future technique for examining emissions of SO 2 (or other gases) from volcanic sources and other point sources.

Description: Middle-atmospheric zonal and meridional wind profiles from polar, tropical and midlatitudes with the ground-based microwave Doppler wind radiometer WIRA Atmospheric Measurement Techniques, 7, 4491-4505, 2014 Author(s): R. Rüfenacht, A. Murk, N. Kämpfer, P. Eriksson, and S. A. Buehler WIRA is a ground-based microwave Doppler spectroradiometer specifically designed for the measurement of profiles of horizontal wind in the upper stratosphere and lower mesosphere region where no other continuously running measurement technique exists. A proof of principle has been delivered in a previous publication. A technical upgrade including a new high-frequency amplifier and sideband filter has improved the signal to noise ratio by a factor of 2.4. Since this upgrade the full horizontal wind field comprising zonal and meridional wind profiles is continuously measured. A completely new retrieval based on optimal estimation has been set up. Its characteristics are detailed in the present paper. Since the start of the routine operation of the first prototype in September 2010, WIRA has been measuring at four different locations at polar, mid- and tropical latitudes (67°22' N/26°38' E, 46°57' N/7°26' E, 43°56' N/5°43' E and 21°04' S/55°23' E) for time periods between 5.5 and 11 months. The data presented in this paper are daily average wind profiles with typical uncertainties and resolutions of 10 to 20 m s −1 and 10 to 16 km, respectively. A comparison between the data series from WIRA and European Centre for Medium-Range Weather Forecasts (ECMWF) model data revealed agreement within 10% in the stratospheric zonal wind. The meridional wind profiles agree within their error bars over the entire sensitive altitude range of WIRA. However, significant differences in the mesospheric zonal wind speed of up to 50% have been found.

Description: Interpretative modelling of a geological cross section from boreholes: sources of uncertainty and their quantification Solid Earth, 5, 1189-1203, 2014 Author(s): R. M. Lark, S. Thorpe, H. Kessler, and S. J. Mathers We conducted a designed experiment to quantify sources of uncertainty in geologists' interpretations of a geological cross section. A group of 28 geologists participated in the experiment. Each interpreted borehole record included up to three Palaeogene bedrock units, including the target unit for the experiment: the London Clay. The set of boreholes was divided into batches from which validation boreholes had been withheld; as a result, we obtained 129 point comparisons between the interpreted elevation of the base of the London Clay and its observed elevation in a borehole not used for that particular interpretation. Analysis of the results showed good general agreement between the observed and interpreted elevations, with no evidence of systematic bias. Between-site variation of the interpretation error was spatially correlated, and the variance appeared to be stationary. The between-geologist component of variance was smaller overall, and depended on the distance to the nearest borehole. There was also evidence that the between-geologist variance depends on the degree of experience of the individual. We used the statistical model of interpretation error to compute confidence intervals for any one interpretation of the base of the London Clay on the cross section, and to provide uncertainty measures for decision support in a hypothetical route-planning process. The statistical model could also be used to quantify error propagation in a full 3-D geological model produced from interpreted cross sections.

Description: Tunable diode laser measurements of hydrothermal/volcanic CO 2 and implications for the global CO 2 budget Solid Earth, 5, 1209-1221, 2014 Author(s): M. Pedone, A. Aiuppa, G. Giudice, F. Grassa, V. Francofonte, B. Bergsson, and E. Ilyinskaya Quantifying the CO 2 flux sustained by low-temperature fumarolic fields in hydrothermal/volcanic environments has remained a challenge, to date. Here, we explored the potential of a commercial infrared tunable laser unit for quantifying such fumarolic volcanic/hydrothermal CO 2 fluxes. Our field tests were conducted between April 2013 and March 2014 at Nea Kameni (Santorini, Greece), Hekla and Krýsuvík (Iceland) and Vulcano (Aeolian Islands, Italy). At these sites, the tunable laser was used to measure the path-integrated CO 2 mixing ratios along cross sections of the fumaroles' atmospheric plumes. By using a tomographic post-processing routine, we then obtained, for each manifestation, the contour maps of CO 2 mixing ratios in the plumes and, from their integration, the CO 2 fluxes. The calculated CO 2 fluxes range from low (5.7 ± 0.9 t d −1 ; Krýsuvík) to moderate (524 ± 108 t d −1 ; La Fossa crater, Vulcano). Overall, we suggest that the cumulative CO 2 contribution from weakly degassing volcanoes in the hydrothermal stage of activity may be significant at the global scale.

Description: Analysis of internal gravity waves with GPS RO density profiles Atmospheric Measurement Techniques, 7, 4123-4132, 2014 Author(s): P. Šácha, U. Foelsche, and P. Pišoft GPS radio occultation (RO) data have proved to be a great tool for atmospheric monitoring and studies. In the past decade, they were frequently used for analyses of the internal gravity waves in the upper troposphere and lower stratosphere region. Atmospheric density is the first quantity of state gained in the retrieval process and is not burdened by additional assumptions. However, there are no studies elaborating in detail the utilization of GPS RO density profiles for gravity wave analyses. In this paper, we introduce a method for density background separation and a methodology for internal gravity wave analysis using the density profiles. Various background choices are discussed and the correspondence between analytical forms of the density and temperature background profiles is examined. In the stratosphere, a comparison between the power spectrum of normalized density and normalized dry temperature fluctuations confirms the suitability of the density profiles' utilization. In the height range of 8–40 km, results of the continuous wavelet transform are presented and discussed. Finally, the limits of our approach are discussed and the advantages of the density usage are listed.

Description: Characterization of model errors in the calculation of tangent heights for atmospheric infrared limb measurements Atmospheric Measurement Techniques, 7, 4117-4122, 2014 Author(s): M. Ridolfi and L. Sgheri We review the main factors driving the calculation of the tangent height of spaceborne limb measurements: the ray-tracing method, the refractive index model and the assumed atmosphere. We find that commonly used ray tracing and refraction models are very accurate, at least in the mid-infrared. The factor with largest effect in the tangent height calculation is the assumed atmosphere. Using a climatological model in place of the real atmosphere may cause tangent height errors up to ± 200 m. Depending on the adopted retrieval scheme, these errors may have a significant impact on the derived profiles.

Description: Recovering long-term aerosol optical depth series (1976–2012) from an astronomical potassium-based resonance scattering spectrometer Atmospheric Measurement Techniques, 7, 4103-4116, 2014 Author(s): A. Barreto, E. Cuevas, P. Pallé, P. M. Romero, C. Guirado, C. J. Wehrli, and F. Almansa A 37-year long-term series of monochromatic aerosol optical depth (AOD) has been recovered from solar irradiance measurements performed with the solar spectrometer Mark-I, deployed at Izaña mountain since 1976. The instrument operation is based on the method of resonant scattering, which affords wavelength absolute reference and stability (long-term stability and high precision) in comparison to other instruments based purely on interference filters. However, it has been specifically designed as a reference instrument for helioseismology, and its ability to determine AOD from transmitted and scattered monochromatic radiation at 769.9 nm inside a potassium vapour cell in the presence of a permanent magnetic field is evaluated in this paper. Particularly, the use of an exposed mirror arrangement to collect sunlight as well as the Sun–laboratory velocity dependence of the scattered component introduces some important inconveniences to overcome when we perform the instrument's calibration. We have solved this problem using a quasi-continuous Langley calibration technique and a refinement procedure to correct for calibration errors as well as for the fictitious diurnal cycle on AOD data. Our results showed similar calibration errors retrieved by means of this quasi-continuous Langley technique applied in different aerosol load events (from 0.04 to 0.3), provided aerosol concentration remains constant throughout the calibration interval. It assures the validity of this technique when it is applied in those periods with relatively high aerosol content. The comparative analysis between the recovered AOD data set from the Mark-I and collocated quasi-simultaneous data from the Cimel-AErosol RObotic NETwork (AERONET) and Precision Filter Radiometer (PFR) instruments showed an absolute mean bias ≤ 0.01 in the 10- and 12-year comparison, respectively. High correlation coefficients between AERONET and Mark-I and PFR/Mark-I pairs confirmed a very good linear relationship between instruments, proving that recovered AOD data series from Mark-I can be used together with PFR and AERONET AOD data to build a long-term AOD data series at the Izaña site (1976–now), suitable for future analysis of aerosols trends and inter-annual variability. Finally, the AOD preliminary trend analysis in the 29-year period from 1984 to 2012 with Mark-I AOD revealed no significant trends.