FracKfinder: A MATLAB Toolbox for Computing Three-Dimensional Hydraulic Conductivity Tensors for Fractured Porous Media

Fractures in porous media have been documented extensively. However, they are often omitted from groundwater flow and mass transport models due to a lack of data on fracture hydraulic properties and the computational burden of simulating fractures explicitly in large model domains. We present a MATLAB toolbox, FracKfinder, that automates HydroGeoSphere (HGS), a variably saturated, control volume finite-element model, to simulate an ensemble of discrete fracture network (DFN) flow experiments on a single cubic model mesh containing a stochastically generated fracture network. Because DFN simulations in HGS can simulate flow in both a porous media and a fracture domain, this toolbox computes tensors for both the matrix and fractures of a porous medium. Each model in the ensemble represents a different orientation of the hydraulic gradient, thus minimizing the likelihood that a single hydraulic gradient orientation will dominate the tensor computation. Linear regression on matrices containing the computed three-dimensional hydraulic conductivity (K) values from each rotation of the hydraulic gradient is used to compute the K tensors. This approach shows that the hydraulic behavior of fracture networks can be simulated where fracture hydraulic data are limited. Simulation of a bromide tracer experiment using K tensors computed with FracKfinder in HGS demonstrates good agreement with a previous large-column, laboratory study. The toolbox provides a potential pathway to upscale groundwater flow and mass transport processes in fractured media to larger scales.     

Effects of Wormhole Configurations on Performance of Cyclic Solvent Injection in Heavy Oil Systems

Natural Resources Research - In this study, the effect of three different wormhole configurations: single-linear, single-branched, and double-branched wormholes on the performance of cyclic solvent...     

Groundwater Depletion and Agricultural Land Use Change in the High Plains: A Case Study from Wichita County,Kansas*

Shifts in resource availability and resulting land use changes are important research themes for coupled human and natural systems and human dimensions of global environmental change. This study examines the relationship between groundwater depletion and agricultural land use change within a single county in the High Plains, where agro-economic systems and populations are dependent on access to reliable groundwater supplies. Results indicate a significant relationship between high proportionate decreases in groundwater and land being removed from irrigation. This case study strengthens the links among local, regional, and global environmental change, and contributes to a greater understanding of the High Plains in the context of regional sustainability and vulnerability concerns.     

全球化与城市国际化:国际城市的一项实证研究

在当今全球化时代，许多城市被确认为世界城市或全球城市，但本文认为世界上只有少数一些城市真正达到“世界城市”的标准，具有控制世界经济的功能。而许多其他城市，无论是支配或被支配的城市，只显示出国际化的若干特征。本文利用全球264个城市与69家高级生产者服务公司的资料，从国际城市的角度研究世界上的重要城市。研究发现不同的城市有不同的国际化水平，除了少数高度国际化的一流城市达到世界城市的标准，许多其他城市只达到局部国际化，且他们的分布在地域上并无集中的趋势。同时研究发现世界上只有少数高级服务公司拥有全球城市服务网络，大部分公司的服务网络只集中在部分地区。     

Vertical-axis rotation controlled by upper crustal stress based on force balance analysis: A case study of the western Transverse Ranges of California

This paper evaluates driving mechanisms of vertical-axis rotation using data from the western Transverse Ranges in southern California. Simple force balance considerations and comparison of torque applied to a rotating block indicate that shear forces applied to the base of the block are not strong enough to produce the motions and deformation observed at the surface. For the measured dimensions of the crustal blocks and crustal viscosities in southern California, stresses transmitted through the upper crust are one to three orders of magnitude stronger than forces generated in the ductile lower crust. These results suggest that the kinematics of crustal blocks in continental deformation zones are primarily controlled by forces within the upper crust rather than a flow field beneath.     

Extension of a fluctuating plume model of tracer dispersion to a sheared boundary layer and to a large array of obstacles

Fluctuating plume models provide a useful conceptual paradigm in the understanding of plume dispersion in a turbulent flow. In particular, these models have enabled analytical predictions of higher-order concentration moments, and the form of the one-point concentration probability density function (PDF). In this paper, we extend the traditional formalism of these models, grounded in the theory of homogeneous and isotropic turbulent flow, to two cases: namely, a simple sheared boundary layer and a large array of regular obstacles. Some very high-resolution measurements of plume dispersion in a water channel, obtained using laser-induced fluorescence (LIF) line-scan techniques are utilised. These data enable us to extract time series of plume centroid position (plume meander) and dispersion in the relative frame of reference in unprecedented detail. Consequently, experimentally extracted PDFs are able to be directly compared with various theoretical forms proposed in the literature. This includes the PDF of plume centroid motion, the PDF of concentration in the relative frame, and a variety of concentration moments in the absolute and relative frames of reference. The analysis confirms the accuracy of some previously proposed functional forms of model components used in fluctuating plume models, as well as suggesting some new forms necessary to deal with the complex boundary conditions in the spatial domain.     

Heat tracing of heterogeneous hyporheic exchange adjacent to in‐stream geomorphic features

There are many field techniques used to quantify rates of hyporheic exchange, which can vary in magnitude and direction spatially over distances of only a few metres, both within and between morphological features. We used in‐stream mini‐piezometers and heat transport modelling of stream and streambed temperatures to quantify the rates and directions of water flux across the streambed interface upstream and downstream of three types of in‐stream geomorphic features: a permanent dam, a beaver dam remnant and a stream meander. We derived hyporheic flux estimates at three different depths at six different sites for a month and then paired those flux rates with measurements of gradient to derive hydraulic conductivity (K) of the streambed sediments. Heat transport modelling provided consistent daily flux estimates that were in agreement directionally with hydraulic gradient measurements and also identified vertical heterogeneities in hydraulic conductivity that led to variable hyporheic exchange. Streambed K varied over an order of magnitude (1·9 × 10^?6 to 5·7 × 10^?5 m/s). Average rates of hyporheic flux ranged from static (q < ±0·02 m/day) to 0·42 m/day. Heat transport modelling results suggest three kinds of flow around the dams and the meander. Copyright © 2010 John Wiley & Sons, Ltd.