双探针型海底热流计的结构优化

本文在现有海底热流探针制作技术条件下,首先建立了脉冲式双探针海底测量单元的有限元数值模型,模拟获得多组参数下的温度-时间数据,作为“实测”数据,再用脉冲加热有限长线热源(PFLS)模型求解待测介质热导率及其相对误差上限(REλ-UL),并以REλ-UL最小为原则,对双探针热流计的结构进行优化.结果表明:(1)在不同探针脉冲强度(q)、温度测量误差(ΔTm)和探针长度(L)组合下,都存在最佳探针间距(Best_r),使得REλ-UL降到最低;(2)随着q增大或ΔTm减小,Best_r逐渐增大;(3)当q、ΔTm及探针半径(a)都给定时,Best_r与探针长度(L)呈线性正相关;(4)当a=1.0 mm,且q、ΔTm分别取为628.0~1100.0 J·m-1、0.5~1.0 mK,若L在20.0~42.0 mm之间时,则Best_r在18.0~30.0 mm之间,此时介质热导率相对误差上限可控制在5.5%以内,同时测量温度可在6 min内达到最大值,即脉冲加热开始后,温度测量只需约7 min,便可满足介质热导率的求解,这比目前常用的Lister型热流计所需海底测量时间缩短8 min左右.     

A non-linear perturbation model considering catchment wetness and its application in river flow forecasting

A non-linear perturbation model for river flow forecasting is developed, based on consideration of catchment wetness using an antecedent precipitation index (API). Catchment seasonality, of the form accounted for in the linear perturbation model (the LPM), and non-linear behaviour both in the runoff generation mechanism and in the flow routing processes are represented by a constrained non-linear model, the NLPM-API. A total of ten catchments, across a range of climatic conditions and catchment area magnitudes, located in China and in other countries, were selected for testing daily rainfall-runoff forecasting with this model. It was found that the NLPM-API model was significantly more efficient than the original linear perturbation model (the LPM). However, restriction of explicit non-linearity to the runoff generation process, in the simpler LMP-API form of the model, did not produce a significantly lower value of the efficiency in flood forecasting, in terms of the model efficiency index R².     

郯庐断裂带是热异常带吗:来自断裂带南段热流的约束

郯庐断裂带是东亚最重要的断裂之一,在构造地质、岩石矿物、地震分布等领域得到了广泛的关注,然而少有针对郯庐断裂带地热学的研究.本文作者对郯庐断裂带南段7口钻井进行了系统测温,采用光学扫描法测试了6口钻井的142块岩心和13块露头样品的热导率,获得了6个高质量的热流数据.结合已发表的热流和地震剖面资料,我们获得了郯庐断裂带南段的岩石圈结构.研究表明,郯庐断裂带南段现今地温梯度从南到北有增大的趋势,介于21.8~30.3℃·km^-1.大地热流值介于44.0~81.7 mW·m^-2,平均61.4±10.8 mW·m^-2,表现为正常大地热流的特征,并非热异常带.最大热流值出现在庐枞盆地的大陆科学钻探井ZK03处,上地壳浅部极高的生热层很可能是高热流的主要因素.沿郯庐断裂南段地震带浅源地震发震带的底界和350℃等温面耦合较好,指示了深部结构的差异.

     

Mantle heat flow and geotherms for major tectonic units in central Europe

The results of seismic measurements along the deep seismic sounding profile VII and terrestrial heat flow measurements used for construction of heat generation models for the crust in the Paleozoic Platform region, the Sudetic Mountains (Variscan Internides) and the European Precambrian Platform show considerable differences in mantle heat flow and temperatures. At the base of the crust variations from 440–510°C in the models of Precambrian Platform to 700–820°C for the Paleozoic Platform and the Variscan Internides (Sudets) are found. These differences are associated with considerable mantle heat flow variations.The calculated models show mantle heat flow of about 8.4–12.6 mW m^–2 for the Precambrian Platform and 31 mW m^–2 to 40.2 mW m^–2 for Paleozoic orogenic areas. The heat flow contribution originating from crustal radioactivity is almost the same for the different tectonic units (from 33.5 mW m^–2 to 37.6 mW m^–2). Considerable physical differences in the lower crust and upper mantle between the Precambrian Platform and the adjacent areas, produced by lateral temperature variations, could be expected. On the basis of carbon ratio data it can be concluded that the Carboniferous paleogeothermal gradient was much lower in the Precambrian Platform area than in the Paleozoic Platform region.     

南海东北部973剖面BSR及其热流特征

对南海东北部973项目采集的地震测线进行了处理,阐明了恒春海脊处天然气水合物似海底反射（BSR）的特征.不同类型剖面的分析表明,单道地震剖面可以揭示南海东北部地区BSR一定的分布,但地震处理对揭示BSR分布全貌起重要作用.南海东北部的BSR具有世界大陆边缘BSR典型的特征,切穿了沉积层理反射,与海底起伏大体平行,为一强振幅的负极性反射.该BSR特征明显,意味着南海东北部地区存在含天然气水合物沉积.利用甲烷水合物与多组分天然气水合物相平衡曲线,计算了稳定带底界的埋深,并与BSR深度进行了对比分析.无论是甲烷水合物稳定带底界还是多组分天然气水合物稳定带的底界,单一地温梯度计算的结果不可能与BSR深度在整个剖面上对应.可知本区横向上地温梯度变化较大.利用BSR资料估算了地温梯度并求得热流值.计算表明,地温梯度与热流值由西向东,随着离海沟距离的增大、离岛弧距离的减小而减小.BSR计算得到的热流值为28～64 mW/m2,与台湾西南实测热流值的结果基本可以对比.     

A model for heat and water balance estimation and its application to land use and climate variation

This paper presents a method, which has been recently developed in the Department of Agrometeorology of the Poznań Agricultural University, for estimating heat balance components based on standard meteorological data, plant development stage, and land-use conditions. Estimates of latent heat flux components made it possible to obtain values of areal evapotranspiration and surface runoff. The method was applied in four catchment areas in Central Poland under present climatic conditions as well as under new, expected climatic conditions caused by an increase of ‘greenhouse gases’ in the atmosphere. Some changes in land use were also considered. The results show significant changes in the water balance when comparing present and future conditions of precipitation and evapotranspiration.     

A heat pulse technique for the determination of small‐scale flow directions and flow velocities in the streambed of sand‐bed streams

The hyporheic zone (HZ) has the capability to eliminate and attenuate nutrients and contaminants in riverine systems. Biogeochemical reactions and the potential elimination of contaminants are strongly controlled by the flow paths and dynamics in the HZ. Nevertheless, an easily applicable method for the field determination of flow patterns in the HZ is still lacking. Therefore, a heat pulse technique, which traces the movement of a short heat pulse in the upper part of the HZ and other sand beds, was developed. Five rods are vertically driven into the sediment of the streambed; one rod with a heater as point source located in about 10‐cm sediment depth and four rods with four temperature sensors in 3 cm distance, arranged concentrically with 7 cm diameter around the heating rod. Subsequently, a heat pulse is applied and the resulting breakthrough curves are indicative of flow velocities and flow directions in the streambed. A rough data analysis procedure is also suggested. In addition, laboratory experiments were performed to test the heat pulse technique. These experiments were validated based on coupled numerical modelling of flow and heat transport. First field tests of the method prove that the method is easily applicable under field conditions. These first field tests showed highly complex flow patterns with flow velocities from 1·8 to 4·9 cm min^?1 and flow directions from parallel to surface flow to opposite to surface flow. This suggests the need for a robust method to quantify hyporheic flow patterns in situ. Copyright © 2011 John Wiley & Sons, Ltd.     

A quadratic element method for evaluating groundwater flow by the inversion of surface tilt with application to the Tono Area, Japan

By modifying a previous method with constant elements, we developed a quadratic element method for more accurately estimating groundwater flow by the inversion of tilt data. In this method: (1) a region of groundwater flow is divided into quadratic elements in which the change in groundwater volume per unit volume of rock (Δv) and the Skempton coefficient (B) vary in a quadratic manner with the coordinates, (2) the values of Δv are set to zero at the boundaries of the region of groundwater flow and (3) the sum of the squared second derivatives of Δv is adopted as a constraining condition that is weighted and added to the sum of the squared errors in tilt. First, analyses were performed for a flow model to determine the accuracy of this method for estimating groundwater flow and also to clarify the effect of the assumed size of a region of groundwater flow. These analyses showed that the quadratic element method proposed in this study gives a much better estimation of Δv than the constant element method and that a large region of groundwater flow should be assumed, rather than a small region, since the values of Δv at points outside of the actual region of groundwater flow are estimated to be nearly zero when a large region is assumed while these values are greatly overestimated when an excessively small region is assumed. Finally, the quadratic element method was applied to the site of the Mizunami Underground Research Laboratory in the Tono area, Japan. Inverse analyses were performed for tilt data measured by four tiltmeters with a resolution of 10⁻⁹ radians during the excavation of two shafts under the assumption that the rock mass is an isotropic and homogeneous half- space. The results showed that the method proposed in this study reproduced the tilt data very accurately. Thus, the distribution of Δv was estimated without sacrificing the reproducibility of the tilt data. The contour maps of B(1 + ν)Δv (ν: Poisson’s ratio) showed that the heterogeneous flow of groundwater occurred at the site and that groundwater volume decreased mainly in the area surrounded by two faults. The latter result is consistent with the finding obtained by previous investigations that these faults have low permeability in the direction perpendicular to the strike and may act as a flow barrier.     

A sensor array system for profiling moisture in unsaturated rock and soil

As the scope of hydrologic investigations extend deeper into the subsurface profile, and increasingly include fractured rock, there is a growing need for techniques that can accurately monitor saturation changes at a high spatial and temporal resolution in this environment. We have developed a technique, the Electrical Resistance Sensor Array System (ERSAS), to track moisture dynamics in vadose zone regions that include both fractured rock and soil. The performance of ERSAS was compared with the time domain reflectometry (TDR) technique under controlled and field conditions. We found that ERSAS was effective in determining patterns of saturation changes along vertical soil/rock profiles. Because of the small size of individual sensors, it was able to resolve travel times associated with a wetting front and peak saturation better than TDR. In addition, ERSAS is significantly cheaper than the TDR system, and the sensor arrays are relatively easier to install in the subsurface profile. Published in 2011 by John Wiley & Sons, Ltd.     

A stochastic model for 3-dimensional flow patterns in infiltration experiments

Modelling the 3-dimensional water flux at field scale is important for the design and the analysis of dye tracer experiments. Furthermore, it enables the estimation of the risk to groundwater by pollutants, and the visualisation and classification of flux patterns. A stochastic model is presented that allows for the modelling of a wide range of flow patterns in soils as they appear in dye tracer experiments. The leading idea is that infiltrating water runs along paths, not necessarily preferential ones, and water spreads into the soil uniformly from the paths into the matrix. The model is essentially based on a Poisson point process and three independent random fields. The point process defines the starting points of the paths at the surface. The values of two random fields determine the course of the paths. The third random field governs the depth of the infiltration front. As an extension of the model, we present two simulated examples for stratified soils.     

Effective elastic thicknesses of the lithosphere in the Central Iberian Peninsula from heat flow: Implications for the rheology of the continental lithospheric mantle

The traditional view of the rheology of the continental lithosphere, sometimes known as the “jelly sandwich model”, consists of a strong upper crust, a weak lower crust, and a strong upper lithospheric mantle. Some authors argue, however, that the lithospheric mantle is weak and contributes little to the total strength and the effective elastic thickness of the lithosphere; this weakness is claimed to be due to the mantle being wet or subjected to temperatures higher than usually believed. This paper uses the relationship between rheology of the lithosphere and heat flow to calculate theoretical effective elastic thicknesses for three regions of the central Iberian Peninsula (the Duero Basin, the Spanish Central System and the Tajo Basin), taking into account the contribution of the crust and the lithospheric mantle, for dry and wet rheologies. We found that a wet peridotite rheology for the lithospheric mantle is generally consistent with independent (based on Bouguer coherence or flexural modeling) estimates of the effective elastic thickness for the study area, whereas a dry peridotite rheology cannot be reconciled with them. Moreover, the contribution of the mantle to the bending moment of the lithosphere, and therefore to both the effective elastic thickness and the total strength of the lithosphere, is important, and it may even be the dominant contribution. Therefore, the jelly sandwich model may be considered valid for the central Iberian Peninsula.     

A novel approach for estimating karst groundwater recharge in mountainous regions and its application in Switzerland

A pragmatic and simple approach for estimating the groundwater recharge of karst aquifers in mountainous regions by extrapolation of the hydrological regimes of gauged and well‐documented systems is presented. Specific discharge rates are derived using annual precipitation and spring measurements by taking into account catchment size and elevation, which are assumed to be the dominant factors. Reference sites with high data reliability are used for calibration and regional extrapolation. This is performed with normalized values employing spatial precipitation deviations and correlation with the elevation of the catchment areas. A tiered step procedure provides minimum and maximum normalized gradients for the relationship between recharge quantity and elevation for karst regions. The normalized recharge can therefore be obtained and extrapolated for any location using the spatial precipitation variability to provide an estimate of annual groundwater recharge. The approach was applied to Switzerland (approximately 7500 km² of karst terrain situated between 200 and over 4000 m a.s.l.) using annual precipitation data from meteorological stations for the years 2000 to 2011. Results show that the average recharge rates of different Swiss karst domains range from 20 to 46 L/km²s, which corresponds to an infiltration ratio between 0.6 and 0.9 of total precipitation. Despite uncertainties inherent in the approach, these results provide a benchmark for renewable karst groundwater resources in Switzerland of about 8.4 km³/year. The approach can be applied to any other mountainous karst region, that is, where a clear relationship between elevation, precipitation and recharge can be assumed. Copyright © 2015 John Wiley & Sons, Ltd.     

A note on benchmarking of numerical models for density dependent flow in porous media

Verification of numerical models for density dependent flow in porous media (DDFPM) by the means of appropriate benchmark problems is a very important step in developing and using these models. Recently, Infinite Horizontal Box (IHB) problem was suggested as a possible benchmark problem for verification of DDFPM codes. IHB is based on Horton–Rogers–Lapwood (HRL) problem. Suitability of this problem for the benchmarking purpose has been investigated in this paper. It is shown that the wavelength of instabilities fails to be a proper criterion to be considered for this problem. However, the threshold of instability formation has been found to be appropriate for benchmarking purpose.     

A 2D finite volume model for bebris flow and its application to events occurred in the Eastern Pyrenees

FLATModel is a 2D finite volume code that contains several original approaches to improve debris-flow simulation. Firstly, FLATModel incorporates a ＂stop-and-go＂ technique in each cell to allow continuous collapses and remobilizations of the debris-flow mass. Secondly, flow velocity and consequently yield stress is directly associated with the type of rheology to improve boundary accuracy. Thirdly, a simple approach for entrainment is also included in the model to analyse the effect of basal erosion of debris flows. FLATMODEL was tested at several events that occurred in the Eastern Pyrenees and simulation results indicated that the model can represent rather well the different characteristics observed in the field.     

一种气泡释放过程连续监测方法及其在三峡水库香溪河CH4通量监测的应用

气泡释放是天然水体释放CH₄的主要途径之一,准确量化水体气泡释放量对于辨析其“汇、源”特性至关重要。自然水体释放气泡的不连续性、不确定性使得监测其过程较为困难。本研究针对水体气泡释放监测难题,通过改进倒置漏斗型气泡通量监测装置提出了一种气泡释放过程连续监测方法。本方法测量对象为定长时间监测水域释放气泡的体积,经室内外实验验证,其理论量程为3.6～132 mL/(m²·min),测量结果能够较好的表征10～40 m水深缓流水体气泡体积通量变化特征。运用该方法于2021年6—11月对三峡水库支流香溪河库湾开展CH4气泡通量连续监测,并分析不同环境因子对其产生的影响。结果表明:监测期间,研究水域CH₄气泡通量变化范围为0.02～8.13 mg/(m²·d),且各采样点间CH₄气泡通量呈现较高的时空变异性;CH₄气泡通量与水温、水体pH呈显著正相关关系,与水深及水体电导率呈显著负相关关系。其中,水深可能是决定水体是否通过气泡形式释放CH₄的重...     

ERORUN-STAFOR: A collaborative observatory for the multidisciplinary study of the critical zone processes in a tropical volcanic watershed including a Tropical Montane Cloud Forest

Tropical volcanic islands are biodiversity hotspots where the Critical Zone (CZ) still remains poorly studied. In such steep topographic environments associated with extreme climatic events (cyclones), deployment and maintenance of monitoring equipment is highly challenging. While a few Critical Zone Observatories (CZOS) are located in tropical volcanic regions, none of them includes a Tropical Montane Cloud Forest (TMCF) at the watershed scale. We present here the dataset of the first observatory from the French network of critical zone observatories (OZCAR) located in an insular tropical and volcanic context, integrating a ‘Tropical Montane Cloud Forest’: The ERORUN-STAFOR observatory. This collaborative observatory is located in the northern part of La Réunion island (Indian Ocean) within the 45.0 km² watershed of Rivière des Pluies (i.e., Rainfall river) which hosts the TMCF of Plaines des Fougères, one of the best preserved natural habitats in La Réunion Island. Since 2014, the ERORUN-STAFOR monitoring in collaboration with local partners collected a multidisciplinary dataset with a constant improvement of the instrumentation over time. At the watershed scale and in its vicinity, the ERORUN-STAFOR observatory includes 10 measurement stations covering the upstream, midstream and downstream part of the watershed. The stations record a total of 48 different variables through continuous (sensors) or periodic (sampling) monitoring. The dataset consists of continuous time series variables related to (i) meteorology, including precipitation, air temperature, relative humidity, wind speed and direction, net radiation, atmospheric pressure, cloud water flux, irradiance, leaf wetness and soil temperature, (ii) hydrology, including water level and temperature, discharge and electrical conductivity (EC) of stream, (iii) hydrogeology, including (ground)water level, water temperature and EC in two piezometers and one horizontally drilled groundwater gallery completed by soil moisture measurements under the canopy. The dataset is completed by periodic time series variables related to (iv) hydrogeochemistry, including field parameters and water analysis results. The periodic sampling survey provides chemical and isotopic compositions of rainfall, groundwater, and stream water at different locations of this watershed. The ERORUN-STAFOR monitoring dataset extends from 2014 to 2022 with an acquisition frequency from 10 min to hourly for the sensor variables and from weekly to monthly frequency for the sampling. Despite the frequent maintenance of the monitoring sites, several data gaps exist due to the remote location of some sites and instrument destruction by cyclones. Preliminary results show that the Rivière des Pluies watershed is characterized by high annual precipitation (>3000 mm y⁻¹) and a fast hydrologic response to precipitation (≈2 h basin lag time). The long-term evolution of the deep groundwater recharge is mainly driven by the occurrence of cyclone events with a seasonal groundwater response. The water chemical results support existing hydrogeological conceptual models suggesting a deep infiltration of the upstream infiltrated rainfall. The TMCF of Plaine des Fougères shows a high water storage capacity (>2000% for the Bryophytes) that makes this one a significant input of water to groundwater recharge which still needs to be quantified. This observatory is a unique research site in an insular volcanic tropical environment offering three windows of observation for the study of critical zone processes through upstream-midstream-downstream measurements sites. This high-resolution dataset is valuable to assess the response of volcanic tropical watersheds and aquifers at both event and long-term scales (i.e., global change). It will also provide insights in the hydrogeological conceptual model of volcanic islands, including the significant role of the TMCFs in the recharge processes as well as the watershed hydrosedimentary responses to extreme climatic events and their respective evolution under changing climatic conditions. All data sets are available at https://doi.org/10.5281/zenodo.7983138 .