Introduction to hydromechanical well tests in fractured rock aquifers

This article introduces hydromechanical well tests as a viable field method for characterizing fractured rock aquifers. These tests involve measuring and analyzing small displacements along with pressure transients. Recent developments in equipment and analyses have simplified hydromechanical well tests, and this article describes initial field results and interpretations during slug and constant-rate pumping tests conducted at a site underlain by fractured biotite gneiss in South Carolina. The field data are characterized by displacements of 0.3 μm to more than 10 μm during head changes up to 10 m. Displacements are a hysteretic function of hydraulic head in the wellbore, with displacements late in a well test always exceeding those at similar wellbore pressures early in the test. Displacement measurements show that hydraulic aperture changes during well tests, and both scaling analyses and field data suggest that T changed by a few percent per meter of drawdown during slug and pumping tests at our field site. Preliminary analyses suggest that displacement data can be used to improve estimates of storativity and to reduce nonuniqueness during hydraulic well tests involving single wells.     

Buruanga peninsula and Antique Range: Two contrasting terranes in Northwest Panay,Philippines featuring an arc–continent collision zone

Alternating chert–clastic sequences juxtaposed with limestone blocks, which are units typical of accretionary complexes, constitute the Buruanga peninsula. New lithostratigraphic units are proposed in this study: the Unidos Formation (Jurassic chert sequence), the Saboncogon Formation (Jurassic siliceous mudstone–terrigenous mudstone and quartz‐rich sandstone), the Gibon Formation (Jurassic(?) bedded pelagic limestone), the Libertad Metamorphics (Jurassic–Cretaceous slate, phyllite, and schist) and the Buruanga Formation (Pliocene–Pleistocene reefal limestone). The first three sedimentary sequences in the Buruanga peninsula show close affinity with the ocean plate stratigraphy of the North Palawan terrane in Busuanga Island: Lower–Middle Jurassic chert sequences overlain by Middle–Upper Jurassic clastics, juxtaposed with pelagic limestone. Moreover, the JR5–JR6 (Callovian to Oxfordian) siliceous mudstone of the Saboncogon Formation in the Buruanga peninsula correlates with the JR5–JR6 siliceous mudstone of the Guinlo Formation in the Middle Busuanga Belt. These findings suggest that the Buruanga peninsula may be part of the North Palawan terrane. The rocks of the Buruanga peninsula completely differ from the Middle Miocene basaltic to andesitic pyroclastic and lava flow deposits with reefal limestone and arkosic sandstone of the Antique Range. Thus, the previously suggested boundary between the Palawan microcontinental block and the Philippine Mobile Belt in the central Philippines, which is the suture zone between the Buruanga peninsula and the Antique Range, is confirmed. This boundary is similarly considered as the collision zone between them.     

The influence of dimensionality on simulations of mass recovery from nonuniform dense non-aqueous phase liquid (DNAPL) source zones

The influence of model dimensionality on predictions of mass recovery from dense non-aqueous phase liquid (DNAPL) source zones in nonuniform permeability fields was investigated using a modified version of the modular three-dimensional transport simulator (MT3DMS). Thirty-two initial two- (2D) and three-dimensional (3D) tetrachloroethene–DNAPL source zone architectures, taken from a recent modeling study, were used as initial conditions for this analysis. Commonly employed source zone metrics were analyzed to determine differences between 2D and 3D predictions: (i) down-gradient flux-averaged contaminant concentration, (ii) reductions in contaminant mass flux through a down-gradient boundary, (iii) source zone ganglia-to-pool (GTP) ratio, and (iv) time required to achieve a remediation objective. 3D flux-averaged contaminant concentrations were approximately 3.5 times lower than concentrations simulated in 2D. This difference was attributed to dilution of the contaminant concentrations down gradient of the source zone. Contaminant flux reduction predictions for a given mass recovery were generally 5% higher in 3D simulations than in 2D simulations. The GTP ratio declined over time as mass was recovered in both 2D and 3D simulations. Although the source longevity (i.e., time required to achieve 99.99% mass recovery) differed between individual 2D and 3D realizations, the mean source longevity for the 2D and 3D simulation ensembles was within 2%. 2D simulations tended to over-predict the time required to achieve lower mass recovery levels (e.g. 50% mass recovery) due to a smaller contaminated area exposed to uncontaminated water. These findings suggest that ensemble averages of 2D numerical simulations of DNAPL migration, entrapment, dissolution, and mass recovery in statistically homogenous, nonuniform media may provide reasonable approximations to average behavior obtained using simulations conducted in fully three-dimensional domains.     

Modelling Watersheds as Spatial Object Hierarchies: Structure and Dynamics

The generation, transport and fate of non-point source pollutants in surface water systems is recognized as a major threat to water supplies, aquatic and coastal ecosystems. The transformation and movement of water, carbon and nutrients through watersheds integrates a set of ecosystem processes along hydrologic flowpaths. Human individual and institutional interactions with these processes involve direct addition or abstraction of these substances, or the alteration of land cover and drainage systems. In natural and developed catchments, these processes often vary at granularities ranging from below the level of a hillslope, up through regional watersheds. This suggests the need for the development of hierarchical analysis tools that can address the integration of a set of biophysical, biogeochemical and socioeconomic processes over a spectrum of scales. We describe and illustrate the use of a watershed model implemented as a spatial object hierarchy, representing successively contained landform classes associated with class specific processes as member functions. The model has been linked in a range of looser and tighter couplings with GRASS and ArcView, supplemented by specific terrain analytical functions. We illustrate the data and model system for an instrumented catchment monitored as part of the Baltimore Ecosystem Study (BES), a Long Term Ecological Research (LTER) site centering on integrated carbon, water and nutrient cycling.     

Profiles of Fecal Porphyrins in River Otters Following the Exxon Valdez Oil Spill

Median levels of Coproporphyrin III (Copro III) in fecal samples of river otters (Lontra canadensis) collected from an oiled area in Prince William Sound, Alaska, USA, during 1990 were significantly higher than in samples collected from the same oiled area during 1996 (p=0.011, one way analysis of variance), a nonoiled reference area in Prince William Sound during 1996 (p=0.002) and a reference area in southeast Alaska during 1998 (p=0.004). An overall test of significance that combined probabilities from the statistical analysis of this porphyrin study with those from other biomarker studies revealed a significant difference in physiological response of river otters between oiled and nonoiled areas of the Sound for 1990 (p < 0.01). We demonstrated that changes in levels of fecal porphyrins may serve as a biomarker that may contribute to a health assessment of wild river otters.     

Velocity predictions for shoaling and breaking waves with a Boussinesq-type model

In this paper, the performance of a 1-D Boussinesq model is evaluated against laboratory data for its ability to predict surf zone velocity moments. Wave evolution over a plane beach and a complex bathymetry both extending into the surf-zone is examined for six cases. For the plane beach, these comprise two cases, a spilling and a plunging cnoidal wave. For the complex bathymetry, these comprise four cases of longer and short wavelengths (spilling and plunging breakers), with regular and irregular periodicity. The model evaluation places emphasis on parameters of the wave field that could be used for the prediction of sediment transport; orbital velocity, undertow, velocity skewness, kurtosis and asymmetry. It is found that, despite an overestimation of the depth-averaged horizontal velocity in the regular waves cases, the predicted higher order velocity moments and undertow are in good agreement with the laboratory data. A bispectral analysis demonstrates that the nonlinear transfers of energy amongst the low order harmonics are well reproduced, but energy exchanges with the higher harmonics are less well predicted. As a result, the model handles velocity moments better in the shorter wave tests than in the long wave cases where triad interactions are stronger. Of the four parameters describing wave breaking, the model behaviour is most sensitive to the critical wave front slope φ_B, especially with regard to velocity skewness and kurtosis predictions. It is also found that increasing the thickness of the surface roller for the case of plunging breakers improves the model's performance.     

Summer Monsoon Rainfall Variability in Central China over the Past 4700 Years and Its Possible Link to Solar Activity

Based on 467 pairs of δ~(18) O and δ~(13) C records and 8 230 Th dates from a stalagmite(BF4) from Xiniu Cave, central China, we present a reconstruction of ~9-yr resolution monsoon rainfall record for the past 4700 years. Our δ~(18) O record shows good coherence with East Asian summer monsoon(EASM) rainfall proxies from adjacent regions during the overlapping intervals, suggesting that δ~(18) O signal in BF4 can be interpreted as a monsoon rainfall proxy. The δ~(13) C variations are related to changes in local processes at the cave site, and regional rainfall and temperature changes.Based on the δ~(18) O record, a series of dry periods can be identified at 4500–4200, 3500–3200, 2800–2500,1900–1600, 1400–1200, 700–500, and 400–200 yr BP, while a series of wet periods can be identified at 4200–3600,2400–2200, 3200–2800, 1100–900, 600–400, and 200–100 yr BP. Power spectrum analysis on our δ~(18) O record reveals significant cycles at ~470 and ~80 yr, coinciding with the typical solar periodic variations. This result suggests that changes in solar activity play a dominant role in driving centennial–decadal monsoon rainfall variation in central China. Due to minor changes in solar irradiance(less than 1.5 W m~(-2)) over the past 4700 years, our record was further compared to the El Ni?o–Southern Oscillation(ENSO) and Pacific Decadal Oscillation(PDO) proxies,confirming that solar forcing on monsoon rainfall changes might be amplified by the ENSO and PDO variations.From 600 to 150 yr BP(the Little Ice Age, LIA), a positive shift of ~2‰ can be revealed in both the δ~(18) O and δ~(13) C records, indicating a cold/dry climatic pattern. By comparing our δ~(18) O and δ~(13) C records with historical documents,we suggest that the climatic deteriorations between 450 and 250 yr BP may have caused serious social unrest at the end of the Ming Dynasty.     

南京直立人的U/Th和U/Pa年代

中国直立人化石的准确定年对于研究人类演化有着极为重要的意义。1993～1994年在南京汤山葫芦洞发现的两个直立人头盖骨化石和一枚牙化石被称为“南京直立人”。其中1号头盖骨化石之上的方解石钙板的U/Th年龄为53.3_(-1.2)~(+1.5)万年;但考虑到定年的准确性,则为53.3_(-3.0)~(+3.5)万年。其~(231)Pa/~(235)U活度比值为0.998±0.006。这表明“南京直立人”的年代应该大于50万年。与”南京直立人”伴生的动物牙化石U/Th年代为18.5～29.0万年;U/Pa年代为13.7～17.2万年。此外,对于同一颗牙化石,牙釉的年龄小于牙本质的年龄。同一样品的U/Pa年龄也显著小于其U/Th年龄。因此,牙化石的U摄取过程并不符合U早期摄取模式。多数牙化石分析点在~(234)U/~(238)U-~(230)Th/~(238)U图上落在位于U早期摄取和线性摄取模式曲线之间,指示牙化石的U摄取过程很可能介于上述两种模式之间。如果这一假设成立,那么牙化石的U/Th和U/Pa线性摄取模式年龄则为其年代的上限。因为不受U摄取过程~(234)U/~(238)U变化的影响,U/Pa线性摄取模式年龄比U/Th较为可靠。最小的U/Pa线性摄取模式年龄为1Ma,这是”南京直立人”上限年龄的估计。从定年结果看,”南京直立人”可能生活在海洋同位素(MIS)16阶段,但这不是最终结论。     

石笋氧同位素和微量元素记录的陕南地区4200～2000a B.P.高分辨率季风降雨变化

基于陕南祥龙洞石笋XL2的19个230Th年龄、218个氧同位素分析以及896个Sr/Ca分析数据,高分辨率重建了4200～1972a B.P.期间陕南地区季风降雨变化.重建结果显示陕南地区这一时期季风降雨有显著的127～105a和57a周期,可能分别受控于太阳活动、PDO和/或AMO的变化.重建时段有3次百年尺度的干旱事件,分别发生于2200～2100a B.P.,2900～2700a B.P.和3600～3400a B.P.,其中2900～2700a B.P.干旱事件对应于北大西洋地区2.8ka冷事件.对比研究显示,尽管祥龙洞石笋和董哥洞石笋δ18O记录整体一致,但除了2900～ 2700aB.P.干旱事件之外,其他两次干旱事件在董哥洞石笋记录中并不明显.而尽管总体上祥龙洞和和尚洞石笋δ18O记录的差异要大,但XL2的3次干旱事件在和尚洞记录都有明显体现.有精确年代控制的祥龙洞、董哥洞及和尚洞石笋氧同位素记录的差异,揭示晚全新世我国季风降雨在十一百年尺度存在区域差异.     

Accretionary rims on inclusions in the Allende meteorite

Many inclusions in Allende, particularly those with irregular shapes, are surrounded by a sequence of thin layers which differ from one another in texture, mineralogy and mineral-chemistry. The layer underlying all others contains either: IA, pyroxene needles + olivine + clumps of hedenbergite and andradite; IB, olivine doughnuts; or IC, rectangular olivine crystals. The next layer outward, II, contains tiny (<5 μm) olivine plates and Layer III large (5–10 μm) olivine laths. The final layer, IV, occurs as clumps of andradite + hedenbergite surrounded by magnesium-rich pyroxene needles. It separates Layer III from the Allende matrix which is more poorly sorted and more sulfide-rich than Layer III. Nepheline and iron sulfide are common constituents of most layers, the latter being particularly fine-grained and abundant in Layer II. Although not every layer is present on every inclusion, the sequence of layers is constant. Evidence that the rims are accretionary aggregates includes the presence of highly disequilibrium mineral assemblages and the fact that they are highly porous masses consisting of many euhedral crystals with few intergrowths. In addition, the layers are thickest in topographic hollows on the surfaces of inclusions and the inner layers are absent or discontinuous beyond such irregularities, suggesting that the probability of accretion of crystals was low initially, except in pockets, and became greater later, after a soft cushion of accreted condensate crystals had already formed. Separation of assemblages of different mineralogy, mineral-chemistry and texture into different rim layers seems best explained by nebular models in which long, slow cooling histories allow differentiation during condensation by grain/gas separation processes.