遂昌麻车─后山垵糜棱岩的发现及其意义

作者在野外工作中，于遂昌─大柘变质岩断块（下称遂昌变质岩断块）北西侧，麻车─后山垵等地的地表或探矿巷道内，相继发现了多处糜棱岩化岩石露头，认为这一带存在一条较大的韧性剪切带。经对该带的地质，地球化学及金银矿分布的特征的综合分析研究，指出该韧性剪切带有关的金银矿化活动可分为两个期次，即早期金在韧性剪切活动中初始富集，晚期的脆性变形在局部地段形成金矿体。提出在该韧性剪切带找矿工作中，应重视继承韧性剪切带的脆性断裂和断裂带内的浅色蚀变组合及矿化特征。     

精密水准与GNSS垂向形变关系研究——以天山构造带为例

以天山构造带区域为例,分别用精密水准和GNSS给出地壳形变垂直运动速率矢量图,对其差异及差异形成的原因进行深入分析和探讨,并基于优势互补原则,给出二者融合的形变图像。根据多年研究经验认为,该方法得到的垂直形变图在强震中长期预测中具有较为重要的应用价值。     

仙女山-九畹溪断裂带附近地区地震活动性与静态库仑应力变化研究

基于仙女山-九畹溪断裂带附近地区地质构造,讨论三峡水库蓄水前后该断裂带附近地区的地震活动特征,同时分析该区域2014-03 M4.5与M4.7地震间的触发关系及2次地震对后续小震的影响。结果表明：1）M4.7地震可能由M4.5地震与库水渗流产生的较大孔隙压力共同触发。2）M4.5与M4.7两次地震的应力扰动导致研究区内地震活动增加,后续地震中约66.9%的地震发生于库仑应力增强区。3）2次地震产生的静态库仑应力变化导致仙女山断裂带附近地区地震活动增加,后续地震活动水平将处于背景地震活动水平之上;九畹溪断裂带活动性相对较弱,略小于背景地震活动水平,且处于应力积累阶段。     

Numerical modeling of fluid‐driven fractures in cohesive poroelastoplastic continuum

In this article, we investigate the main parameters that influence the propagation of a fluid‐driven fracture in a poroelastoplastic continuum. These parameters include the cohesive zone, the stress anisotropy, and the pore pressure field. The fracture is driven in a permeable porous domain that corresponds to weak formation by pumping of an incompressible viscous fluid at the fracture inlet under plane strain conditions. Rock deformation is modeled with the Mohr–Coulomb yield criterion with associative flow rule. Fluid flow in the fracture is modeled by the lubrication theory. The movement of the pore fluid in the surrounding medium is assumed to obey the Darcy law and is of the same nature as the fracturing fluid. The cohesive zone approach is used as the fracture propagation criterion. The problem is modeled numerically with the finite element method to obtain the solution for the fracture length, the fracture opening, and the propagation pressure as a function of the time and distance from the pumping inlet. It is demonstrated that the plastic yielding that is associated with the rock dilation in an elastoplastic saturated porous continuum is significantly affected by the cohesive zone characteristics, the stress anisotropy, and the pore pressure field. These influences result in larger fracture profiles and propagation pressures due to the larger plastic zones that are developing during the fracture propagation. Furthermore, it is also found that the diffusion process that is a major mechanism in hydraulic fracture operations influences further the obtained results on the fracture dimensions, plastic yielding, and fluid pressures. These findings may explain partially the discrepancies in net pressures between field measurements and conventional model predictions. Copyright © 2012 John Wiley & Sons, Ltd.     

The role of lithology,catchment size and the alluvial zone on the hydrogeochemistry of two intermittent Mediterranean streams

Hydrology and solute concentrations of two intermittent Mediterranean streams draining two nested catchments were compared. The two catchments were mainly underlain by granitic rocks and different types of sericitic schists. Only the lowland catchment had an alluvial zone and a well‐developed riparian forest. The rainfall–runoff relationship and the correlation between daily flow concentrations showed that hydrological behaviour was similar at both sites during most of the year. However, reverse fluxes were detected during the wetting and drying up periods only in the stream with an alluvial zone. The intermittence in stream flow also had effects on absolute solute concentrations, temporal solute dynamics and streamwater stoichiometry. Streamwater chemistry was not affected by drainage area, except for cations produced mainly by bedrock dissolution (i.e. calcium and magnesium) that increased with increasing catchment size. Differences in the relationship among cations and anions were detected between the two streams, which could be attributed to biogeochemical processes occurring in the alluvial zone. The multivariate model used in this study showed that stoichiometry was more useful than absolute concentrations when analyzing the influence of different lithologies on streamwater chemistry. Such differences were amplified in autumn, likely due to a low hydrological connectivity between the two nested catchments. Copyright © 2007 John Wiley & Sons, Ltd.     

Spatio‐temporal distribution of near‐surface and root zone soil moisture at the catchment scale

Soil moisture is highly variable both spatially and temporally. It is widely recognized that improving the knowledge and understanding of soil moisture and the processes underpinning its spatial and temporal distribution is critical. This paper addresses the relationship between near‐surface and root zone soil moisture, the way in which they vary spatially and temporally, and the effect of sampling design for determining catchment scale soil moisture dynamics. In this study, catchment scale near‐surface (0–50 mm) and root zone (0–300 mm) soil moisture were monitored over a four‐week period. Measurements of near‐surface soil moisture were recorded at various resolutions, and near‐surface and root zone soil moisture data were also monitored continuously within a network of recording sensors. Catchment average near‐surface soil moisture derived from detailed spatial measurements and continuous observations at fixed points were found to be significantly correlated (r² = 0·96; P = 0·0063; n = 4). Root zone soil moisture was also found to be highly correlated with catchment average near‐surface, continuously monitored (r² = 0·81; P < 0·0001; n = 26) and with detailed spatial measurements of near‐surface soil moisture (r² = 0·84). The weaker relationship observed between near‐surface and root zone soil moisture is considered to be caused by the different responses to rainfall and the different factors controlling soil moisture for the soil depths of 0–50 mm and 0–300 mm. Aspect is considered to be the main factor influencing the spatial and temporal distribution of near‐surface soil moisture, while topography and soil type are considered important for root zone soil moisture. The ability of a limited number of monitoring stations to provide accurate estimates of catchment scale average soil moisture for both near‐surface and root zone is thus demonstrated, as opposed to high resolution spatial measurements. Similarly, the use of near‐surface soil moisture measurements to obtain a reliable estimate of deeper soil moisture levels at the small catchment scale was demonstrated. Copyright © 2007 John Wiley & Sons, Ltd.     

Assess hydrological responses to a warming climate at the Lysina Critical Zone Observatory in Central Europe

Climate warming is having profound effects on the hydrological cycle by increasing atmospheric demand, changing water availability, and snow seasonality. Europe suffered three distinct heat waves in 2019, and 11 of the 12 hottest years ever recorded took place in the past two decades, which will potentially change seasonal streamflow patterns and long-term trends. Central Europe exhibited six dry years in a row since 2014. This study uses data from a well-documented headwater catchment in Central Europe (Lysina) to explore hydrological responses to a warming climate. We applied a lumped parameter hydrologic model Brook90 and a distributed model Penn State Integrated Hydrologic Model (PIHM) to simulate long-term hydrological change under future climate scenarios. Both models performed well on historic streamflow and in agreement with each other according to the catchment water budget. In addition, PIHM was able to simulate lateral groundwater redistribution within the catchment validated by the groundwater table dynamics. The long-term trends in runoff and low flow were captured by PIHM only. We applied different EURO-CORDEX models with two emission scenarios (Representative Concentration Pathways RCP 4.5, 8.5) and found significant impacts on runoff and evapotranspiration (ET) for the period of 2071–2100. Results from both models suggested reduced runoff and increased ET, while the monthly distribution of runoff was different. We used this catchment study to understand the importance of subsurface processes in projection of hydrologic response to a warming climate.     

Using a H∞ filter assimilation procedure to estimate root zone soil water content

Root zone soil water content impacts plant water availability, land energy and water balances. Because of unknown hydrological model error, observation errors and the statistical characteristics of the errors, the widely used Kalman filter (KF) and its extensions are challenged to retrieve the root zone soil water content using the surface soil water content. If the soil hydraulic parameters are poorly estimated, the KF and its extensions fail to accurately estimate the root zone soil water. The H‐infinity filter (HF) represents a robust version of the KF. The HF is widely used in data assimilation and is superior to the KF, especially when the performance of the model is not well understood. The objective of this study is to study the impact of uncertain soil hydraulic parameters, initial soil moisture content and observation period on the ability of HF assimilation to predict in situ soil water content. In this article, we study seven cases. The results show that the soil hydraulic parameters hold a critical role in the course of assimilation. When the soil hydraulic parameters are poorly estimated, an accurate estimation of root soil water content cannot be retrieved by the HF assimilation approach. When the estimated soil hydraulic parameters are similar to actual values, the soil water content at various depths can be accurately retrieved by the HF assimilation. The HF assimilation is not very sensitive to the initial soil water content, and the impact of the initial soil water content on the assimilation scheme can be eliminated after about 5–7 days. The observation interval is important for soil water profile distribution retrieval with the HF, and the shorter the observation interval, the shorter the time required to achieve actual soil water content. However, the retrieval results are not very accurate at a depth of 100 cm. Also it is complex to determine the weighting coefficient and the error attenuation parameter in the HF assimilation. In this article, the trial‐and‐error method was used to determine the weighting coefficient and the error attenuation parameter. After the first establishment of limited range of the parameters, ‘the best parameter set’ was selected from the range of values. For the soil conditions investigated, the HF assimilation results are better than the open‐loop results. Copyright © 2010 John Wiley & Sons, Ltd.     

Copper‐gold mineralisation in New Guinea: Numerical modelling of collision,fluid flow and intrusion‐related hydrothermal systems

Two‐ and three‐dimensional numerical modelling techniques, constrained by key geodynamic data, provide insights into the controls on development of porphyry‐related Cu–Au mineralisation in the Tertiary collision zone of New Guinea. Modelling shows that the creation of local dilation to facilitate magma emplacement can be caused by reactivation of arc‐normal transfer faults, where they cut the weakened fold belt. Additionally, dilation occurs where fluid overpressuring is caused by collision‐related, south‐directed fluid flow being localised into the more permeable units of the Mesozoic passive‐margin sedimentary succession. Rapid uplift and erosion, which may be a mechanism for magmatic fluid release in these systems, is shown to be greatest in the west of West Papua, where the stronger Australian crust acts as a buttress. Within the Papuan Fold Belt, uplift is greatest near the margins, where the weaker fold belt abuts the stronger crust and/or major faults have been reactivated. Increased orographically induced precipitation and erosion exposes the lower parts of the stratigraphy within or on the margins of these uplifted zones. On a smaller scale, 2–D coupled fluid‐flow ‐ thermal‐chemical modelling uses a scenario of fluid mixing to calculate metal precipitation distribution and magnitude around an individual intrusive complex. Modelling highlights the interdependence of the spatial permeability structure, the regional temperature gradient, and the geometry of the convection cells and how this impacts on the distribution of metal precipitation.