Chandrayaan-1 observations of backscattered solar wind protons from the lunar regolith: Dependence on the solar wind speed

We study the backscattering of solar wind protons from the lunar regolith using the Solar Wind Monitor of the Sub-keV Atom Reflecting Analyzer on Chandrayaan-1. Our study focuses on the component of the backscattered particles that leaves the regolith with a positive charge. We find that the fraction of the incident solar wind protons that backscatter as protons, i.e., the proton-backscattering efficiency, has an exponential dependence on the solar wind speed that varies from ~0.01% to ~1% for solar wind speeds of 250 km/s to 550 km/s. We also study the speed distribution of the backscattered protons in the fast (~550 km/s) solar wind case and find both a peak speed at ~80% of the solar wind speed and a spread of ~85 km/s. The observed flux variations and speed distribution of the backscattered protons can be explained by a speed-dependent charge state of the backscattered particles.

Variability of daily winter wind speed distribution over Northern Europe during the past millennium in regional and global climate simulations

We analyse the variability of the probability distribution of daily wind speed in wintertime over Northern and Central Europe in a series of global and regional climate simulations covering the last centuries, and in reanalysis products covering approximately the last 60 years. The focus of the study lies on identifying the link of the variations in the wind speed distribution to the regional near-surface temperature, to the meridional temperature gradient and to the North Atlantic Oscillation. Our main result is that the link between the daily wind distribution and the regional climate drivers is strongly model dependent. The global models tend to behave similarly, although they show some discrepancies. The two regional models also tend to behave similarly to each other, but surprisingly the results derived from each regional model strongly deviates from the results derived from its driving global model. In addition, considering multi-centennial timescales, we find in two global simulations a long-term tendency for the probability distribution of daily wind speed to widen through the last centuries. The cause for this widening is likely the effect of the deforestation prescribed in these simulations. We conclude that no clear systematic relationship between the mean temperature, the temperature gradient and/or the North Atlantic Oscillation, with the daily wind speed statistics can be inferred from these simulations. The understand- ing of past and future changes in the distribution of wind speeds, and thus of wind speed extremes, will require a detailed analysis of the representation of the interaction between large-scale and small-scale dynamics.

Sensitivity of the WRF model to PBL parametrisations and nesting techniques: evaluation of wind storms over complex terrain

Simulating surface wind over complex terrain is a challenge in regional climate modelling. Therefore, this study aims at identifying a set-up of the Weather Research and Forecasting Model (WRF) model that minimises system- atic errors of surface winds in hindcast simulations. Major factors of the model configuration are tested to find a suitable set-up: the horizontal resolution, the planetary boundary layer (PBL) parameterisation scheme and the way the WRF is nested to the driving data set. Hence, a number of sensitivity simulations at a spatial resolution of 2 km are carried out and compared to observations. Given the importance of wind storms, the analysis is based on case studies of 24 historical wind storms that caused great economic damage in Switzerland. Each of these events is downscaled using eight different model set-ups, but sharing the same driving data set. The results show that the lack of representation of the unresolved topography leads to a general overestimation of wind speed in WRF. However, this bias can be substantially reduced by using a PBL scheme that explicitly considers the effects of non-resolved topography, which also improves the spatial structure of wind speed over Switzerland. The wind direction, although generally well reproduced, is not very sensitive to the PBL scheme. Further sensitivity tests include four types of nesting methods: nesting only at the boundaries of the outermost domain, analysis nudging, spectral nudging, and the so-called re-forecast method, where the simulation is frequently restarted. These simulations show that restricting the freedom of the model to develop large-scale disturbances slightly increases the temporal agreement with the observations, at the same time that it further reduces the overestimation of wind speed, especially for maximum wind peaks. The model performance is also evaluated in the outermost domains, where the resolution is coarser. The results demonstrate the important role of horizontal resolution, where the step from 6 to 2 km significantly improves model performance. In summary, the combination of a grid size of 2 km, the non-local PBL scheme modified to explicitly account for non-resolved orography, as well as analysis or spectral nudging, is a superior combination when dynamical downscaling is aimed at reproducing real wind fields.

Reflection of solar wind hydrogen from the lunar surface

The solar wind continuously flows out from the Sun and directly interacts with the surfaces of dust and airless planetary bodies throughout the solar system. A significant fraction of solar wind ions reflect from an object's surface as energetic neutral atoms (ENAs). ENA emission from the Moon was first observed during commissioning of the Interstellar Boundary Explorer (IBEX) mission on 3 December 2008. We present the analysis of 10 additional IBEX observations of the Moon while it was illuminated by the solar wind. For the viewing geometry and energy range (> 250 eV) of the IBEX-Hi ENA imager, we find that the spectral shape of the ENA emission from the Moon is well-represented by a linearly decreasing flux with increasing energy. The fraction of the incident solar wind ions reflected as ENAs, which is the ENA albedo and defined quantitatively as the ENA reflection coefficient RN, depends on the incident solar wind speed, ranging from ~0.2 for slow solar wind to ~0.08 for fast solar wind. The average energy per incident solar wind ion that is reflected to space is 30 eV for slow solar wind and 45 eV for fast solar wind. Once ionized, these ENAs can become pickup ions in the solar wind with a unique spectral signature that reaches 3vSW. These results apply beyond the solar system; the reflection process heats plasmas that have significant bulk flow relative to interstellar dust and cools plasmas having no net bulk flow relative to the dust.

Plasticity to wind is modular and genetically variable in Arabidopsis thaliana

Thigmomorphogenesis, the characteristic phenotypic changes by which plants react to mechanical stress, is a widespread and probably adaptive type of phenotypic plasticity. However, little is known about its genetic basis and population variation. Here, we examine genetic variation for thigmomorphogenesis within and among natural populations of the model system Arabidopsis thaliana. Offspring from 17 field-collected European populations was subjected to three levels of mechanical stress exerted by wind. Overall, plants were remarkably tolerant to mechanical stress. Even high wind speed did not significantly alter the correlation structure among phenotypic traits. However, wind significantly affected plant growth and phenology, and there was genetic variation for some aspects of plasticity to wind among A. thaliana populations. Our most interesting finding was that phenotypic traits were organized into three distinct and to a large degree statistically independent covariance modules associated with plant size, phenology, and growth form, respectively. These phenotypic modules differed in their responsiveness to wind, in the degree of genetic variability for plasticity, and in the extent to which plasticity affected fitness. It is likely, therefore, that thigmomorphogenesis in this species evolves quasi-independently in different phenotypic modules.

Pollen dispersal properties of Poaceae and Cyperaceae: First estimates of their absolute pollen productivities

Pollen-trap results from the Swiss Alps 1996–2009 were used to assess the pollen dispersal–deposition properties of Poaceae (grasses) and Cyperaceae (sedges). Dispersal parameter values were investigated for a modified version of the Prentice–Sugita pollen dispersal–deposition model. Appropriate values (i.e. realistic in the field and allowing realistic modelling results) for wind speed are suggested to be in the range of 3–7 m s− 1 and for pollen an injection height of 0.03–0.1 m above the ground. The appropriate range of pollen injection height values for grasses and sedges differs from that of trees in the same area, suggesting different pollen dispersal properties between herbs and trees. In addition, logarithmic weighting of the vegetation was tested as an alternative to the modified Prentice–Sugita model. This yielded very similar results, suggesting that the use of such much simpler approximations of the pollen–vegetation relationship is a plausible alternative. Based on the modified Prentice–Sugita model, absolute pollen productivity for Poaceae was estimated to 7300 ± 400 grains cm− 2 year− 1 (1 SE). The data basis for Cyperaceae is smaller than for Poaceae, but the dispersal parameter values determined as appropriate for Poaceae yield good results. Absolute pollen productivity for Cyperaceae was estimated to 6300 ± 1100 grains cm− 2 year− 1 (1 SE).

Observations of middle atmospheric H2O and O3 during the 2010 major sudden stratospheric warming by a network of microwave radiometers

In this study, we present middle atmospheric water vapor (H2O) and ozone (O3) measurements obtained by ground-based microwave radiometers at three European locations in Bern (47° N), Onsala (57° N) and Sodankylä (67° N) during Northern winter 2009/2010. In January 2010, a major sudden stratospheric warming (SSW) occurred in the Northern Hemisphere whose signatures are evident in the ground-based observations of H2O and O3. The observed anomalies in H2O and O3 are mostly explained by the relative location of the polar vortex with respect to the measurement locations. The SSW started on 26 January 2010 and was most pronounced by the end of January. The zonal mean temperature in the middle stratosphere (10 hPa) increased by approximately 25 Kelvin within a few days. The stratospheric vortex weakened during the SSW and shifted towards Europe. In the mesosphere, the vortex broke down, which lead to large scale mixing of polar and midlatitudinal air. After the warming, the polar vortex in the stratosphere split into two weaker vortices and in the mesosphere, a new, pole-centered vortex formed with maximum wind speed of 70 m s−1 at approximately 40° N. The shift of the stratospheric vortex towards Europe was observed in Bern as an increase in stratospheric H2O and a decrease in O3. The breakdown of the mesospheric vortex during the SSW was observed at Onsala and Sodankylä as a sudden increase in mesospheric H2O. The following large-scale descent inside the newly formed mesospheric vortex was well captured by the H2O observations in Sodankylä. In order to combine the H2O observations from the three different locations, we applied the trajectory mapping technique on our H2O observations to derive synoptic scale maps of the H2O distribution. Based on our observations and the 3-D wind field, this method allows determining the approximate development of the stratospheric and mesospheric polar vortex and demonstrates the potential of a network of ground-based instruments.

A physics-based correction model for homogenizing sub-daily temperature series

A new physics-based technique for correcting inhomogeneities present in sub-daily temperature records is proposed. The approach accounts for changes in the sensor-shield characteristics that affect the energy balance dependent on ambient weather conditions (radiation, wind). An empirical model is formulated that reflects the main atmospheric processes and can be used in the correction step of a homogenization procedure. The model accounts for short- and long-wave radiation fluxes (including a snow cover component for albedo calculation) of a measurement system, such as a radiation shield. One part of the flux is further modulated by ventilation. The model requires only cloud cover and wind speed for each day, but detailed site-specific information is necessary. The final model has three free parameters, one of which is a constant offset. The three parameters can be determined, e.g., using the mean offsets for three observation times. The model is developed using the example of the change from the Wild screen to the Stevenson screen in the temperature record of Basel, Switzerland, in 1966. It is evaluated based on parallel measurements of both systems during a sub-period at this location, which were discovered during the writing of this paper. The model can be used in the correction step of homogenization to distribute a known mean step-size to every single measurement, thus providing a reasonable alternative correction procedure for high-resolution historical climate series. It also constitutes an error model, which may be applied, e.g., in data assimilation approaches.

Methodical study of nitrous oxide eddy covariance measurements using quantum cascade laser spectrometery over a Swiss forest

Nitrous oxide fluxes were measured at the Lägeren CarboEurope IP flux site over the multi-species mixed forest dominated by European beech and Norway spruce. Measurements were carried out during a four-week period in October–November 2005 during leaf senescence. Fluxes were measured with a standard ultrasonic anemometer in combination with a quantum cascade laser absorption spectrometer that measured N2O, CO2, and H2O mixing ratios simultaneously at 5 Hz time resolution. To distinguish insignificant fluxes from significant ones it is proposed to use a new approach based on the significance of the correlation coefficient between vertical wind speed and mixing ratio fluctuations. This procedure eliminated roughly 56% of our half-hourly fluxes. Based on the remaining, quality checked N2O fluxes we quantified the mean efflux at 0.8±0.4 μmol m−2 h−1 (mean ± standard error). Most of the contribution to the N2O flux occurred during a 6.5-h period starting 4.5 h before each precipitation event. No relation with precipitation amount could be found. Visibility data representing fog density and duration at the site indicate that wetting of the canopy may have as strong an effect on N2O effluxes as does below-ground microbial activity. It is speculated that above-ground N2O production from the senescing leaves at high moisture (fog, drizzle, onset of precipitation event) may be responsible for part of the measured flux.

Strong Down-Valley Low-Level Jets over the Atacama Desert: Observational Characterization

The near-surface wind and temperature regime at three points in the Atacama Desert of northern Chile is described using two-year multi-level measurements from 80-m towers located in an altitude range between 2100 and 2700 m ASL. The data reveal the frequent development of strong nocturnal drainage flows at all sites. Down-valley nose-shaped wind speed profiles are observed with maximum values occurring at heights between 20 m and 60 m AGL. The flow intensity shows considerable inter-daily variability and a seasonal modulation of maximum speeds, which in the cold season can attain hourly average values larger than 20 m s−1. Turbulent mixing appears significant over the full tower layer, affecting the curvature of the nighttime temperature profile and possibly explaining the observed increase of surface temperatures in the down-valley direction. Nocturnal valley winds and temperatures are weakly controlled by upper-air conditions observed at the nearest aerological station. Estimates of terms in the momentum budget for the development and the quasi-stationary phases of the down-valley flows suggest that the pressure gradient force due to the near-surface cooling along the sloping valley axes plays an important role in these drainage flows. A scale for the jet nose height of equilibrium turbulent down-slope jets is proposed, based on surface friction velocity and surface inversion intensity. At one of the sites this scale explains about 70% of the case-to-case observed variance of jet nose heights. Further modeling and observational work is needed, however, in order to better define the dynamics, extent and turbulence structure of this flow system, which has significant wind-energy, climatic and environmental implications.

Spatial heterogeneity and lake morphology affect diffusive greenhouse gas emission estimates of lakes

Most estimates of diffusive flux (F) of methane (CH4) and carbon dioxide (CO2) from lakes are based on single-point flux chamber measurements or on piston velocity (k) modeled from wind speed and single-point measurements of surface water gas concentrations (C-aq). We analyzed spatial variability of F of CH4 and CO2, as well as C-aq and k in 22 European lakes during late summer. F and k were higher in the lake centers, leading to considerable bias when extrapolating single-point chamber measurements to whole-lake estimates. The ratio of our empirical k estimates to wind speed-modeled k was related to lake size and shape, suggesting a lake morphology effect on the relationship between wind speed and k. This indicates that the error inherent to established wind speed models can be reduced by determining k and C-aq at multiple sites on lakes to calibrate wind speed-modeled k to the local system.

Topoclimatological case-study of Alpine pastures near the Albula Pass in the eastern Swiss Alps

Alpine grasslands are an important source of fodder for the cattle of Alpine farmers. Only during the short summer season can these pastures be used for grazing. With the anticipated climate change, it is likely that plant production – and thus the fodder basis for the cattle – will be influenced. Investigating the dependence of biomass production on topoclimatic factors will allow us to better understand how anticipated climate change may influence this traditional Alpine farming system. Because small-scale topoclimatological variations of the main meteorological variables: temperature, humidity, precipitation, shortwave incoming radiation and wind speed are not easily derived from available long-term climate stations in mountainous terrain, it was our goal to investigate the topoclimatic variations over the pastures belonging to the Alp Weissenstein research station north of the Albula Pass in the eastern Swiss Alps. We present a basic assessment of current topoclimatic conditions as a site characterization for ongoing ecological climate change studies. To be able to link short-term studies with long-term climate records, we related agrometeorological measurements with those of surrounding long-term sites run by MeteoSwiss, both on valley bottoms (Davos, Samedan), and on mountain tops (Weissfluhjoch, Piz Corvatsch). We found that the Davos climate station north of the study area is most closely correlated with the local climate of Alp Weissenstein, although a much closer site (Samedan) exists on the other side of the Albula Pass. Mountain top stations, however, did not provide a convincing approximation for the climate at Alp Weissenstein. Direct comparisons of near-surface measurements from a set of 11 small weather stations distributed over the domain where cattle and sheep are grazed indicate that nocturnal minimum air temperature and minimum vapor pressure deficit are mostly governed by the altitudinal gradient, whereas daily maxima – including also wind speed – are more strongly depending on vegetation cover and less on the altitude.

Observation errors in early historical upper-air observations

Upper-air observations are a fundamental data source for global atmospheric data products, but uncertainties, particularly in the early years, are not well known. Most of the early observations, which have now been digitized, are prone to a large variety of undocumented uncertainties (errors) that need to be quantified, e.g., for their assimilation in reanalysis projects. We apply a novel approach to estimate errors in upper-air temperature, geopotential height, and wind observations from the Comprehensive Historical Upper-Air Network for the time period from 1923 to 1966. We distinguish between random errors, biases, and a term that quantifies the representativity of the observations. The method is based on a comparison of neighboring observations and is hence independent of metadata, making it applicable to a wide scope of observational data sets. The estimated mean random errors for all observations within the study period are 1.5 K for air temperature, 1.3 hPa for pressure, 3.0 ms−1for wind speed, and 21.4° for wind direction. The estimates are compared to results of previous studies and analyzed with respect to their spatial and temporal variability.