A note of caution when interpreting parameters of the distribution of excesses

In climatology and hydrology, univariate Extreme Value Theory has become a powerful tool to model the distribution of extreme events. The Generalized Pareto Distribution (GPD) is routinely applied to model excesses in space or time by letting the two GPD parameters depend on appropriate covariates. Two possible pitfalls of this strategy are the modeling and the interpretation of the scale and shape GPD parameters estimates which are often and incorrectly viewed as independent variables. In this note we first recall a statistical technique that makes the GPD estimates less correlated within a Maximum Likelihood (ML) estimation approach. In a second step we propose novel reparametrizations for two method-of-moments particularly popular in hydrology: the Probability Weighted Moment (PWM) method and its generalized version (GPWM). Finally these three inference methods (ML, PWM and GPWM) are compared and discussed with respect to the issue of correlations.     

How does micropetrography help us to understand the permeability and poromechanical behaviour of a rock?

Many poromechanical studies of rocks give little regard to the micropetrographic observation of their constituents. The aim of this study is to show how such a detailed study can help understand and provide more precision on the poromechanical behaviour of the material. The rock studied here is an oolitic limestone. Micropetrography reveals the nature of the constituents (calcite, mostly in the form of fine‐grained micrite), their abundance, their organization (micrite organized in ooids linked by a calcitic cement, small‐sized highly connected isotropic porosity, scarce unconnected porosity) and the structure of the rock (global isotropy, small heterogeneities due to sedimentological phenomena, no microcracks). Micropetrography completes the information obtained by poromechanics in that it allows the observation of elements deduced from macroscopical tests. The aim of this paper is to show the effect of microporosity on permeability and flow in a rock where microporous grains are in contact with each other.     

Rb-Sr and Sm-Nd ages and isotopic geochemistry of Archaean granodioritic gneisses from Eastern Finland

Twenty granodioritic rocks and one amphibolitic enclave of the “basement” of the Suomussalmi-Kuhmo Archaean (2.65 Ga) greenstone belts (central-eastern Finland), have been chosen together with one greenstone sample for Rb-Sr and Sm-Nd geochronological and isotopic studies.The granitoïd rocks are subdivided into three groups: two generations of grey gneisses and a post-belt augen gneiss. The Rb-Sr ages of the first and second generation of grey gneisses are 2.86 ± 0.09 and 2.62 ± 0.07 Ga, respectively. These results are corroborated by Sm-Nd data. The post-belt augen gneiss gives an age of 2.51 ± 0.11 Ga. The results show that the two generations of grey gneisses, the greenstone belts and the post-greenstone augen gneiss, were developed over a period > 350 Ma. The two generations of grey gneisses show identical I_Sr values (0.7023 ± 8 and 0.7024 ± 6) which contrast with that of the augen gneiss (0.7049 ± 8). The low I_Sr and the near-chondritic ?_T^CHUR values indicate that the grey gneisses cannot derived from much older continental materials. Trace element studies suggest that these grey gneisses have had a multi-stage development. The augen gneiss with a moderately high I_Sr is likely to be derived from a granodiorite originated by partial melting of older sialic crust. The more probable parent rock seems to be the first generation grey gneisses. The I_Sr and average Rb/Sr values preclude the greenstone belt and the second generation of grey gneisses as the protolith.     

Impact of Terrain Heterogeneity on Coherent Structure Properties: Numerical Approach

A three-dimensional large-eddy simulation (LES) model, which includes the effects of plant–atmosphere interactions, is used to study the effects of surface inhomogeneities on near-surface coherent structures over an open field and behind a forest canopy. These simulated conditions are representative of two wind sectors of the Site Instrumental de Recherche par Télédétection Atmosphérique (SIRTA) experimental site at the Institut Pierre Simon Laplace, Palaiseau, France. Coherent structure properties deduced from wavelet transforms of the simulated near-surface vertical velocity time series are not modified by upstream terrain heterogeneities, in agreement with site measurements. This feature is related to the nature of structures detected from the vertical velocity time series. The turbulence close to the surface seems composed of both local coherent structures and large coherent structures reflecting outer-layer properties, which depend on the overall surface heterogeneity or upstream heterogeneity. It is argued that the streamwise velocity is representative of these large outer-layer structures that impinge onto the ground through a top-down mechanism as identified through the space–time correlation of the wind velocity components. In contrast, the vertical velocity is more representative of small structures resulting from the impingement of the large outer-layer structures. These small structures represent locally-generated, active turbulence, which adjusts rapidly to local surface conditions, and consequently they are only weakly dependent on upstream heterogeneities.     

Lander radioscience for obtaining the rotation and orientation of Mars

The paper presents the concept, the objectives, the approach used, and the expected performances and accuracies of a radioscience experiment based on a radio link between the Earth and the surface of Mars. This experiment involves radioscience equipment installed on a lander at the surface of Mars. The experiment with the generic name lander radioscience (LaRa) consists of an X-band transponder that has been designed to obtain, over at least one Martian year, two-way Doppler measurements from the radio link between the ExoMars lander and the Earth (ExoMars is an ESA mission to Mars due to launch in 2013). These Doppler measurements will be used to obtain Mars’ orientation in space and rotation (precession and nutations, and length-of-day variations). More specifically, the relative position of the lander on the surface of Mars with respect to the Earth ground stations allows reconstructing Mars’ time varying orientation and rotation in space.Precession will be determined with an accuracy better by a factor of 4 (better than the 0.1% level) with respect to the present-day accuracy after only a few months at the Martian surface. This precession determination will, in turn, improve the determination of the moment of inertia of the whole planet (mantle plus core) and the radius of the core: for a specific interior composition or even for a range of possible compositions, the core radius is expected to be determined with a precision decreasing to a few tens of kilometers.A fairly precise measurement of variations in the orientation of Mars’ spin axis will enable, in addition to the determination of the moment of inertia of the core, an even better determination of the size of the core via the core resonance in the nutation amplitudes. When the core is liquid, the free core nutation (FCN) resonance induces a change in the nutation amplitudes, with respect to their values for a solid planet, at the percent level in the large semi-annual prograde nutation amplitude and even more (a few percent, a few tens of percent or more, depending on the FCN period) for the retrograde ter-annual nutation amplitude. The resonance amplification depends on the size, moment of inertia, and flattening of the core. For a large core, the amplification can be very large, ensuring the detection of the FCN, and determination of the core moment of inertia.The measurement of variations in Mars’ rotation also determines variations of the angular momentum due to seasonal mass transfer between the atmosphere and ice caps. Observations even for a short period of 180 days at the surface of Mars will decrease the uncertainty by a factor of two with respect to the present knowledge of these quantities (at the 10% level).The ultimate objectives of the proposed experiment are to obtain information on Mars’ interior and on the sublimation/condensation of CO₂ in Mars’ atmosphere. Improved knowledge of the interior will help us to better understand the formation and evolution of Mars. Improved knowledge of the CO₂ sublimation/condensation cycle will enable better understanding of the circulation and dynamics of Mars’ atmosphere.     

Intercomparison of the LASCO-C2, SECCHI-COR1, SECCHI-COR2, and Mk4 Coronagraphs

In order to assess the reliability and consistency of white-light coronagraph measurements, we report on quantitative comparisons between polarized brightness [pB] and total brightness [B] images taken by the following white-light coronagraphs: LASCO-C2 on SOHO, SECCHI-COR1 and -COR2 on STEREO, and the ground-based MLSO-Mk4. The data for this comparison were taken on 16?April 2007, when both STEREO spacecraft were within 3.1^° of Earth??s heliographic longitude, affording essentially the same view of the Sun for all of the instruments. Due to the difficulties of estimating stray-light backgrounds in COR1 and COR2, only Mk4 and C2 produce reliable coronal-hole values (but not at overlapping heights), and these cannot be validated without rocket flights or ground-based eclipse measurements. Generally, the agreement between all of the instruments?? pB values is within the uncertainties in bright streamer structures, implying that measurements of bright CMEs also should be trustworthy. Dominant sources of uncertainty and stray light are discussed, as is the design of future coronagraphs from the perspective of the experiences with these instruments.     

Mapping the mean monthly precipitation of a small island using kriging with external drifts

This study focuses on the spatial distribution of mean annual and monthly precipitation in a small island (1128 km²) named Martinique, located in the Lesser Antilles. Only 35 meteorological stations are available on the territory, which has a complex topography. With a digital elevation model (DEM), 17 covariates that are likely to explain precipitation were built. Several interpolation methods, such as regression-kriging (????????, ????????, and ????????) and external drift kriging (??????) were tested using a cross-validation procedure. For the regression methods, predictors were chosen by established techniques whereas a new approach is proposed to select external drifts in a kriging which is based on a stepwise model selection by the Akaike Information Criterion (AIC). The prediction accuracy was assessed at validation sites with three different skill scores. Results show that using methods with no predictors such as inverse distance weighting (??????) or universal kriging (????) is inappropriate in such a territory. ?????? appears to outperform regression methods for any criteria, and selecting predictors by our approach improves the prediction of mean annual precipitation compared to kriging with only elevation as drift. Finally, the predicting performance was also studied by varying the size of the training set leading to less conclusive results for ?????? and its performance. Nevertheless, the proposed method seems to be a good way to improve the mapping of climatic variables in a small island.     

Chemistry of rainwater in the Massif Central (France): a strontium isotope and major element study

Atmospheric aerosols (sea salt, crustal dust, and biogenic aerosols) are the primary source of dissolved species in rainwater as well as one of the sources of dissolved species in river water. Chemical weathering studies require quantification of this atmospheric input. The crustal component of atmospheric input can have various origins, both distant and local. The proportions of the various inputs (marine, distant or local) are determined in this study.Strontium isotope ratios and Ca, Na, K, Mg, Al, Cl, SO₄, NO₃ and Sr concentrations were measured in rainwater samples collected in the Massif Central (France) over a period of one year. Each sample, collected automatically, represents a monthly series of rain events. Chemical composition of the rainwater samples varied considerably and the ⁸⁷Sr/⁸⁶Sr ratios ranged between 0.709198 and 0.713143.Using Na as an indicator of marine origin, and Al for the crustal input in rain samples, the proportion of marine and crustal elements was estimated from elemental ratios. A marine origin of 4 to 100% of Cl, of 0.6 to 20% of the SO₄, of <1 to 10% of Ca, <1 to 40% of K, 4 to 100% of Mg and 1 to 44% of Sr was determined.Strontium isotopes were used to characterize the crustal sources. The ⁸⁷Sr/⁸⁶Sr ratios of the crustal sources varied considerably from 0.7092 to 0.71625 and indicate the occurrence of multiple sources for the crustal component in the analysed rainwaters.     

Multiple-point statistics using multi-resolution images

Stochastic Environmental Research and Risk Assessment - Multiple-point statistics (MPS) is a simulation technique allowing to generate images that reproduce the spatial features present in a...     

Nearby Protoclusters as Laboratories for Understanding Star Formation on Galactic Scales

Detailed studies of nearby cluster-forming molecular clouds can help us understand the physical processes by which most stars form in galaxies. I review recent advances made on this subject. Submillimeter observations of nearby protoclusters suggest that stars are generally built from finite, detached reservoirs of mass inside molecular cloud cores, and point to a cloud fragmentation origin for the IMF. Much progress in this field will come from future large submillimeter instruments such as Herschel and ALMA. This revised version was published online in September 2006 with corrections to the Cover Date.