Flexural deformation controls on Late Quaternary sediment dispersal in the Garo-Rajmahal Gap,NW Bengal Basin

Subsurface deformation is a driver for river path selection when deformation rates become comparable to the autogenic mobility rate of rivers. Here we combine geomorphology, soil and sediment facies analyses, and geophysical data of the Late Quaternary sediments of the central Garo-Rajmahal Gap in Northwest Bengal to link subsurface deformation with surface processes. We show variable sedimentation characteristics, from slow rates (<0.8 mm/year) in the Tista megafan at the foot of the Himalaya to nondeposition at the exposed surface of the Barind Tract to the south, enabling the development of mature soils. Combined subsidence in the Tista fan and uplift of the Barind Tract are consistent with a N-S flexural response of the Indian plate to loading of the Himalaya Mountains given a low value of elastic thickness (15–25 km). Provenance analysis based on bulk strontium concentration suggests a dispersal of sediment consistent with this flexural deformation—in particular the abandonment of the Barind Tract by a Pleistocene Brahmaputra River and the current extents of the Tista megafan lobes. Overall, these results highlight the control by deeply rooted deformation patterns on the routing of sediment by large rivers in foreland settings.     

Feasibility study of a multiple flyby mission of main belt asteroids

Specific information on the surface morphology, composition, mean density, and internal structure of asteroids, which is necessary to advance our understanding of asteroids, can be obtained only by a detailed investigation of individual bodies: this will require space missions to individual targets. Since an essential characteristic of the asteroids is their variety, several objects must be visited. The Ariane launcher developed presently in Europe makes a multiple flyby mission possible. The first results of our feasibility study are particularly encouraging: during one revolution, five to six preselected main belt asteroids may be approached to within 1000 km with relative velocities which lie between 3 and 14 km/sec using a total impulse correction on the order of 1 to 2 km/sec. The weight of the spacecraft, excluding the engine and the propellant, would be at least 250 kg. This allows a scientific payload of 50 to 60 kg, in which priority will be given to an imaging system and radar altimeter.     

Extreme value theory based on the r largest annual events: a robust approach

In Smith (1986, J. Hydrol. 86, 27–43), a family of statistical distributions and estimators for extreme values based on a fixed number r > = 1 of the largest annual events are presented. The method of estimation was numerical maximum likelihood. In this paper, we consider the robust estimation of parameters in such families of distributions. The estimation technique, which is based on optimal B-robust estimates, will assign weights to each observation and give estimates of the parameters based on the data which are well modeled by the distribution. Thus, observations which are not consistent with the proposed distribution can be identified and the validity of the model can be assessed. The method is illustrated on Venice sea level data.     

Turbulent heat fluxes by profile and inertial dissipation methods: analysis of the atmospheric surface layer from shipboard measurements during the SOFIA/ASTEX and SEMAPHORE experiments

Heat flux estimates obtained using the inertial dissipation method, and the profile method applied to radiosonde soundings, are assessed with emphasis on the parameterization of the roughness lengths for temperature and specific humidity. Results from the inertial dissipation method show a decrease of the temperature and humidity roughness lengths for increasing neutral wind speed, in agreement with previous studies. The sensible heat flux estimates were obtained using the temperature estimated from the speed of sound determined by a sonic anemometer. This method seems very attractive for estimating heat fluxes over the ocean. However allowance must be made in the inertial dissipation method for non-neutral stratification. The SOFIA/ASTEX and SEMAPHORE results show that, in unstable stratification, a term due to the transport terms in the turbulent kinetic energy budget, has to be included in order to determine the friction velocity with better accuracy. Using the profile method with radiosonde data, the roughness length values showed large scatter. A reliable estimate of the temperature roughness length could not be obtained. The humidity roughness length values were compatible with those found using the inertial dissipation method.     

Study of the air-sea interactions at the mesoscale: the SEMAPHORE experiment

The SEMAPHORE (Structure des Echanges Mer-Atmosphère, Propriétés des Hétérogénéités Océaniques: Recherche Expérimentale) experiment has been conducted from June to November 1993 in the Northeast Atlantic between the Azores and Madeira. It was centered on the study of the mesoscale ocean circulation and air-sea interactions. The experimental investigation was achieved at the mesoscale using moorings, floats, and ship hydrological survey, and at a smaller scale by one dedicated ship, two instrumented aircraft, and surface drifting buoys, for one and a half month in October-November (IOP: intense observing period). Observations from meteorological operational satellites as well as spaceborne microwave sensors were used in complement. The main studies undertaken concern the mesoscale ocean, the upper ocean, the atmospheric boundary layer, and the sea surface, and first results are presented for the various topics. From data analysis and model simulations, the main characteristics of the ocean circulation were deduced, showing the close relationship between the Azores front meander and the occurrence of Mediterranean water lenses (meddies), and the shift between the Azores current frontal signature at the surface and within the thermocline. Using drifting buoys and ship data in the upper ocean, the gap between the scales of the atmospheric forcing and the oceanic variability was made evident. A 2 °C decrease and a 40-m deepening of the mixed layer were measured within the IOP, associated with a heating loss of about 100 W m⁻². This evolution was shown to be strongly connected to the occurrence of storms at the beginning and the end of October. Above the surface, turbulent measurements from ship and aircraft were analyzed across the surface thermal front, showing a 30% difference in heat fluxes between both sides during a 4-day period, and the respective contributions of the wind and the surface temperature were evaluated. The classical momentum flux bulk parameterization was found to fail in low wind and unstable conditions. Finally, the sea surface was investigated using airborne and satellite radars and wave buoys. A wave model, operationally used, was found to get better results compared with radar and wave-buoy measurements, when initialized using an improved wind field, obtained by assimilating satellite and buoy wind data in a meteorological model. A detailed analysis of a 2-day period showed that the swell component, propagating from a far source area, is underestimated in the wave model. A data base has been created, containing all experimental measurements. It will allow us to pursue the interpretation of observations and to test model simulations in the ocean, at the surface and in the atmospheric boundary layer, and to investigate the ocean-atmosphere coupling at the local and mesoscales.     

Modelling peak accelerations from earthquakes

This paper deals with the prediction of peak horizontal accelerations with emphasis on seismic risk and insurance concerns. Non‐linear mixed effects models are used to analyse well‐known earthquake data and the consequences of mis‐specifying assumptions on the error term are quantified. A robust fit of the usual model, using recently developed robust weighted maximum likelihood estimators, is presented. Outlying data are automatically identified and subsequently investigated. A more appropriate model accounting for the extreme value nature of the responses, is also developed and implemented. The implication on acceleration predictions is demonstrated. Copyright © 2006 John Wiley & Sons, Ltd.     

Upper Limit of Iodine Oxide in the Lower Stratosphere

It has been suggested that iodine oxides, IOx, could play a significant role in the ozone destruction in the lower stratosphere. To investigate this suggestion, spectra from nine SAOZ uv-visible spectrometer balloon flights were examined for the IO absorption signature between 405 and 450 nm. IO was not detected, either at mid- or high latitude, in the morning or the evening, in summer or winter. An upper limit of 0.2 parts per trillion by volume (pptv) at 20 km and 0.1 pptv at 15 km at the 95% confidence level (2), was derived from the best measurements at 90° SZA at sunset and sunrise. Since a photochemical model shows that 70% of inorganic iodine should be in the form IO at that time, it is concluded that unless iodine chemistry is different from that assumed at the moment, the role of iodine in stratospheric ozone depletion is small.     

FUSE Observations of EC14026 Stars

We present a FUSE abundance analysis of EC14026 stars. We compare the abundances of heavy elements in the atmospheres of EC14026 stars to non-pulsating stars with similar atmospheric parameters, and investigate whether weak stellar winds could explain the coexistence of variable and non-variable sdB stars in the log g – T_eff diagram. We also present preliminary results on time-dependent diffusion calculations of iron in presence of radiative levitation and mass loss, and show how weak stellar winds can affect the diffusive equilibrium between gravitational settling and radiative support.     

Preliminary hydrodynamic results of a field experiment on a barred beach, Truc Vert beach on October 2001

A field experiment conducted on a sandy barred beach, situated on the southern part of the French Atlantic coastline, allowed us to investigate the impact of the intertidal bar on the wave-energy dissipation on the beach face in presence of a high-energy long-incoming swell (significant wave height of about 1.7 to 3.0 m in 56 m water depth and significant wave period about 12 s). Data were collected along three parallel cross-shore transects deployed along an intertidal ridge and runnel system. Wave heights in the inner surf zone are depth-limited, consistent with previous works, and the wave-energy dissipation in the inner surf zone appears to be relatively independent of the offshore energy level. On the other hand, the presence of the bar seems to scatter the data. In models of surf-zone hydrodynamics, wave-energy dissipation is often parameterized in terms of , the ratio of the sea-swell significant wave height to the local mean water depth. The observed values of are not constant along a cross-shore transect, and increase onshore. Furthermore, the observed values observed onshore the intertidal bar are higher than those observed outside the influence of the intertidal bar, and this cannot be fully explained by the different local beach slope.Responsible Editor: Iris Grabemann     

Direct forecasting of reservoir performance using production data without history matching

The conventional paradigm for predicting future reservoir performance from existing production data involves the construction of reservoir models that match the historical data through iterative history matching. This is generally an expensive and difficult task and often results in models that do not accurately assess the uncertainty of the forecast. We propose an alternative re-formulation of the problem, in which the role of the reservoir model is reconsidered. Instead of using the model to match the historical production, and then forecasting, the model is used in combination with Monte Carlo sampling to establish a statistical relationship between the historical and forecast variables. The estimated relationship is then used in conjunction with the actual production data to produce a statistical forecast. This allows quantifying posterior uncertainty on the forecast variable without explicit inversion or history matching. The main rationale behind this is that the reservoir model is highly complex and even so, still remains a simplified representation of the actual subsurface. As statistical relationships can generally only be constructed in low dimensions, compression and dimension reduction of the reservoir models themselves would result in further oversimplification. Conversely, production data and forecast variables are time series data, which are simpler and much more applicable for dimension reduction techniques. We present a dimension reduction approach based on functional data analysis (FDA), and mixed principal component analysis (mixed PCA), followed by canonical correlation analysis (CCA) to maximize the linear correlation between the forecast and production variables. Using these transformed variables, it is then possible to apply linear Gaussian regression and estimate the statistical relationship between the forecast and historical variables. This relationship is used in combination with the actual observed historical data to estimate the posterior distribution of the forecast variable. Sampling from this posterior and reconstructing the corresponding forecast time series, allows assessing uncertainty on the forecast. This workflow will be demonstrated on a case based on a Libyan reservoir and compared with traditional history matching.