Evolution of an organic‐rich lake basin – stratigraphy,climate and tectonics: Piceance Creek basin,Eocene Green River Formation

The Piceance Creek basin formed as a continental foreland basin ca 53 to 48 Ma in the early to middle Eocene. On a global basis, the basin contains one of the richest oil shale resources known, where the profundal oil shale deposits, kerogen‐rich mudstones (clay and carbonate), exist over most of the basin. Despite its economic importance, the evolution of the Piceance Creek basin is still somewhat unclear. Based on facies association analysis, depositional trends, and gamma ray and Fischer assay data, six evolutionary lake stages are recognized: (i) fresh lake; (ii) transitional lake; (iii) highly fluctuating lake; (iv) rising lake; (v) high lake; and (vi) closing lake. Lake stages are composed of depositional units and characterize large‐scale changes in sedimentological patterns, depositional trends and fluctuations in the oil shale richness related to changes in climate and tectonics. Lake stage evolution is also consistent with the global Eocene climate trend. Stage 1 formed prior to the Eocene climate optimum. At the beginning of the Eocene climate optimum, a saline‐restricted lake formed (Stage 2) and evolved into the highly fluctuating lake (Stage 3) indicating rapid climate changes during the peak of the Eocene climate optimum. This stage was followed by the rising and high lakes (Stages 4 and 5) after the climate optimum and during a change to a more humid climate. The closing of the lake (Stage 6) was caused by increased sand input from the north, indicating the influence of both tectonics and climate. Based on depositional trends and climate evolution, it is suggested that, during the arid climate, laterally heterogeneous highly cyclic depositional units dominate, whereas, during the humid climate, depositional units form laterally continuous sediments that can be traced over long distances.     

Response of unconfined turbidity current to deep‐water fold and thrust belt topography: Orthogonal incidence on solitary and segmented folds

Sea‐floor topography of deep‐water folds is widely considered to have a major impact on turbidity currents and their depositional systems, but understanding the flow response to such features was limited mainly to conceptual notions inspired by small‐scale laboratory experiments. High‐resolution three‐dimensional numerical experiments can compensate for the lack of natural‐scale flow observations. The present study combines numerical modelling of thrusts with fault‐propagation folds by Trishear3D software with computational fluid dynamics simulations of a natural‐scale unconfined turbidity current by MassFlow‐3D? software. The study reveals the hydraulic and depositional responses of a turbidity current (ca 50 m thick) to typical topographic features that it might encounter in an orthogonal incidence on a sea‐floor deep‐water fold and thrust belt. The supercritical current (ca 10 m sec^?1) decelerated and thickened due to the hydraulic jump on the fold backlimb counter‐slope, where a reverse overflow formed through current self‐reflection and a reverse underflow was issued by backward squeezing of a dense near‐bed sediment load. The reverse flows were re‐feeding sediment to the parental current, reducing its waning rate and extending its runout. The low‐efficiency current, carrying sand and silt, outran a downslope distance of >17 km with only modest deposition (<0·2 m) beyond the fold. Most of the flow volume diverted sideways along the backlimb to surround the fold and spread further downslope, with some overspill across the fold and another hydraulic jump at the forelimb toe. In the case of a segmented fold, a large part of the flow went downslope through the segment boundary. Preferential deposition (0·2 to 1·8 m) occurred on the fold backlimb and directly upslope, and on the forelimb slope in the case of a smaller fold. The spatial patterns of sand entrapment revealed by the study may serve as guidelines for assessing the influence of substrate folds on turbiditic sedimentation in a basin.     

Electrical Resistivity Characterization of Anisotropy in the Biscayne Aquifer

Electrical anisotropy occurs when electric current flow varies with azimuth. In porous media, this may correspond to anisotropy in the hydraulic conductivity resulting from sedimentary fabric, fractures, or dissolution. In this study, a 28‐electrode resistivity imaging system was used to investigate electrical anisotropy at 13 sites in the Biscayne Aquifer of SE Florida using the rotated square array method. The measured coefficient of electrical anisotropy generally ranged from 1.01 to 1.12 with values as high as 1.36 found at one site. The observed electrical anisotropy was used to estimate hydraulic anisotropy (ratio of maximum to minimum hydraulic conductivity) which ranged from 1.18 to 2.83. The largest values generally were located on the Atlantic Coastal Ridge while the lowest values were in low elevation areas on the margin of the Everglades to the west. The higher values of anisotropy found on the ridge may be due to increased dissolution rates of the oolitic facies of the Miami formation limestone compared with the bryozoan facies to the west. The predominate trend of minimum resistivity and maximum hydraulic conductivity was E‐W/SE‐NW beneath the ridge and E‐W/SW‐NE farther west. The anisotropy directions are similar to the predevelopment groundwater flow direction as indicated in published studies. This suggests that the observed anisotropy is related to the paleo‐groundwater flow in the Biscayne Aquifer.     

Rating curve uncertainty and change detection in discharge time series: case study with 44‐year historic data from the Nyangores River,Kenya

The intersection of the developing topic of rating curve and discharge series uncertainty with the topic of hydrological change detection (e.g., in response to land cover or climatic change) has not yet been well studied. The work herein explores this intersection, with consideration of a long‐term discharge response (1964–2007) for a ~650‐km² headwater basin of the Mara River in west Kenya, starting with stream rating and daily gauge height data. A rating model was calibrated using Bayesian methods to quantify uncertainty intervals in model parameters and predictions. There was an unknown balance of random and systemic error in rating data scatter (a scenario not likely unique to this basin), which led to an unknown balance of noise and information in the calibrated statistical error model. This had implications on testing for hydrological change. Overall, indications were that shifts in basin's discharge response were rather subtle over the 44‐year period. A null hypothesis for change using flow duration curves (FDCs) from four different 8‐year data intervals could be either accepted or rejected over much of the net flow domain depending on different applications of the statistical error model (each with precedence in the literature). The only unambiguous indication of change in FDC comparisons appeared to be a reduction in lowest baseflow in recent years (flows with >98% exceedance probability). We defined a subjective uncertainty interval based on an intermediate balance of random and systematic error in the rating model that suggested a possibility of more prevalent impacts. These results have relevance to management in the Mara basin and to future studies that might establish linkages to historic land use and climatic factors. The concern about uncertain uncertainty intervals (uncertainty²) extends beyond the Mara and is relevant to testing change where non‐random rating errors may be important and subtle responses are investigated. Copyright © 2013 John Wiley & Sons, Ltd.     

Conceptualization of the hydrogeological system of some sedimentary aquifers in Savelugu–Nanton and surrounding areas,Northern Ghana

Groundwater resources play a pivotal role in the rural water delivery system in Ghana. The hydrogeological system of Middle Voltaian terrain was simulated using available data on hydraulic heads and boundary conditions. The objective was to characterize the general groundwater flow pattern and provide local estimates of the distribution of hydraulic conductivity and recharge fields. The results suggest a predominant NE–SW flow direction, which ties in with the general regional structural trend and indicates that the hydrogeological conditions of the rocks are controlled by structural entities created in the wake of fracturing and/or weathering of the rocks whose primary permeabilities are considerably reduced because of high compaction and low‐grade metamorphism. Calibrated hydraulic conductivities range between 1.90 and 10.81 m/d. The spatial distribution appears to reflect the intensity of fracturing and/or weathering of the rock and the proportion of the clay fraction of the weathered zone. Vertical groundwater recharge has been estimated to range between 0.3% and 4.1% of the annual rainfall. This recharge rate is quite low and reflects the imperviousness of the thick overburden because of high clay content in some places and high compaction in others. Despite this apparently low recharge rate, groundwater resources potential in the area appear to be high, and increased abstraction from existing abstraction wells by up to 50% does not appear to register significant effects on groundwater budgets at the simulated recharge rates. This suggests that the well yields are much lower than the potential of the aquifer system. The apparently low yields might be associated with poor well development and the choice of inappropriate well completion materials. This study recommends a monitoring system to be developed for a much more regional groundwater flow simulation under transient conditions. Copyright © 2012 John Wiley & Sons, Ltd.     

Characterizing the effect of elastic interactions on the effective elastic properties of porous,cracked rocks

Elastic interactions between pores and cracks reflect how they are organized or spatially distributed in porous rocks. The principle goal of this paper is to understand and characterize the effect of elastic interactions on the effective elastic properties. We perform finite element modelling to quantitatively study how the spatial arrangement of inclusions affects stress distribution and the resulting overall elasticity. It is found that the stress field can be significantly altered by elastic interactions. Compared with a non‐interacting situation, stress shielding considerably stiffens the effective media, while stress amplification appreciably reduces the effective elasticity. We also demonstrate that the T‐matrix approach, which takes into account the ellipsoid distribution of pores or cracks, can successfully characterize the competing effects between stress shielding and stress amplification. Numerical results suggest that, when the concentrations of cracks increase beyond the dilute limit, the single parameter crack density is not sufficient to characterize the contribution of the cracks to the effective elasticity. In order to obtain more reliable and accurate predictions for the effective elastic responses and seismic anisotropies, the spatial distribution of pores and cracks should be included. Additionally, such elastic interaction effects are also dependent on both the pore shapes and the fluid infill.     

Retrieval of suspended sediment concentrations in large turbid rivers using Landsat ETM+: an example from the Yangtze River,China

Suspended sediment concentration (SSC) is a critical parameter in the study of river sediment transport and water quality variation, but traditional measurement methods are costly and time‐consuming. This paper is focused on presenting a methodology that may be useful in estimating SSC which is of key importance in process geomorphology and hydrology. In previous studies, remote sensing has been applied to estimate the SSC of sea waters as well as low turbid inland waters like lakes, reservoirs and short river reaches visible within a single Landsat satellite image coverage. Rivers, especially highly turbid large rivers, have largely been ignored. The dataset used in this paper includes measured SSC and multi‐temporal Landsat ETM+ images covering most part of the Yangtze River. Using an effective easy‐to‐use atmospheric correction method that does not require in situ atmospheric conditions, retrieved water reflectance of Band 4 was found to be a good SSC indicator within the large SSC range 22–2610 mg l^–1. The newly developed regression relation between SSC and water reflectance of Band 4 appears to be able to provide a relatively accurate SSC estimate directly from Landsat ETM+ images for the Yangtze River from the upper, the middle to the lower reaches. With the relation it is possible to estimate or map out SSC dynamics of large rivers which lack SSC data due to constraints of conventional measurements. Copyright © 2009 John Wiley & Sons, Ltd.     

Comparison of groundwater age models for assessing nitrate loading,transport pathways,and management options in a complex aquifer system

In an aquifer system with complex hydrogeology, mixing of groundwater with different ages could occur associated with various flow pathways. In this study, we applied different groundwater age‐estimation techniques (lumped parameter model and numerical model) to characterize groundwater age distributions and the major pathways of nitrate contamination in the Gosan agricultural field, Jeju Island. According to the lumped parameter model, groundwater age in the study area could be explained by the binary mixing of the young groundwater (4–33 years) and the old water component (>60 years). The complex hydrogeologic regimes and local heterogeneity observed in the study area (multilayered aquifer, well leakage hydraulics) were particularly well reflected in the numerical model. The numerical model predicted that the regional aquifer of Gosan responded to the fertilizer applications more rapidly (mean age: 9.7–22.3 years) than as estimated by other models. Our study results demonstrated that application and comparison of multiple age‐estimation methods can be useful to understand better the flow regimes and the mixing characteristics of groundwater with different ages (pathways), and accordingly, to reduce the risk of improper groundwater management plans arising from the aquifer heterogeneity.