An extended finite volume method and fixed-stress approach for modeling fluid injection–induced tensile opening in fractured reservoirs

Fluid injection–induced tensile opening is modeled using an extended finite volume method (XFVM). An embedded fracture strategy is used for the flow problem, that is, the fractures are discretized using finite volume segments without resolving the grid around them. Further, the discontinuities across fractures are modeled using special basis functions. The fracture openings due to enhanced fluid pressure and the associated shear slip due to traction free boundary condition on the fracture segments are both modeled using these special discontinuity basis functions. Mass transfer between fractures and matrix is modeled using the pressure difference. The enhancement of fracture storativity due to tensile opening leads to stronger coupling between flow and mechanics. An iterative scheme relying on the fixed-stress approach for fractures, which conserves the stress dependent terms over each iteration of the flow problem, has been introduced. Tensile opening has been simulated for single fractures embedded in two- and three-dimensional matrices. The convergence criterion for sequentially implicit fixed-stress scheme for fractures embedded in elastic media is established and has been validated numerically. Further, for 2D simulations, the effect of the matrix permeability for fracture propagation due to tensile opening has been studied.     

A geomorphology‐based integrated stream–aquifer interaction model for semi‐gauged catchments

Many of the existing stream–aquifer interaction models available in the literature are very complex with limited applicability in semi‐gauged and ungauged catchments. In this study, to estimate the influent and effluent subsurface water fluxes under limited geo‐hydrometeorological data availability conditions, a simple stream–aquifer interaction model, namely, the variable parameter McCarthy–Muskingum (VPMM) hillslope‐storage Boussinesq (hsB) model, has been developed. This novel model couples the VPMM streamflow transport with the hsB groundwater flow transport modules in online mode. In this integrated model, the surface water–groundwater flux exchange process is modelled by the Darcian approach with the variable hydraulic heads between the river stage and groundwater table accounting for the rainfall forcing. Considering the exchange fluxes in the hyporheic zone and lateral overland flow contribution, this approach is field tested in a typical 48‐km stretch of the Brahmani River in eastern India to simulate the streamflow and its depth with the minimum Nash–Sutcliffe efficiency of 94% and 88%; the maximum root mean square error of 134 m³/s and 0.35 m; and the minimum index of agreement of 98% and 97%, respectively. This modelling approach could be very well utilized in data‐scarce world‐river basins to estimate the stream–aquifer exchange flux due to rainfall forcings.     

Near‐surface hydraulic conductivity of northern hemisphere glaciers

The hydrology of near‐surface glacier ice remains a neglected aspect of glacier hydrology despite its role in modulating meltwater delivery to downstream environments. To elucidate the hydrological characteristics of this near‐surface glacial weathering crust, we describe the design and operation of a capacitance‐based piezometer that enables rapid, economical deployment across multiple sites and provides an accurate, high‐resolution record of near‐surface water‐level fluctuations. Piezometers were employed at 10 northern hemisphere glaciers, and through the application of standard bail–recharge techniques, we derive hydraulic conductivity (K) values from 0.003 to 3.519 m day^?1, with a mean of 0.185 ± 0.019 m day^?1. These results are comparable to those obtained in other discrete studies of glacier near‐surface ice, and for firn, and indicate that the weathering crust represents a hydrologically inefficient aquifer. Hydraulic conductivity correlated positively with water table height but negatively with altitude and cumulative short‐wave radiation since the last synoptic period of either negative air temperatures or turbulent energy flux dominance. The large range of K observed suggests complex interactions between meteorological influences and differences arising from variability in ice structure and crystallography. Our data demonstrate a greater complexity of near‐surface ice hydrology than hitherto appreciated and support the notion that the weathering crust can regulate the supraglacial discharge response to melt production. The conductivities reported here, coupled with typical supraglacial channel spacing, suggest that meltwater can be retained within the weathering crust for at least several days. Not only does this have implications for the accuracy of predictive meltwater run‐off models, but we also argue for biogeochemical processes and transfers that are strongly conditioned by water residence time and the efficacy of the cascade of sediments, impurities, microbes, and nutrients to downstream ecosystems. Because continued atmospheric warming will incur rising snowline elevations and glacier thinning, the supraglacial hydrological system may assume greater importance in many mountainous regions, and consequently, detailing weathering crust hydraulics represents a research priority because the flow path it represents remains poorly constrained.     

Long‐term and seasonal hydrologic performance of an extensive green roof

Sustainable strategies such as green roofs have been implemented as stormwater management tools to mitigate disturbance of the hydrologic cycle resulting from urbanization. Green roofs, also referred to as vegetated roofs, can improve the urban landscape by reducing heat island effects, providing ecosystem services, and facilitating the retention and treatment of stormwater. Green roofs have received particular attention because they do not require acquisition and development of land and represent an application of biomimicry in design and construction. In this paper, we evaluate the effects of precipitation, evapotranspiration (ET), antecedent dry period (ADP), and seasonal variation on the run‐off quantity and distribution of an extensive, sedum covered, green roof on a commercial building in Syracuse, NY, USA. The green roof greatly facilitated retention of precipitation events without significant changes over the 4‐year study. The green roof retained on average 95.9 ± 3.6% (6.5 ± 5.6 mm) per rainfall event, with a range from 75% to 99.6% (33.2 to 3.3 mm). However, as precipitation quantity increased, the retention of water decreased. This high water retention capacity was the result of the combined effects of ET, stormwater storage (plants, growth media, and stormwater retention layer), and limited surface run‐off from the roof deck due to variation in the sloping of the green roof and the tapered insulation to the deck drains. The water retention capacity of the green roof did not change significantly between growing and nongrowing seasons. Slightly greater precipitation during the growing season coincided with increased ET. Average potential ET during the growing season was approximately 3 times greater than during the nongrowing season. The hydrologic performance of the green roof was not significantly impacted by an ADP greater than 2 days.     

Spatially distributed tracer‐aided modelling to explore water and isotope transport,storage and mixing in a pristine,humid tropical catchment

Rapidly transforming headwater catchments in the humid tropics provide important resources for drinking water, irrigation, hydropower, and ecosystem connectivity. However, such resources for downstream use remain unstudied. To improve understanding of the behaviour and influence of pristine rainforests on water and tracer fluxes, we adapted the relatively parsimonious, spatially distributed tracer‐aided rainfall–runoff (STARR) model using event‐based stable isotope data for the 3.2‐km² San Lorencito catchment in Costa Rica. STARR was used to simulate rainforest interception of water and stable isotopes, which showed a significant isotopic enrichment in throughfall compared with gross rainfall. Acceptable concurrent simulations of discharge (Kling–Gupta efficiency [KGE] ~0.8) and stable isotopes in stream water (KGE ~0.6) at high spatial (10 m) and temporal (hourly) resolution indicated a rapidly responding system. Around 90% of average annual streamflow (2,099 mm) was composed of quick, near‐surface runoff components, whereas only ~10% originated from groundwater in deeper layers. Simulated actual evapotranspiration (ET) from interception and soil storage were low (~420 mm/year) due to high relative humidity (average 96%) and cloud cover limiting radiation inputs. Modelling suggested a highly variable groundwater storage (~10 to 500 mm) in this steep, fractured volcanic catchment that sustains dry season baseflows. This groundwater is concentrated in riparian areas as an alluvial–colluvial aquifer connected to the stream. This was supported by rainfall–runoff isotope simulations, showing a “flashy” stream response to rainfall with only a moderate damping effect and a constant isotope signature from deeper groundwater (~400‐mm additional mixing volume) during baseflow. The work serves as a first attempt to apply a spatially distributed tracer‐aided model to a tropical rainforest environment exploring the hydrological functioning of a steep, fractured‐volcanic catchment. We also highlight limitations and propose a roadmap for future data collection and spatially distributed tracer‐aided model development in tropical headwater catchments.     

High‐T and low‐P metamorphism in the Xilingol Complex of central Inner Mongolia,China: An indicator of extension in a previous orogeny

The Xilingol Complex comprises biotite gneisses and amphibolite interlayers with extensive migmatization. Four representative samples were documented and found to record either two or three metamorphic stages. Phase modelling using thermocalc suggests that the observed assemblages represent the final stages that underwent cooling from temperature peaks, and are consistent with a fluid‐absent solidus in P–T pseudosections. Their P–T conditions are further constrained to be 5–6 kbar/680–725°C and 4–5 kbar/650–680°C for two garnet‐bearing gneiss samples, 4–5 kbar/660–730°C for a cordierite‐bearing gneiss sample, and 4–5 kbar/680–710°C for an amphibolite sample based on mineral composition isopleths, involving measured Mg content in biotite, anorthite in plagioclase, grossular and pyrope in garnet and Ti content in amphibole. The peak temperature conditions recovered are 760–790°C or >760°C at 5–6 kbar based on the composition isopleths of plagioclase, biotite, garnet and especially the comparison of melt contents between the calculated and observed. A pre‐peak heating process with slight decompression can be suggested for some samples on the basis of the core–rim increase in the plagioclase anorthite, and the stability of ilmenite. Zircon U–Pb dating using the LA‐ICP‐MS method provides systemic constraints on the metamorphic ages of the Xilingol Complex to be 348–305 Ma, interpreted to represent the post‐peak cooling stages. Moreover, metagabbroic dykes that intruded into the Xilingol Complex yield 317 ± 3 Ma from magmatic zircon, and are considered to have played a significant role for heat advection triggering the high‐T and low‐P metamorphism. Thus, the clockwise P–T paths involving pre‐peak heating, peak and post‐peak cooling recovered for the Xilingol Complex are consistent with an extensional setting in the Carboniferous that developed on a previous orogen in response to addition of mantle‐derived materials probably together with upwelling of the asthenospheric mantle.     

Catchment properties,function, and conceptual model representation: is there a correspondence?

This study investigates the possible correspondence between catchment structure, as represented by perceptual hydrological models developed from fieldwork investigations, and mathematical model structures, selected on the basis of reproducing observed catchment hydrographs. Three Luxembourgish headwater catchments are considered, where previous fieldwork suggested distinct flow‐generating mechanisms and hydrological dynamics. A set of lumped conceptual model structures are hypothesized and implemented using the SUPERFLEX framework. Following parameter calibration, the model performance is examined in terms of predictive accuracy, quantification of uncertainty, and the ability to reproduce the flow–duration curve signature. Our key research question is whether differences in the performance of the conceptual model structures can be interpreted based on the dominant catchment processes suggested from fieldwork investigations. For example, we propose that the permeable bedrock and the presence of multiple aquifers in the Huewelerbach catchment may explain the superior performance of model structures with storage elements connected in parallel. Conversely, model structures with serial connections perform better in the Weierbach and Wollefsbach catchments, which are characterized by impermeable bedrock and dominated by lateral flow. The presence of threshold dynamics in the Weierbach and Wollefsbach catchments may favour nonlinear models, while the smoother dynamics of the larger Huewelerbach catchment were suitably reproduced by linear models. It is also shown how hydrologically distinct processes can be effectively described by the same mathematical model components. Major research questions are reviewed, including the correspondence between hydrological processes at different levels of scale and how best to synthesize the experimentalist's and modeller's perspectives. Copyright © 2013 John Wiley & Sons, Ltd.     

A rainfall‐runoff model parameterized from GIS and runoff data

A guiding principle in hydrological modelling should be to keep the number of calibration parameters to a minimum. A reduced number of parameters to be calibrated, while maintaining the accuracy and detail required by modern hydrological models, will reduce parameter and model structure uncertainty and improve model diagnostics. In this study, the dynamics of runoff are derived from the distribution of distances from points in the catchments to the nearest stream. This distribution is unique for each catchment and can be determined from a geographical information system. The distribution of distances, will, when a celerity of (subsurface) flow is introduced, provide a distribution of travel times, or a unit hydrograph (UH). For spatially varying levels of saturation deficit, we have different celerities and, hence, different UHs. Runoff is derived from the superposition of the different UHs. This study shows how celerities can be estimated if we assume that recession events represent the combined UHs for different levels of saturation deficit. A new soil moisture routine which estimates saturated and unsaturated volumes of subsurface water and with only one parameter to calibrate is included in the new model. The performance of the new model is compared with that of the Swedish HBV model and is found to perform equally well for eight Norwegian catchments although the number of parameters to be calibrated in the module concerning soil moisture and runoff dynamics is reduced from seven in the HBV model to one in the new model. It is also shown that the new model has a more realistic representation of the subsurface hydrology. Copyright © 2013 John Wiley & Sons, Ltd.     

GIS research to address tensions in geography

The rise of Geographic Information Systems (GIS) technologies has been producing powerful tools for spatial data processing, management, analysis, modeling, and visualization. While supporting many tasks, GIS technologies have evoked new thinking and advanced intellectual inquiries in geography. Technological advances in other fields often stimulate new research questions and lead to revolutionary discoveries. The Hubble telescope revolutionizes our understanding of the universe, and 3D digital microscopes transform our knowledge of the coordination among biological, neurological, and physiological systems in living organisms. Can GIS claim similar revolutionary effects on geography? The answer is much up for debate. With GIS technologies, geographic studies can explore a broader extent across multiple scales in space and time and tackle problems through increasingly complex spatial statistics, visual analytics, computation, simulation, machine learning, and artificial intelligence. Both the Hubble telescope and 3D digital microscope were built based on scientific research that subserves the technological advances for inquiries into space and life systems. Likewise, GIScience research innovates GIS methods that enable novel geographic investigations and therefore contribute to geographic knowledge production. For the sake of simplicity, the term GIS used here represents both GIScience and GIS technologies. This essay attempts to clarify the intellectual contributions of GIS to geography on the following two questions: (1) What novel geographic thinking is driven by GIS? (2) How may GIS provoke new geographic inquiries and knowledge? Building on Nystuen's notion of four tensions that trigger geographic questions, the essay discusses how GIS innovations mediate historical tension, space‐time tension, dimensional tension and scale tensions.