Constraining a Flow Model with Field Measurements to Assess Water Transit Time Through a Vadose Zone

The modeling of thick vadose zones is particularly challenging because of difficulties in collecting a variety of measured sediment properties, which are required for parameterizing the model. Some models rely on synthetic data, whereas others are simplified by running as homogeneous sediment domains and relying on a single set of sediment properties. Few studies have simulated flow processes through a thick vadose zone using real and comprehensive data sets comprising multiple measurements. Here, we develop a flow model for a 7-m-thick vadose zone. This model, combining the numerical codes CTRAN/W with SEEP/W, includes the measured sediment hydraulic properties of the investigated vadose zone and incorporates the actual climate and subsurface conditions of the study site (precipitations, water-table elevations, and stable isotope data). The model is calibrated by fitting the simulated and measured vertical profiles of water content. Our flow model calculates a transit time of 1 year for the travel of water through the 7-m vadose zone; this estimate matches stable isotope-based results obtained previously for this site. A homogeneous sediment domain flow model, which considers only a single set of sediment properties, produces a transit time that is approximately half the duration of that of the heterogeneous flow model. This difference highlights the importance of assuming heterogeneous material within models of thick vadose zones and testifies to the advantage gained when using real sediment hydraulic properties to parametrize a flow model.     

Uncertainties in simulating regional climate of Southern Africa: sensitivity to physical parameterizations using WRF

This study aims at quantifying seasonal biases of regional climate model outputs during southern African summer, against a dense in situ measurement network (daily rain-gauge and surface air temperature records, and 12?h UTC radiosondes), and uncertainties associated with some physical parameterizations. Using the non-hydrostatic Advanced Research Weather Forecast (WRF) laterally forced by ERA40 reanalysis, twenty-seven experiments configured with three schemes of cumulus (CU), planetary boundary layer (PBL) and microphysics (MP), are performed at 35?km horizontal resolution during the core of a summer rainy season (December 1993 to February 1994 season) representative of the South African rainfall climatology. WRF simulates accurately seasonal large-scale rainfall patterns, as well as seasonal gradients of South African rainfall and 2-m temperature, and seasonal vertical profiles of the air temperature and humidity. However seasonal biases fluctuate strongly from an experiment to another, denoting considerable uncertainties generated by the physical package. Rainfall amounts are the most sensitive parameter to the tested schemes. Their geography, intensity, and intraseasonal characteristics are predominantly sensitive to CU schemes, and much less to PBL and MP schemes. Some CU-PBL combinations produce additive effects, which can dramatically either reduce or increase biases. Satisfactory configurations are found for South African climate, which would not have been possible without testing numerous physical parameterizations.     

The crustal structure of the Guayana Shield, Venezuela, from seismic refraction and gravity data

We present results from a seismic refraction experiment on the northern margin of the Guayana Shield performed during June 1998, along nine profiles of up to 320 km length, using the daily blasts of the Cerro Bolívar mines as energy source, as well as from gravimetric measurements. Clear Moho arrivals can be observed on the main E–W profile on the shield, whereas the profiles entering the Oriental Basin to the north are more noisy. The crustal thickness of the shield is unusually high with up to 46 km on the Archean segment in the west and 43 km on the Proterozoic segment in the east. A 20 km thick upper crust with P-wave velocities between 6.0 and 6.3 km/s can be separated from a lower crust with velocities ranging from 6.5 to 7.2 km/s. A lower crustal low velocity zone with a velocity reduction to 6.3 km/s is observed between 25 and 25 km depth. The average crustal velocity is 6.5 km/s. The changes in the Bouguer Anomaly, positive (30 mGal) in the west and negative (−20 mGal) in the east, cannot be explained by the observed seismic crustal features alone. Lateral variations in the crust or in the upper mantle must be responsible for these observations.     

Amplitude versus offset analyses of the deep sedimentary structures at the northern flank of the Porcupine Basin,SW of Ireland

This paper represents an attempt of analysing the amplitude versus offset (AVO) behaviour and specific seismic attributes of sedimentary structures from the Porcupine Basin, SW of Ireland. During the last decade, a huge number of carbonate mounds were investigated in this region in water depths of 600–1,000 m, but the genesis and growth of these mounds are still not clearly identified. The aim of this paper is to give a better understanding of the connection between fluid migration pathways in the deeper underground and surface expressions of their fluid expulsions like gas chimneys and pockmarks through which the mounds may generate themselves. The data used in this study to determine boundary conditions for the physical properties of the underlaying strata were gathered from the northern flank of the Porcupine Basin, where a huge amount of fluid and/or gas chimneys covers the seabed. Marine seismic reflection data contain information about the elastic properties of the underlying earth, mainly based on the observed variations in the seismic reflection amplitude at different shot–receiver offsets. To extract elastic parameters from the data, inversion techniques were used, which presume that input amplitudes are proportional to reflection coefficients for plane wave reflection. To calibrate the AVO analyses with the existing stratigraphy in the working area we have used the well logs from several bore holes in the region. The results of this study show clearly that the investigated and identified pockmarks on the seafloor are the surface expression of hydrocarbon seepage in the deeper sedimentary underground.     

Impact of partial steps and momentum advection schemes in a global ocean circulation model at eddy-permitting resolution

Series of sensitivity tests were performed with a z-coordinate, global eddy-permitting (1/4°) ocean/sea-ice model (the ORCA-R025 model configuration developed for the DRAKKAR project) to carefully evaluate the impact of recent state-of-the-art numerical schemes on model solutions. The combination of an energy–enstrophy conserving (EEN) scheme for momentum advection with a partial step (PS) representation of the bottom topography yields significant improvements in the mean circulation. Well known biases in the representation of western boundary currents, such as in the Atlantic the detachment of the Gulf Stream, the path of the North Atlantic Current, the location of the Confluence, and the strength of the Zapiola Eddy in the south Atlantic, are partly corrected. Similar improvements are found in the Pacific, Indian, and Southern Oceans, and characteristics of the mean flow are generally much closer to observations. Comparisons with other state-of-the-art models show that the ORCA-R025 configuration generally performs better at similar resolution. In addition, the model solution is often comparable to solutions obtained at 1/6 or 1/10° resolution in some aspects concerning mean flow patterns and distribution of eddy kinetic energy. Although the reasons for these improvements are not analyzed in detail in this paper, evidence is shown that the combination of EEN with PS reduces numerical noise near the bottom, which is likely to affect current–topography interactions in a systematic way. We conclude that significant corrections of the mean biases presently seen in general circulation model solutions at eddy-permitting resolution can still be expected from the development of numerical methods, which represent an alternative to increasing resolution.     

Climatological Ertel’s potential-vorticity flux and mean meridional circulation in the extratropical troposphere – lower stratosphere

We investigated to what extent the isentropic, non-geostrophic formulation of zonally averaged circulation derived for stratospheric conditions is applicable to climatological transport in the extratropical troposphere and lower stratosphere. The study is based on 10 years of daily data of ECMWF analysis and on the ECHAM3 climate model of the German Climate Computing Centre. The main result is a scalar isentropic mixing coefficient, K_yy, and a mean meridional transport circulation consistently derived from the same data base. For both data sources, isentropic mean meridional circulation is derived from horizontal mass flow rate for 4 representative months. Alternatively, a mean meridional circulation is calculated from total diabatic heating rates of the ECHAM3 model. It is shown that only the latter is in good agreement with the ECMWF mean meridional circulation. Isentropic analysis also comprises the seasonal cycle of the climatological meridional gradient and flux of Ertels potential vorticity (PV). Application of Tungs flux-gradient relation yields that for all seasons K_yy is positive in height-latitude regions where statistical significance is reached. Large K_yy values, marking regions of more efficient mixing, have been found in the subtropical vertical band of weak westerly wind and in mid-latitudes in regions of upward-propagating baroclinic wave activity in the middle and upper troposphere. Based on the ECMWF data and results of baroclinic-wave behaviour, strong indications are presented that positive zonally averaged PV flux polewards of the jet core in the NH is strengthened by stationary waves and nonlinear effects. Reduced eddy transport is apparent in winter and spring slightly below the subtropical tropopause jet. The seasonal cycle of K_yy from ECHAM3 data is to a great extent in agreement with the result based on ECMWF analysis. In the model, reduced interannual variability enlarges the height-latitude range where sign of K_yy is significant.     

Contribution à l'étude de la vitesse critique d'érosion des sols cohésifs

The phenomena of erosion and sedimentation in rivers are treated qualitatively and quantitatively. An analytical solution of the bedload critical velocity for erosion and sedimentation is proposed, depending on material property, geometry and on flow characteristics. The critical velocity of erosion of the river bed, as defined by Hjulström in 1935, can be explained by a cohesive behaviour of the grains. Whereas the cohesive force is dominated by an 1/r² (r is the radius of grains) force when grains are small, it can be modelled by the cohesion C and the friction angle φ for larger grains. To cite this article: J. Gargani, C. R. Geoscience 336 (2004).     

Quality Assurance of Hydrocarbon Measurements for the German Tropospheric Research Focus (TFS)

In the German Focus on Tropospheric Research (TFS) independent quality assurance procedures were implemented in order to obtain information on data quality and comparability of the different measurements made in the different field campaigns. This paper describes the results for measurements of hydrocarbons using in-situ gas chromatographictechniques and off line analysis of samples collected in canisters (analysed by two central laboratories CL-1 and CL-2) and samples collected on adsorption tubes. The QA-procedures included comparisons with synthetic standards (prepared by EN 45001 certified laboratories), absolute calibration with a diffusion source, and an instrument/methodology comparison in ambient air. Harmonisation of the ambient measurements was achieved with a complex mixture (NCAR/BERLIOZ) containing 70 commonly observed hydrocarbons at mixing ratios of 0.2 to 10 ppb (mole fraction) in nitrogen, which was calibrated by referencing to hydrocarbon standards of the National Institute for Standards and Technology (NIST).For the certified synthetic standards, the experienced groups agreed to within ± 20% for most compounds. Much larger discrepancies were observed for the new Airmotec HC2010 instruments due to problems with identification, co-elution and blanks/memory effects. The results in ambient air were similar: Reasonable agreement was found for the results from the experienced groups with well characterised in situ instruments andfor the charcoal tubes, whereas larger discrepancies were observed for the results from new groups and instruments. For the latter, only selected compounds met the data quality objectives (DQO). The canister samples analysed by CL-2 were in good agreement with the reference instrument, whereas large deviations were found for a number of compounds in the analysis of the same canisters by CL-1. The results of the comparison provided the final basis for flagging and harmonising the data from all participants prior to their submission to the TFS central data bank.