Temporal variation of soil erosion resistance on sloping farmland during the growth stages of maize (Zea mays L.)

Maize growth has great effects on soil properties and thus likely induces the changes in soil erosion resistance on sloping farmland. However, temporal variation of soil erosion resistance during the growth stages of maize is still unclear in the mountainous yellow soil area where maize is the dominant crop. In this study, four maize plots (MP) and four bare land plots (CK) were conducted to investigate soil erosion resistance, and multiple indicators of soil erosion resistance were measured including the total soil anti-scourability (TAS), mean weight diameter (MWD), soil erodibility K factor and soil shear strength (SH). A comprehensive soil erosion resistance index (CSERI) was employed to quantify the temporal variation of soil erosion resistance during the growth stages of maize (seedling stage, SS; jointing stage, JS; tasselling stage, TS; maturing stage, MS). The results showed that TAS, MWD, SH increased significantly with maize growth and SH decreased when at MS. But K factor decreased significantly over time. CSERI increased significantly during the growth stages of maize and the CSERI of JS, TS, MS increased on average by 74.72, 180.68 and 234.57% than that of SS. Compared to CK, CSERI of MP increased by 49.90, 66.82, 55.60 and 38.61% during the growth stages of maize. The temporal variation of soil erosion resistance was closely related to the changes in maize cover, maize roots and soil organic carbon. The findings demonstrated that it is necessary to consider the temporal variation of soil erosion resistance in the mountainous yellow soil area.     

Watershed delineation using hydrographic features and a DEM in plain river network region

Watershed delineation is a required step when conducting any spatially distributed hydrological modelling. Automated approaches are often proposed to delineate a watershed based on a river network extracted from the digital elevation model (DEM) using the deterministic eight‐neighbour (D8) method. However, a realistic river network cannot be derived from conventional DEM processing methods for a large flat area with a complex network of rivers, lakes, reservoirs, and polders, referred to as a plain river network region (PRNR). In this study, a new approach, which uses both hydrographic features and DEM, has been developed to address the problems of watershed delineation in PRNR. It extracts the river nodes and determines the flow directions of the river network based on a vector‐based hydrographic feature data model. The river network, lakes, reservoirs, and polders are then used to modify the flow directions of grid cells determined by D8 approach. The watershed is eventually delineated into four types of catchments including lakes, reservoirs, polders, and overland catchments based on the flow direction matrix and the location of river nodes. Multiple flow directions of grid cells are represented using a multi‐direction encoding method, and multiple outflows of catchments are also reflected in the topology of catchments. The proposed approach is applied to the western Taihu watershed in China. Comparisons between the results obtained from the D8 approach, the ‘stream burning’ approach, and those from the proposed approach clearly demonstrate an improvement of the new approach over the conventional approaches. This approach will benefit the development of distributed hydrological models in PRNR for the consideration of different types and multiple inlets and outlets of catchments. Copyright © 2015 John Wiley & Sons, Ltd.     

Carboniferous eclogite and garnet–omphacite granulite from northeastern Hainan Island,South China: Implications for the evolution of the eastern Palaeo‐Tethys

High‐P (HP) eclogite and associated garnet–omphacite granulite have recently been discovered in the Mulantou area, northeastern Hainan Island, South China. These rocks consist mainly of garnet, omphacite, hornblende, quartz and rutile/ilmenite, with or without zoisite and plagioclase. Textural relationships, mineral compositions and thermobarometric calculations demonstrate that the eclogite and garnet–omphacite granulite share the same three‐stage metamorphic evolution, with prograde, peak and retrograde P?T conditions of 620–680°C and 8.7–11.1 kbar, 820–860°C and 17.0–18.2 kbar, and 700–730°C and 7.1–8.5 kbar respectively. Sensitive high‐resolution ion microprobe U–Pb zircon dating, coupled with the identification of mineral inclusions in zircon, reveals the formation of mafic protoliths before 355 Ma, prograde metamorphism at c. 340–330 Ma, peak to retrograde metamorphism at c. 310–300 Ma, and subsequent pegmatite intrusion at 295 Ma. Trace element geochemistry shows that most of the rocks have a MORB affinity, with initial ε_Nd values of +2.4 to +6.7. As with similar transitional eclogite–HP granulite facies rocks in the thickened root in the European Variscan orogen, the occurrence of relatively high P?T metamorphic rocks of oceanic origin in northeastern Hainan Island suggests Carboniferous oceanic subduction leading to collision of the Hainan continental block, or at least part of it, with the South China Block in the eastern Palaeo‐Tethyan tectonic domain.     

Generation of continuous surface soil moisture dataset using combined optical and thermal infrared images

Surface soil moisture (SSM) is a critical variable for understanding water and energy flux between the atmosphere and the Earth's surface. An easy to apply algorithm for deriving SSM time series that primarily uses temporal parameters derived from simulated and in situ datasets has recently been reported. This algorithm must be assessed for different biophysical and atmospheric conditions by using actual geostationary satellite images. In this study, two currently available coarse‐scale SSM datasets (microwave and reanalysis product) and aggregated in situ SSM measurements were implemented to calibrate the time‐invariable coefficients of the SSM retrieval algorithm for conditions in which conventional observations are rare. These coefficients were subsequently used to obtain SSM time series directly from Meteosat Second Generation (MSG) images over the study area of a well‐organized soil moisture network named REMEDHUS in Spain. The results show a high degree of consistency between the estimated and actual SSM time series values when using the three SSM dataset‐calibrated time‐invariable coefficients to retrieve SSM, with coefficients of determination (R²) varying from 0.304 to 0.534 and root mean square errors ranging from 0.020 m³/m³ to 0.029 m³/m³. Further evaluation with different land use types results in acceptable debiased root mean square errors between 0.021 m³/m³ and 0.048 m³/m³ when comparing the estimated MSG pixel‐scale SSM with in situ measurements. These results indicate that the investigated method is practical for deriving time‐invariable coefficients when using publicly accessed coarse‐scale SSM datasets, which is beneficial for generating continuous SSM dataset at the MSG pixel scale.     

Revised runoff curve number for runoff prediction in the Loess Plateau of China

The soil conservation service (now Natural Resources Conservation Service) Curve Number (SCS-CN), one of the most commonly used methods for surface runoff prediction. The runoff calculated by this method was very sensitive to CN values. In this study, CN values were calculated by both arithmetic mean (CN_C) and least square fit method (CN_F) using observed rainfall-runoff data from 43 sites in the Loess Plateau region, which are considerably different from the CN₂ values obtained from the USDA-SCS handbook table (CN_T). The results showed that using CN_C instead of CN_T for each watershed produce little improvement, while replacing CN_T with CN_F improves the performance of the original SCS-CN method, but still performs poorly in most study sites. This is mainly due to the SCS-CN method using a constant CN value and discounting of the temporal variation in rainfall-runoff process. Therefore, three factors—soil moisture, rainfall depth and intensity—affecting the surface runoff variability are considered to reflect the variation of CN in each watershed, and a new CN value was developed. The reliability of the proposed method was tested with data from 38 watersheds, and then applied to the remaining five typical watersheds using the optimized parameters. The results indicated that the proposed method, which boosted the model efficiencies to 81.83% and 74.23% during calibration and validation cases, respectively, performed better than the original SCS-CN and the Shi and Wang (2020b) method, a modified SCS-CN method based on tabulated CN value. Thus, the proposed method incorporating the influence of the temporal variability of soil moisture, rainfall depth, and intensity factors suggests an accurate runoff prediction for general applications under different hydrological and climatic conditions on the Loess Plateau region.     

An analytical study on three‐dimensional versus two‐dimensional water table‐induced flow patterns in a Tóthian basin

Although it has been increasingly acknowledged that groundwater flow pattern is complicated in the three‐dimensional (3‐D) domain, two‐dimensional (2‐D) water table‐induced flow models are still widely used to delineate basin‐scale groundwater circulation. However, the validity of 2‐D cross‐sectional flow field induced by water table has been seldom examined. Here, we derive the analytical solution of 3‐D water table‐induced hydraulic head in a Tóthian basin and then examine the validity of 2‐D cross‐sectional models by comparing the flow fields of selected cross sections calculated by the 2‐D cross‐sectional model with those by the 3‐D model, which represents the “true” cases. For cross sections in the recharge or discharge area of the 3‐D basin, even if head difference is not significant, the 2‐D cross‐sectional models result in flow patterns absolutely different from the true ones. For the cross section following the principal direction of groundwater flow, although 2‐D cross‐sectional models would overestimate the penetrating depth of local flow systems and underestimate the recharge/discharge flux, the flow pattern from the cross‐sectional model is similar to the true one and could be close enough to the true one by adjusting the decay exponent and anisotropy ratio of permeability. Consequently, to determine whether a 2‐D cross‐sectional model is applicable, a comparison of hydraulic head difference between 2‐D and 3‐D solutions is not enough. Instead, the similarity of flow pattern should be considered to determine whether a cross‐sectional model is applicable. This study improves understanding of groundwater flow induced by more natural water table undulations in the 3‐D domain and the limitations of 2‐D models accounting for cross‐sectional water table undulation only.     

Steady‐state analytical solutions of flow and deformation coupling due to a point sink within a finite fluid‐saturated poroelastic layer

An exact steady‐state closed‐form solution is presented for coupled flow and deformation of an axisymmetric isotropic homogeneous fluid‐saturated poroelastic layer with a finite radius due to a point sink. The hydromechanical behavior of the poroelastic layer is governed by Biot's consolidation theory. Boundary conditions on the lateral surface are specifically chosen to match the appropriate finite Hankel transforms and simplify the transforms of the governing equations. Ordinary differential equations in the transformed domain are solved, and then the analytical solutions in the physical space for the pore pressure and the displacements are finally obtained by using finite Hankel inversions. The analytical solutions at some special locations such as the top and bottom surfaces, lateral surface, and the symmetrical axis are given and analyzed. And a case study for the consolidation of a water‐saturated soft clay layer due to pumping is conducted. The analytical solution is verified against the finite element solution. Meanwhile, an analysis of coupled hydromechanical behavior is carried out herein. The presented analytical solution is an exact solution to the practical poroelastic problem within an axisymmetric finite layer. It can provide us a better understanding of the poroelastic behavior of the finite layer due to fluid extraction. Besides, it can be applied to calibrate numerical schemes of axisymmetric poroelasticity within finite domains. Copyright © 2017 John Wiley & Sons, Ltd.     

Theoretical and experimental study of breaking wave on sloping bottoms

This paper studies the continuous evolution of breaking wave for the surface water waves propagating on a sloping beach. A Lagrangian asymptotic solution is derived. According to the solution coupled with the wave breaking criteria and the equations of water particles motion, the wave deformation and the continuous wave breaking processes for the progressive water waves propagating on a sloping bottom can be derived. A series of experiments are also conducted to compare with the theoretical solution. The results show that the present solution can reasonably describe the plunging or spilling wave breaking phenomenon.     

Multiple overpressure systems and their relationship with hydrocarbon accumulation in the Qikou Depression of the Bohai Bay Basin,China

Pressure measurements using drill stem tests and estimates from log data calculation indicate that three vertically stacked regional pressure compartments exist in the Qikou Depression of Bohai Bay Basin, N. China. The compartments comprise hydrostatic, upper weak, and lower overpressure systems. Laterally, overpressure (pressure coefficient > 1.2) occurs in the deeper areas and weakens gradually from the centre to the margin of the depression. The accumulation of oil and gas exhibits the interesting characteristics of oil‐bearing layers above gas‐bearing layers in the Qikou Depression. The pattern can be accounted for by the evolution of overpressure system, the maturity process of the source rock and the main fault activity. In the late Dongying Formation (Ed, 30 Ma), the lower overpressure system began to form shape, and the hydrocarbon sources generated a large volume of oil. However, because there was no migration pathway, the oil only accumulated in the original strata. In the late Guantao Formation (Ng, 12 Ma), the gas was generated, the upper overpressure system formed gradually, and the activity of the main fault gradually increased. Then, the overpressure pushed the early gathered oil to flow from the lower overpressure system into the upper overpressure system. Afterwards, the activity of the main fault decreased again and remains weak until now. Thus, later generated natural gas cannot keep migrating along the main fault and can only accumulate in the lower overpressure system. Copyright © 2014 John Wiley & Sons, Ltd.