Rainfall intensity affects runoff responses in a semi-arid catchment

The source and hydrochemical makeup of a stream reflects the connectivity between rainfall, groundwater, the stream, and is reflected to water quantity and quality of the catchment. However, in a semi-arid, thick, loess covered catchment, temporal variation of stream source and event associated behaviours are lesser known. Thus, the isotopic and chemical hydrographs in a widely distributed, deep loess, semi-arid catchment of the northern Chinese Loess Plateau were characterized to determine the source and hydrochemical behaviours of the stream during intra-rainfall events. Rainfall and streamflow were sampled during six hydrologic events coupled with measurements of stream baseflow and groundwater. The deuterium isotope (²H), major ions (Cl⁻, SO₄²⁻, NO₃⁻, Ca²⁺, K⁺, Mg²⁺, and Na⁺) were evaluated in water samples obtained during rainfall events. Temporal variation of ²H and Cl⁻ measured in the groundwater and stream baseflow prior to rainfall was similar; however, the isotope compositions of the streamflow fluctuated significantly and responded quickly to rainfall events, likely due to an infiltration excess, overland dominated surface runoff during torrential rainfall events. Time source separation using ²H demonstrated greater than 72% on average, the stream composition was event water during torrential rainfall events, with the proportion increasing with rainfall intensity. Solutes concentrations in the stream had loglinear relationships with stream discharge, with an outling anomaly with an example of an intra-rainfall event on Oct. 24, 2015. Stream Cl⁻ behaved nonconservative during rainfall events, temporal variation of Cl⁻ indicated a flush and washout at the onset of small rainfall events, a dilution but still high concentration pattern in high discharge and old water dominated in regression flow period. This study indicates rainfall intensity affects runoff responses in a semi-arid catchment, and the stored water in the thick, loess covered areas was less connected with stream runoff. Solute transport may threaten water quality in the area, requiring further analysis of the performance of the eco-restoration project.     

Beaver dam analogues drive heterogeneous groundwater–surface water interactions

Beaver dam analogues (BDAs) are a cost-effective stream restoration approach that leverages the recognized environmental benefits of natural beaver dams on channel stability and local hydrology. Although natural beaver dams are known to exert considerable influence on the hydrologic conditions of a stream system by mediating geomorphic processes, nutrient cycling, and groundwater–surface water interactions, the impacts of beaver-derived restoration methods on groundwater–surface water exchange are poorly characterized. To address this deficit, we monitored hyporheic exchange fluxes and streambed porewater biogeochemistry across a sequence of BDAs installed along a central Wyoming stream during the summer of 2019. Streambed fluxes were quantified by heat tracing methods and vertical hydraulic gradients. Biogeochemical activity was evaluated using major ion porewater chemistry and principal component analysis. Vertical fluxes of approximately 1.0 m/day were observed around the BDAs, as was the development of spatially heterogeneous zones of nitrate production, groundwater upwelling, and anaerobic reduction. Strong contrasts in hyporheic zone processes were observed across BDAs of differing sizes. This suggests that structures may function with size-dependent behaviour, only altering groundwater–surface water interactions after a threshold hydraulic step height is exceeded. Patterns of hyporheic exchange and biogeochemical cycling around the studied BDAs resemble those around natural beaver dams, suggesting that BDAs may provide comparable benefits to channel complexity and near-stream function over a 1-year period.     

Fully coupled elastoplastic hydro-mechanical analysis of unsaturated porous media using a meshfree method

This paper presents the first application of an advanced meshfree method, ie, the edge-based smoothed point interpolation method (ESPIM), in simulation of the coupled hydro-mechanical behaviour of unsaturated porous media. In the proposed technique, the problem domain is spatially discretised using a triangular background mesh, and the polynomial point interpolation method combined with a simple node selection scheme is adopted for creating nodal shape functions. Smoothing domains are formed on top of the background mesh, and a constant smoothed strain, created by applying the smoothing operation over the smoothing domains, is assigned to each smoothing domain. The deformation and flow models are developed based on the equilibrium equation of the mixture, and linear momentum and mass balance equations of the fluid phases, respectively. The effective stress approach is followed to account for the coupling between the flow and deformation models. Further coupling among the phases is captured through a hysteretic soil water retention model that evolves with changes in void ratio. An advanced elastoplastic constitutive model within the context of the bounding surface plasticity theory is employed for predicting the nonlinear behaviour of soil skeleton. Time discretisation is performed by adopting a three-point discretisation method with growing time steps to avoid temporal instabilities. A modified Newton-Raphson framework is designed for dealing with nonlinearities of the discretised system of equations. The performance of the numerical model is examined through a number of numerical examples. The state-of-the-art computational scheme developed is useful for simulation of geotechnical engineering problems involving unsaturated soils.     

High-frequency monitoring of hydrological and biogeochemical fluxes in forested catchments of southern Chile

The variability of rainfall-dependent streamflow at catchment scale modulates many ecosystem processes in wet temperate forests. Runoff in small mountain catchments is characterized by a quick response to rainfall pulses which affects biogeochemical fluxes to all downstream systems. In wet-temperate climates, water erosion is the most important natural factor driving downstream soil and nutrient losses from upland ecosystems. Most hydrochemical studies have focused on water flux measurements at hourly scales, along with weekly or monthly samples for water chemistry. Here, we assessed how water and element flows from broad-leaved, evergreen forested catchments in southwestern South America, are influenced by different successional stages, quantifying runoff, sediment transport and nutrient fluxes during hourly rainfall events of different intensities. Hydrograph comparisons among different successional stages indicated that forested catchments differed in their responses to high intensity rainfall, with greater runoff in areas covered by secondary forests (SF), compared to old-growth forest cover (OG) and dense scrub vegetation (CH). Further, throughfall water was greatly nutrient enriched for all forest types. Suspended sediment loads varied between successional stages. SF catchments exported 455 kg of sediments per ha, followed by OG with 91 kg/ha and CH with 14 kg/ha, corresponding to 11 rainfall events measured from December 2013 to April 2014. Total nitrogen (TN) and phosphorus (TP) concentrations in stream water also varied with rainfall intensity. In seven rainfall events sampled during the study period, CH catchments exported less nutrients (46 kg/ha TN and 7 kg/ha TP) than SF catchments (718 kg/ha TN and 107 kg/ha TP), while OG catchments exported intermediate sediment loads (201 kg/ha TN and 23 kg/ha TP). Further, we found significant effects of successional stage attributes (vegetation structure and soil physical properties) and catchment morphometry on runoff and sediment concentrations, and greater nutrients retention in OG and CH catchments. We conclude that in these southern hemisphere, broad-leaved evergreen temperate forests, hydrological processes are driven by multiple interacting phenomena, including climate, vegetation, soils, topography, and disturbance history.     

Assessment of alternative adsorption models and global sensitivity analysis to characterize hexavalent chromium loss from soil to surface runoff

We investigate our ability to assess transfer of hexavalent chromium, Cr(VI), from the soil to surface runoff by considering the effect of coupling diverse adsorption models with a two‐layer solute transfer model. Our analyses are grounded on a set of two experiments associated with soils characterized by diverse particle size distributions. Our study is motivated by the observation that Cr(VI) is receiving much attention for the assessment of environmental risks due to its high solubility, mobility, and toxicological significance. Adsorption of Cr(VI) is considered to be at equilibrium in the mixing layer under our experimental conditions. Four adsorption models, that is, the Langmuir, Freundlich, Temkin, and linear models, constitute our set of alternative (competing) mathematical formulations. Experimental results reveal that the soil samples characterized by the finest grain sizes are associated with the highest release of Cr(VI) to runoff. We compare the relative abilities of the four models to interpret experimental results through maximum likelihood model calibration and four model identification criteria (i.e., the Akaike information criteria [AIC and AIC_C] and the Bayesian and Kashyap information criteria). Our study results enable us to rank the tested models on the basis of a set of posterior weights assigned to each of them. A classical variance‐based global sensitivity analysis is then performed to assess the relative importance of the uncertain parameters associated with each of the models considered, within subregions of the parameter space. In this context, the modelling strategy resulting from coupling the Langmuir isotherm with a two‐layer solute transfer model is then evaluated as the most skilful for the overall interpretation of both sets of experiments. Our results document that (a) the depth of the mixing layer is the most influential factor for all models tested, with the exception of the Freundlich isotherm, and (b) the total sensitivity of the adsorption parameters varies in time, with a trend to increase as time progresses for all of the models. These results suggest that adsorption has a significant effect on the uncertainty associated with the release of Cr(VI) from the soil to the surface runoff component.     

Predicting quantiles of water quality from catchment characteristics

Water quality is often highly variable both in space and time, which poses challenges for modelling the more extreme concentrations. This study developed an alternative approach to predicting water quality quantiles at individual locations. We focused on river water quality data that were collected over 25 years, at 102 catchments across the State of Victoria, Australia. We analysed and modelled spatial patterns of the 10th, 25th, 50th, 75th and 90th percentiles of the concentrations of sediments, nutrients and salt, with six common constituents: total suspended solids (TSS), total phosphorus (TP), filterable reactive phosphorus (FRP), total Kjeldahl nitrogen (TKN), nitrate-nitrite (NO_x), and electrical conductivity (EC). To predict the spatial variation of each quantile for each constituent, we developed statistical regression models and exhaustively searched through 50 catchment characteristics to identify the best set of predictors for that quantile. The models predict the spatial variation in individual quantiles of TSS, TKN and EC well (66%–96% spatial variation explained), while those for TP, FRP and NO_x have lower performance (37%–73% spatial variation explained). The most common factors that influence the spatial variations of the different constituents and quantiles are: annual temperature, percentage of cropping land area in catchment and channel slope. The statistical models developed can be used to predict how low- and high-concentration quantiles change with landscape characteristics, and thus provide a useful tool for catchment managers to inform planning and policy making with changing climate and land use conditions.     

A COMPARISON OF EASTERN AND WESTERN HONG KONG PHYTOPLANKTON FROM WEEKLY SAMPLES (1997-1999)

ＩＮＴＲＯＤＵＣＴＩＯＮＴｈｅｍａｉｎｐｕｒｐｏｓｅｏｆｔｈｉｓｓｔｕｄｙｗａｓｔｏｄｅｔｅｒｍｉｎｅｔｈｅｍａｇｎｉｔｕｄｅａｎｄｆｒｅｑｕｅｎｃｙｏｆｈａｒｍｆｕｌｄｉａｔｏｍａｎｄｄｉｎｏｆｌａｇｅｌｌａｔｅｏｕｔｂｒｅａｋｓｉｎＨｏｎｇＫｏｎｇ’ｓｗｅｓｔｅｒｎａｎｄｅａｓｔｅｒｎｗａｔｅｒｓ (ｅ .ｇ.ＬａｍｍａＳｔｒａｉｔｓａｎｄＰｏｒｔＳｈｅｌｔｅｒｒｅｓｐｅｃｔｉｖｅｌｙ)ｆｏｒｔｈｅｐｅｒｉｏｄＪａｎｕａｒｙ 1 997ｔｏＤｅｃｅｍｂｅｒ1 999.Ｔｈｉｓｗａｓｄｏｎｅｉｎｏｒｄｅｒｔｏｄ…     

用双三次样条函数和GPS资料反演现今中国大陆构造形变场

将中国大陆现今构造变动视为一种连续的地壳变形,利用双三次样条函数模拟了近期GPS测定的大陆内部及周边地区412个测站速率,反演大陆地区自洽的构造变动速度场和应变率场.模拟结果显示:印度板块与欧亚板块的碰撞、挤压是构成中国大陆内部岩石层水平形变的主要驱动力.印度板块在东喜马拉雅构造结深深插入青藏高原,造成地壳大规模的缩短和抬升.青藏高原东南部的喜马拉雅带、拉萨和羌塘地块以及青藏高原东南边的川滇地区,内部构造活动强烈,其内部的构造变形包含地壳碎片的冲断、褶皱和侧向逃逸.大陆地壳(或岩石圈)的增厚,尤其是喜马拉雅山脉南北向的快速缩短和青藏高原东西向的缓慢拉张,大约吸收了印欧板块会聚量的85%,西藏中东地区东西向的拉张速率达到了(16±2.0)mm/a,且顺时针方向扭转明显.印度板块相对欧亚板块运动的欧拉极为(29.7°N, 19.3°E, 0.392°/Ma);华南地块相对于欧亚大陆向东(102°±7.4°)南的运动速率是(11±1.54)mm/a,华南块体相对欧亚板块运动的欧拉极为(62.25°N, 126.56°E, 0.141°/Ma);塔里木地块相对较稳定,其西部运动速度高于东部运动速度,作顺时针方向旋转.总体上讲,中国大陆运动方向为北偏东呈辐射状,从西部近南北方向的运动转向东部地区东南方向的运动,绕东喜马拉雅构造结有一顺时针方向的旋转.横穿喜马拉雅构造带及青藏内部的南北向压缩速率为(19±2.0)mm/a,横穿西天山构造带的南北向压缩平均速率为(13±1.5)mm/a,横穿东天山构造带的南北向压缩平均速率为(6.0±1.4)mm/a.阿尔金断裂带的左旋走滑速率为(6±1.2)mm/a.     

低轨道人造卫星(CHAMP、GRACE、GOCE)与高精度地球重力场——卫星重力大地测量的最新发展及其对地球科学的重大影响

评述了卫星重力大地测量的最新发展及其对地球科学的重大影响。为了更好地理解地球内部物理构造与海洋动力学，以及大陆，冰川和海洋的相互作用，改善现有地球重力场模型（包括精度和空间解析度）是非常重要的。IUGG等国际组织对此已经强调了很多年。最近，由德国的GFZ（GeoForschungsZentrum)，美国的NASA（National Aeronautics and Space Adminitration)以及欧洲宇航局ESA（European Space Agency)开发研制了最先进的地球监测技术－SST（Satellite-to-Sateilite Tracking)。其主要特点是利用现有的GPS连续追踪新发射低轨道卫星，并由低轨道卫星对地球重力场作精密观测。已经发射和即将发射的卫星有3颗：GHAMP（Challenging Mini-Satellite Payload for Geophysical Research an Application)已经于2000年发射；GRACE（Gravity Recovery and Climate Experimert)定于2002年发射；GOCE（Gravity Field and Steady-state Ocean Cirulation Explorer)计划2004年发射，它们可以统称为重力卫星。载有SST技术的人造卫星的主要目的是获得具有前所未有的高精度和高空间解析度的全球重力场和大地水准面模型，加强人们对地球内部构造的理解并为海洋和气象研究提供更好地参考。上述3个重力卫星工作在有明显区别的不同波谱内，它们有不同的科学应用，仅有一小部分重合。所以，就应用而言它们是完全互补的。它们在地球科学中的应用将是广泛的，特别对于固体地球物理学，海洋学以及大地测量学等领域，它们将会带来革命性的变化，其意义不亚于GPS。     

THE POSITION AND ROLES OF PORTS AND RAILWAYS IN THE EURASIA LAND BRIDGE TRANSPORTATION IN NORTHEAST CHINA

This paper systematically analyzes and proves the favorable factors of utilizing the ports and the railways in Northeast China for the Asia-Europe land bridge transportation. It will be more beneficial to Dalian Port and HaDa (Harbin-Dalian) and Binzhou (Harbin-Manzhouli) railways in Northeast China to the Siberia railway in Russia than 1) to other ports and the corresponding railways in China, 2) to the Nakhodka Port in Russia and the corresponding railways and 3) to the Chongjin Port in North Korea and the corresponding railways. This paper also puts forward the reform measures to adopt the ports and the railways in the northeast region in China for the land bridge transportation and the problems in transportation policies and management systems.