Discrete fracture model for coupled flow and geomechanics

We present a fully implicit formulation of coupled flow and geomechanics for fractured three-dimensional subsurface formations. The Reservoir Characterization Model (RCM) consists of a computational grid, in which the fractures are represented explicitly. The Discrete Fracture Model (DFM) has been widely used to model the flow and transport in natural geological porous formations. Here, we extend the DFM approach to model deformation. The flow equations are discretized using a finite-volume method, and the poroelasticity equations are discretized using a Galerkin finite-element approximation. The two discretizations—flow and mechanics—share the same three-dimensional unstructured grid. The mechanical behavior of the fractures is modeled as a contact problem between two computational planes. The set of fully coupled nonlinear equations is solved implicitly. The implementation is validated for two problems with analytical solutions. The methodology is then applied to a shale-gas production scenario where a synthetic reservoir with 100 natural fractures is produced using a hydraulically fractured horizontal well.     

Internal variability of regional climate models

 Two regional climate models have been applied to the task of generating an ensemble of realizations of the year 1982 with observed boundary conditions in areas covering parts of the Mediterranean countries. These realizations were generated by applying boundary conditions from the ECMWF ERA reanalysis project consecutively, carrying over the soil variables from the regional models from one iteration to the next. Monthly mean fields for six iterations of each model have been used as statistical ensembles in order to investigate the internal variability of the regional model dynamics. This internal variability is a necessary consequence of the non-linear physical feedback mechanisms of the RCM being active. A small value of internal variability will give better statistics for climate sensitivity signals, but will make these results less credible. The internal variability is important for the quantitative assessment of a climate sensitivity signal. With the present choice of models and integration domains the internal variabilities of surface fields and precipitation do reach levels that are less than, but in summer of comparable order of magnitude to, corresponding atmospheric variabilities of an atmospheric general circulation model. Received: 26 October 1999 / Accepted: 18 December 2000     

Uncertainties associated with Ra and Ra measurements in water via a Delayed Coincidence Counter (RaDeCC)

The short-lived radium isotopes, ²²³Ra (T_1/2 = 11.4 days) and ²²⁴Ra (T_1/2 = 3.66 days), have been successfully used as tracers of several environmental processes, e.g., submarine groundwater discharge, coastal mixing processes, and water residence times. In this paper, the uncertainties associated with ²²³Ra and ²²⁴Ra measurements using a Radium Delayed Coincidence Counter are determined on a detailed error propagation basis with a confidence interval of 1σ. From the data analyses of several groups of coastal water samples, the calculated relative uncertainties averaged 12% for the ²²³Ra and 7% for the ²²⁴Ra. These uncertainties can decrease for radium-enriched groundwater samples although asymptotic limits have been found at 7% relative uncertainty for ²²³Ra and 4% for ²²⁴Ra. In this paper, the influence of sampling and measurement parameters on the final radium uncertainties is evaluated in order to provide guidance to optimize these factors and obtain more reliable results.     

Evaluation of food competition and resource partitioning of recruiting fish with permanent residents in a seagrass habitat

ABSTRACT

Recruiting King George whiting Sillaginodes punctatus were studied to assess the potential for food competition with permanent resident fish species in a nursery habitat. Marine migrant post larval S. punctatus (<60?mm TL) consumed primarily harpacticoid copepods and had high (>0.6) diet overlap with permanent resident fish species Favonigobius leteralis and Stigmatopora nigra. Food electivity index indicated that S. punctatus juveniles preferred harpacticoid copepods and amphipods, while juvenile Heteroclinus adelaide and Gymnapistes marmoratus migrating to the nursery habitat targeted larger prey such as amphipods. Preference for larger prey by H. adelaide and G. marmoratus species coupled with differences in prey composition in the stomach was due to mouth size and feeding habits, resulting in different food preferences to S. punctatus. The whiting showed an ontogenetic shift in diet with early settlers (>60?mm TL) consuming less copepods and more amphipods, while previous year recruits (>120?mm TL) consumed polychaete worms. This study indicates that competition for food resources between the new recruiting S. punctatus juveniles and permanent resident juveniles is reduced through differences in temporal and spatial feeding behaviours, mouth morphology, and ontogenetic shift in prey consumption.     

A methodology for regionalizing 3‐D effective porosity at watershed scale in crystalline aquifers

An innovative approach for regionalizing the 3‐D effective porosity field is presented and applied to two large, overexploited, and deeply weathered crystalline aquifers located in southern India. The method derives from earlier work on regionalizing a 2‐D effective porosity field in that part of an aquifer where the water table fluctuates, which is now extended over the entire aquifer using a 3‐D approach. A method based on geological and geophysical surveys has also been developed for mapping the weathering profile layers (saprolite and fractured layers). The method for regionalizing 3‐D effective porosity combines water table fluctuation and groundwater budget techniques at various cell sizes with the use of satellite‐based data (for groundwater abstraction), the structure of the weathering profile, and geostatistical techniques. The approach is presented in detail for the Kudaliar watershed (983 km²) and tested on the 730 km² Anantapur watershed. At watershed scale, the effective porosity of the aquifer ranges from 0.5% to 2% in Kudaliar and between 0.3% and 1% in Anantapur, which agrees with earlier works. Results show that (a) depending on the geology and on the structure of the weathering profile, the vertical distribution of effective porosity can be very different and that the fractured layers in crystalline aquifers are not necessarily characterized by a rapid decrease in effective porosity and (b) that the lateral variations in effective porosity can be larger than the vertical ones. These variations suggest that within a same weathering profile, the density of open fractures and/or degree of weathering in the fractured zone may significantly vary from a place to another. The proposed method provides information on the spatial distribution of effective porosity that is of prime interest in terms of flux and contaminant transport in crystalline aquifers. Implications for mapping groundwater storage and scarcity are also discussed, which should help in improving groundwater resource management strategies.     

The 1-mm brightness temperature of titan

Titan has been observed with the 5-m Hale telescope at an effective wavelength of 1 mm. Adopting a value of 2700 km for the radius of Titan, we find a brightness temperature of 86±12°K at 1 mm. Comparing our results with previous measurements at longer wavelengths, we conclude that the satellite surface is the source of the 1-mm radiation. Since our measured brightness temperature is close to the equilibrium temperature of a blackbody at the distance of Saturn, we believe there is no significant greenhouse effect on Titan.     

Slab behaviour and its surface expression: new insights from gravity modelling in the SE-Carpathians

We use lithosphere-scale gravity models to calculate gravity anomalies resulting from oceanic subduction, continental collision, slab steepening, delamination, and break-off. Local isostasy was assumed for determining vertical movements caused by mass changes related to these tectonic processes. Our results show that subduction is accompanied by basin subsidence on the upper plate caused by the heavy lithospheric root of the subducting slab. The basin evolution goes parallel with the slab evolution and shows considerable modifications when the processes at depth change (slab steepening, delamination, break-off). Characteristic gravity anomaly curves were acquired for the different tectonic scenarios. These curves together with other data (e.g. basin evolution on the upper and the lower plate) were used for the reconstruction of the tectonic evolution of the SE-Carpathians which includes Tertiary subduction and collision followed by slab steepening and delamination.