共查询到20条相似文献，搜索用时 15 毫秒
1.
Heat as a groundwater tracer in shallow and deep heterogeneous media: Analytical solution,spreadsheet tool,and field applications 
下载免费PDF全文
下载免费PDF全文 Barret L. Kurylyk Dylan J. Irvine Sean K. Carey Martin A. Briggs Dale D. Werkema Mariah Bonham 《水文研究》2017,31(14):2648-2661
Groundwater flow advects heat, and thus, the deviation of subsurface temperatures from an expected conduction‐dominated regime can be analysed to estimate vertical water fluxes. A number of analytical approaches have been proposed for using heat as a groundwater tracer, and these have typically assumed a homogeneous medium. However, heterogeneous thermal properties are ubiquitous in subsurface environments, both at the scale of geologic strata and at finer scales in streambeds. Herein, we apply the analytical solution of Shan and Bodvarsson ( 2004 ), developed for estimating vertical water fluxes in layered systems, in 2 new environments distinct from previous vadose zone applications. The utility of the solution for studying groundwater‐surface water exchange is demonstrated using temperature data collected from an upwelling streambed with sediment layers, and a simple sensitivity analysis using these data indicates the solution is relatively robust. Also, a deeper temperature profile recorded in a borehole in South Australia is analysed to estimate deeper water fluxes. The analytical solution is able to match observed thermal gradients, including the change in slope at sediment interfaces. Results indicate that not accounting for layering can yield errors in the magnitude and even direction of the inferred Darcy fluxes. A simple automated spreadsheet tool (Flux‐LM) is presented to allow users to input temperature and layer data and solve the inverse problem to estimate groundwater flux rates from shallow (e.g., <1 m) or deep (e.g., up to 100 m) profiles. The solution is not transient, and thus, it should be cautiously applied where diel signals propagate or in deeper zones where multi‐decadal surface signals have disturbed subsurface thermal regimes. 相似文献
2.
A simple numerical model is presented for estimating vertical groundwater flux from transient subsurface temperature profiles obtained from field measurements. The model developed utilizes the MacCormack scheme, which is based on the Finite Difference Method (FDM), for solving the governing partial differential equation of convection–diffusion heat transport with appropriate initial and boundary conditions within the subsurface. In order to validate the model, numerical solutions obtained for the study area located in the Nagoka plain, Japan are compared with the published measured data and results obtained by others. Results obtained show good agreement and fit the observed data with a correlation coefficient, R2, of 0·88. The estimated groundwater flux is 1·85 × 10−7 m s−1. Sensitivity analyses were also carried out to investigate the effect of variations in groundwater fluxes, thermal properties and the annual thermal variability due to climatic changes on the transient subsurface temperature profiles and to have a better understanding of the subsurface thermal dynamics. A substantial effect of annual climatic variability is observed on the temporal distributions of temperature depth profiles, and a better estimate of thermal parameters is required to estimate vertical groundwater flux. The largest change in subsurface temperature depth profiles due to groundwater flux over a year is within ± 4 °C. The influence of groundwater flux on subsurface temperature distributions in space and time may be more pronounced in areas where the top of the saturated layer fluctuates considerably. Variation in thermal diffusivity results in temperature change up to ± 1·5% and may cause change in groundwater flux estimate by ± 18%. The model presented has merits over analytical solutions (type curve matching techniques) in terms of suitability and applicability to real field problems, and can be a good asset to hydrological models as quantifying groundwater recharge or deducing it from other quantities, such as rainfall, evapotranspiration and runoff, is often complicated. Copyright © 2007 John Wiley & Sons, Ltd. 相似文献
3.
Estimation of vertical water fluxes from temperature time series by the inverse numerical computer program FLUX‐BOT 下载免费PDF全文
下载免费PDF全文 The application of heat as a hydrological tracer has become a standard method for quantifying water fluxes between groundwater and surface water. The typical application is to estimate vertical water fluxes in the shallow subsurface beneath streams or lakes. For this purpose, time series of temperatures in the surface water and in the sediment are measured and evaluated by a vertical 1D representation of heat transport by advection and conduction. Several analytical solutions exist to calculate the vertical water flux from the measured temperatures. Although analytical solutions can be easily implemented, they are restricted to specific boundary conditions such as a sinusoidal upper temperature boundary. Numerical solutions offer higher flexibility in the selection of the boundary conditions. This, in turn, reduces the effort of data preprocessing, such as the extraction of the diurnal temperature variation from the raw data. Here, we present software to estimate water fluxes based on temperatures—FLUX‐BOT. FLUX‐BOT is a numerical code written in MATLAB that calculates vertical water fluxes in saturated sediments based on the inversion of measured temperature time series observed at multiple depths. FLUX‐BOT applies a centred Crank–Nicolson implicit finite difference scheme to solve the one‐dimensional heat advection–conduction equation. FLUX‐BOT includes functions for the inverse numerical routines, functions for visualizing the results, and a function for performing uncertainty analysis. We present applications of FLUX‐BOT to synthetic and to real temperature data to demonstrate its performance. 相似文献
4.
Use of shallow groundwater temperature profiles to infer climate and land use change: interpretation and measurement challenges 下载免费PDF全文
下载免费PDF全文 Temperature is often used to infer the effect of land use and climate conditions on aquifers. Reliable data are needed to examine the temperature behaviour in the subsurface; thus, the use of robust acquisition techniques is unavoidable. Three temperature measurement techniques were applied to assess the sources of bias that could occur during temperature logging in a shallow Quaternary coastal aquifer in Ferrara (Northern Italy). Open borehole temperature logging, multilevel sampling straddle packers isolated temperature measurements within a flow cell above ground and multilevel sampling straddle packers isolated temperature measurements via an in‐well level logger (MLS‐IW) were compared for several coastal monitoring wells to gain insights on the limitations of each technique. Results show that the source of bias between the three applied techniques are different: (i) the open borehole temperature logging method tends to record heat convection through the open borehole and is not representative of the aquifer temperature distribution; (ii) the multilevel sampling straddle packers isolated temperature measurements within a flow cell above ground method is swayed by the air temperature and the heating of the submersible pump used to lift groundwater above ground; and (iii) the MLS‐IW provides the most reliable vertical thermal profiling both in summer and winter, because groundwater temperature is directly measured at the selected monitoring depth. The implementation of a 1D flow model demonstrates that if precise temperature profiles are needed to infer the influence that land use and climate changes have on groundwater, the MLS‐IW method is a reliable method that could be applied to existing monitoring wells. Copyright © 2016 John Wiley & Sons, Ltd. 相似文献
5.
V. Čermák J. Šafanda L. Bodri M. Yamano E. Gordeev 《Studia Geophysica et Geodaetica》2006,50(4):675-695
To reconstruct the recent climate history in Kamchatka, a series of repeated precise temperature logs were performed in a
number of boreholes located in a broad east-west strip (between 52 and 54°N) in the central part of Kamchatka west of Petropavlovsk-Kamchatski.
Within three years more than 30 temperature logs were performed in 10 holes (one up to six logs per hole) to the depth of
up to 400 metres. Measured temperature gradients varied in a broad interval 0 to 60 mK/m and in some holes a sizeable variation
in the subsurface temperatures due to advective heat transport by underground water was observed. Measured data were compared
with older temperature profiles obtained in the early eighties by Sugrobov and Yanovsky (1993). Even when older data are of
poorer precision (accuracy of about 0.1 K), they presented valuable information of the subsurface temperature conditions existing
20–25 years ago. Borehole observations and the inverted ground surface temperature histories (GSTHs) used for the paleoclimate
reconstruction were complemented with a detailed survey of meteorological data. Namely, the long-term surface air temperature
(SAT) and precipitation records from Petropavlovsk station (in operation since 1890) were used together with similar data
from a number of local subsidiary meteo-stations operating in Central Kamchatka since 1950. Regardless of extreme complexity
of the local meteorological/climate conditions, diversity of borehole sites and calibration of measuring devices used during
the whole campaign, the results of the climate reconstruction supported a general warming of about 1 K characteristic for
the 20th century, which followed an inexpressive cooler period typical for the most of the 19th century. In the last three
to four decades the warming rate has been locally increasing up to 0.02 K/year. It was also shown that the snow cover played
a dominant role in the penetration of the climate “signal” to depth and could considerably smooth down the subsurface response
to the changes occurred on the surface. 相似文献
6.
We estimated the effects of hydrogeological and surface temperature warming on subsurface thermal regime from the temperature-depth profiles and hydrological data of groundwater quality both in the quaternary and tertiary systems in the Sendai Plain as a preliminary step toward reconstruction of climate changes.Annual mean air temperature in the plain has increased about 1.5 °C in the last 70 years and this surface warming resulted in low or negative thermal gradient. However, anomaly of thermal gradient was recognized in not all temperature-depth profiles. Groundwater chemical compositions and stable isotope data (δD and δ18O) show that the groundwater flow system has marked difference between those of tertiary and quaternary systems. Calculated results of three dimensional groundwater flow and heat transport model ensure the above hypothesis and shows that thermal gradient changes at close to basement of the quaternary system. The differences in groundwater flow systems are expressed as subsurface thermal gradient anomalies in the temperature-depth profiles in the Sendai Plain. Furthermore, one-dimensional numerical analyze including the effect of surface warming indicates that calculated profile has departure from steady state line at depths in 60-80 m agree well with observed one. 相似文献
7.
Aufeis (also known as icings) are large sheet-like masses of layered ice that form in river channels in arctic environments in the winter as groundwater discharges to the land surface and subsequently freezes. Aufeis are important sources of water for Arctic river ecosystems, bolstering late summer river discharge and providing habitat for caribou escaping insect harassment. The aim of this research is to use numerical simulations to evaluate a conceptual model of subsurface hydrogeothermal conditions that can lead to the formation of aufeis. We used a conceptual model based on geophysical data from the Kuparuk aufeis field on the North Slope of Alaska to develop a two-dimensional heterogeneous vertical profile model of groundwater flow, heat transport, and freeze/thaw dynamics. Modelling results showed that groundwater can flow to the land surface through subvertical high permeability pathways during winter months when the lower permeability soils near the land surface are frozen. The groundwater discharge can freeze on the surface, contributing to aufeis formation throughout the winter. We performed sensitivity analyses on subsurface properties and surface temperature and found that aufeis formation is most sensitive to the volume of unfrozen water available in the subsurface and the rate at which the subsurface water travels to the land surface. Although a trend of warming air temperatures will lead to a greater volume of unfrozen subsurface water, the aufeis volume can be reduced under warming conditions if the period of time for which air temperatures are below freezing is reduced. 相似文献
8.
We developed a numerical model, RFLUX, which uses the heat tracer method for vertical groundwater flux estimation, and applied it to the Leizhou Peninsula, South China, to provide information to inform local groundwater resource utilization and management. The temperature–depth (TD) profiles of 24 boreholes, along with the observed ground surface temperature (GST) and surface air temperature (SAT) series in recent decades, were collected in this area. Underground TD data demonstrated the capacity to identify groundwater flow patterns, and local GST and SAT data demonstrated a strong correlation with each other over monthly, seasonal, and annual scales. In the RFLUX model, the average GST and SAT series were applied as an upper boundary condition, and a nonlinear initial condition was set using an analytical solution from the literature. The model results of selected TD profiles demonstrated that the annual vertical groundwater flux was about 0.15 m a?1, which tended to be overestimated if a linear initial condition was used. This model can be easily applied with minor modifications, considering its clear purpose and simplicity. 相似文献
9.
Subsurface temperatures were analyzed to evaluate the effects of urbanization and groundwater flow in Osaka, the second largest metropolitan in Japan. The temperature-depth profiles were classified into three (recharge, intermediate and discharge) types based on the thermal gradient, shape of the profiles, and temperature itself. The locations of the three types of profiles represent the places of regional and local groundwater flow systems in Osaka plain. The analyses of subsurface temperature using heat conduction-convection theory under the condition of surface warming showed that the depth of minimum temperature in the profile increased with the groundwater recharge rate and time lapse from the beginning of surface warming due to urbanization. Comparisons between observed and estimated temperature profiles showed that the surface warming due to urbanization is larger and occurred earlier in the middle of the city than those in surrounding area, which agreed with the air temperature records. 相似文献
10.
Čermák Vladimír Šafanda Jan Krešl Milan Dědeček Petr Bodri Louise 《Studia Geophysica et Geodaetica》2000,44(3):430-441
Long-term (1961 – 1996) meteorological air temperature series together with the reconstructed ground surface temperature histories, obtained by inverting borehole temperature-depth profiles, were used to project regional patterns of the recent (climate) warming rate on the territory of the Czech Republic. The characteristic magnitude of the warming rate of 0.02 –0.03 K/yr was confirmed by the results of several years of monitoring the temperature in two experimental boreholes. The monitoring of shallow temperatures at depths of about 30 –40 m, i.e. below the reach of the seasonal surface temperature variations, can serve as an alternative tool of direct quantitative assessment of the present warming rate. The data also seem to sustain a potential man-made component contributing to the more pronounced recent warming rate observed in the areas of large agglomeration. 相似文献
11.
《Advances in water resources》1998,21(4):297-313
Diffusion of dissolved contaminants in multilayer porous media is an important phenomenon affecting both contamination and remediation in natural aqueous environments, including diffusion in groundwater aquitards and contaminated bed sediments. This study presents a new analytical solution for solute diffusion in a semi-infinite two-layer porous medium for arbitrary boundary and initial conditions. The solution was obtained by using the Green's function approach in the Laplace domain with the application of the binomial theorem to facilitate inversion back to the real time domain. Results based on this solution were found to be simple both in form and ease of calculation and to be in good agreement with those obtained using numerical calculations based on the Crank-Nicolson finite difference method. Applications of the solution are presented in the context of a contaminated groundwater aquitard to demonstrate how different boundary and initial conditions can greatly affect the contamination and decontamination of porous media, and to illustrate how diffusion modeling might be used in a forensic sense. 相似文献
12.
Long‐term heating of shallow urban aquifers is observed worldwide. Our measurements in the city of Cologne, Germany revealed that the groundwater temperatures found in the city centre are more than 5 K higher than the undisturbed background. To explore the role of groundwater flow for the development of subsurface urban heat islands, a numerical flow and heat transport model is set up, which describes the hydraulic conditions of Cologne and simulates the transient evolution of thermal anomalies in the urban ground. A main focus is on the influence of horizontal groundwater flow, groundwater recharge and trends in local ground warming. To examine heat transport in groundwater, a scenario consisting of a local hot spot with a length of 1 km of long‐term ground heating was set up in the centre of the city. Groundwater temperature‐depth profiles at upstream, central and downstream locations of this hot spot are inspected. The simulation results indicate that the main thermal transport mechanisms are long‐term vertical conductive heat input, horizontal advection and transverse dispersion. Groundwater recharge rates in the city are low (<100 mm a?1) and thus do not significantly contribute to heat transport into the urban aquifer. With groundwater flow, local vertical temperature profiles become very complex and are hard to interpret, if local flow conditions and heat sources are not thoroughly known. Copyright © 2014 John Wiley & Sons, Ltd. 相似文献
13.
Tamara Kolbe Jean Marçais Jean-Raynald de Dreuzy Thierry Labasque Kevin Bishop 《水文研究》2020,34(10):2176-2189
Large proportions of rainwater and snowmelt infiltrate into the subsurface before contributing to stream flow and stream water quality. Subsurface flow dynamics steer the transport and transformation of contaminants, carbon, weathering products and other biogeochemistry. The distribution of groundwater ages with depth is a key feature of these flow dynamics. Predicting these ages are a strong test of hypotheses about subsurface structures and time-varying processes. Chlorofluorocarbon (CFC)-based groundwater ages revealed an unexpected groundwater age stratification in a 0.47 km2 forested catchment called Svartberget in northern Sweden. An overall groundwater age stratification, representative for the Svartberget site, was derived by measuring CFCs from nine different wells with depths of 2–18 m close to the stream network. Immediately below the water table, CFC-based groundwater ages of already 30 years that increased with depth were found. Using complementary groundwater flow models, we could reproduce the observed groundwater age stratification and show that the 30 year lag in rejuvenation comes from return flow of groundwater at a subsurface discharge zone that evolves along the interface between two soil types. By comparing the observed groundwater age stratification with a simple analytical approximation, we show that the observed lag in rejuvenation can be a powerful indicator of the extent and structure of the subsurface discharge zone, while the vertical gradient of the age-depth-relationship can still be used as a proxy of the overall aquifer recharge even when sampled in the discharge zone. The single age stratification profile measured in the discharge zone, close to the aquifer outlet, can reveal the main structure of the groundwater flow pattern from recharge to discharge. This groundwater flow pattern provides information on the participation of groundwater in the hydrological cycle and indicates the lower boundary of hydrological connectivity. 相似文献
14.
15.
Water temperatures in mountain streams are likely to rise under future climate change, with negative impacts on ecosystems and water quality. However, it is difficult to predict which streams are most vulnerable due to sparse historical records of mountain stream temperatures as well as complex interactions between snowpack, groundwater, streamflow and water temperature. Minimum flow volumes are a potentially useful proxy for stream temperature, since daily streamflow records are much more common. We confirmed that there is a strong inverse relationship between annual low flows and peak water temperature using observed data from unimpaired streams throughout the montane regions of the United States' west coast. We then used linear models to explore the relationships between snowpack, potential evapotranspiration and other climate-related variables with annual low flow volumes and peak water temperatures. We also incorporated previous years' flow volumes into these models to account for groundwater carryover from year to year. We found that annual peak snowpack water storage is a strong predictor of summer low flows in the more arid watersheds studied. This relationship is mediated by atmospheric water demand and carryover subsurface water storage from previous years, such that multi-year droughts with high evapotranspiration lead to especially low flow volumes. We conclude that watershed management to help retain snow and increase baseflows may help counteract some of the streamflow temperature rises expected from a warming climate, especially in arid watersheds. 相似文献
16.
Jacek A. Majorowicz 《Pure and Applied Geophysics》1993,140(4):655-666
High quality temperature measurements have been made to depths of 30 to 220 m at 42 sites in 62 observational hydrogeological wells in Alberta. The temperature profiles commonly show near-surface inversions with a minimum temperature at depths of 30 to 50 m. Thermal modelling suggests a surface temperature history with warming reaching 2°C over the past 30 to 60 years. Recent climate warming evident from the analysis of the air temperature data in the region seems to provide at least a partial explanation of the increased ground temperatures. A sudden increase of the surface ground temperature caused by land clearing may be the other explanation, although modelling of such a sudden increase can only explain the observed temperature-depth data if the onset of such warming is 20–30 years old, which is in disagreement with the history of land development in the studied area. The effect of near-surface inversions of the temperature profiles also has been observed in the forested areas. The above support the climate based effect. The superposition of the climatic effect and man-made activity effect upon the ground warming is a very complicated process calling for considerably more research. 相似文献
17.
Water exchange across the sediment–water interface of streams impresses a characteristic thermal pattern at the interface. The use of fibre optic distributed temperature sensing at the sediment–water interface in a small sand‐bed stream identifies such temperature patterns. Groundwater and interflow can be differentiated based on the temporal evolution of temperature patterns. Additionally, sudden temperature changes at the sediment–water interface observed during the transit of floods enable spatial identification of local up and downwelling. Electromagnetic induction geophysics can detect subsurface texture structures that support groundwater–surface water exchange. Our results show that areas of permanent temperature anomalies observed with fibre optic distributed temperature sensing match areas of comparatively homogeneous electrical conductivity. This indicates groundwater discharge and enables differentiating groundwater discharge from interflow and local downwelling. 相似文献
18.
Thermal and isotopic evidence for surface and subsurface water contributions to baseflow in a high Arctic river 下载免费PDF全文
下载免费PDF全文 Seven longitudinal water temperature tow surveys were conducted to attempt to identify the location of surface and subsurface river water exchanges along the length of the West River at the Cape Bounty Arctic Watershed Observatory, Melville Island, Nunavut, Canada (74°55′ N, 109°35′ W). Water temperature data were collected using a calibrated thermistor with an accuracy of ±0.002 °C (resolution <0.00005 °C) along the river during July 2014 in conjunction with stable water isotope sampling to support the thermal data and to determine the extent of surface water mixing from different sources such as precipitation, snowmelt, and surface/subsurface water contributions to the river. Atmospheric conditions were found to be the main contributor to seasonal temperature variance in the river, whereas tributary inflows and residual channel snow also had important thermal effects to river temperatures. Residual channel snow was a sustained source of cold water during much of the 2014 summer season (June–August) and substantially offset downstream warming. The longitudinal temperature profiles indicate notable changes to the thermal state of the river, which are interpreted to be indicative of subsurface and surface water exchange through inputs of relatively cold or warm water. Broadly, surface inflows were found to provide warmer water relative to the West River, and contributed to downstream warming of the river, along with downstream enrichment of δD and δ18O. Subsurface inflows provided cooler water relative to the river, and contributed to downstream depletion of δD and δ18O and downstream cooling of river temperatures. These results demonstrate that localized changes in river temperature, in conjunction with isotopic tracers, can be used to track channel–slope water interactions in Arctic hydrological systems, work previously limited to alpine and temperate settings. 相似文献
19.
In contaminant hydrogeology, investigations at fractured rock sites are typically undertaken to improve understanding of the fracture networks and associated groundwater flow that govern past and/or future contaminant transport. Conventional hydrogeologic, geophysical, and hydrophysical techniques used to develop a conceptual model are often implemented in open boreholes under conditions of cross-connected flow. A new approach using high-resolution temperature (±0.001°C) profiles measured within static water columns of boreholes sealed using continuous, water-inflated, flexible liners (FLUTe™) identifies hydraulically active fractures under ambient (natural) groundwater flow conditions. The value of this approach is assessed by comparisons of temperature profiles from holes (100 to 200 m deep) with and without liners at four contaminated sites with distinctly different hydrogeologic conditions. The results from the lined holes consistently show many more hydraulically active fractures than the open-hole profiles, in which the influence of vertical flow through the borehole between a few fractures masks important intermediary flow zones. Temperature measurements in temporarily sealed boreholes not only improve the sensitivity and accuracy of identifying hydraulically active fractures under ambient conditions but also offer new insights regarding previously unresolvable flow distributions in fractured rock systems, while leaving the borehole available for other forms of testing and monitoring device installation. 相似文献
20.
Freezing characteristics were investigated for a sedge covered floating fen and spruce covered swamp located beside a shallow lake in the Western Boreal Forest of Canada. Thermal properties were measured in situ for one freeze‐thaw cycle, and for two freeze‐thaw cycles in laboratory columns. Thermal conductivity and liquid water content were related to a range of subsurface temperatures above and below the freezing thresholds, and clearly illustrate hysteresis between the freezing and thawing process. Thermal hysteresis occurs because of the large change in thermal conductivity between water and ice, high water content of the peat, and wide variation in pore sizes that govern ice formation. Field and laboratory results were combined to develop linear freezing functions, which were tested in a heat transfer model. For surface temperature boundary conditions, subsurface temperatures were simulated for the over‐winter period and compared with field measurements. Replication of the transient subsurface thermal regime required that freezing functions transition gradually from thawed to frozen state (spanning the ?0·25 to ?2 °C range) as opposed to a more abrupt step function. Subsurface temperatures indicate that the floating fen underwent complete phase change (from water to ice) and froze to approximately the same depth as lake ice thickness. Therefore, the floating fen peatland froze as a ‘shelf’ adjacent to the lake, whereas the spruce covered swamp had a higher capacity for thermal buffering, and subsurface freezing was both more gradual and limited in depth. These thermal properties, and the timing and duration of frozen state, are expected to control the interaction of water and nutrients between surface water and groundwater, which will be affected by changes in air temperature associated with global climate change. Copyright © 2009 John Wiley & Sons, Ltd. 相似文献