Faulting, fracturing and in situ stress prediction in the Ahnet Basin, Algeria — a finite element approach

Many low-efficiency hydrocarbon reservoirs are productive largely because effective reservoir permeability is controlled by faults and natural fractures. Accurate and low-cost information on basic fault and fracture properties, orientation in particular, is critical in reducing well costs and increasing well recoveries. This paper describes how we used an advanced numerical modelling technique, the finite element method (FEM), to compute site-specific in situ stresses and rock deformation and to predict fracture attributes as a function of material properties, structural position and tectonic stress. Presented are the numerical results of two-dimensional, plane-strain end-member FEM models of a hydrocarbon-bearing fault-propagation-fold structure. Interpretation of the modelling results remains qualitative because of the intrinsic limitations of numerical modelling; however, it still allows comparisons with (the little available) geological and geophysical data.

In all models, the weak mechanical strength and flow properties of a thick shale layer (the main seal) leads to a decoupling of the structural deformation of the shallower sediments from the underlying sediments and basement, and results in flexural slip across the shale layer. All models predict rock fracturing to initiate at the surface and to expand with depth under increasing horizontal tectonic compression. The stress regime for the formation of new fractures changes from compressional to shear with depth. If pre-existing fractures exist, only (sub)horizontal fractures are predicted to open, thus defining the principal orientation of effective reservoir permeability. In models that do not include a blind thrust fault in the basement, flexural amplification of the initial fold structure generates additional fracturing in the crest of the anticline controlled by the material properties of the rocks. The folding-induced fracturing expands laterally along the stratigraphic boundaries under enhanced tectonic loading. Models incorporating a blind thrust fault correctly predict the formation of secondary syn- and anti-thetic mesoscale faults in the basement and sediments of the hanging wall. Some of these faults cut reservoir and/or seal layers, and thus may influence effective reservoir permeability and affect seal integrity. The predicted faults divide the sediments across the anticline in several compartments with different stress levels and different rock failure (and proximity to failure). These numerical model outcomes can assist classic interpretation of seismic and well bore data in search of fractured and overpressured hydrocarbon reservoirs.     

Mammoth tooth enamel growth rates inferred from stable isotope analysis and histology

Mammoth (Mammuthus sp.) teeth are relatively abundant in Quaternary deposits from Eurasia and North America, and their isotopic compositions can be used to reconstruct past seasonal patterns in precipitation, diet, and migration. Strategies for collecting and interpreting such data, however, are strongly dependent on growth rates, which can vary among species, individuals, and within teeth. In this study, we use histological and isotopic measurements to determine enamel growth rates for a Columbian mammoth (Mammuthus columbi) tooth in two directions. Using histology, the growth rate through the enamel thickness (ET; perpendicular to the height of the tooth) is estimated at 0.8 to 1.5 mm/yr. Isotopic sampling through the innermost 0.36 mm of the ET recovered less than half a period of variation (i.e., half an inferred year of growth), which is consistent with the histological estimate for ET growth rate. A combination of histological and isotopic measurements suggests that the enamel extension rate (growth in the height of the tooth) is 13–14 mm/yr. Knowledge of enamel growth rates should improve the design and interpretation of future isotopic studies of mammoth teeth. The combination of histological and isotopic measurements may also prove useful in determining growth rates for other extinct taxa.     

Use of composite data sets for source-tracking enterococci in the water column and shoreline interstitial waters on Pensacola Beach, Florida

Sources of Enterococcus faecalis isolates from Pensacola Beach, FL. were identified using a library-based approach by applying the statistical method of average similarity to single and composite data sets generated from separate analyses. Data sets included antibiotic resistance analysis (ARA), rep-fingerprints, and fatty acid methyl ester (FAME) profiles. Use of a composite data set composed of ARA and rep-fingerprints, added to the confidence of the identifications. The addition of FAME data to composite data sets did not add to the confidence of identifications. Source identification was performed to better understand risk associated with higher densities of enterococci found in swash zone interstitial water (SZIW) as compared to adjacent bathing water on Pensacola Beach, FL. The "swash zone" is that area of the beach continually washed over by waves. As the potential sources of enterococci were limited in this environment, only two library units, sea gull and human, were constructed. Identification of the beach isolates using a composite data set indicated a sea gull origin. The clonality of the beach isolates suggested that the beach environment selects certain subspecies of E. faecalis.     

Scalar-Flux Similarity in the Layer Near the Surface Over Mountainous Terrain

The scaled standard deviations of temperature and humidity are investigated in complex terrain. The study area is a steep Alpine valley, with six measurement sites of different slope, orientation and roughness (i-Box experimental site, Inn Valley, Austria). Examined here are several assumptions forming the basis of Monin–Obukhov similarity theory (MOST), including constant turbulence fluxes with height and the degree of self-correlation between the involved turbulence variables. Since the basic assumptions for the applicability of the MOST approach—horizontally homogeneous and flat conditions—are violated, the analysis is performed based on a local similarity hypothesis. The scaled standard deviations as a function of local stability are compared with previous studies from horizontally homogeneous and flat terrain, horizontally inhomogeneous and flat terrain, weakly inhomogeneous and flat terrain, as well as complex terrain. As a reference, similarity relations for unstable and stable conditions are evaluated using turbulence data from the weakly inhomogeneous and flat terrain of the Cabauw experimental site in the Netherlands, and assessed with the same post-processing method as the i-Box data. Significant differences from the reference curve and also among the i-Box sites are noted, especially for data derived from the i-Box sites with steep slopes. These differences concern the slope and the magnitude of the best-fit curves, illustrating the site dependence of any similarity theory.     

Derivation of Structure Parameters of Temperature and Humidity in the Convective Boundary Layer from Large-Eddy Simulations and Implications for the Interpretation of Scintillometer Observations

We derive the turbulent structure parameters of temperature $C_{T}^2$ and humidity $C_q^2$ from high-resolution large-eddy simulations (LES) of a homogeneously-heated convective boundary layer. Boundary conditions and model forcing were derived from measurements at Cabauw in The Netherlands. Three different methods to obtain the structure-parameters from LES are investigated. The shape of the vertical structure-parameter profiles from all three methods compare well with former experimental and LES results. Depending on the method, deviations in the magnitude up to a factor of two are found and traced back to the effects of discretization and numerical dissipation of the advection scheme. Furthermore, we validate the LES data with airborne and large-aperture scintillometer (LAS) measurements at Cabauw. Virtual path measurements are used to study the variability of $C_{T}^2$ in the mixed layer and surface layer and its implications for airborne and LAS measurements. A high variability of $C_{T}^2$ along a given horizontal path in the LES data is associated with plumes (high values) and downdrafts (low values). The path average of $C_{T}^2$ varies rapidly in time due to the limited path length. The LES results suggest that measured path averages require sufficient temporal averaging and an adequate ratio of path length to height above the ground for the LAS in order to approach the domain average of $C_{T}^2$ .     

Simulated impacts of the South Atlantic Ocean Dipole on summer precipitation at the Guinea Coast

An intermediate complexity atmospheric general circulation model has been used to investigate the influence of the South Atlantic Ocean (SAO) dipole (SAOD) on summer precipitation over the Guinea Coast of West Africa. Two ensemble integrations in which idealized but realistic SAOD-type sea surface temperature (SST) anomaly is prescribed only in the SAO, and then globally are performed and inter-compared. Consistently, above (below) the average precipitation is simulated over the Guinea Coast during the positive (negative) phase of the SAOD. Comparison of the two set of experiments reveal that in its active years, the SAOD is a dominant mechanism that shapes the spatial character of summer precipitation at the Guinea coast, the global SST variability merely slightly moderate its effects. During the SAOD, cool SST anomaly in the extra-tropical SAO off the Brazil–Uruguay–Argentina coast gives rise to suppressed convection and mass divergence. In turn, the subsidence tends to amplify the sub-tropical arm of anomalous Hadley-type circulation and consequently large scale convection and mass flux convergence in the equatorial Atlantic Ocean/Gulf of Guinea region bordering on the coastal fringes of West Africa. Precipitation is therefore increased at the Guinea Coast.     

Are peatlands in different states with respect to their thermodynamic behaviour? A simple test of peatland energy and entropy budgets

Whilst all ecosystems must obey the second law of thermodynamics, these physical bounds and controls on ecosystem evolution and development are largely ignored across the ecohydrological literature. To unravel the importance of these underlying restraints on ecosystem form and function, and their power to inform our scientific understanding, we have calculated the entropy budget of a range of peat ecosystems. We hypothesize that less disturbed peatlands are ‘near equilibrium’ with respect to the second law of thermodynamics and thus respond to change by minimizing entropy production. This ‘near equilibrium’ state is best achieved by limiting evaporative losses. Alternatively, peatlands ‘far-from-equilibrium’ respond to a change in energy inputs by maximizing entropy production which is best achieved by increasing evapotranspiration. To test these alternatives this study examined the energy balance time series from seven peatlands across a disturbance gradient. We estimate the entropy budgets for each and determine how a change in net radiation (ΔR_n) was transferred to a change in latent heat flux (ΔλE). The study showed that: (i) The transfer of net radiation to latent heat differed significantly between peatlands. One group transferred up to 64% of the change in net radiation to a change in latent heat flux, while the second transferred as little as 27%. (ii) Sites that transferred the most energy to latent heat flux were those that produced the greatest entropy. The study shows that an ecosystem could be ‘near equilibrium’ rather than ‘far from equilibrium’.     

Implications of the relationship between catchment vegetation type and the variability of annual runoff

The impact of changing catchment vegetation type on mean annual runoff has been known for some time, however, the impact on the variability of annual runoff has been established only recently. Differences in annual actual evapotranspiration between vegetation types and the potential effect of changing vegetation type on mean annual runoff and the variability of annual runoff are briefly reviewed. The magnitude of any change in the variability of annual runoff owing to a change in catchment vegetation type is related to the pre‐ and post‐change vegetation types and the seasonality of precipitation, assuming that the variability of annual precipitation remains constant throughout. Significant implications of the relationship between vegetation type and the variability of annual runoff are presented and discussed for water resource management, stream ecology and fluvial geomorphology. Copyright © 2002 John Wiley & Sons, Ltd.     

Assessment of potential transport of pollutants into the Barents Sea via sea ice--an observational approach

The present estimates of ice drift in the Arctic include utilization of satellite imagery data (special sensor microwave/imager) and a reconstruction of air pressure for the period 1899-1998. A significant part of the sea ice in the Arctic Ocean has its origin in the Kara Sea and melts in the Greenland and the Barents Sea (BS). Consequently there may be a particular risk of pollutants in the Kara Sea entering the food webs of the Greenland and BS. The ice export from the Kara Sea between 1988 and 1994 was about 208,000 km2 (154 km3) per year. The import of ice into the BS was during the same period 161,000 km2 (183 km3) per year while the ice drift through the Fram Strait into the Greenland Sea was 583,000 km2 (1859 km3) per year. Ice which formed adjacent to the Ob and Yenisey rivers in early January, drifted into the BS within two years (with a probability of about 50%.