A set-membership approach for high integrity height-aided satellite positioning

A Robust Set-Inversion via Interval Analysis method in a bounded-error framework is used to compute three-dimensional location zones in real time, at a given confidence level. This approach differs significantly from the usual Gaussian error model paradigm, since the satellite positions and the pseudorange measurements are represented by intervals encompassing the true value with a particular level of confidence. The method computes a location zone recursively, using contractions and bisections of an arbitrarily large initial location box. Such an approach can also handle an arbitrary number of erroneous measurements using a q-relaxed solver and allows the integration of geographic and cartographic information such as digital elevation models or three-dimensional maps. With enough data redundancy, inconsistent measurements can be detected and even rejected. The integrity risk of the location zone comes only from the measurement bounds settings, since the solver is guaranteed. A method for setting these bounds for a particular location zone confidence level is proposed. An experimental validation using real L1 code measurements and a digital elevation model is also reported in order to illustrate the performance of the method on real data.     

A mass‐balance approach to estimate in‐stream processes in a large river

A mass‐balance approach was used to estimate in‐stream processes related to inorganic nitrogen species (NH₄⁺, NO₂^? and NO₃^?) in a large river characterized by highly variable hydrological conditions, the Garonne River (south‐west France). Studies were conducted in two consecutive reaches of 30 km located downstream of the Toulouse agglomeration (population 760 000, seventh order), impacted by modification of discharge regime and high nitrogen concentrations. The mass‐balance was calculated by two methods: the first is based on a variable residence time (VRT) simulated by a one‐dimensional (1‐D) hydraulic model; the second is a based on a calculation using constant residence time (CRT) evaluated according to hydrographic peaks. In the context of the study, removal of dissolved inorganic nitrogen (DIN) for a reach of 30 km is underestimated by 11% with the CRT method. In sub‐reaches, the discrepancy between the two methods led to a 50% overestimation of DIN removal in the upper reach (13 km) and a 43% underestimation in the lower reach (17 km) using the CRT method. The study highlights the importance of residence time determination when using modelling approaches in the assessment of whole stream processes in short‐duration mass‐balance for a large river under variable hydrological conditions. Copyright © 2007 John Wiley & Sons, Ltd.     

Influence of a hedge surrounding bottomland on seasonal soil‐water movement

The influence of a hedge surrounding bottomland on soil‐water movement along the hillslope was studied on a plot scale for 28 months. The study was based on the comparison of two transects, one with a hedge, the other without, using mainly a dense grid of tensiometers. The influence of the bottomland hedge was located in the area where tree roots were developed, several metres upslope from the hedge, and could be observed both in the saturated and non‐saturated zone, from May to December. The hedge induced a high rate of soil drying, because of the high evaporative capacity of the trees. We evaluated that water uptake by the hedge during the growing season was at least 100 mm higher than without a hedge. This increased drying rate led to a delayed rewetting of the soils upslope from the hedge in autumn, of about 1 month compared with the situation without a hedge. Several consequences of this delayed rewetting are expected: a delay in the return of subsurface transfer from the hillslope to the riparian zone, a buffering effect of hedges on floods, already observed at the catchment scale, and an increased residence time of pollutants. Copyright © 2003 John Wiley & Sons, Ltd.     

Estimating Long-Term Solar Irradiance Variability: A New Approach

The detection of solar irradiance variations (both bolometric and at various wavelengths) by satellite-based experiments during the last one-and-a-half decades stimulated modeling efforts to help identify their causes and to provide estimates of irradiance data for those time intervals when no satellite observations exist. In this paper we present estimates of the long-term irradiance changes developed with Fourier and wavelet transforms. The month-to-month irradiance variations, after removing the solar cycle related long-term changes, are studied with the cross-correlation technique. Results of the analysis reveal a significant phase shift at 3 months between the full-disk magnetic field strength and total solar and UV irradiance, with irradiance leading the magnetic field variability. In addition to this time delay between the changes in solar irradiance and the magnetic field, a 10-month phase shift has been found between the UV flux at 280 nm and total solar irradiance corrected for sunspot darkening. The existence of these phase shifts suggests the possibility of a coupling between the physical processes taking place below, in, and above the photosphere.     

The Greenland ice sheet and greenhouse warming

Increased melting on glaciers and ice sheets and rising sea level are often mentioned as important aspects of the anticipated greenhouse warming of the earth's atmosphere. This paper deals with the sensitivity of Greenland's ice mass budget and presents a tentative projection of the Greenland component of future sea level rise for the next few hundred years. To do this, the ‘Villach II temperature scenario’ is prescribed,output from a comprehensive mass balance model is used to drive a high-resolution 3-D thermomechanic model of the ice sheet.The mass balance model consists of two parts: the accumulation part is based on presently observed values and is forced by changes in mean anr tempeerature. The ablation model is based on the degree-day method and accounts for daily and annual temperature cycle, a different degree-day factor for ice and snow melting and superimposed ice formation. Under present-day climatic conditions, the following total mass balance results (in ice equivalent per years): 599.3 × 10⁹ m³ of accumulation, 281.7 × 10⁹m³ of runoff assuming a balanced budget, 317.6 × 10⁹m³ of iceberg calving. A 1K uniform warming is then calculated to increase the runoff by 119.5 × 10⁹ m³. Since accumulation also increases by 32 × 10⁹ m³, this leads to reduction of the total mass balance by 887.5 × 10⁹ m³ of ice, corresponding to a sea level rise of 0.22 mm/yr. For temperature increase larger than 2.7 K, runoff, exceeds accumulation, and if ice sheet dynamics were to remain unchanged, this would add an extra amount of 0.8 mmyr to the worl's oceans.Imposing the Villach II scenario (warming up to 4.23 K) and accumulating mass balance changes forward in time (static response) would then result in a global sea level rise of 7.1 cm by 2100 AD, but this figure may go up to as much as 40 cm per century in case the warming is doubled. In a subsequent dynamic model involving the ice flow, the ice sheet is found to produce a counteracting effect by dynamically producing steeper slopes at the margin, thereby reducing the area over which runoff can take place. This effect is particularly apparent in the northeastern part of the ice sheet, and is also more pronounced for the smaller temperature perturbations. Nevertheless, all these experiments certainly highlight the vulnerability of the Greenland ice sheet with respect to a climatic warming.     

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.