A numerical model for estimating groundwater flux from subsurface temperature profiles

A simple numerical model is presented for estimating vertical groundwater flux from transient subsurface temperature profiles obtained from field measurements. The model developed utilizes the MacCormack scheme, which is based on the Finite Difference Method (FDM), for solving the governing partial differential equation of convection–diffusion heat transport with appropriate initial and boundary conditions within the subsurface. In order to validate the model, numerical solutions obtained for the study area located in the Nagoka plain, Japan are compared with the published measured data and results obtained by others. Results obtained show good agreement and fit the observed data with a correlation coefficient, R², of 0·88. The estimated groundwater flux is 1·85 × 10⁻⁷ m s⁻¹. Sensitivity analyses were also carried out to investigate the effect of variations in groundwater fluxes, thermal properties and the annual thermal variability due to climatic changes on the transient subsurface temperature profiles and to have a better understanding of the subsurface thermal dynamics. A substantial effect of annual climatic variability is observed on the temporal distributions of temperature depth profiles, and a better estimate of thermal parameters is required to estimate vertical groundwater flux. The largest change in subsurface temperature depth profiles due to groundwater flux over a year is within ± 4 °C. The influence of groundwater flux on subsurface temperature distributions in space and time may be more pronounced in areas where the top of the saturated layer fluctuates considerably. Variation in thermal diffusivity results in temperature change up to ± 1·5% and may cause change in groundwater flux estimate by ± 18%. The model presented has merits over analytical solutions (type curve matching techniques) in terms of suitability and applicability to real field problems, and can be a good asset to hydrological models as quantifying groundwater recharge or deducing it from other quantities, such as rainfall, evapotranspiration and runoff, is often complicated. Copyright © 2007 John Wiley & Sons, Ltd.     

Lithofacies modeling by multipoint statistics and economic evaluation by NPV volume for the early Cretaceous Wabiskaw Member in Athabasca oilsands area,Canada

The static modeling and dynamic simulation are essential and critical processes in petroleum exploration and development. In this study, lithofacies models for Wabiskaw Member in Athabasca, Canada are generated by multipoint statistics(MPS) and then compared with the models built by sequential indicator simulation(SIS). Three training images(Tls) are selected from modern depositional environments;the Orinoco River Delta estuary, Cobequid bay-Salmon River estuary, and Danube River delta environment. In order to validate lithofacies models, average and variance of similarity in lithofacies are calculated through random and zonal blind-well tests.In random six-blind-well test, similarity average of MPS models is higher than that of SIS model. The Salmon MPS model closely resembles facies pattern of Wabiskaw Member in subsurface. Zonal blind-well tests show that successful lithofacies modeling for transitional depositional setting requires additional or proper zonation information on horizontal variation, vertical proportion, and secondary data.As Wabiskaw Member is frontier oilsands lease, it is impossible to evaluate the economics from production data or dynamic simulation. In this study, a dynamic steam assisted gravity drainage(SAGD)performance indicator(SPIDER) on the basis of reservoir characteristics is calculated to build 3 D reservoir model for the evaluation of the SAGD feasibility in Wabiskaw Member. SPIDER depends on reservoir properties, economic limit of steam-oil ratio, and bitumen price. Reservoir properties like porosity,permeability, and water saturation are measured from 13 cores and calculated from 201 well-logs. Three dimensional volumes of reservoir properties are constructed mostly based on relationships among properties. Finally, net present value(NPV) volume can be built by equation relating NPV and SPIDER. The economic area exceeding criterion of US$ 10,000 is identified, and the ranges of reservoir properties are estimated. NPV-volume-generation workflow from reservoir parameter to static model provides costand time-effective method to evaluate the oilsands SAGD project.     

Understanding the impact of droughts in the Yarmouk Basin,Jordan: monitoring droughts through meteorological and hydrological drought indices

This article assesses drought status in the Yarmouk Basin (YB), in northern Jordan, using the Standardized Precipitation Index (SPI), the Standardized Water-Level Index (SWI), and the Percent Departure from Normal rainfall (PDNimd) during the years 1993–2014. The results showed that the YB suffers from frequent and irregular periods of drought as variations in drought intensity and frequency have been observed. The SPI results revealed that the highest drought magnitude of ??2.34 appeared at Nuaimeh rainfall station in 1991. This station has also experienced severe drought particularly in years 1995, 1999, 2005, and 2012 with SPI values ranging from ??1.51 to ??1.59. Some other rainfall stations such as Baqura, Ibbin, Khanasiri, Kharja, Mafraq police, Ramtha, Turra, and Umm Qais have also suffered several periods of drought mostly in 1993. The SWI results show the highest extreme drought events in 2001 in Souf well while other extreme drought periods were observed at Wadi Elyabis well in 1994 and at Mafraq well in 1995. As compared to SPI maps, our SWI maps reflect severe and extreme drought events in most years, negatively impacting the groundwater levels in the study area.     

Deep flow and transport through the Ulleung Interplain Gap in the Southwestern East/Japan Sea

Deep circulation in the southwestern East/Japan Sea through the Ulleung Interplain Gap (UIG), a possible pathway for deep-water exchange, was directly measured for the first time. Five concurrent current meter moorings were positioned to effectively span the UIG between the islands of Ulleungdo to the west and Dokdo to the east. They provided a 495-day time series of deep currents below 1800 m depth spanning the full breadth of the East Sea Deep and Bottom Water flowing from the Japan Basin into the Ulleung Basin. The UIG circulation is found to be mainly a two-way flow with relatively weak southward flows directed into the Ulleung Basin over about two-thirds of the western UIG. A strong, persistent, and narrow compensating northward outflow occurs in the eastern UIG near Dokdo and is first referred to here as the Dokdo Abyssal Current. The width of the abyssal current is about 20 km below 1800 m depth. The low-frequency variability of the transports is dominated by fluctuations with a period of about 40 days for inflow and outflow transports. The 40-day fluctuations of both transports are statistically coherent, and occur almost concurrently. The overall mean transport of the deep water below 1800 m into the Ulleung Basin over the 16.5 months is about 0.005 Sv (1 Sv=10⁶ m³ s^?1), with an uncertainty of 0.025 Sv indicating net transport is negligible below 1800 m through the UIG.     

Infiltration pattern in a regolith–fractured bedrock profile: field observation of a dye stain pattern

We examined the infiltration pattern of water in a regolith–bedrock profile consisting of two overburdens (OB1 and OB2), a buried rice paddy soil (PS), two texturally distinctive weathered materials (WM1 and WM2) and a fractured sedimentary rock (BR), using a Brilliant Blue FCF dye tracer. A black‐coloured coating in conducting fractures in WM1, WM2 and BR was analysed by X‐ray diffraction and scanning electron microscopy. The dye tracer penetrated to greater than 2 m depth in the profile. The macropore flow and saturated interflow were the major infiltration patterns in the profile. Macropore flow and saturated interflow were observed along fractures in WM1, WM2 and BR and at the dipping interfaces of PS–WM1, PS–WM2 and PS–BR respectively. Heterogeneous matrix flow occurred in upper overburden (OB1) and PS. Compared with OB1, the coarser textured OB2 acted as a physical barrier for vertical flow of water. The PS with low bulk density and many fine roots was another major conducting route of water in the profile. Manganese oxide and iron oxide were positively identified in the black coating material and had low crystallinity and high surface area, indicating their high reactivity with conducting contaminants. Copyright © 2006 John Wiley & Sons, Ltd.