Horizontal Mass Flux Research Articles

Difference in wind erosion characteristics between loamy and sandy farmlands and the implications for soil dust emission potential

AbstractDust emission depends directly and indirectly on soil particle distribution, but little information is available that examines dust emission potential and its relationship with wind erosion for different soils. This study aimed to elucidate the difference in wind erosion characteristics between loamy farmland (LF) and sandy farmland (SF) and discuss the influence of this difference on dust emission potential. Field observations were conducted during 2012–2016 in Bashang district, North China. With increasing friction wind speed, the horizontal mass flux increased exponentially in LF but obeyed the power function in SF. The horizontal mass flux was higher in SF than in LF, but the gap narrowed from 766 times at friction wind speed of 0.3 m s−1 to 50 times at 0.6 m s−1. The dry aggregate size distribution of sediment in LF was finer than in SF. Fine particle content (diameter < 0.05 mm) in windblown sediment was much higher in LF (31.39%) than in SF (5.92%), which was influenced by parent soils. With height increasing, the mass percentage of sediment decreased more rapidly in SF than in LF, with 86% and 66% of sediment clustering at height 0–20 cm, respectively. Because the horizontal mass flux and dust content in sediment both increased more quickly in LF than in SF with increasing friction wind speed, the dust emission potential was larger in LF than in SF under strong winds. The adverse effect of dust emission in loamy soils on atmospheric pollution should receive more attention in strong wind erosion events.

Local and non-local effects of building arrangements on pollutant fluxes within the urban canopy

This work investigates the vertical and horizontal mass (scalar) flux of a contaminant emitted from an area source located in an array of blocks representing an urban environment. Arrays consisting of buildings with random and uniform heights and staggered and aligned arrangements were tested. Results shows that the vertical scalar flux close to the source can affect downwind clean zones. It is also shown that taller buildings increase the vertical scalar flux and the fluctuations of the vertical velocity above the smaller buildings. The vertical advective scalar flux was found to have an effect on dispersion in the vicinity of the building (a local effect), while the vertical turbulent fluxes are associated with pollutant transportation downwind above the smaller buildings (a non-local effect).

Unseen Dust Emission and Global Dust Abundance: Documenting Dust Emission from the Mojave Desert (USA) by Daily Remote Camera Imagery and Wind‐Erosion Measurements

AbstractA large gap in understanding the effects of atmospheric dust at all spatial scales is uncertainty about how much and whence dust is emitted annually. Direct digital recording of dust emission at high spatial and temporal resolution would together with periodic flux measurements support improved estimates of local‐scale dust flux. Remote camera recording of dust‐emitting settings on and around Soda Lake (Mojave Desert) was conducted every 15 min during daylight between 10 November 2010 and 31 December 2016. Examination of 135,000 images revealed frequent dust events characterized as dust days, 68 per year on average, when dust plumes obscured mountains beyond areas of dust generation. Dust was not observed from examination of satellite retrievals (Moderate Resolution Imaging Spectroradiometer and Geostationary Operational Environmental Satellite System) during six cloudless days of extreme emission. Estimated masses of emitted dust fractions were made from measured horizontal mass flux, sediment particle sizes, and areas of dust generation. Between April 2000 and June 2013, nearly 4 Tg/year of dust (particles < 63 μm) were emitted across the study area. Higher rates (about 7 Tg/year) were estimated following the December 2010 Mojave River flood that deposited sediment across the lake basin. Within the Mojave and Great Basin deserts of western North America, many settings akin to those at Soda Lake similarly emit dust that is rarely detected in satellite retrievals, implying that local and regional dust emissions from western North America are far underestimated and that, by extension to relatively small dust‐source areas across all drylands, global dust emissions may also be greatly underestimated.

Elevated aeolian sediment transport on the Colorado Plateau, USA: The role of grazing, vehicle disturbance, and increasing aridity

AbstractDryland wind transport of sediment can accelerate soil erosion, degrade air quality, mobilize dunes, decrease water supply, and damage infrastructure. We measured aeolian sediment horizontal mass flux (q) at 100 cm height using passive aspirated sediment traps to better understand q variability on the Colorado Plateau. Measured q ‘hot spots’ rival the highest ever recorded including 7,460 g m−2 day−1 in an off‐highway vehicle (OHV) area, but were more commonly 50‐2,000 g m−2 day−1. Overall mean q on rangeland sites was 5.14 g m−2 day−1, considerably lower than areas with concentrated livestock use (9‐19 g m−2 day−1), OHV use (414 g m−2 day−1), and downwind of unpaved roads (13.14 g m−2 day−1), but were higher than areas with minimal soil disturbance (1.60 g m−2 day−1). Rangeland q increased with increasing annual temperature, increased winds, and decreasing precipitation. Spatial modeling suggests that ~92‐93% of regional q occurs in rangelands versus ~7‐8% along unpaved roads. Four of the five largest road q values (n=33) measured were along roads used primarily for oil or gas wells. Our findings indicate that predicted future mega‐droughts will increase q disproportionately in disturbed rangelands, and potentially further compromise air quality, hydrologic cycles, and other ecosystem services. Published 2018. This article is a U.S. Government work and is in the public domain in the USA.

Modeling natural dust emissions in the central Middle East: Parameterizations and sensitivity

Middle East Area and especially the central part is a dust dominated domain and therefore estimation of natural dust emissions is critical. Aim of this study is to assess the sensitivity of state-of-the art dust modules, utilized in multiple modelling schemes, and especially of their three main components; namely the horizontal mass flux, the drag partition and the sandblasting efficiency. To accomplish this, several simulations with the Natural Emission MOdel (NEMO), driven by the Weather Research and Forecasting (WRF) model, were made for the period April–June 2015 focusing in the central Middle East (CME). First, the meteorological model was evaluated over the study period, as well as during a Shamal event, showing a satisfactory performance on wind speed which is the most important meteorological variable inserted in the dust modules. Second, dust emissions were found to be a few tenths of μg/m2s on average, while the daily amount exceeds 1000 mg/m2day. Over the studied period, the fluxes range from 15 to 750 g/m2. The total emissions in CME were estimated with the NEMO's basic configuration at 139 Tg. The most active areas are associated with clay content ranging between 17 and 25%, highlighting the importance of soil texture on the definition of the dust sources. During a Shamal event in the second half of June 2015, the mean dust emission fluxes reached up to several hundreds of μg/m2s, contributing 21 Tg in CME. Few such events could cover a large portion of the dust emissions. Finally, concerning the sensitivity of the total emissions, the drag partition scheme has the largest effect, followed by sandblasting efficiency and horizontal mass flux. The total dust emissions budget, among the 16 simulations was found to be in the range of 21–139 Tg, with the choice of the drag partition scheme to be mandatory for the narrowing of the total emissions in CME.

Exploring dust emission responses to land cover change using an ecological land classification

Abstract Despite efforts to quantify the impacts of land cover change on wind erosion, assessment uncertainty remains large. We address this uncertainty by evaluating the application of ecological site concepts and state-and-transition models (STMs) for detecting and quantitatively describing the impacts of land cover change on wind erosion. We apply a dust emission model over a rangeland study area in the northern Chihuahuan Desert, New Mexico, USA, and evaluate spatiotemporal patterns of modelled horizontal sediment mass flux and dust emission in the context of ecological sites and their vegetation states; representing a diversity of land cover types. Our results demonstrate how the impacts of land cover change on dust emission can be quantified, compared across land cover classes, and interpreted in the context of an ecological model that encapsulates land management intensity and change. Results also reveal the importance of established weaknesses in the dust model soil characterisation and drag partition scheme, which appeared generally insensitive to the impacts of land cover change. New models that address these weaknesses, coupled with ecological site concepts and field measurements across land cover types, could significantly reduce assessment uncertainties and provide opportunities for identifying land management options.

Using rare earth elements to trace wind-driven dispersion of sediments from a point source

Abstract The entrainment and movement of aeolian sediments is determined by the direction and intensity of erosive winds. Although erosive winds may blow from all directions, in most regions there is a predominant direction. Dust emission causes the removal preferentially of soil nutrients and contaminants which may be transported tens to even thousands of kilometers from the source and deposited into other ecosystems. It would be beneficial to understand spatially and temporally how the soil source may be degraded and depositional zones enriched. A stable chemical tracer not found in the soil but applied to the surface of all particles in the surface soil would facilitate this endeavor. This study examined whether solution-applied rare earth elements (REEs) could be used to trace aeolian sediment movement from a point source through space and time at the field scale. We applied erbium nitrate solution to a 5 m2 area in the center of a 100 m diameter field 7854 m2 on the Southern High Plains of Texas. The solution application resulted in a soil-borne concentration three orders of magnitude greater than natively found in the field soil. We installed BSNE sampler masts in circular configurations and collected the trapped sediment weekly. We found that REE-tagged sediment was blown into every sampler mast during the course of the study but that there was a predominant direction of transport during the spring. This preliminary investigation suggests that the REEs provide a viable and incisive technique to study spatial and temporal variation of aeolian sediment movement from specific sources to identifiable locations of deposition or locations through which the sediments were transported as horizontal mass flux and the relative contribution of the specific source to the total mass flux.

Quantifying Postfire Aeolian Sediment Transport Using Rare Earth Element Tracers

AbstractGrasslands, which provide fundamental ecosystem services in many arid and semiarid regions of the world, are undergoing rapid increases in fire activity and are highly susceptible to postfire‐accelerated soil erosion by wind. A quantitative assessment of physical processes that integrates fire‐wind erosion feedbacks is therefore needed relative to vegetation change, soil biogeochemical cycling, air quality, and landscape evolution. We investigated the applicability of a novel tracer technique—the use of multiple rare earth elements (REE)—to quantify soil transport by wind and to identify sources and sinks of wind‐blown sediments in both burned and unburned shrub‐grass transition zone in the Chihuahuan Desert, NM, USA. Results indicate that the horizontal mass flux of wind‐borne sediment increased approximately threefold following the fire. The REE tracer analysis of wind‐borne sediments shows that the source of the horizontal mass flux in the unburned site was derived from bare microsites (88.5%), while in the burned site it was primarily sourced from shrub (42.3%) and bare (39.1%) microsites. Vegetated microsites which were predominantly sinks of aeolian sediments in the unburned areas became sediment sources following the fire. The burned areas showed a spatial homogenization of sediment tracers, highlighting a potential negative feedback on landscape heterogeneity induced by shrub encroachment into grasslands. Though fires are known to increase aeolian sediment transport, accompanying changes in the sources and sinks of wind‐borne sediments may influence biogeochemical cycling and land degradation dynamics. Furthermore, our experiment demonstrated that REEs can be used as reliable tracers for field‐scale aeolian studies.

Mass flux analysis of 137Cs plumes emitted from the Fukushima Daiichi nuclear power plant

The flow vectors of radioactive cesium-137 (137Cs) plume emitted from the Fukushima Daiichi nuclear power plant in March 2011 were quantitatively depicted by a mass flux analysis in this study. 137Cs plumes were calculated by an Eulerian dispersion model with a 3-km horizontal resolution. The vertically column-integrated mass flux was consistent with the flow approximation based on ground surface 137Cs observations, even though there were some discrepancies that were caused by differences in the wind direction between the ground surface and the dominant plume layer. These discrepancies were explained by combining the use of the ground surface horizontal mass flux with the column-integrated mass flux. The mass flux analysis clearly provided an illustration of 137Cs dominant stream locations, directions, and depositions by reducing high-dimensional model outputs into a lower-dimensional plot. Mass flux (i.e. the product of the mass density and wind velocity) has often been used in dynamic meteorology but has not been used as frequently in atmospheric chemistry or pollutant dispersion studies. However, the concept of mass flux is a robust alternative for conventional validation approaches that only utilize a time series of pollutant concentrations. Mass flux analyses can be used further in atmospheric chemistry as a quantitative visualization tool to track the emission, advection, dispersion, and deposition of atmospheric constituents.

Wind tunnel observations of weak and strong snow saltation dynamics

AbstractTheoretical considerations suggest that saltation dynamics is dominated by either aerodynamic entrainment or by a combination of ejection and rebound at a given time and location. Calling these two regimes “weak” and “strong” saltation, respectively, we have investigated high‐resolution snow mass flux measurements by shadowgraphy in a cold wind tunnel to determine whether these two regimes can be experimentally reproduced. In this contribution, we first suggest that aerodynamic entrainment should lead to a negative correlation in the growth rates of horizontal versus vertical mass fluxes, while for the ejection regime, a positive correlation is plausible. Based on this criterion, we find evidence of weak and strong saltation in our data. Weak saltation is characterized by smaller peaks of horizontal mass flux; however, these peaks have a higher growth rate and therefore are typically narrow. The larger peaks are associated with strong saltation, are wider, and their frequency of occurrence increases with increasing overall saltation mass flux.

Modeling and Process Features of Plug Flow Reactor with Internal Recirculation for Biomass Pyrolysis

A new fluidized bed design is considered for the pyrolysis of biomass, in a plug-flow with internal-recirculation (PFIR) arrangement, where the movement of solids from feed zone to product zone and the circulation within a closed loop are achieved, by the use of directional high-speed tuyeres that induce both the suspension and horizontal movement of the charge. Separation between feed and product zones is provided through one or more underflow weirs. Wet and hard-to-fluidize biomass is targeted, where fluidization and pyrolysis of the feed is achieved by the internal hot sand recirculation. Plug flow can allow for both stages (pyrolysis and heating) to occur inside a single vessel, but at separate zones. Computer simulation of the PFIR reactor is demonstrated to verify that plug-flow with internal- recirculation can be achieved under reasonable operating conditions. A commercial software (CPFD-Barracuda) is used, which shows promising results, including: solids circulation due to directional tuyeres; large solids mass transport rates relative to reactor size; and controllability of the circulation rate via tuyere velocity. It is shown that the underflow weir(s) offer little resistance compared to the frictional losses and, therefore, multiple zones within a single reactor can be achieved with little impact on the solids recirculation rate. The horizontal mass flux is shown to occur predominantly near the bottom Tuyeres (inclined jet zone), hence enabling the underflow weir to effectively separate the feed and product zones, with a small bottom passage for the solids.

Quantifying horizontal and vertical tracer mass fluxes in an idealized valley during daytime

Abstract. The transport and mixing of pollution during the daytime evolution of a valley boundary layer is studied in an idealized way. The goal is to quantify horizontal and vertical tracer mass fluxes between four different valley volumes: the convective boundary layer, the slope wind layer, the stable core, and the atmosphere above the valley. For this purpose, large eddy simulations (LES) are conducted with the Weather Research and Forecasting (WRF) model for a quasi-two-dimensional valley. The valley geometry consists of two slopes with constant slope angle and is homogeneous in the along-valley direction. The surface sensible heat flux is horizontally homogeneous and prescribed by a sine function. The initial sounding is characterized by an atmosphere at rest and a constant Brunt–Väisälä frequency. Various experiments are conducted for different combinations of surface heating amplitudes and initial stability conditions. A passive tracer is released with an arbitrary but constant rate at the valley floor and resulting tracer mass fluxes are evaluated between the aforementioned volumes.As a result of the surface heating, a convective boundary layer is established in the lower part of the valley with a stable layer on top – the so-called stable core. The height of the slope wind layer, as well as the wind speed within, decreases with height due to the vertically increasing stability. Hence, the mass flux within the slope wind layer decreases with height as well. Due to mass continuity, this along-slope mass flux convergence leads to a partial redirection of the flow from the slope wind layer towards the valley centre and the formation of a horizontal intrusion above the convective boundary layer. This intrusion is associated with a transport of tracer mass from the slope wind layer towards the valley centre. A strong static stability and/or weak forcing lead to large tracer mass fluxes associated with this phenomenon. The total export of tracer mass out of the valley atmosphere increases with decreasing stability and increasing forcing. The effects of initial stability and forcing can be combined to a single parameter, the breakup parameter B. An analytical function is presented that describes the exponential decrease of the percentage of exported tracer mass with increasing B. This study is limited by the idealization of the terrain shape, stratification, and forcing, but quantifies transport processes for a large range of forcing amplitudes and atmospheric stability.

Using albedo to reform wind erosion modelling, mapping and monitoring

Wind erosion and dust emission models are used to assess the impacts of dust on radiative forcing in the atmosphere, cloud formation, nutrient fertilisation and human health. The models are underpinned by a two-dimensional geometric property (lateral cover; L) used to characterise the three-dimensional aerodynamic roughness (sheltered area or wakes) of the Earth’s surface and calibrate the momentum it extracts from the wind. We reveal a fundamental weakness in L and demonstrate that values are an order of magnitude too small and significant aerodynamic interactions between roughness elements and their sheltered areas have been omitted, particularly under sparse surface roughness. We describe a solution which develops published work to establish a relation between sheltered area and the proportion of shadow over a given area; the inverse of direct beam directional hemispherical reflectance (black sky albedo; BSA). We show direct relations between shadow and wind tunnel measurements and thereby provide direct calibrations of key aerodynamic properties. Estimation of the aerodynamic parameters from albedo enables wind erosion assessments over areas, across platforms from the field to airborne and readily available satellite data. Our new approach demonstrated redundancy in existing wind erosion models and thereby reduced model complexity and improved fidelity. We found that the use of albedo enabled an adequate description of aerodynamic sheltering to characterise fluid dynamics and predict sediment transport without the use of a drag partition scheme (Rt) or threshold friction velocity (u∗t). We applied the calibrations to produce global maps of aerodynamic properties which showed very similar spatial patterns to each other and confirmed the redundancy in the traditional parameters of wind erosion modelling. We evaluated temporal patterns of predicted horizontal mass flux at locations across Australia which revealed variation between land cover types that would not be detected using traditional models. Our new approach provided new opportunities to investigate the dynamics of wind erosion in space and time and elucidate aeolian processes across scales.

Effects of self-pulsation on the spray characteristics of gas–liquid swirl coaxial injector

To understand the influence of self-pulsation on the spray characteristics of gas–liquid swirl coaxial injector, a back-lighting photography technique has been employed to capture the instantaneous self-pulsated spray and stable spray images with a high speed camera. The diameter and velocity of the droplets in the spray have been characterized with a Dantec Phase Doppler Anemometry (PDA) system. The effects of self-pulsation on the spray pattern, primary breakup, spray angle, diameter and velocity distribution and mass flow rate distribution are analyzed and discussed. The results show that the spray morphology is greatly influenced by self-pulsation. The stable spray has a cone shape, while the self-pulsated spray looks like a Christmas tree. The main difference of these two sprays is the primary breakup. The liquid film of stable spray keeps stable while that of self-pulsated spray oscillates periodically. The film width of self-pulsated spray varies in a large range with ‘neck’ and 'shoulder’ features existing. The liquid film of self-pulsated spray breaks up at the second neck, and then the second shoulder begins to breakup into ligaments. The self-pulsated spray produces droplet clusters periodically, varies horizontal spray width and mass flux periodically. From the point of spatial distribution, self-pulsation is good for the spray, it uniformizes the mass flux along radius and increases the spray angle. However, when self-pulsation occurs, the SMD distribution varies from an inverted V shape to a hollow cone shape, and SMD increases at all the measuring points. Namely, from the point of atomization performance, self-pulsation has negative effects even when the breakup length is smaller. The effects of self-pulsation on the diameter and velocity distributions of the spray are mainly in the center part of the spray. The periphery of stable and self-pulsated spray has similar diameter and velocity distribution.

Wind tunnel investigation of horizontal and vertical sand fluxes of ascending and descending sand particles in aeolian sand transport

AbstractThe horizontal and vertical sand mass fluxes in aeolian sand transport are investigated in a wind tunnel by PTV (particle tracking velocimetry). According to the particle velocity and volume fraction of each individual particle from PTV images, the total horizontal sand mass flux, the horizontal mass fluxes of ascending and descending sand particles, and upward and downward vertical sand fluxes are analyzed. The results show that the horizontal mass fluxes of ascending and descending sand particles generally decrease with the increase of height and can be described by an exponential function above about 0.03 m height. At the same friction velocity, the decay heights of the total horizontal sand mass flux and the horizontal mass fluxes of ascending and descending sand particles are very similar. The proportion of horizontal mass flux of ascending sand particles is generally about 0.3–0.42, this means the horizontal mass flux of descending sand particles makes an important contribution to the total horizontal sand mass flux. Both the upward and downward vertical sand mass fluxes generally decrease with height and they are approximately equal at the same height and friction velocity. The relation between upward (or downward) vertical sand mass flux and horizontal sand mass flux can be described by a power function. The present study is used to help understand the transport of ascending and descending sand particles. Copyright © 2016 John Wiley & Sons, Ltd.

Effects of Recess on the Self-Pulsation Characteristics of Liquid-Centered Swirl Coaxial Injectors

The self-pulsation characteristics of liquid-centered swirl coaxial injectors have been investigated using a Photron Fastcam SA-1.1 camera and a PULSE system (B&K Corp., 3650 type). The effects of both recess length and injection condition were analyzed and discussed. The flow pattern in the recess chamber of a liquid-centered swirl coaxial injector can be classified into three types, namely, outer, critical, and inner mixing flow. Once the critical flow approaches, self-pulsation is about to occur. When self-pulsation occurs, the spray oscillates with a Christmas tree shape. The Christmas-tree-like spray periodically produces droplet clusters, making horizontal spray width and mass flux oscillate periodically. Self-pulsation is the strongest when the flow is around the critical mixing flow. Both variations of recess length and injection condition can change the flow pattern in the recess chamber. When the flow pattern moves away from the critical mixing flow, self-pulsation is suppressed. All the frequency behaviors of self-pulsation indicate that self-pulsation is strongly related to Kelvin–Helmholtz instability. A possible mechanism of self-pulsation is also presented. This mechanism interprets the experiments of the present study and the results of literature well.

A wind tunnel study of sand-cemented bodies on wind erosion intensity and sand transport

Wind tunnel experiments were used to test the capacity of sand-cemented bodies (SCB) on mulch beds. The total sand transport rate decreased as the level of SCB coverage increased. At higher SCB coverage (more than 40 %), the sand transport was basically unaffected by further increases in SCB coverage. While at low SCB coverage (less than 10 %), wind velocity played an important role in sand transport. Under the same SCB coverage, the sand transport depends on the increasing SCB size, due to the decrease in SCB density. The wind erosion intensity exponentially decreased with increasing SCB coverage (less than 40 %). The vertical profiles of horizontal mass flux from the SCB mulch–sand surface were also described by an exponential relationship. The vertical sand movement of particles was more sensitive to changes in SCB coverage at 20–40 %, compared with at less than 10 %. When the SCB coverage was more than 40 %, the decay rate of sand transport with height was nearly invariable. In summary, increases in SCB coverage had anti-erosion benefits for the underlying sand surface and could be considered for the development of a new type of sand fixation technology.

Threshold wind velocity dynamics as a driver of aeolian sediment mass flux

Abstract Horizontal (saltation) mass flux is a key driver of aeolian dust emission. Estimates of the horizontal mass flux underpin assessments of the global dust budget and influence our understanding of the dust cycle and its interactions. Current equations for predicting horizontal mass flux are based on limited field data and are constrained to representing transport-limited equilibrium saltation, driven by the wind momentum flux in excess of an entrainment threshold. This can result in large overestimation of the sediment mass flux. Here we compare measurements of the soil entrainment threshold, horizontal mass flux, and their temporal variability for five undisturbed dryland soils to explore the role of threshold in controlling the magnitude of mass flux. Average and median entrainment threshold showed relatively small variability among sites and relatively small variability between seasons, despite significant differences in soil surface conditions. Physical and biological soil crusts had little effect on the threshold value, and threshold appeared to play a minor role in determining the magnitude of sediment transport. Our results suggest that horizontal mass flux was controlled more by the supply limitation and abrasion efficiency of saltators present as loose erodible material or originating from neighboring soil sources. The omission of sediment supply and explicit representation of saltation bombardment from horizontal flux equations is inconsistent with the process representation in dust emission schemes and contributes to uncertainty in model predictions. This uncertainty can be reduced by developing greater process fidelity in models to predict horizontal mass flux under both supply- and transport-limited conditions.

Consistency of wind erosion assessments across land use and land cover types: A critical analysis

In recent decades, large areas of rangeland have been converted to cropland or vice versa in the western United States and elsewhere in the world, driven largely by increased crop prices, loss of access to irrigation water, and agricultural expansion/contraction. Wind erosion and dust emissions are key processes that have not been well studied during land use and associated land cover changes. This assessment is challenging because currently no model is available that can provide field- to landscape-scale estimates of wind erosion on both rangeland and cropland, and account for soil, vegetation and management changes. In this paper, we compare aeolian sediment transport estimates from available cropland models and a number of mass flux equations developed for rangelands, for a bare soil surface with different levels of crust and surface roughness under different wind speeds. Our results show that the simulated horizontal sediment mass fluxes are similar for cropland and rangeland models at large surface crust coverage and aerodynamic roughness. In situations of small to moderate crust cover and soil roughness, horizontal mass fluxes varied by over three orders of magnitude among the tested models. A correlation analysis shows that horizontal mass fluxes simulated by cropland and rangeland models are correlated, with correlation R 2 of 0.37–0.99 across different models. Finally, we propose an approach to estimate changes in

Numerical assessment of potential impacts of hydraulically fractured Bowland Shale on overlying aquifers

AbstractNatural gas extracted from hydraulically fractured shale formations potentially has a big impact on the global energy landscape. However, there are concerns of potential environmental impacts of hydraulic fracturing of the shale formations, particularly those related to water quality. To evaluate the potential impact of hydraulically fractured shale on overlying aquifers, we conduct realizations of numerical modeling simulations to assess fluid flow and chloride transport from a synthetic Bowland Shale over a period of 11,000 years. The synthetic fractured shale was represented by a three‐dimensional discrete fracture model that was developed by using the data from a Bowland Shale gas exploration in Lancashire, UK. Chloride mass exchange between fractures and the rock matrix was fully accounted for in the model. The assessment was carried out to investigate fluid and chloride mass fluxes before, during, and after hydraulic fracturing of the Bowland Shale. Impacts of the upward fracture height and aperture, as well as hydraulic conductivity of the multilayered bedrock system, are also included this assessment. This modeling revealed that the hydraulically fractured Bowland Shale is unlikely to pose a risk to its overlying groundwater quality when the induced fracture aperture is ≤200 µm. With the fracture aperture ≥1000 µm, the upward chloride flux becomes very sensitive to the upward fracture height growth and hydraulic conductivity of the multilayered bedrock system. In the extremely unlikely event of the upward fracture growth directly connecting the shale formation to the overlying Sherwood Sandstone aquifer with the fracture aperture ≥1000 µm, the upward chloride mass flux could potentially pose risks to the overlying aquifer in 100 years. The model study also revealed that the upward mass flux is significantly intercepted by the horizontal mass flux within a high permeable layer between the Bowland Shale and its overlying aquifers, reducing further upward flux toward the overlying aquifers.