A New GIS Approach for Estimating Shallow Groundwater Recharge and Discharge

Groundwater recharge and discharge define the relationships among groundwater storage, precipitation, evapotranspiration, and surface water. Understanding and quantifying groundwater recharge and discharge is crucial for developing effective management options for water supply. Traditional recharge and discharge studies usually require time-consuming (months or years) and labor-intensive field characterization. Therefore, a fast initial recharge and discharge estimation will help scientists and decision makers design more accurate and cost-effective research plans and management strategies. A new GIS based methodology along with an ArcMap plug-in package, PRO-GRADE, was developed to help hydrogeologists estimate recharge and discharge in a more efficient way than conventional methods. The Groundwater Recharge And Discharge Estimator (GRADE-GIS) in the PRO-GRADE package uses a mass balance approach for steady state two-dimensional unconfined aquifers. GRADE-GIS only requires relatively short preparation time (hours) using readily available data for water table and bedrock elevations, and hydraulic conductivities. The Pattern Recognition Organizer (PRO-GIS) in the PRO-GRADE package uses multiple image processing algorithms to estimate and visualize shallow recharge and discharge patterns and rates with GIS. The output of PRO-GRADE can be used as an initial guideline for conventional field study planning and decision making, or an initial condition for numerical simulations.     

Geostatistical improvements of evapotranspiration spatial information using satellite land surface and weather stations data

The objective of the present study was to use the simple cokriging methodology to characterize the spatial variability of Penman–Monteith reference evapotranspiration and Thornthwaite potential evapotranspiration methods based on Moderate Resolution Imaging Spetroradiometer (MODIS) global evapotranspiration products and high-resolution surfaces of WordClim temperature and precipitation data. The climatic element data referred to 39 National Institute of Meteorology climatic stations located in Minas Gerais state, Brazil and surrounding states. The use of geostatistics and simple cokriging technique enabled the characterization of the spatial variability of the evapotranspiration providing uncertainty information on the spatial prediction pattern. Evapotranspiration and precipitation surfaces were implemented for the climatic classification in Minas Gerais. Multivariate geostatistical determined improvements of evapotranspiration spatial information. The regions in the south of Minas Gerais derived from the moisture index estimated with the MODIS evapotranspiration (2000–2010), presented divergence of humid conditions when compared to the moisture index derived from the simple kriged and cokriged evapotranspiration (1961–1990), indicating climate change in this region. There was stronger pattern of crossed covariance between evapotranspiration and precipitation rather than temperature, indicating that trends in precipitation could be one of the main external drivers of the evapotranspiration in Minas Gerais state, Brazil.     

Monte Carlo simulation of permeability fields and reservoir performance predictions with SVD parameterization in RML compared with EnKF

In a previous paper, we developed a theoretical basis for parameterization of reservoir model parameters based on truncated singular value decomposition (SVD) of the dimensionless sensitivity matrix. Two gradient-based algorithms based on truncated SVD were developed for history matching. In general, the best of these “SVD” algorithms requires on the order of 1/2 the number of equivalent reservoir simulation runs that are required by the limited memory Broyden–Fletcher–Goldfarb–Shanno (LBFGS) algorithm. In this work, we show that when combining SVD parameterization with the randomized maximum likelihood method, we can achieve significant additional computational savings by history matching all models simultaneously using a SVD parameterization based on a particular sensitivity matrix at each iteration. We present two new algorithms based on this idea, one which relies only on updating the SVD parameterization at each iteration and one which combines an inner iteration based on an adjoint gradient where during the inner iteration the truncated SVD parameterization does not vary. Results generated with our algorithms are compared with results obtained from the ensemble Kalman filter (EnKF). Finally, we show that by combining EnKF with the SVD-algorithm, we can improve the reliability of EnKF estimates.     

Combining sensitivities and prior information for covariance localization in the ensemble Kalman filter for petroleum reservoir applications

Sampling errors can severely degrade the reliability of estimates of conditional means and uncertainty quantification obtained by the application of the ensemble Kalman filter (EnKF) for data assimilation. A standard recommendation for reducing the spurious correlations and loss of variance due to sampling errors is to use covariance localization. In distance-based localization, the prior (forecast) covariance matrix at each data assimilation step is replaced with the Schur product of a correlation matrix with compact support and the forecast covariance matrix. The most important decision to be made in this localization procedure is the choice of the critical length(s) used to generate this correlation matrix. Here, we give a simple argument that the appropriate choice of critical length(s) should be based both on the underlying principal correlation length(s) of the geological model and the range of the sensitivity matrices. Based on this result, we implement a procedure for covariance localization and demonstrate with a set of distinctive reservoir history-matching examples that this procedure yields improved results over the standard EnKF implementation and over covariance localization with other choices of critical length.     

Plastic Pollution at a Sea Turtle Conservation Area in NE Brazil: Contrasting Developed and Undeveloped Beaches

Sea turtles are highly susceptible to plastic ingestion and entanglement. Beach debris were surveyed along the most important sea turtle nesting beaches in Brazil (Costa dos Coqueiros, Bahia State). No significant differences among developed and undeveloped beaches were observed in terms of total number of items. Local sources (tourism activities) represented 70% of debris on developed beaches, where cigarette butts, straws, paper fragments, soft plastic fragments, and food packaging were the most abundant items. Non-local sources (domestic and fishing activities) accounted for about 70% of debris on undeveloped beaches, where the most abundant items were rigid plastic fragments, ropes, soft plastic fragments, caps, and polystyrene. The projected surface area of beach debris did not vary among developed and undeveloped beaches. Overseas containers accounted for about 25% of regional plastic pollution, implying that international pollution prevention agreements are not being respected off the Brazilian coast.     

Significance of the clay mineral distribution in fluvial sediments of the Neogene to Recent Himalayan Foreland Basin (west‐central Nepal)

Clay mineral assemblages of the Neogene Himalayan foreland basin are studied to decipher their significance with respect to tectonic and climate processes. Fluvial deposits of the Siwalik Group (west‐central Nepal), and sediment of the Ganga River drainage system were analysed for clay mineralogy. The observed clay mineral assemblages are mainly composed of illite (dominant), chlorite, smectite and kaolinite. Illite and chlorite are chiefly of detrital origin, derived from Himalayan sources. Kaolinite and smectite are authigenic, and mainly developed within pore space and as coating of detrital particles. With increasing burial, diagenetic processes affected the original clay mineral signature. Illitisation of smectite and kaolinite occurred below 2500 and 3500 m depth, respectively. Therefore, illite in the lower parts of the Siwalik Group consists of a mixture of inherited illite and illitised smectite and kaolinite, as suggested by illite crystallinity. Detrital grains that make up the framework of the Siwalik Group sandstones mainly consist of quartz, feldspar and lithic fragments, which are principally of sedimentary and metamorphic origin. Lithoclast content increases over time at the expense of quartz and K‐feldspar in response to uplift and erosion of the Lesser Himalaya Series since about 11–10 Ma. Despite mainly felsic source rocks, dominantly physical erosion processes in the Himalayan belt, and high‐energy fluvial depositional systems, smectite is abundant in the <7 Ma Siwalik Group deposits. Analyses of the Siwalik deposits and comparison with the clay mineralogy of the modern drainage system suggest that smectite preferentially formed in floodplains and intermontane valleys during early diagenesis because of downward percolating fluids rich in cations from weathering and soil development. In general, increasing seasonality and aridity linked to variability of the Asian monsoon from about 8 Ma enhanced clay mineral formation and development of authigenic smectite in paleo‐plains on the southern side of the Himalaya.     

Late Quaternary history of the coastal Wahiba Sands,Sultanate of Oman

Continental sediments and geomorphological features of the coastal Wahiba Sands, Sultanate of Oman, reflect environmental variability in southeastern Arabia during the late Quaternary. Weakly cemented dune sands, interdune deposits and coastal sediments were dated by luminescence methods to establish an absolute chronology of changes in sedimentary dynamics. The dating results confirm previous assumptions that during times of low global sea level sand was transported by southerly winds from the exposed shelf onto the Arabian Peninsula. Two prominent phases of sand accumulation in the coastal area took place just before and after the last glacial maximum (LGM). A final significant period of dune consolidation is recognised during the early Holocene. However, no major consolidation of dunes appears to have occurred during the LGM and the Younger Dryas. In the northern part of the Wahiba Sands, these two periods are characterised by substantial sand deposition. This discrepancy is explained by the lack of conservation potential for dunes in the coastal area, probably caused by a low groundwater table due to low sea level and decreased precipitation. While the times of aeolian activity reflect arid to hyper‐arid conditions, lacustrine and pedogenically altered interdune deposits indicate wetter conditions than today caused by increased monsoonal circulation during the Holocene climatic optimum. Copyright © 2005 John Wiley & Sons, Ltd.     

II. Snow hydrology processes and remote sensing

In the last 20 years remote sensing research has led to significant progress in monitoring and measuring certain snow hydrology processes. Snow distribution in a drainage basin can be adequately assessed by visible sensors. Although there are still some interpretation problems, the NOAA-AVHRR sensor can provide frequent views of the areal snow cover in a basin, and snow cover maps are produced operationally by the National Weather Service on about 3000 drainage basins in North America. Measurement of snow accumulation or snow water equivalent with microwave remote sensing has great potential because of the capabilities for depth penetration, all-weather observation and night-time viewing. Several critical areas of research remain, namely, the acquisition of snow grain size information for input to microwave models and improvement in passive microwave resolution from space. Methods that combine both airborne gamma ray and visible satellite remote sensing of the snowpack with field measurements also hold promise for determining areal snow water equivalent. Some remote sensing techniques can also be used to detect different stages of snow metamorphism. Various aspects of snowpack ripening can be detected using microwave and thermal infra-red capabilities. The capabilities for measurement of snow albedo and surface temperature have direct application in both snow metamorphism and snowpack energy balance studies. The potentially most profitable research area here is the study of the bidirectional reflectance distribution function to improve snow albedo measurements. Most of the remote sensing capabilities in snow hydrology have been developed for improving snowmelt-run-off forecasting. Most applications have used the input of snow cover extent to deterministic models, both of the degree day and energy balance types. Snowmelt-run-off forecasts using satellite derived snow cover depletion curves and the models have been successfully made. As the extraction of additional snow cover characteristics becomes possible, remote sensing will have an even greater impact on snow hydrology. Important remote sensing capabilities will become available in the next 20 years through space platform observing systems that will improve our capability to observe the snowpack on an operational basis.     

Paleoclimate implications of high latitude precession-scale mineralogic fluctuations during early Oligocene Antarctic glaciation: the Great Australian Bight record

Sediments from ODP Site 1128 in the Great Australian Bight record isotopic and mineralogic variations corresponding to orbital parameters and regional climate change during the early Oligocene climate transition and Oi1 glacial event. Bulk carbonate stable isotope analyses reveal prominent positive oxygen and carbon isotope shifts related to the inferred major increase in glaciation at approximately 33.6 to 33.48 Ma. The oxygen isotope excursion corresponds to a prolonged period of low eccentricity, suggesting ice-sheet growth during low seasonality conditions. The clay mineralogy is dominated by smectite throughout. The exclusive occurrence of highly crystalline smectite from 33.6 to 33.5 Ma suggests the occurrence of explosive volcanism that correlates with the positive oxygen isotope shift. The dominance of mixed-layer smectite from 33.5 to 33.4 Ma and an increase in illite following 33.4 Ma indicates a transition from cool, wet conditions to cool, dry conditions over Australia during the Oi1 glaciation. Clay mineralogy and carbonate percentages reveal precession-scale oscillations during the Oi1 event. Kaolinite varies inversely with smectite and percent carbonate. Variations in precipitation and runoff, and wind velocities during southern hemisphere summer perihelion and high eccentricity intervals may account for the precession-scale oscillations.