Groundwater-flow modeling in the Yucatan karstic aquifer, Mexico

The current conceptual model of the unconfined karstic aquifer in the Yucatan Peninsula, Mexico, is that a fresh-water lens floats above denser saline water that penetrates more than 40 km inland. The transmissivity of the aquifer is very high so the hydraulic gradient is very low, ranging from 7–10 mm/km through most of the northern part of the peninsula. The computer modeling program AQUIFER was used to investigate the regional groundwater flow in the aquifer. The karstified zone was modeled using the assumption that it acts hydraulically similar to a granular, porous medium. As part of the calibration, the following hypotheses were tested: (1) karstic features play an important role in the groundwater-flow system; (2) a ring or belt of sinkholes in the area is a manifestation of a zone of high transmissivity that facilitates the channeling of groundwater toward the Gulf of Mexico; and (3) the geologic features in the southern part of Yucatan influence the groundwater-flow system. The model shows that the Sierrita de Ticul fault, in the southwestern part of the study area, acts as a flow barrier and head values decline toward the northeast. The modeling also shows that the regional flow-system dynamics have not been altered despite the large number of pumping wells because the volume of water pumped is small compared with the volume of recharge, and the well-developed karst system of the region has a very high hydraulic conductivity. Electronic Publication     

县市级政府决策信息系统的开发应用:专业GIS 软件SK set V4 .0功能介绍———以贵州沿河土家族自治县资源遥感信息系统为例

SK setV4.0是中国地质科学院岩溶地质研究所GIS中心2001年4月开发成功的、具有自主板权的专业GIS平台软件。它集成了现今较为先进的组件、控件重用，面向对象编程等技术，并以美国Inprise公司的BorladnC^ Builder编译软件为开发平台，提供了全中文菜单环境，不仅具有基本信息和多图层查询、相关分析动态监测、多因子分析、图斑编辑、自动演示、系统维护等功能，而且还具有快速响应的工具条、友好的对话框和较强的纠错机制，很适合缺乏专业GIS人才的县市级政府管理部门使用。文中还以贵州省沿河县资源遥感信息系统为例，对SK setV4.0的功能和使用方法进行了具体的介绍。     

Analogue study of particle segregation in pyroclastic density currents, with implications for the emplacement mechanisms of large ignimbrites

Abstract Analogue flume experiments were conducted to investigate the transport and sedimentation behaviour of turbulent pyroclastic density currents. The experimental currents were scaled approximately to the natural environment in three ways: (1) they were fully turbulent; (2) they had a very wide range of particle sizes and associated Rouse numbers (the ratio of particle settling velocity to effective turbulent eddy velocity in the current); and (3) they contained particles of two different densities. Two sets of surge‐type experiments were conducted in a 5 m long, water‐filled lock‐exchange flume at five different volumetric particle concentrations from 0·6% to 23%. In one set (one‐component experiments), the currents contained just dense particles; in the other set (two‐component experiments), they contained both light and dense particles in equal volume proportions. In both sets of experiments, the population of each component had a log‐normal size distribution. In the two‐component experiments, the size range of the light particle population was selected in order to be in hydrodynamic equivalence with that of the dense particles. Dense particles were normally graded, both vertically and downstream, in the deposits from both sets of experiments. The mass loading (normalized to the initial mass of the suspension) and grain size of the dense component in the deposits decreased with distance from the reservoir and were insensitive to initial total particle concentration in the currents. On the other hand, in the two‐component experiments, the light particles were extremely sensitive to concentration. They were deposited in hydrodynamic equivalence with the dense particles from dilute currents, but were segregated efficiently at concentrations higher than a few per cent. With increasing particle concentration, the large, light particles were carried progressively further down the flume because of buoyancy effects. Deposits from the high‐concentration currents exhibited reverse vertical grading of the large, light particles. Efficient segregation of the light component was observed even if the bulk density of the current was less than that of the light particles. In both sets of experiments, marked inflexions in the rate of downstream decline in mass loading and maximum grain size of the dense component can be attributed to the presence of two different particle settling regimes in the flow: (1) particles with Rouse numbers >2·5, which did not respond to the turbulence and settled rapidly; and (2) particles with Rouse numbers <2·5, which followed the turbulent eddies and settled slowly. The results are applied to the transport and sedimentation dynamics of pyroclastic density currents that generate large, widespread ignimbrites. Field data fail to reveal significant departures from aerodynamic equivalence between pumice and lithic clasts in three such ignimbrites: the particulate loads of some large ignimbrites are transported principally in turbulent suspensions of low concentration. In some ignimbrites, the well‐developed inflexions in curves of maximum lithic (ML) size vs. distance can be attributed to the existence of distinct high and low Rouse number particle settling regimes that mark the transition from an overcharged state to one in which the residual particulate load is transported more effectively by turbulence.     

Anatomy of shelf deltas at the edge of a prograding Eocene shelf margin, Spitsbergen

Abstract Although shelf‐edge deltas are well‐imaged seismic features of Holocene and Pleistocene shelf margins, documented outcrop analogues of these important sand‐prone reservoirs are rare. The facies and stratigraphic architecture of an outcropping shelf‐edge delta system in the Eocene Battfjellet Formation, Spitsbergen, is presented here, as well as the implications of this delta system for the generation of sand‐prone, shelf‐margin clinoforms. The shelf‐edge deltas of the Battfjellet Formation on Litledalsfjellet and Høgsnyta produced a 3–5 × 15 km, shelf edge‐attached, slope apron (70 m of sandstones proximally, tapering to zero on the lower slope). The slope apron consists of distributary channel and mouth‐bar deposits in its shelf‐edge reaches, passing downslope to slope channels/chutes that fed turbiditic lobes and spillover sheets. In the transgressive phase of the slope apron, estuaries developed at the shelf edge, and these also produced minor lobes on the slope. The short‐headed mountainous rivers that drained the adjacent orogenic belt and fed the narrow shelf, and the shelf‐edge position of the discharging deltas, made an appropriate setting for the generation of hyperpycnal turbidity currents on the slope of the shelf margin. The abundance of organic matter and of coal fragments in the slope turbidites is consistent with this notion. Evidence that many of the slope turbidites were generated by sustained turbidity currents that waxed then waned includes the presence of scour surfaces and thick intervals of plane‐parallel laminae within turbidite beds in the slope channels, and thick spillover lobes with repetitive alternations of massive and flat‐laminated intervals. The examined shelf‐edge to slope system, now preserved mainly below the shelf break and dominated by sediment gravity‐flow deposits, has a threefold stratigraphic architecture: a lower, progradational part, in which the clinoforms have a slight downward‐directed trajectory; a thin aggradational zone; and an upper part in which clinoforms backstep up onto the shelf edge. A greatly increased density of erosional channels and chutes marks the regressive‐to‐transgressive turnaround within the slope apron, and this zone becomes an angular unconformity up near the shelf edge. This unconformity, with both subaerial and subaqueous components, is interpreted as a sequence boundary and developed by vigorous sand delivery and bypass across the shelf edge during the time interval of falling relative sea level. The studied shelf‐margin clinoforms accreted mostly during falling stage (sea level below the shelf edge), but the outer shelf later became estuarine as sea level became re‐established above the shelf edge.     

A consideration of the dune:antidune transition in fine gravel

Hydraulic data defining the dune:antidune transition in fine gravel are compared with potential flow theory, and information is drawn from published experiments and field‐based studies. Attention is given to both transitional bedforms and the development of downstream‐migrating antidunes. In the latter case, most data pertain to sand beds and not to gravel. Empirical data provide some weak support for the theoretical notion that the transition occurs at progressively lower Froude numbers at greater relative depths. Although a critical Froude number of 0·84 may reasonably be applied for the beginning of the dune to antidune transformation, lag effects (and a possible depth limitation) ensure that transitional bedforms may persist across a broad range of Froude numbers from 0·5 to 1·8. This latter observation has great relevance for palaeohydraulic estimates derived from outcrop data. Whereas the application of theoretical bedform existence fields, based upon potential flow theory, to fine gravel was previously purely speculative, the addition of experimental and field data to these plots provides a degree of confidence in applying stability theory to practical geological problems. For the first time, laboratory data pertaining to downstream‐migrating gravel antidunes are compared with theory. These bedforms have been reported from certain experimental near‐critical flows above sand or gravel beds, but have been observed infrequently in natural streams. However, there are no detailed studies from natural rivers and only a few contentious identifications from outcrops. Nevertheless, the limited hydraulic data conform to theoretical expectations.     

Heat flow and thermal modeling of the Yinggehai Basin, South China Sea

Geothermal gradients are estimated to vary from 31 to 43 °C/km in the Yinggehai Basin based on 99 temperature data sets compiled from oil well data. Thirty-seven thermal conductivity measurements on core samples were made and the effects of porosity and water saturation were corrected. Thermal conductivities of mudstone and sandstone range from 1.2 to 2.7 W/m K, with a mean of 2.0±0.5 W/m K after approximate correction. Heat flow at six sites in the Yinggehai Basin range from 69 to 86 mW/m², with a mean value of 79±7 mW/m². Thick sediments and high sedimentation rates resulted in a considerable radiogenic contribution, but also depressed the heat flow. Measurements indicate the radiogenic heat production in the sediment is 1.28 μW/m³, which contributes 20% to the surface heat flow. After subtracting radiogenic heat contribution of the sediment, and sedimentation correction, the average basal heat flow from basement is about 86 mW/m².Three stages of extension are recognized in the subsidence history, and a kinematic model is used to study the thermal evolution of the basin since the Cenozoic era. Model results show that the peak value of basal heat flow was getting higher and higher through the Cenozoic. The maximum basal heat flow increased from 65 mW/m² in the first stage to 75 mW/m² in the second stage, and then 90 mW/m² in the third stage. The present temperature field of the lithosphere of the Yinggehai Basin, which is still transient, is the result of the multistage extension, but was primarily associated with the Pliocene extension.     

Physical properties of ultrahigh-pressure metamorphic rocks from the Sulu terrain, eastern central China: implications for the seismic structure at the Donghai (CCSD) drilling site

About 30 samples representing major lithologies of Sulu ultrahigh-pressure (UHP) metamorphic rocks were collected from surface exposures and exploration wells, and compressional (Vp) and shear wave (Vs) velocities and their directional dependence (anisotropy) were determined over a range of constant confining pressures up to 600 MPa and temperatures ranging from 20 to 600 °C. Samples range in composition from acidic to ultramafic. P- and S-wave velocities measured at 600 MPa vary from 5.08 to 8.64 km/s and 2.34 to 4.93 km/s, respectively. Densities are in the range from 2.60 to 3.68 g/cm³. To make a direct tie between seismic measurements (refraction and reflection) and subsurface lithologies, the experimental velocity data (corresponding to shallow depths) were used to calculate velocity profiles for the different lithologies and profiles of reflection coefficients at possible lithologic interfaces across the projected 5000-m Chinese Continental Scientific Drilling Program (CCSD) crustal segment. Comparison of calculated in situ velocities with respective intrinsic velocities suggests that the in situ velocities at shallow depths are lowered by an increased abundance of open microcracks. The strongly reflective zone beneath the Donghai drill site can be explained by the impedance contrasts between the different lithologies. Contacts between eclogite/peridotite and felsic rocks (gt-gneiss, granitic gneiss), in particular, may give rise to strong seismic reflections. In addition, shear-induced (lattice preferred orientation (LPO)-related) seismic anisotropy can increase reflectivity. For the explanation of the high velocity bodies (>6.4 km/s) around 1000 m and below 3200-m depth, large proportions of eclogite/peridotite (about 40 and 30 vol.%, respectively) are needed.     

Strike-slip fault bridge fluid pumping mechanism: insights from field-based palaeostress analysis and numerical modelling

We present a finite-element study of stress perturbation in evolving compressive and extensional strike-slip fault bridges. The results are compared with a fracture study of a compressive bridge at St Donats, UK. Horizontally interbedded calcareous mudstone and bioclastic calcilutite at St Donats have a distinct vertical permeability anisotropy. This sedimentary sequence behaves as a set of horizontal aquifers. The fluid flow in these aquifers is sensitive to mean stress gradients. Paleostress analysis of field fracture data, verified by finite-element modelling, indicates a rotation of σ₁ towards parallelism with boundary faults inside the growing compressive bridge. Boundary faults and bridge faults recorded numerous fluid flow events. The modelled mean stress pattern shows a regional maximum within the bridge and local maxima/minima pairs at boundary fault tips.Finite-element modelling of an extensional bridge indicates that σ₃ rotates towards parallelism with boundary faults. The mean stress pattern is similar to the pattern in compressive bridge but with maxima and minima locations interchanged. The stress patterns are reestablished by each stress build-up preceding the rupturation of the boundary faults throughout the development stages of strike-slip fault bridges. Mean stress gradients developed pre-failure control the fluid flow in fractures of the strike-slip fault system at and after the end of each stress build-up and the fluid flow in boundary faults post-failure. Fracture reactivation and new fracture generation within an evolving bridge is a process consisting of multiple successive events that retain the storage capacity of the bridge. Rupture and sealing of the main bounding-faults is a step-wise process that opens and closes fluid conduits between areas with different pressures.     

Analyzing CO2 emissions mitigation by technology improvement in Central and Eastern Europe

As a legacy of the centrally planned economy, the economies in transition of Central and Eastern Europe (CEE) have a unique potential to reduce their greenhouse gas emissions through the improvement in their high energy intensities. Since much of this `low-hanging fruit' in energy-efficiency improvements can be highly cost-effective, many developed countries facing difficulties in meeting their greenhouse gas (GHG) emission targets domestically are eager to find such opportunities in the CEE region. Therefore, studies analysing the potentials and costs of carbon dioxide reduction through technology improvement in the region have come into the limelight. While there are a few excellent studies in the region aimed at analysing climate change abatement potentials, they all embark on different assumptions, methodologies and boundary conditions. It is hence difficult, if not impossible, to compare and analyse the results of these studies across different authors, countries or time horizons. Consequently, the purpose of this paper is to place four leading studies on GHG mitigation through technology improvement from the CEE region into an internationally comparable framework. Four studies were selected from three countries, Poland, Hungary and Estonia, which are all the results of major national and international efforts to assess costs and potentials of GHG reduction. The paper places their assumptions, methods and final results into a framework which enables policy-makers and project designers to compare these across geographical and technological boundaries. Since other studies from around the globe have been analysed in this framework in the literature, this paper provides a vehicle for the findings of these four studies to be compared to others worldwide. In addition, the paper highlights a few areas where similar studies to be completed in the future in the region may be enhanced by incorporating features used in GHG mitigation research in other parts of the world. This revised version was published online in July 2006 with corrections to the Cover Date.     

Electrons in feldspar II: a consideration of the influence of conduction band-tail states on luminescence processes

 Most natural feldspars contain many charged impurities, and display a range of bond angles, distributed about the ideal. These effects can lead to complications in the structure of the conduction band, giving rise to a tail of energy states (below the high-mobility conduction band) through which electrons can travel, but with reduced mobility: transport through these states is expected to be thermally activated. The purpose of this article is twofold. Firstly, we consider what kind of lattice perturbations could give rise to both localized and extended conduction band-tail states. Secondly, we consider what influence the band tails have on the luminescence properties of feldspar, where electrons travel through the sample prior to recombination. The work highlights the dominant role that 0.04–0.05-eV phonons play in both the luminescence excitation and emission processes of these materials. It also has relevance in the dating of feldspar sediments at elevated temperatures. Received: 11 May 2001 / Accepted: 6 September 2001