The Dunhuang Basin, a typical inland basin in northwestern China, suffers a net loss of groundwater and the occasional disappearance of the Crescent Lake. Within this region, the groundwater/surface-water interactions are important for the sustainability of the groundwater resources. A three-dimensional transient groundwater flow model was established and calibrated using MODFLOW 2000, which was used to predict changes to these interactions once a water diversion project is completed. The simulated results indicate that introducing water from outside of the basin into the Shule and Danghe rivers could reverse the negative groundwater balance in the Basin. River-water/groundwater interactions control the groundwater hydrology, where river leakage to the groundwater in the Basin will increase from 3,114?×?104 m3/year in 2017 to 11,875?×?104 m3/year in 2021, and to 17,039?×?104 m3/year in 2036. In comparison, groundwater discharge to the rivers will decrease from 3277?×?104 m3/year in 2017 to 1857?×?104 m3/year in 2021, and to 510?×?104 m3/year by 2036; thus, the hydrology will switch from groundwater discharge to groundwater recharge after implementing the water diversion project. The simulation indicates that the increased net river infiltration due to the water diversion project will raise the water table and then effectively increasing the water level of the Crescent Lake, as the lake level is contiguous with the water table. However, the regional phreatic evaporation will be enhanced, which may intensify soil salinization in the Dunhuang Basin. These results can guide the water allocation scheme for the water diversion project to alleviate groundwater depletion and mitigate geo-environmental problem. 相似文献
Submarine groundwater discharge (SGD) is an important pathway for groundwater and associated chemicals to discharge to the sea. Groundwater levels monitored along a transect perpendicular to the shoreline are used to calculate SGD flux from the nearshore aquifer to Tolo Harbor, Hong Kong (China). The calculated SGD flux—recharge/discharge measured with Darcy’s Law methods—agrees well with estimates based on geo-tracer techniques and seepage meter in Tolo Harbor during previous studies. The estimated freshwater SGD is 1.69–2.0 m2/d at the study site and 0.3?±?0.04 cm/d for the whole of Tolo Harbor, which is comparable to the river discharge (0.25?±?0.07 cm/d) and precipitation (0.45?±?0.15 cm/d). The tide-driven SGD in the intertidal zone is 13.98–17.59 m2/d at the study site and 2.42?±?0.56 cm/d for the whole of Tolo Harbor. The SGD occurring in the subtidal zone and the bottom of Tolo Harbor is 3.12?±?4.63 cm/d. Fresh SGD accounts for ~5% of the total SGD, while the rest (~95%) is contributed by saline SGD driven by various forces. About 96% of the tide-driven SGD in the intertidal zone occurs in the ebbing tide period because the head difference between the groundwater level and sea level is great during this period. Tide-driven SGD in the spring tide is ~1.2 times that during neap tide. The tidal fluctuation amplitude and tide-driven SGD in the intertidal zone are positively correlated to each other; thus, a spring neap variation of the tide-driven SGD is observed. 相似文献
This study investigates the impact of model complexity and multi-scale prior hydrogeological data on the interpretation of pumping test data in a dual-porosity aquifer (the Chalk aquifer in England, UK). In order to characterize the hydrogeological properties, different approaches ranging from a traditional analytical solution (Theis approach) to more sophisticated numerical models with automatically calibrated input parameters are applied. Comparisons of results from the different approaches show that neither traditional analytical solutions nor a numerical model assuming a homogenous and isotropic aquifer can adequately explain the observed drawdowns. A better reproduction of the observed drawdowns in all seven monitoring locations is instead achieved when medium and local-scale prior information about the vertical hydraulic conductivity (K) distribution is used to constrain the model calibration process. In particular, the integration of medium-scale vertical K variations based on flowmeter measurements lead to an improvement in the goodness-of-fit of the simulated drawdowns of about 30%. Further improvements (up to 70%) were observed when a simple upscaling approach was used to integrate small-scale K data to constrain the automatic calibration process of the numerical model. Although the analysis focuses on a specific case study, these results provide insights about the representativeness of the estimates of hydrogeological properties based on different interpretations of pumping test data, and promote the integration of multi-scale data for the characterization of heterogeneous aquifers in complex hydrogeological settings. 相似文献
Scattered data interpolation schemes using kriging and radial basis functions (RBFs) have the advantage of being meshless and dimensional independent; however, for the datasets having insufficient observations, RBFs have the advantage over geostatistical methods as the latter requires variogram study and statistical expertise. Moreover, RBFs can be used for scattered data interpolation with very good convergence, which makes them desirable for shape function interpolation in meshless methods for numerical solution of partial differential equations. For interpolation of large datasets, however, RBFs in their usual form, lead to solving an ill-conditioned system of equations, for which, a small error in the data can cause a significantly large error in the interpolated solution. In order to reduce this limitation, we propose a hybrid kernel by using the conventional Gaussian and a shape parameter independent cubic kernel. Global particle swarm optimization method has been used to analyze the optimal values of the shape parameter as well as the weight coefficients controlling the Gaussian and the cubic part in the hybridization. Through a series of numerical tests, we demonstrate that such hybridization stabilizes the interpolation scheme by yielding a far superior implementation compared to those obtained by using only the Gaussian or cubic kernels. The proposed kernel maintains the accuracy and stability at small shape parameter as well as relatively large degrees of freedom, which exhibit its potential for scattered data interpolation and intrigues its application in global as well as local meshless methods for numerical solution of PDEs. 相似文献
Steelmaking-coal waste rock placed in mountain catchments in the Elk Valley, British Columbia, Canada, drain constituents of interest (CIs) to surface water downgradient of the waste rock dumps. The role of groundwater in transporting CIs in the headwaters of mountain catchments is not well understood. This study characterizes the physical hydrogeology of a portion of a 10-km2 headwater catchment (West Line Creek) downgradient of a 2.7-km2 waste rock dump placed over a natural headwater valley-bottom groundwater system. The study site was instrumented with 13 monitoring wells. Drill core samples were collected to determine subsurface lithology and geotechnical properties. The groundwater system was characterized using field testing and water-level monitoring. The valley-bottom sediments were composed of unconsolidated glacial and meltwater successions (<64 m thick) deposited as a series of cut and fill structures overlying shale bedrock. An unconfined basal alluvial aquifer located above fractured bedrock was identified as the primary conduit for groundwater flow toward Line Creek (650 m from the toe of the dump). Discharge through the basal alluvial aquifer was estimated using the geometric mean hydraulic conductivity (±1 standard deviation). These calculations suggest groundwater discharge could account for approximately 15% (ranging from 2 to 60%) of the total water discharged from the watershed. The residence time from the base of the waste rock dump to Line Creek was estimated at <3 years. The groundwater system was defined as a snowmelt (i.e., nival) regime dominated by direct recharge (percolation of precipitation) across the catchment. 相似文献
An accurate estimate of the groundwater inflow to a tunnel is one of the most challenging but essential tasks in tunnel design and construction. Most of the numerical or analytical solutions that have been developed ignore tunnel seepage conditions, material properties and hydraulic-head changes along the tunnel route during the excavation process, leading to inaccurate prediction of inflow rates. A method is introduced that uses MODFLOW code of GMS software to predict inflow rate as the tunnel boring machine (TBM) gradually advances. In this method, the tunnel boundary condition is conceptualized and defined using Drain package, which is simulated by dividing the drilling process into a series of successive intervals based on the tunnel excavation rates. In addition, the drain elevations are specified as the respective tunnel elevations, and the conductance parameters are assigned to intervals, depending on the TBM type and the tunnel seepage condition. The Qomroud water conveyance tunnel, located in Lorestan province of Iran, is 36 km in length. Since the Qomroud tunnel involved groundwater inrush during excavating, it is considered as a good case study to evaluate the presented method. The groundwater inflow to this tunnel during the TBM advance is simulated using the proposed method and the predicted rates are compared with observed rates. The results show that the presented method can satisfactorily predict the inflow rates as the TBM advances. 相似文献
Intensive irrigated agriculture substantially modifies the hydrological cycle and often has major environmental impacts. The article focuses upon a specific concern—the tendency for progressive long-term increases in the salinity of groundwater recharge derived from irrigated permeable soils and replenishment of unconfined aquifers in more arid regions. This process has received only scant attention in the water-resource literature and has not been considered by agricultural science. This work makes an original contribution by analysing, from scientific principles, how the salinisation of groundwater recharge arises and identifies the factors affecting its severity. If not proactively managed, the process eventually will impact irrigation waterwell salinity, the productivity of agriculture itself, and can even lead to land abandonment. The types of management measure required for mitigation are discussed through three detailed case histories of areas with high-value groundwater-irrigated agriculture (in Spain, Argentina and Pakistan), which provide a long-term perspective on the evolution of the problem over various decades. 相似文献
Serpentine soils and ultramafic laterites develop over ultramafic bedrock and are important geological materials from environmental, geochemical, and industrial standpoints. They have naturally elevated concentrations of trace metals, such as Ni, Cr, and Co, and also high levels of Fe and Mg. Minerals host these trace metals and influence metal mobility. Ni in particular is an important trace metal in these soils, and the objective of this research was to use microscale (µ) techniques to identify naturally occurring minerals that contain Ni and Ni correlations with other trace metals, such as Fe, Mn, and Cr. Synchrotron based µ-XRF, µ-XRD, and µ-XAS were used. Ni was often located in the octahedral layer of serpentine minerals, such as lizardite, and in other layered phyllosilicate minerals with similar octahedral structure, such as chlorite group minerals including clinochlore and chamosite. Ni was also present in goethite, hematite, magnetite, and ferrihydrite. Goethite was present with lizardite and antigorite on the micrometer scale. Lizardite integrated both Ni and Mn simultaneously in its octahedral layer. Enstatite, pargasite, chamosite, phlogopite, and forsterite incorporated various amounts of Ni and Fe over the micrometer spatial scale. Ni content increased six to seven times within the same 500 µm µ-XRD transect on chamosite and phlogopite. Data are shown down to an 8 µm spatial scale. Ni was not associated with chromite or zincochromite particles. Ni often correlated with Fe and Mn, and generally did not correlate with Cr, Zn, Ca, or K in µ-XRF maps. A split shoulder feature in the µ-XAS data at 8400 eV (3.7 Å?1 in k-space) is highly correlated (94% of averaged LCF results) to Ni located in the octahedral sheet of layered phyllosilicate minerals, such as serpentine and chlorite-group minerals. A comparison of bulk-XAS LCF to averaged µ-XAS LCF results showed good representation of the bulk soil via the µ-XAS technique for two of the three soils. In the locations analyzed by µ-XAS, average Ni speciation was dominated by layered phyllosilicate and serpentine minerals (76%), iron oxides (18%), and manganese oxides (9%). In the locations analyzed by µ-XRD, average Ni speciation was dominated by layered phyllosilicate, serpentine, and ultramafic-related minerals (71%) and iron oxides (17%), illustrating the complementary nature of these two methods.
Estuaries act as an organic matter and nutrient filter in the transition between the land, rivers and the ocean. In the past, high nutrient and organic carbon load and low oxygen concentration made the Elbe River estuary (NW Europe) a sink for dissolved inorganic nitrogen. A recent reduction in loads and subsequent recovery of the estuary changed its biogeochemical function, so that nitrate is no longer removed on its transition towards the coastal North Sea. Nowadays in the estuary, nitrification appears to be a significant nitrate source. To quantify nitrification and determine actively nitrifying regions in the estuary, we measured the concentrations of ammonium, nitrite and nitrate, the dual stable isotopes of nitrate and net nitrification rates in the estuary on five cruises from August 2012 to August 2013. The nitrate concentration increased markedly downstream of the port of Hamburg in summer and spring, accompanied by a decrease of nitrate isotope values that was clearest in summer exactly at the location where nitrate concentration started to increase. Ammonium and nitrite peaked in the Hamburg port region (up to 18 and 8 μmol L?1, respectively), and nitrification rates in this region were up to 7 μmol L?1 day?1. Our data show that coupled re-mineralization and nitrification are significant internal nitrate sources that almost double the estuary’s summer nitrate concentration. Furthermore, we find that the port of Hamburg is a hot spot of nitrification, whereas the maximum turbidity zone (MTZ) only plays a subordinate role in turnover of nitrate. 相似文献
We developed a spatially explicit, individual-based model to analyze how hypoxia effects on reproduction, growth, and mortality of Atlantic croaker in the northwestern Gulf of Mexico lead to population-level responses. The model follows the hourly growth, mortality, reproduction, and movement of individuals on a 300 × 800 spatial grid of 1-km2 cells for 140 years. Chlorophyll-a concentration, water temperature, and dissolved oxygen (DO) were specified daily for each grid cell and repeated for each year of the simulation. A bioenergetics model was used to represent growth, mortality was assumed stage- and age-dependent, and the movement behavior of juveniles and adults was modeled based on temperature and avoidance of low DO. Hypoxia effects were imposed using exposure effect submodels that converted time-varying exposures to low DO to reduced hourly growth, increased hourly mortality, and reduced annual fecundity. Results showed that 100 years of either mild or intermediate hypoxia produced small reductions in population abundance, while repeated severe hypoxia caused a 19% reduction in long-term population abundance. Relatively few individuals were exposed to low DO each hour, but many individuals experienced some exposure. The response was dominated by a 5% average reduction in annual fecundity of individuals. Under conditions of random sequences of mild, intermediate, and severe hypoxia years occurring in proportion to their historical frequency, the model predicted a 10% decrease in the long-term population abundance of croaker. A companion paper substitutes hourly DO values from a three-dimensional water quality model for the idealized hypoxia and results in a more realistic population reduction of about 25%. 相似文献