ALBERT

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

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

Description: Magnetic field and dynamic pressure ULF fluctuations in coronal-mass-ejection-driven sheath regions Annales Geophysicae, 31, 1559-1567, 2013 Author(s): E. K. J. Kilpua, H. Hietala, H. E. J. Koskinen, D. Fontaine, and L. Turc Compressed sheath regions form ahead of interplanetary coronal mass ejections (ICMEs) that are sufficiently faster than the preceding solar wind. The turbulent sheath regions are important drivers of magnetospheric activity, but due to their complex internal structure, relatively little is known on the distribution of the magnetic field and plasma variations in them. In this paper we investigate ultra low frequency (ULF) fluctuations in the interplanetary magnetic field (IMF) and in dynamic pressure ( P dyn ) using a superposed epoch analysis of 41 sheath regions observed during solar cycle 23. We find strongest fluctuation power near the shock and in the vicinity of the ICME leading edge. The IMF and P dyn ULF power have different profiles within the sheath; the former is enhanced in the leading part of the sheath, while the latter is increased in the trailing part of the sheath. We also find that the ICME properties affect the level and distribution of the ULF power in sheath regions. For example, sheath regions associated with strong or fast ICMEs, or those that are crossed at intermediate distances from the center, have strongest ULF power and large variation in the power throughout the sheath region. The weaker or slower ICMEs, or those that are crossed centrally, have in general considerably weaker ULF power with relatively smooth profiles. The strong and abrupt decrease of the IMF ULF power at the ICME leading edge could be used to distinguish the ICME from the preceding sheath plasma.

Description: Cluster observation of few-hour-scale evolution of structured plasma in the inner magnetosphere Annales Geophysicae, 31, 1569-1578, 2013 Author(s): M. Yamauchi, I. Dandouras, H. Rème, R. Lundin, and L. M. Kistler Using Cluster Ion Spectrometry (CIS) data from the spacecraft-4 perigee traversals during the 2001–2006 period (nearly 500 traversals after removing those that are highly contaminated by radiation belt particles), we statistically examined the local time distribution of structured trapped ions at sub- to few-keV range as well as inbound–outbound differences of these ion signatures in intensities and energy–latitude dispersion directions. Since the Cluster orbit during this period was almost constant and approximately north–south symmetric at nearly constant local time near the perigee, inbound–outbound differences are attributed to temporal developments in a 1–2 h timescale. Three types of structured ions at sub- to few keV range that are commonly found in the inner magnetosphere are examined: – Energy–latitude dispersed structured ions at less than a few keV, – Short-lived dispersionless ion stripes at wide energy range extending 0.1–10 keV, – Short-lived low-energy ion bursts at less than a few hundred eV. The statistics revealed that the wedge-like dispersed ions are most often observed in the dawn sector (60% of traversals), and a large portion of them show significant enhancement during the traversals at all local times. The short-lived ion stripes are predominantly found near midnight, where most stripes are significantly enhanced during the traversals and are associated with substorm activities with geomagnetic AL 〈 −300 nT. The low-energy bursts are observed at all local times and under all geomagnetic conditions, with moderate peak of the occurrence rate in the afternoon sector. A large portion of them again show significant enhancement or decay during the traversals.

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

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

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

Description: Influence of the terrestrial magnetic field geometry on the cutoff rigidity of cosmic ray particles Annales Geophysicae, 31, 1637-1643, 2013 Author(s): K. Herbst, A. Kopp, and B. Heber Studies of the propagation of charged energetic particles in the Earth's magnetic field go back to Carl Størmer. In the end, his investigations finally lead to the definition of the so-called cutoff rigidity R C ; that is, the minimum momentum per charge a particle must have in order to reach a certain geographical location. Employing Monte Carlo simulations with the PLANETOCOSMICS code we investigate the correlation between the geomagnetic field structure and the cutoff rigidity. We show that the geometry of the magnetic field has a considerable influence on the resulting cutoff rigidity distribution. Furthermore, we will present a simple geometry-based parameter, δ B , which is able to reflect the location-dependent cutoff rigidity. We show that this correlation is also visible in the temporal evolution of the Earth's magnetic field, at least over the last 100 yr. Using latitude scans with neutron monitors, changes of the relative counting rates at different positions are calculated, showing small variations for, e.g., Kiel and Moscow, while large ones occur at Mexico City as well as on the British Virgin Islands.

Description: Electron acceleration at Jupiter: input from cyclotron-resonant interaction with whistler-mode chorus waves Annales Geophysicae, 31, 1619-1630, 2013 Author(s): E. E. Woodfield, R. B. Horne, S. A. Glauert, J. D. Menietti, and Y. Y. Shprits Jupiter has the most intense radiation belts of all the outer planets. It is not yet known how electrons can be accelerated to energies of 10 MeV or more. It has been suggested that cyclotron-resonant wave-particle interactions by chorus waves could accelerate electrons to a few MeV near the orbit of Io. Here we use the chorus wave intensities observed by the Galileo spacecraft to calculate the changes in electron flux as a result of pitch angle and energy diffusion. We show that, when the bandwidth of the waves and its variation with L are taken into account, pitch angle and energy diffusion due to chorus waves is a factor of 8 larger at L -shells greater than 10 than previously shown. We have used the latitudinal wave intensity profile from Galileo data to model the time evolution of the electron flux using the British Antarctic Survey Radiation Belt (BAS) model. This profile confines intense chorus waves near the magnetic equator with a peak intensity at ∼5° latitude. Electron fluxes in the BAS model increase by an order of magnitude for energies around 3 MeV. Extending our results to L = 14 shows that cyclotron-resonant interactions with chorus waves are equally important for electron acceleration beyond L = 10. These results suggest that there is significant electron acceleration by cyclotron-resonant interactions at Jupiter contributing to the creation of Jupiter's radiation belts and also increasing the range of L -shells over which this mechanism should be considered.

Description: The coupling between the solar wind and proton fluxes at GEO Annales Geophysicae, 31, 1631-1636, 2013 Author(s): R. J. Boynton, S. A. Billings, O. A. Amariutei, and I. Moiseenko The relationship between the solar wind and the proton flux at geosynchronous Earth orbit (GEO) is investigated using the error reduction ratio (ERR) analysis. The ERR analysis is able to search for the most appropriate inputs that control the evolution of the system. This approach is a black box method and is able to derive a mathematical model of a system from input-output data. This method is used to analyse eight energy ranges of the proton flux at GEO from 80 keV to 14.5 MeV. The inputs to the algorithm were solar wind velocity, density and pressure; the Dst index; the solar energetic proton (SEP) flux; and a function of the interplanetary magnetic field (IMF) tangential magnitude and clock angle. The results show that for lowest five energy channels (80 to 800 keV) the GEO proton fluxes are controlled by the solar wind velocity with a lag of two to three days. However, above 350 keV, the SEP fluxes, accounts for a significant portion of the GEO proton flux variance. For the highest three energy channels (0.74 to 14.5 MeV), the SEPs account for the majority of the ERR. The results also show an anisotropy of protons with gyrocenters inside GEO and outside GEO, where the protons inside GEO are controlled partly by the Dst index and also an IMF-clock angle function.

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

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

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

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

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

Description: A statistical study on O + flux in the dayside magnetosheath Annales Geophysicae, 31, 1005-1010, 2013 Author(s): R. Slapak, H. Nilsson, and L. G. Westerberg Studies on terrestrial oxygen ion (O + ) escape into the interplanetary space have considered a number of different escape paths. Recent observations however suggest a yet insufficiently investigated additional escape route for hot O + : along open magnetic field lines in the high altitude cusp and mantle. Here we present a statistical study on O + flux in the high-latitude dayside magnetosheath. The O + is generally seen relatively close to the magnetopause, consistent with observations of O + flowing primarily tangentially to the magnetopause. We estimate the total escape flux in this region to be ~ 7 × 10 24 s −1 , implying this escape route to significantly contribute to the overall total O + loss into interplanetary space.

Description: Simulation of the influence of historical land cover changes on the global climate Annales Geophysicae, 31, 995-1004, 2013 Author(s): Y. Wang, X. Yan, and Z. Wang In order to estimate biogeophysical effects of historical land cover change on climate during last three centuries, a set of experiments with a climate system model of intermediate complexity (MPM-2) is performed. In response to historical deforestation, the model simulates a decrease in annual mean global temperature in the range of 0.07–0.14 °C based on different grassland albedos. The effect of land cover changes is most pronounced in the middle northern latitudes with maximum cooling reaching approximately 0.6 °C during northern summer. The cooling reaches 0.57 °C during northern spring owing to the large effects of land surface albedo. These results suggest that land cover forcing is important for study on historical climate change and that more research is necessary in the assessment of land management options for climate change mitigation.

Description: IMF effect on the polar cap contraction and expansion during a period of substorms Annales Geophysicae, 31, 1021-1034, 2013 Author(s): A. T. Aikio, T. Pitkänen, I. Honkonen, M. Palmroth, and O. Amm The polar cap boundary (PCB) location and motion in the nightside ionosphere has been studied by using measurements from the EISCAT radars and the MIRACLE magnetometers during a period of four substorms on 18 February 2004. The OMNI database has been used for observations of the solar wind and the Geotail satellite for magnetospheric measurements. In addition, the event was modelled by the GUMICS-4 MHD simulation. The simulation of the PCB location was in a rather good agreement with the experimental estimates at the EISCAT longitude. During the first three substorm expansion phases, neither the local observations nor the global simulation showed any poleward motions of the PCB, even though the electrojets intensified. Rapid poleward motions of the PCB took place only in the early recovery phases of the substorms. Hence, in these cases the nightside reconnection rate was locally higher in the recovery phase than in the expansion phase. In addition, we suggest that the IMF B z component correlated with the nightside tail inclination angle and the PCB location with about a 17-min delay from the bow shock. By taking the delay into account, the IMF northward turnings were associated with dipolarizations of the magnetotail and poleward motions of the PCB in the recovery phase. The mechanism behind this effect should be studied further.

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

Description: A comparison of bow shock models with Cluster observations during low Alfvén Mach number magnetic clouds Annales Geophysicae, 31, 1011-1019, 2013 Author(s): L. Turc, D. Fontaine, P. Savoini, H. Hietala, and E. K. J. Kilpua Magnetic clouds (MCs) are very geoeffective solar wind structures. Their properties in the interplanetary medium have been extensively studied, yet little is known about their characteristics in the Earth's magnetosheath. The Cluster spacecraft offer the opportunity to observe MCs in the magnetosheath, but before MCs reach the magnetosphere, their structure is altered when they interact with the terrestrial bow shock (BS). The physics taking place at the BS strongly depends on Θ Bn , the angle between the shock normal and the interplanetary magnetic field. However, in situ observations of the BS during an MC's crossing are seldom available. In order to relate magnetosheath observations to solar wind conditions, we need to rely on a model to determine the shock's position and normal direction. Yet during MCs, the models tend to be less accurate, because the Alfvén Mach number ( M A ) is often significantly lower than in regular solar wind. On the contrary, the models are generally optimised for high M A conditions. In this study, we compare the predictions of four widely used models available in the literature (Wu et al., 2000; Chapman and Cairns, 2003; Jeřáb et al., 2005; Měrka et al., 2005b) to Cluster's dayside BS crossings observed during five MC events. Our analysis shows that the Θ Bn angle is well predicted by all four models. On the other hand, the Jeřáb et al. (2005) model yields the best estimates of the BS position during low M A MCs. The other models locate the BS either too far from or too close to Earth. The results of this paper can be directly used to estimate the BS parameters in all studies of MC interaction with Earth's magnetosphere.

Description: Investigation of radiative effects of the optically thick dust layer over the Indian tropical region Annales Geophysicae, 31, 647-663, 2013 Author(s): S. K. Das, J.-P. Chen, M. Venkat Ratnam, and A. Jayaraman Optical and physical properties of aerosols derived from multi-satellite observations (MODIS-Aqua, OMI-Aura, MISR-Terra, CALIOP-CALIPSO) have been used to estimate radiative effects of the dust layer over southern India. The vertical distribution of aerosol radiative forcing and heating rates are calculated with 100 m resolution in the lower atmosphere, using temperature and relative humidity data from balloon-borne radiosonde observations. The present study investigates the optically thick dust layer of optical thickness 0.18 ± 0.06 at an altitude of 2.5 ± 0.7 km over Gadanki, transported from the Thar Desert, producing radiative forcing and heating rate of 11.5 ± 3.3 W m −2 and 0.6 ± 0.26 K day −1 , respectively, with a forcing efficiency of 43 W m −2 and an effective heating rate of 4 K day −1 per unit dust optical depth. Presence of the dust layer increases radiative forcing by 60% and heating rate by 60 times at that altitude compared to non-dusty cloud-free days. Calculation shows that the radiative effects of the dust layer strongly depend on the boundary layer aerosol type and mass loading. An increase of 25% of heating by the dust layer is found over relatively cleaner regions than urban regions in southern India and further 15% of heating increases over the marine region. Such heating differences in free troposphere may have significant consequences in the atmospheric circulation and hydrological cycle over the tropical Indian region.

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

Description: Low level jet intensification by mineral dust aerosols Annales Geophysicae, 31, 625-632, 2013 Author(s): O. Alizadeh Choobari, P. Zawar-Reza, and A. Sturman Modification of the intensity of a low level jet (LLJ) and near-surface wind speed by mineral dust is important as it has implications for dust emission and its long-range transport. Using the Weather Research and Forecasting with Chemistry (WRF/Chem) regional model, it is shown that direct radiative forcing by mineral dust reduces temperature in the lower atmosphere, but increases it in the layers aloft. The surface cooling is shown to be associated with a reduction of turbulent kinetic energy (TKE) and hence vertical mixing of horizontal momentum. Changes in the vertical profile of temperature over the regions that are under the influence of a LLJ are shown to result in an intensification of the LLJ and near-surface wind speed, but a decrease of winds aloft. These changes in the wind speed profile differ from results of previous research which suggested a decrease of wind speed in the lower atmosphere and its increase in the upper boundary layer.

Description: Parametric validations of analytical lifetime estimates for radiation belt electron diffusion by whistler waves Annales Geophysicae, 31, 599-624, 2013 Author(s): A. V. Artemyev, D. Mourenas, O. V. Agapitov, and V. V. Krasnoselskikh The lifetimes of electrons trapped in Earth's radiation belts can be calculated from quasi-linear pitch-angle diffusion by whistler-mode waves, provided that their frequency spectrum is broad enough and/or their average amplitude is not too large. Extensive comparisons between improved analytical lifetime estimates and full numerical calculations have been performed in a broad parameter range representative of a large part of the magnetosphere from L ~ 2 to 6. The effects of observed very oblique whistler waves are taken into account in both numerical and analytical calculations. Analytical lifetimes (and pitch-angle diffusion coefficients) are found to be in good agreement with full numerical calculations based on CRRES and Cluster hiss and lightning-generated wave measurements inside the plasmasphere and Cluster lower-band chorus waves measurements in the outer belt for electron energies ranging from 100 keV to 5 MeV. Comparisons with lifetimes recently obtained from electron flux measurements on SAMPEX, SCATHA, SAC-C and DEMETER also show reasonable agreement.

Description: Coherent amplitude modulation of electron-beam-driven Langmuir waves Annales Geophysicae, 31, 633-638, 2013 Author(s): K. Baumgärtel A linear approach to the phenomenon of irregular amplitude modulation of beam-driven Langmuir waves, developed in a previous paper, is extended to explain periodic modulation as well. It comes about by beating of the fastest growing mode of the instability with beam-aligned plasma oscillations. They are naturally generated in a uniform domain of beam–plasma interaction prior to the onset of the instability. Particle-in-cell (PIC) simulations support the results of the linear analysis.

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

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

Description: Seasonal variation of the ion upflow in the topside ionosphere during SAPS (subauroral polarization stream) periods Annales Geophysicae, 31, 1521-1534, 2013 Author(s): H. Wang and H. Lühr A statistical study has been performed by using two years of DMSP (Defense Meteorological Satellite Program) plasma observations to investigate the seasonal effect of SAPS (subauroral polarization stream) on the ion upflow in the duskside ionosphere of the Northern Hemisphere. There are obvious upflows occurring in the topside ionosphere around the SAPS region, exceeding 200 m s −1 at winter solstice, indicating an important relationship between SAPS and the local plasma upward motion. Both SAPS and ion upward velocities show similar seasonal variations, largest in winter and smallest in summer, irrespective of geomagnetic activity. A good correlation is found and a linear relationship is derived between SAPS and the ion upflow velocities. During December solstice the average upflow flux can reach about 2 × 10 8 cm −2 s −1 for more disturbed periods, which is comparable to the typical upflow flux in the dayside cusp region. The depression of the ion temperatures around the peak SAPS region can be understood in terms of the adiabatic cooling. The hot ion cools down when expanding into the low ion concentration region. The electron temperature elevates around the SAPS region because of the reduced Coulomb cooling in the low ion density region. Both the changes of ion and electron temperatures are larger in winter than in summer, however, for Kp 〈 4 the electron temperatures are almost seasonably independent. The present work highlights the important role of the SAPS-related frictional heating at mid-latitudes on the local formation of the strong upward flow, which might provide a direct ionospheric ion source for the ring current and plasmasphere in the duskside sector.

Description: Magnetic field generation in a jet-sheath plasma via the kinetic Kelvin-Helmholtz instability Annales Geophysicae, 31, 1535-1541, 2013 Author(s): K.-I. Nishikawa, P. Hardee, B. Zhang, I. Duţan, M. Medvedev, E. J. Choi, K. W. Min, J. Niemiec, Y. Mizuno, A. Nordlund, J. T. Frederiksen, H. Sol, M. Pohl, and D. H. Hartmann We have investigated the generation of magnetic fields associated with velocity shear between an unmagnetized relativistic jet and an unmagnetized sheath plasma. We have examined the strong magnetic fields generated by kinetic shear (Kelvin–Helmholtz) instabilities. Compared to the previous studies using counter-streaming performed by Alves et al. (2012), the structure of the kinetic Kelvin–Helmholtz instability (KKHI) of our jet-sheath configuration is slightly different, even for the global evolution of the strong transverse magnetic field. In our simulations the major components of growing modes are the electric field E z , perpendicular to the flow boundary, and the magnetic field B y , transverse to the flow direction. After the B y component is excited, an induced electric field E x , parallel to the flow direction, becomes significant. However, other field components remain small. We find that the structure and growth rate of KKHI with mass ratios m i / m e = 1836 and m i / m e = 20 are similar. In our simulations saturation in the nonlinear stage is not as clear as in counter-streaming cases. The growth rate for a mildly-relativistic jet case (γ j = 1.5) is larger than for a relativistic jet case (γ j = 15).

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

Description: Ionospheric response to total solar eclipse of 22 July 2009 in different Indian regions Annales Geophysicae, 31, 1549-1558, 2013 Author(s): S. Kumar, A. K. Singh, and R. P. Singh The variability of ionospheric response to the total solar eclipse of 22 July 2009 has been studied analyzing the GPS data recorded at the four Indian low-latitude stations Varanasi (100% obscuration), Kanpur (95% obscuration), Hyderabad (84% obscuration) and Bangalore (72% obscuration). The retrieved ionospheric vertical total electron content (VTEC) shows a significant reduction (reflected by all PRNs (satellites) at all stations) with a maximum of 48% at Varanasi (PRN 14), which decreases to 30% at Bangalore (PRN 14). Data from PRN 31 show a maximum of 54% at Kanpur and 26% at Hyderabad. The maximum decrement in VTEC occurs some time (2–15 min) after the maximum obscuration. The reduction in VTEC compared to the quiet mean VTEC depends on latitude as well as longitude, which also depends on the location of the satellite with respect to the solar eclipse path. The amount of reduction in VTEC decreases as the present obscuration decreases, which is directly related to the electron production by the photoionization process. The analysis of electron density height profile derived from the COSMIC (Constellation Observing System for Meteorology, Ionosphere & Climate) satellite over the Indian region shows significant reduction from 100 km altitude up to 800 km altitude with a maximum of 48% at 360 km altitude. The oscillatory nature in total electron content data at all stations is observed with different wave periods lying between 40 and 120 min, which are attributed to gravity wave effects generated in the lower atmosphere during the total solar eclipse.

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

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

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

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

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

Description: Application of nonlinear autoregressive moving average exogenous input models to geospace: advances in understanding and space weather forecasts Annales Geophysicae, 31, 1579-1589, 2013 Author(s): R. J. Boynton, M. A. Balikhin, S. A. Billings, and O. A. Amariutei The nonlinear autoregressive moving average with exogenous inputs (NARMAX) system identification technique is applied to various aspects of the magnetospheres dynamics. It is shown, from an example system, how the inputs to a system can be found from the error reduction ratio (ERR) analysis, a key concept of the NARMAX approach. The application of the NARMAX approach to the Dst (disturbance storm time) index and the electron fluxes at geostationary Earth orbit (GEO) are reviewed, revealing new insight into the physics of the system. The review of studies into the Dst index illustrate how the NARMAX approach is able to find a coupling function for the Dst index from data, which was then analytically justified from first principles. While the review of the electron flux demonstrates how NARMAX is able to reveal new insight into the physics of the acceleration and loss processes within the radiation belt.

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

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

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

Description: Energy exchange and wave action conservation for magnetohydrodynamic (MHD) waves in a general, slowly varying medium Annales Geophysicae, 32, 1495-1510, 2014 Author(s): A. D. M. Walker Magnetohydrodynamic (MHD) waves in the solar wind and magnetosphere are propagated in a medium whose velocity is comparable to or greater than the wave velocity and which varies in both space and time. In the approximation where the scales of the time and space variation are long compared with the period and wavelength, the ray-tracing equations can be generalized and then include an additional first-order differential equation that determines the variation of frequency. In such circumstances the wave can exchange energy with the background: wave energy is not conserved. In such processes the wave action theorem shows that the wave action, defined as the ratio of the wave energy to the frequency in the local rest frame, is conserved. In this paper we discuss ray-tracing techniques and the energy exchange relation for MHD waves. We then provide a unified account of how to deal with energy transport by MHD waves in non-uniform media. The wave action theorem is derived directly from the basic MHD equations for sound waves, transverse Alfvén waves, and the fast and slow magnetosonic waves. The techniques described are applied to a number of illustrative cases. These include a sound wave in a medium undergoing a uniform compression, an isotropic Alfvén wave in a steady-state shear layer, and a transverse Alfvén wave in a simple model of the magnetotail undergoing compression. In each case the nature and magnitude of the energy exchange between wave and background is found.

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

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

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

Description: On the influence of CMEs on the global 3-D coronal electron density Annales Geophysicae, 29, 1019-1028, 2011 Author(s): M. Kramar, J. Davila, H. Xie, and S. Antiochos In order to analize the influence of a Coronal Mass Ejection (CME) on the coronal streamer belt, we made 3-D reconstructions of the electron density in the corona at heliospheric distances from 1.5 to 4 R ⊙ for periods before and after a CME occured. The reconstructions were performed using a tomography technique. We studied two CME cases: (i) a slow CME on 1 June 2008; (ii) two fast CMEs on 31 December 2007 and 2 January 2008. For the first case of slow CME, it was found: (i) the potential magnetic field configuration in the CME initiation region before the CME does not agree with the coronal density structure while after the CME the agreement between the field and density is much better. This could be manifistation of that that the field was non-potential before the CME and after the CME the field relaxes towards a more potential state. (ii) It was shown that the dimming caused by the slow CME is not due to rotation of the corona and a line-of-sight (LOS) effect but a streamer blow out effect took place.

Description: Longitudinal study of the ionospheric response to the geomagnetic storm of 15 May 2005 and manifestation of TADs Annales Geophysicae, 29, 1063-1070, 2011 Author(s): S. Sharma, P. Galav, N. Dashora, and R. Pandey Response of low latitude ionosphere to the geomagnetic storm of 15 May 2005 has been studied using total electron content (TEC) data, obtained from three GPS stations namely, Yibal, Udaipur and Kunming situated near the northern crest of equatorial ionization anomaly at different longitudes. Solar wind parameters, north-south component of the interplanetary magnetic field (IMF B z ) and AE index data have been used to infer the strength of the geomagnetic storm. A large value of eastward interplanetary electric field at 06:15 UT, during the time of maximum southward IMF B z has been used to infer the transmission of an eastward prompt penetration electric field (PPEF) which resulted in a peak in TEC at 07:45 UT due to the local uplift of plasma in the low latitudes near the anomaly crest over a wide range of longitudes. Wave-like modulations superposed over the second enhancement in TEC between 09:15 UT to 10:30 UT have been observed at all the three stations. The second enhancement in TEC along with the modulations of up to 5 TECU have been attributed to the combined effect of super plasma fountain and traveling atmospheric disturbances (TAD). Observed large enhancements in TEC are a cause of concern for satellite based navigation and ground positioning. Increased [O/N 2 ] ratio between 09:15 UT to 10:15 UT when modulations in TEC have been also observed, confirms the presence of TADs over a wide range of longitudes.

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

Description: Volume cross section of auroral radar backscatter and RMS plasma fluctuations inferred from coherent and incoherent scatter data: a response on backscatter volume parameters Annales Geophysicae, 29, 1081-1092, 2011 Author(s): M. V. Uspensky, P. Janhunen, A. V. Koustov, and K. Kauristie Norway and Finland STARE radar measurements in the eastward auroral electrojet are combined with EISCAT CP-1 measurements of the electron density and electric field vector in the common scattering volume to investigate the variation of the auroral radar volume cross section (VCS) with the flow angle of observations (radar look direction with respect to the E × B electron drift). The data set available consists of ~6000 points for flow angles of 40–85° and electron drifts between 500 and 2000 m s −1 . The EISCAT electron density N(h) -profile data are used to estimate the effective electron density, aspect angle and thickness of the backscattering layer. It is shown that the flow angle variation of the VCS is rather weak, only ~5 dB within the range of the considered flow angles. The VCS values themselves respond almost linearly to the square of both the electron drift velocity magnitude and the effective electron density. By adopting the inferred shape of the VCS variation with the flow angle and the VCS dependence upon wavelength, the relative amplitude of electrostatic electron density fluctuations over all scales is estimated. Inferred values of 2–4 percent react nearly linearly to the electron drift velocity in the range of 500–1000 m s −1 but the rate of increase slows down at electron drifts 〉1000 m s −1 and density fluctuations of ~5.5 percent due to, perhaps, progressively growing nonlinear wave losses.

Description: Flux quanta, magnetic field lines, merging – some sub-microscale relations of interest in space plasma physics Annales Geophysicae, 29, 1121-1127, 2011 Author(s): R. A. Treumann, R. Nakamura, and W. Baumjohann We clarify the notion of magnetic field lines in plasma by referring to sub-microscale (quantum mechanical) particle dynamics. It is demonstrated that magnetic field lines in a field of strength B carry single magnetic flux quanta Φ 0 = h/e . The radius of a field line in the given magnetic field B is calculated. It is shown that such field lines can merge and annihilate only over the length ℓ ∥ of their strictly anti-parallel sections, for which case we estimate the power generated. The length ℓ ∥ becomes a function of the inclination angle θ of the two merging magnetic flux tubes (field lines). Merging is possible only in the interval 1 2 πθ≤π. This provides a sub-microscopic basis for "component reconnection" in classical macro-scale reconnection. We also find that the magnetic diffusion coefficient in plasma appears in quanta D 0 m = e Φ 0 / m e = h/m e . This lets us conclude that the bulk perpendicular plasma resistivity is limited and cannot be less than η 0⊥ =μ 0 e Φ 0 / m e =μ 0 h/m e ~10 −9 Ohm m. This resistance is an invariant.

Description: Could periodic patterns in human mortality be sensitive to solar activity? Annales Geophysicae, 29, 1113-1120, 2011 Author(s): R. Díaz-Sandoval, R. Erdélyi, and R. Maheswaran Seasonal behaviour of human diseases have been observed and reported in the literature for years. Although the Sun plays an essential role in the origin and evolution of life on Earth, it is barely taken into account in biological processes for the development of a specific disease. Higher mortality rates occur during the winter season in the Northern Hemisphere for several diseases, particularly diseases of the cardiovascular and respiratory systems. This increment has been associated with seasonal and social causes. However, is there more behind these correlations, in particular in terms of solar variability? In this paper we attempt to make a first step towards answering this question. A detailed wavelet analysis of periodicities for diseases from England and Wales seem to reveal that mortality periodicities (3 days to half a year) could be due to the Earth's position around the Sun. Moreover, crosswavelet and wavelet coherence analysis show common features between medical diseases and solar proxies around solar maximum activity suggesting that this relation, if any, has to be searched in times of high solar activity.

Description: Excitation of zero-frequency magnetic field-aligned currents by ionospheric heating Annales Geophysicae, 29, 1147-1152, 2011 Author(s): A. V. Streltsov and T. R. Pedersen Time-dependent, three-dimensional numerical simulations of the reduced MHD model describing shear Alfvén waves in the magnetosphere provide an interesting prediction superficially similar to results of several ionospheric heating experiments conducted at high altitudes. In these experiments, heating of the ionospheric F-region with a constant/zero-frequency beam of HF waves causes luminous structures in the ionosphere in the form of a ring or a solid spot with a characteristic size comparable to the size of the heated spot. Simulations suggest that spots/rings or similar optical appearance might be associated with a magnetic field-aligned current system produced by the ionospheric heating. Two of the most interesting features of this current system are (1) strong localization across the ambient magnetic field and (2) distinctive non-symmetrical luminous signatures (ring/spot) in magnetically conjugate locations in the ionosphere.

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

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

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

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

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

Description: Comment on "The PC index: review of methods", by McCreadie and Menvielle (2010) Annales Geophysicae, 29, 1137-1146, 2011 Author(s): P. Stauning The Polar Cap (PC) index is a controversial topic within the IAGA scientific community. Since 1997 discussions of the validity of the index to be endorsed as an official IAGA index have ensued. The article: "The PC index: review of methods", written by two members of the IAGA PC index committee, H. McCreadie and M. Menvielle, holds a critical review of some aspects of the methods used to derive PC index values. However, a number of incorrect statements and factual errors have been found and shall be called attention to and discussed in this commentary. Further critical comments concern the Corrigendum issued by the same authors and published in Ann. Geophys., 29, 813–814, 2011.

Description: The role of the bow shock in solar wind-magnetosphere coupling Annales Geophysicae, 29, 1129-1135, 2011 Author(s): R. E. Lopez, V. G. Merkin, and J. G. Lyon In this paper we examine the role of the bow shock in coupling solar wind energy to the magnetosphere using global magnetohydrodynamic simulations of the solar wind-magnetosphere interaction with southward IMF. During typical solar wind conditions, there are two significant dynamo currents in the magnetospheric system, one in the high-latitude mantle region tailward of the cusp and the other in the bow shock. As the magnitude of the (southward) IMF increases and the solar wind becomes a low Mach number flow, there is a significant change in solar wind-magnetosphere coupling. The high-latitude magnetopause dynamo becomes insignificant compared to the bow shock and a large load appears right outside the magnetopause. This leaves the bow shock current as the only substantial dynamo current in the system, and the only place where a significant amount of mechanical energy is extracted from the solar wind. That energy appears primarily as electromagnetic energy, and the Poynting flux generated at the bow shock feeds energy back into the plasma, reaccelerating it to solar wind speeds. Some small fraction of that Poynting flux is directed into the magnetosphere, supplying the energy needed for magnetospheric dynamics. Thus during periods when the solar wind flow has a low Mach number, the main dynamo in the solar wind-magnetosphere system is the bow shock.

Description: Measurements of aerosol intensive properties over Visakhapatnam, India for 2007 Annales Geophysicae, 29, 973-985, 2011 Author(s): K. Niranjan, B. Spandana, T. Anjana Devi, V. Sreekanth, and B. L. Madhavan Information on the aerosol intensive properties like Single Scattering Albedo (SSA) and asymmetry parameter are very limited, particularly over the peninsular India, though extensive reports are available on the aerosol bulk properties. In view of the importance of these parameters in evaluating the aerosol radiative forcing, we present for the first time the temporal variation in SSA with measurements on aerosol absorption and scattering coefficients over Visakhapatnam (17.72° N, 83.32° E; located on the east coast of India) for the year 2007. The inferred SSA was in the range of 0.65 and 0.9 with an annual mean of 0.76 ± 0.013 and with a probable value of 0.80, indicating a marginal atmospheric warming over the region. The mixed layer contribution to column Aerosol Optical depth is found to be 35 % in summer while it is well above 35 % in winter, indicating the confinement of aerosol within the boundary layer during winter. The asymmetry parameter which represents the angular scattering in radiative forcing estimation is found to be around 0.65 ± 0.1 for the location. The day to day variability in SSA is found to be well correlated with the variations in surface BC mass concentrations and/or the relative dominance of the fine/coarse mode aerosol. The results are discussed in light of the aerosol physical and optical properties and the asymmetry parameter.

Description: Corrigendum to "Downward auroral currents from the reconnection Hall-region", published in Ann. Geophys., 29, 679–685, 2011 Annales Geophysicae, 29, 1061-1061, 2011 Author(s): R. A. Treumann, R. Nakamura, and W. Baumjohann No abstract available.

Description: Hybrid models of solar wind plasma heating Annales Geophysicae, 29, 1071-1079, 2011 Author(s): L. Ofman, A.-F. Viñas, and P. S. Moya Remote sensing and in-situ observations show that solar wind ions are often hotter than electrons, and the heavy ions flow faster than the protons by up to an Alfvén speed. Turbulent spectrum of Alfvénic fluctuations and shocks were detected in solar wind plasma. Cross-field inhomogeneities in the corona were observed to extend to several tens of solar radii from the Sun. The acceleration and heating of solar wind plasma is studied via 1-D and 2-D hybrid simulations. The models describe the kinetics of protons and heavy ions, and electrons are treated as neutralizing fluid.The expansion of the solar wind is considered in 1-D hybrid model. A spectrum of Alfvénic fluctuations is injected at the computational boundary, produced by differential streaming instability, or initial ion temperature anisotropy, and the parametric dependence of the perpendicular heating of H + -He ++ solar wind plasma is studied. It is found that He ++ ions are heated efficiently by the Alfvénic wave spectrum below the proton gyroperiod.

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

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

Description: On the effect of IMF turning on ion dynamics at Mercury Annales Geophysicae, 29, 987-996, 2011 Author(s): D. C. Delcourt, T. E. Moore, and M.-C. H. Fok We investigate the effect of a rotation of the Interplanetary Magnetic Field (IMF) on the transport of magnetospheric ion populations at Mercury. We focus on ions of planetary origin and investigate their large-scale circulation using three-dimensional single-particle simulations. We show that a nonzero B X component of the IMF leads to a pronounced asymmetry in the overall circulation pattern. In particular, we demonstrate that the centrifugal acceleration due to curvature of the E × B drift paths is more pronounced in one hemisphere than the other, leading to filling of the magnetospheric lobes and plasma sheet with more or less energetic material depending upon the hemisphere of origin. Using a time-varying electric and magnetic field model, we investigate the response of ions to rapid (a few tens of seconds) re-orientation of the IMF. We show that, for ions with gyroperiods comparable to the field variation time scale, the inductive electric field should lead to significant nonadiabatic energization, up to several hundreds of eVs or a few keVs. It thus appears that IMF turning at Mercury should lead to localized loading of the magnetosphere with energetic material of planetary origin (e.g., Na + ).

Description: Instability of combined gravity-inertial-Rossby waves in atmospheres and oceans Annales Geophysicae, 29, 997-1003, 2011 Author(s): J. F. McKenzie The properties of the instability of combined gravity-inertial-Rossby waves on a β-plane are investigated. The wave-energy exchange equation shows that there is an exchange of energy with the background stratified medium. The energy source driving the instability lies in the background enthalpy released by the gravitational buoyancy force. It is shown that if the phase speed of the westward propagating low frequency-long wavelength Rossby wave exceeds the Poincaré-Kelvin (or "equivalent" shallow water) wave speed, instability arises from the merging of Rossby and Poincaré modes. There are two key parameters in this instability condition; namely, the equatorial/rotational Mach (or Froude) number M and the latitude θ 0 of the β-plane. In general waves equatorward of a critical latitude for given M can be driven unstable, with corresponding growth rates of the order of a day or so. Although these conclusions may only be safely drawn for short wavelengths corresponding to a JWKB wave packet propagating internally and located far from boundaries, nevertheless such a local instability may play a significant role in atmosphere-ocean dynamics.

Description: The correlation between solar and geomagnetic activity – Part 3: An integral response model Annales Geophysicae, 29, 1005-1018, 2011 Author(s): Z. L. Du An integral response model is proposed to describe the relationship between geomagnetic activity ( aa index) and solar activity (represented by sunspot number R z ): The aa at a given time t is the integral of R z at past times ( t '≤ t ) multiplied by an exponential decay factor of the time differences ( e −(t−t')/τ ), where τ is the decay time scale (~40 months). The correlation coefficient of aa with the reconstructed series based on this model ( r f =0.85) is much higher than that of aa with R z ( r 0 =0.61). If this model is applied to each solar cycle, the correlation coefficient will be higher ( r f =0.95). This model can naturally explain some phenomena related to aa and R z , such as (i) the significant increase in the aa index (and its baseline) over the twentieth century; (ii) the longer lag times of aa to R z at solar cycle maxima than at minima; and (iii) the variations in the correlations related to solar and Hale cycles. These results demonstrate that aa depends not only on the present R z but also on past values. The profile of aa can be better predicted from R z by this model than by point-point correspondence.

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

Description: Cluster observations of trapped ions interacting with magnetosheath mirror modes Annales Geophysicae, 29, 1049-1060, 2011 Author(s): J. Soucek and C. P. Escoubet Mirror modes are among the most intense low frequency plasma wave phenomena observed in the magnetosheaths of magnetized planets. They appear as large amplitude non-propagating fluctuations in the magnetic field magnitude and plasma density. These structures are widely accepted to represent a non-linear stage of the mirror instability, dominant in plasmas with large ion beta and a significant ion temperature anisotropy T ⊥ / T ∥ 〉1. It has long been recognized that the mirror instability both in the linear and non-linear stage is a kinetic process and that the behavior of resonant particles at small parallel velocities is crucial for its development and saturation. While the dynamics of the instability and the effect of trapped particles have been studied extensively in theoretical models and numerical simulations, only spurious observations of the trapped ions were published to date. In this work we used data from the Cluster spacecraft to perform the first detailed experimental study of ion velocity distribution associated with mirror mode oscillations. We show a conclusive evidence for the predicted cooling of resonant ions at small parallel velocities and heating of trapped ions at intermediate pitch angles.

Description: MHD waves generated by high-frequency photospheric vortex motions Annales Geophysicae, 29, 1029-1035, 2011 Author(s): V. Fedun, S. Shelyag, G. Verth, M. Mathioudakis, and R. Erdélyi In this paper, we discuss simulations of MHD wave generation and propagation through a three-dimensional open magnetic flux tube in the lower solar atmosphere. By using self-similar analytical solutions for modelling the magnetic field in Cartesian coordinate system, we have constructed a 3-D magnetohydrostatic configuration which is used as the initial condition for non-linear MHD wave simulations. For a driver we have implemented a high-frequency vortex-type motion at the footpoint region of the open magnetic flux tube. It is found that the implemented swirly source is able to excite different types of wave modes, i.e. sausage, kink and torsional Alfvén modes. Analysing these waves by magneto-seismology tools could provide insight into the magnetic structure of the lower solar atmosphere.

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

Description: Average properties and small-scale variations of the mesospheric Na and Fe layers as observed simultaneously by two closely colocated lidars at 30° N Annales Geophysicae, 29, 1037-1048, 2011 Author(s): L. Chen and F. Yi We report the average properties and small-scale variation features of the mesospheric Na and Fe layers at 30° N from extensive simultaneous and common-volume Na and Fe lidar measurements at Wuhan, China. The annual mean Na and Fe density profiles are derived in terms of an averaging method taken from an early literature. The mean Na and Fe profiles preserve the sharp gradients present in most individual density profiles near the layer bottom. Near the bottommost of the layers the mean Na and Fe scale heights are respectively −0.42 and −0.30 km. The mean layer parameters coincide well with the previous report. The Na and Fe densities in the lowest several kilometers of the layers consistently exhibit nearly the same time variations. A clear-cut distinction between the Na and Fe time variations always appears in an altitude range near 90 km. A relatively weak positive correlation between them persistently occurs also in an altitude range near 100 km. The mean increase and decrease rates for both Na and Fe are altitude dependent and have a single-peak structure. The time constant of the layer variation is ~0.07–2.0 h for Na and ~0.02–1.7 h for Fe, suggesting that the variability is dominated by small-scale processes. However, there is also a slow net increase in each of the annual mean column abundances (Na and Fe) during night.

Description: Polarisation in the auroral red line during coordinated EISCAT Svalbard Radar/optical experiments Annales Geophysicae, 29, 1101-1112, 2011 Author(s): M. Barthélémy, J. Lilensten, F. Pitout, C. Simon Wedlund, R. Thissen, D. Lorentzen, F. Sigernes, J. Moen, G. Gronoff, I. McCrea, H. Rothkael, H. Ménager, and A. Aruliah The polarisation of the atomic oxygen red line in the Earth's thermosphere is observed in different configurations with respect to the magnetic field line at high latitude during several coordinated Incoherent Scatter radar/optical experiment campaigns. When pointing northward with a line-of-sight nearly perpendicular to the magnetic field, we show that, as expected, the polarisation is due to precipitated electrons with characteristic energies of a few hundreds of electron Volts. When pointing toward the zenith or southward with a line-of-sight more parallel to the magnetic field, we show that the polarisation practically disappears. This confirms experimentally the predictions deduced from the recent discovery of the red line polarisation. We show that the polarisation direction is parallel to the magnetic field line during geomagnetic activity intensification and that these results are in agreement with theoretical work.

Description: New insights into the stratospheric and mesosphere-lower thermospheric ozone response to the abrupt changes in solar forcing Annales Geophysicae, 29, 1093-1099, 2011 Author(s): Karanam Kishore Kumar, K. V. Subrahmanyam, and Sherine Rachel John Using a unique set of satellite based observations of the vertical distribution of ozone during the recent annular solar eclipse of 15 January 2010, we demonstrate for the first time, a complete picture of the response of stratospheric ozone to abrupt changes in solar forcing. The stratospheric ozone decreased after the maximum obscuration of the Sun and then gradually increased with time. A dramatic increase in stratospheric ozone of up to 4 ppmv is observed 3 h after the maximum obscuration of the Sun. The present study also reports for the first time the mesosphere-lower thermospheric ozone response to solar eclipse. Thus it is envisaged that the present results will have important implications in understanding the ozone response to abrupt changes in solar forcing and time-scales involved in such response.

Description: Corrigendum to "The PC index: review of methods" published in Ann. Geophys., 28, 1887–1903, 2010 Annales Geophysicae, 29, 813-814, 2011 Author(s): H. McCreadie and M. Menvielle No abstract available.

Description: Pre-monsoon aerosol characteristics over the Indo-Gangetic Basin: implications to climatic impact Annales Geophysicae, 29, 789-804, 2011 Author(s): A. K. Srivastava, S. Tiwari, P. C. S. Devara, D. S. Bisht, Manoj K. Srivastava, S. N. Tripathi, P. Goloub, and B. N. Holben Sun/sky radiometer observations over the Indo-Gangetic Basin (IGB) region during pre-monsoon (from April–June 2009) have been processed to analyze various aerosol characteristics in the central and eastern IGB region, represented by Kanpur and Gandhi College, respectively, and their impacts on climate in terms of radiative forcing. Monthly mean aerosol optical depth (AOD at 500 nm) and corresponding Angstrom Exponent (AE at 440–870 nm, given within the brackets) was observed to be about 0.50 (0.49) and 0.51 (0.65) in April, 0.65 (0.74) and 0.67 (0.91) in May and 0.69 (0.45) and 0.77 (0.71) in June at Kanpur and Gandhi College, respectively. Results show a positive gradient in AOD and AE from central to eastern IGB region with the advancement of the pre-monsoon, which may be caused due to diverse geographical location of the stations having different meteorological conditions and emission sources. Relatively lower SSA was observed at the eastern IGB (0.89) than the central IGB (0.92) region during the period, which suggests relative dominance of absorbing aerosols at the eastern IGB as compared to central IGB region. The absorbing aerosol optical properties over the station suggest that the atmospheric absorption over central IGB region is mainly due to dominance of coarse-mode dust particles; however, absorption over eastern IGB region is mainly due to dominance of fine-particle pollution. The derived properties from sun/sky radiometer during pre-monsoon period are used in a radiative-transfer model to estimate aerosol radiative forcing at the top-of-the atmosphere (TOA) and at the surface over the IGB region. Relatively large TOA and surface cooling was observed at the eastern IGB as compared to the central IGB region. This translates into large heating of the atmosphere ranging from 0.45 to 0.55 K day −1 at Kanpur and from 0.45 to 0.59 K day −1 at Gandhi College.

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

Description: Annual and semiannual variations of vertical total electron content during high solar activity based on GPS observations Annales Geophysicae, 29, 865-873, 2011 Author(s): M. P. Natali and A. Meza Annual, semiannual and seasonal variations of the Vertical Total Electron Content (VTEC) have been investigated during high solar activity in 2000. In this work we use Global IGS VTEC maps and Principal Component Analysis to study spatial and temporal ionospheric variability. The behavior of VTEC variations at two-hour periods, at noon and at night is analyzed. Particular characteristics associated with each period and the geomagnetic regions are highlighted. The variations at night are smaller than those obtained at noon. At noon it is possible to see patterns of the seasonal variation at high latitude, and patterns of the semiannual anomaly at low latitudes with a slow decrease towards mid latitudes. At night there is no evidence of seasonal or annual anomaly for any region, but it was possible to see the semiannual anomaly at low latitudes with a sudden decrease towards mid latitudes. In general, the semiannual behavior shows March–April equinox at least 40 % higher than September one. Similarities and differences are analyzed also with regard to the same analysis done for a period of low solar activity.

Description: Modelling of spacecraft spin period during eclipse Annales Geophysicae, 29, 875-882, 2011 Author(s): E. Georgescu, F. Plaschke, U. Auster, K.-H. Fornaçon, and H. U. Frey The majority of scientific satellites investigating the Earth magnetosphere are spin stabilized. The attitude information comes usually from a sun sensor and is missing in the umbra; hence, the accurate experimental determination of vector quantities is not possible during eclipses. The spin period of the spacecraft is generally not constant during these times because the moment of inertia changes due to heat dissipation. The temperature dependence of the moment of inertia for each spacecraft has a specific signature determined by its design and distribution of mass. We developed an "eclipse-spin" model for the spacecraft spin period behaviour using magnetic field vector measurements close to the Earth, where the magnetic field is dominated by the dipole field, and in the magnetospheric lobes, where the magnetic field direction is mostly constant. The modelled spin periods give us extraordinarily good results with accumulated phase deviations over one hour of less than 10 degrees. Using the eclipse spin model satellite experiments depending on correct spin phase information can deliver science data even during eclipses. Two applications for THEMIS B, one in the lobe and the other in the lunar wake, are presented.

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

Description: Apparent temperature anisotropies due to wave activity in the solar wind Annales Geophysicae, 29, 909-917, 2011 Author(s): D. Verscharen and E. Marsch The fast solar wind is a collisionless plasma permeated by plasma waves on many different scales. A plasma wave represents the natural interplay between the periodic changes of the electromagnetic field and the associated coherent motions of the plasma particles. In this paper, a model velocity distribution function is derived for a plasma in a single, coherent, large-amplitude wave. This model allows one to study the kinetic effects of wave motions on particle distributions. They are by in-situ spacecraft measured by counting, over a certain sampling time, the particles coming from various directions and having different energies. We compare our results with the measurements by the Helios spacecraft, and thus find that by assuming high wave activity we are able to explain key observed features of the measured distributions within the framework of our model. We also address the recent discussions on nonresonant wave–particle interactions and apparent heating. The applied time-averaging procedure leads to an apparent ion temperature anisotropy which is connected but not identical to the intrinsic temperature of the underlying distribution function.

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

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

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

Description: The numerical simulation on ionospheric perturbations in electric field before large earthquakes Annales Geophysicae, 32, 1487-1493, 2014 Author(s): S. F. Zhao, X. M. Zhang, Z. Y. Zhao, and X. H. Shen Many observational results have shown electromagnetic abnormality in the ionosphere before large earthquakes. The theoretical simulation can help us to understand the internal mechanism of these anomalous electromagnetic signals resulted from seismic regions. In this paper, the horizontal and vertical components of electric and magnetic field at the topside ionosphere are simulated by using the full wave method that is based on an improved transfer matrix method in the lossy anisotropic horizontally stratified ionosphere. Taken account into two earthquakes with electric field perturbations recorded by the DEMETER satellite, the numerical results reveal that the propagation and penetration of ULF (ultra-low-frequency) electromagnetic waves into the ionosphere is related to the spatial distribution of electron and ion densities at different time and locations, in which the ion density has less effect than electron density on the field intensity. Compared with different frequency signals, the minimum values of electric and magnetic field excited by earthquakes can be detected by satellite in current detection capability have also been calculated, and the lower frequency wave can be detected easier.

Description: A comparison between VEGA 1, 2 and Giotto flybys of comet 1P/Halley: implications for Rosetta Annales Geophysicae, 32, 1441-1453, 2014 Author(s): M. Volwerk, K.-H. Glassmeier, M. Delva, D. Schmid, C. Koenders, I. Richter, and K. Szegö Three flybys of comet 1P/Halley, by VEGA 1, 2 and Giotto, are investigated with respect to the occurrence of mirror mode waves in the cometosheath and field line draping in the magnetic pile-up region around the nucleus. The time interval covered by these flybys is approximately 8 days, which is also the approximate length of an orbit or flyby of Rosetta around comet 67P/Churyumov–Gerasimenko. Thus any significant changes observed around Halley are changes that might occur for Rosetta during one pass of 67P/CG. It is found that the occurrence of mirror mode waves in the cometosheath is strongly influenced by the dynamical pressure of the solar wind and the outgassing rate of the comet. Field line draping happens in the magnetic pile-up region. Changes in nested draping regions (i.e. regions with different B x directions) can occur within a few days, possibly influenced by changes in the outgassing rate of the comet and thereby the conductivity of the cometary ionosphere.

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

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

Description: Magnetic curvature driven Rayleigh-Taylor instability revisited Annales Geophysicae, 29, 411-413, 2011 Author(s): O. A. Pokhotelov and O. G. Onishchenko The problem of incomplete finite ion Larmor radius (FLR) stabilization of the magnetic curvature driven Rayleigh-Taylor instability (RTI) in low beta plasma with homogeneous ion temperature is investigated. For this purpose a model hydrodynamic description of nonlinear flute waves with arbitrary spatial scales compared to the ion Larmor radius is developed. It is shown that the RTI is not stabilized by FLR effects in a plasma with cold electrons when the ratio of characteristic spatial scale of the plasma inhomogeneity to local effective radius of curvature of the magnetic field lines is larger than 1/4. The crucial role in the absence of the complete FLR stabilization plays the contribution of the compressibility of the polarization part of the ion velocity.

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

Description: Forecasting relativistic electron flux using dynamic multiple regression models Annales Geophysicae, 29, 415-420, 2011 Author(s): H.-L. Wei, S. A. Billings, A. Surjalal Sharma, S. Wing, R. J. Boynton, and S. N. Walker The forecast of high energy electron fluxes in the radiation belts is important because the exposure of modern spacecraft to high energy particles can result in significant damage to onboard systems. A comprehensive physical model of processes related to electron energisation that can be used for such a forecast has not yet been developed. In the present paper a systems identification approach is exploited to deduce a dynamic multiple regression model that can be used to predict the daily maximum of high energy electron fluxes at geosynchronous orbit from data. It is shown that the model developed provides reliable predictions.

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

Description: Airglow observations over the equatorial ionization anomaly zone in Taiwan Annales Geophysicae, 29, 749-757, 2011 Author(s): J. Y. Liu, P. K. Rajesh, I. T. Lee, and T. C. Chow Airglow imaging at mid-latitude stations often show intensity modulations associated with medium scale travelling ionospheric disturbances (MSTID), while those carried out near the equatorial regions reveal depletions caused by equatorial plasma bubbles (EPB). Two all sky cameras are used to observe plasma depletions in the 630.0 nm emission over the equatorial ionization anomaly (EIA) region, Taiwan (23° N, 121° E; 13.5° N Magnetic) during 1998–2002 and 2006–2007. The results show EPB and MSTID depletions in different solar activity conditions. Several new features of the EPB depletions such as bifurcation, secondary structure on the walls, westward tilt, etc., are discussed in this paper. Evidence of tilted depletions with secondary structures developing on the eastern wall that later evolve to appear as bifurcations, are presented for the first time. Moreover, detail investigations are carried out using International Reference Ionosphere (IRI) model as well as the electron density from Ionosonde and Global Positioning System (GPS) Occultation Experiment (GOX) onboard FORMOSAT-3/COSMIC satellite, to understand the conditions that favor the propagation of MSTID to the latitude of Taiwan.

Description: 3-D mesoscale MHD simulations of a cusp-like magnetic configuration: method and first results Annales Geophysicae, 29, 759-770, 2011 Author(s): E. Adamson and K. Nykyri We present a local mesoscale model of the magnetospheric cusp region with high resolution (up to 300 km). We discuss the construction and implementation of the initial configuration and give a detailed description of the numerical simulation. An overview of simulation results for the case of strongly northward interplanetary magnetic field (IMF) is then presented and compared with data from Cluster 2 spacecraft from 14 February 2003. Results show a cusp diamagnetic cavity (CDC) with depth normal to the magnetospheric boundary on the order of 1–2 R E and a much larger extent of ~5–9 R E tangential to the boundary, bounded by a gradual inner boundary with the magnetospheric lobe and a more distinct exterior boundary with the magnetosheath. These results are qualitatively consistent with observational data.

Description: Polar cap ion beams during periods of northward IMF: Cluster statistical results Annales Geophysicae, 29, 771-787, 2011 Author(s): R. Maggiolo, M. Echim, J. De Keyser, D. Fontaine, C. Jacquey, and I. Dandouras Above the polar caps and during prolonged periods of northward IMF, the Cluster satellites detect upward accelerated ion beams with energies up to a few keV. They are associated with converging electric field structures indicating that the acceleration is caused by a quasi-static field-aligned electric field that can extend to altitudes higher than 7 R E (Maggiolo et al., 2006; Teste et al., 2007). Using the AMDA science analysis service provided by the Centre de Données de la Physique des Plasmas, we have been able to extract about 200 events of accelerated upgoing ion beams above the polar caps from the Cluster database. Most of these observations are taken at altitudes lower than 7 R E and in the Northern Hemisphere. We investigate the statistical properties of these ion beams. We analyze their geometry, the properties of the plasma populations and of the electric field inside and around the beams, as well as their dependence on solar wind and IMF conditions. We show that ~40 % of the ion beams are collocated with a relatively hot and isotropic plasma population. The density and temperature of the isotropic population are highly variable but suggest that this plasma originates from the plasma sheet. The ion beam properties do not change significantly when the isotropic, hot background population is present. Furthermore, during one single polar cap crossing by Cluster it is possible to detect upgoing ion beams both with and without an accompanying isotropic component. The analysis of the variation of the IMF B Z component prior to the detection of the beams indicates that the delay between a northward/southward turning of IMF and the appearance/disappearance of the beams is respectively ~2 h and 20 min. The observed electrodynamic characteristics of high altitude polar cap ion beams suggest that they are closely connected to polar cap auroral arcs. We discuss the implications of these Cluster observations above the polar cap on the magnetospheric dynamics and configuration during prolonged periods of northward IMF.

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

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

Description: Diurnal variation of short-period (20–120 min) gravity waves in the equatorial Mesosphere and Lower Thermosphere and its relation to deep tropical convection Annales Geophysicae, 29, 623-629, 2011 Author(s): N. Venkateswara Rao, Y. Shibagaki, and T. Tsuda We study short period gravity waves (20–120 min) in the equatorial Mesosphere and Lower Thermosphere (MLT) using a Medium Frequency (MF) radar at Pameungpeuk (7.4° S, 107.4° E), Indonesia. In particular, we study local time and seasonal variation of the gravity wave variance and its relation to tropical convection. The gravity wave variance at 88 km enhances between 20:00 LT and 07:00 LT, with a peak at 02:00–03:00 LT. The enhancement is mainly observed during February–April and September–October and shows inter-annual variability. Convective activity over the same location persists from 16:00–21:00 LT with a peak activity ~18:00 LT and enhances between November–April. Time delay between the peak of convection and that of gravity wave activity ranges 1–15 h, which is consistent with theoretical calculations and previous reports based on reverse ray tracing analysis.

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

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

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