Calculating the effective permeability of sandstone with multiscale lattice Boltzmann/finite element simulations

The lattice Boltzmann (LB) method is an efficient technique for simulating fluid flow through individual pores of complex porous media. The ease with which the LB method handles complex boundary conditions, combined with the algorithm’s inherent parallelism, makes it an elegant approach to solving flow problems at the sub-continuum scale. However, the realities of current computational resources can limit the size and resolution of these simulations. A major research focus is developing methodologies for upscaling microscale techniques for use in macroscale problems of engineering interest. In this paper, we propose a hybrid, multiscale framework for simulating diffusion through porous media. We use the finite element (FE) method to solve the continuum boundary-value problem at the macroscale. Each finite element is treated as a sub-cell and assigned permeabilities calculated from subcontinuum simulations using the LB method. This framework allows us to efficiently find a macroscale solution while still maintaining information about microscale heterogeneities. As input to these simulations, we use synchrotron-computed 3D microtomographic images of a sandstone, with sample resolution of 3.34 μm. We discuss the predictive ability of these simulations, as well as implementation issues. We also quantify the lower limit of the continuum (Darcy) scale, as well as identify the optimal representative elementary volume for the hybrid LB–FE simulations.     

Mainstreaming local perceptions of hurricane risk into policymaking: A case study of community GIS in Mexico

This article suggests a framework for incorporating and communicating local perceptions of hurricane risk into policymaking through a case study conducted at El Zapotito commune in the State of Veracruz, Mexico. The authors constructed a geographical information system (GIS)-based model to quantify and spatially assess specific household-level vulnerabilities from information generated through interviews. This research developed a household vulnerability index applied to a participatory GIS to map vulnerability to hurricane hazard. The results indicate that infrastructural weaknesses are the most important factor contributing to vulnerability, explaining on their own 72.2% of the variation in the vulnerability patterns. These findings are corroborated by a vulnerability and capacity assessment (VCA), which shows that the community lacks strategies to cope with unsafe housing. It is suggested that linking community participation with modern techniques to analyse risk can empower communities and mobilise their capacities to address very specific vulnerabilities.     

Investigating the response of Cladocera to a major saltwater intrusion event in an Arctic lake from the outer Mackenzie Delta (NT, Canada)

An increase in the frequency and intensity of marine storm surges is a predicted consequence of climate warming, and therefore it is important to better understand the biological responses to such events in coastal regions. In late September 1999, a major storm surge resulted in a saltwater intrusion event over a large area of the Mackenzie Delta (NT, Canada) front, causing rapid salinization of lakes on the alluvial plain. Due to a lack of long-term ecological monitoring data in the region, the impacts that the saltwater intrusion event had on the biota of affected lakes were unknown. We used high-resolution paleolimnological approaches to reconstruct past assemblage changes in Cladocera from impacted Lake DZO-29 (unofficial name) in order to determine how different cladoceran species responded to a major increase in lake salinity following the 1999 storm surge. Camptocercus were extirpated from the lake following the saltwater intrusion and have not recovered. We also observed an initial decrease in Alona relative abundance following the marine flooding, likely reflecting a loss of A. quadrangularis, A. barbulata, and A. costata from the lake. A. circumfimbriata, Chydorus biovatus, C. brevilabris, and Bosmina spp. were abundant both before and after the saltwater intrusion, and Paralona pigra was present following the storm surge, but not prior to it. The most notable shift in Cladocera in the recent sedimentary record, however, occurred much earlier, with an increase in pelagic Bosmina taxa and a subsequent decrease in the benthic/littoral taxa Chydorus and Camptocercus, an assemblage shift that is consistent with a response to climate warming in this region, and strongly correlated to other changes in the lake inferred to be as a result of regional warming.     

Tracing metamorphism,exhumation and topographic evolution in orogenic belts by multiple thermochronology: a case study from the Nízke Tatry Mts., Western Carpathians

A combination of four thermochronometers [zircon fission track (ZFT), zircon (U–Th)/He (ZHe), apatite fission track (AFT) and apatite (U–Th–[Sm])/He (AHe) dating methods] applied to a valley to ridge transect is used to resolve the issues of metamorphic, exhumation and topographic evolution of the Nízke Tatry Mts. in the Western Carpathians. The ZFT ages of 132.1 ± 8.3, 155.1 ± 12.9, 146.8 ± 8.6 and 144.9 ± 11.0 Ma show that Variscan crystalline basement of the Nízke Tatry Mts. was heated to temperatures >210°C during the Mesozoic and experienced a low-grade Alpine metamorphic overprint. ZHe and AFT ages, clustering at ~55–40 and ~45–40 Ma, respectively, revealed a rapid Eocene cooling event, documenting erosional and/or tectonic exhumation related to the collapse of the Carpathian orogenic wedge. This is the first evidence that exhumation of crystalline cores in the Western Carpathians took place in the Eocene and not in the Cretaceous as traditionally believed. Bimodal AFT length distributions, Early Miocene AHe ages and thermal modelling results suggest that the samples were heated to temperatures of ~55–90°C during Oligocene–Miocene times. This thermal event may be related either to the Oligocene/Miocene sedimentary burial, or Miocene magmatic activity and increased heat flow. This finding supports the concept of thermal instability of the Carpathian crystalline bodies during the post-Eocene period.     

Identification of groundwater recharge sources and processes in a heterogeneous alluvial aquifer: results from multi‐level monitoring of hydrochemistry and environmental isotopes in a riverside agricultural area in Korea

We evaluated sources and pathways of groundwater recharge for a heterogeneous alluvial aquifer beneath an agricultural field, based on multi‐level monitoring of hydrochemistry and environmental isotopes of a riverside groundwater system at Buyeo, Korea. Two distinct groundwater zones were identified with depth: (1) a shallow oxic groundwater zone, characterized by elevated concentrations of NO₃^? and (2) a deeper (>10–14 m from the ground surface) sub‐oxic groundwater zone with high concentrations of dissolved Fe, silica, and HCO₃^?, but little nitrate. The change of redox zones occurred at a depth where the aquifer sediments change from an upper sandy stratum to a silty stratum with mud caps. The δ¹⁸O and δ²H values of groundwater were also different between the two zones. Hydrochemical and δ¹⁸O? δ²H data of oxic groundwater are similar to those of soil water. This illustrates that recharge of oxic groundwater mainly occurs through direct infiltration of rain and irrigation water in the sandy soil area where vegetable cropping with abundant fertilizer use is predominant. Oxic groundwater is therefore severely contaminated by agrochemical pollutants such as nitrate. In contrast, deeper sub‐oxic groundwater contains only small amounts of dissolved oxygen (DO) and NO₃^?. The ³H contents and elevated silica concentrations in sub‐oxic groundwater indicate a somewhat longer mean residence time of groundwater within this part of the aquifer. Sub‐oxic groundwater was also characterized by higher δ¹⁸O and δ²H values and lower d‐excess values, indicating significant evaporation during recharge. We suggest that recharge of sub‐oxic groundwater occurs in the areas of paddy rice fields where standing irrigation and rain water are affected by strong evaporation, and that reducing conditions develop during subsequent sub‐surface infiltration. This study illustrates the existence of two groundwater bodies with different recharge processes within an alluvial aquifer. Copyright © 2009 John Wiley & Sons, Ltd.     

Bulk density effects on soil hydrologic and thermal characteristics: A numerical investigation

Soil bulk density (ρ_b) is commonly treated as static in studies of land surface dynamics. Magnitudes of errors associated with this assumption are largely unknown. Our objectives were to (a) quantify ρ_b effects on soil hydrologic and thermal properties and (b) evaluate effects of ρ_b on surface energy balance and heat and water transfer. We evaluated 6 soil properties, volumetric heat capacity, thermal conductivity, soil thermal diffusivity, water retention characteristics, hydraulic conductivity, and vapour diffusivity, over a range of ρ_b, using a combination of 6 models. Thermal conductivity, water retention, hydraulic conductivity, and vapour diffusivity were most sensitive to ρ_b, each changing by fractions greater than the associated fractional changes in ρ_b. A 10% change in ρ_b led to 10–11% change in thermal conductivity, 6–11% change in saturated and residual water content, 49–54% change in saturated hydraulic conductivity, and 80% change in vapour diffusivity. Subsequently, 3 field seasons were simulated with a numerical model (HYDRUS‐1D) for a range of ρ_b values. When ρ_b increased 25% (from 1.2 to 1.5 Mg m^?3), soil temperature variation decreased by 2.1 °C in shallow layers and increased by 1 °C in subsurface layers. Surface water content differed by 0.02 m³ m^?3 for various ρ_b values during drying events but differences mostly disappeared in the subsurface. Matric potential varied by >100 m of water. Surface energy balance showed clear trends with ρ_b. Latent heat flux decreased 6%, sensible heat flux increased 9%, and magnitude of ground heat flux varied by 18% (with a 25% ρ_b increase). Transient ρ_b impacted surface conditions and fluxes, and clearly, it warrants consideration in field and modelling investigations.     

Simulated impacts of artificial groundwater recharge and discharge on the source area and source volume of an Atlantic Coastal Plain stream, Delaware, USA

A numerical groundwater-flow model was used to characterize the source area and volume of Phillips Branch, a baseflow-dominated stream incising a highly permeable unconfined aquifer on the low relief Delmarva Peninsula, USA. Particle-tracking analyses indicate that the source area (5.51 km²) is ～20% smaller than the topographically defined watershed (6.85 km²), and recharge entering ～37% of the surface watershed does not discharge to Phillips Branch. Groundwater residence time within the source volume ranges from a few days to almost 100 years, with 95% of the volume “flushing” within 50 years. Artificial discharge from groundwater pumping alters the shape of the source area and reduces baseflow due to the interception of stream flow paths, but has limited impacts on the residence time of groundwater discharged as baseflow. In contrast, artificial recharge from land-based wastewater disposal substantially reduces the source area, lowers the range in residence time due to the elimination of older flow paths to the stream, and leads to increased discharge to adjacent surface-water bodies. This research suggests that, in this and similar hydrogeologic settings, the “watershed” approach to water-resource management may be limited, particularly where anthropogenic stresses alter the transport of soluble contaminants through highly permeable unconfined aquifers.     

Impact of land-atmospheric coupling in CFSv2 on drought prediction

Recent summers in the United States have been plagued by intense droughts that have caused significant damage to crops and have had a large impact on society. The ability to forecasts such events would allow for preparations that could help reduce the impact on society. Coupled land–atmosphere–ocean models were created to provide such forecasts but there are large uncertainties associated with their predictions. The predictive skill of these models is particularly low during the convective season due to the weaker connections with the oceans and an increase in the land–atmosphere interactions. To better understand the degradation of forecasts skill during the summer months and its connection to the land–atmosphere interactions we analyze National Centers for Environmental Prediction’s Climate Forecast System Version 2 (CFSv2) in terms of its climatological land–atmosphere interactions. To do this we use a recently developed classification of land–atmosphere interactions and other diagnostic variables to compare the reanalysis from the Climate Forecast System (CFSR) with CFSv2 re-forecasts (CFSRR) over the period 1982–2009. Coupling in the CFSRR tends toward the wet coupling regime for most areas east of the Rocky Mountains. Although the specific mechanism driving CFSRR to wet coupling state varies by region, the overall cause is enhanced vegetation rooting depth, originally implemented to address a near-surface warm bias in CFSR. The long-term tendency to wet coupling precludes the forecast model from consistently predicting and maintaining drought over the continental US.