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.     

Isotope geochemistry and Re–Os geochronology of the Yanjiagou Mo deposit in the central North China Craton

The Yanjiagou deposit, located in the central North China Craton (NCC), is a newly found porphyry‐type Mo deposit. The Mo mineralization here is spatially associated with the Mapeng batholith. In this study, we identify four stages of ore formation in this deposit: pyrite phyllic stage (I), quartz–pyrite stage (II), quartz–pyrite–molybdenite stage (III), which is the main mineralization stage, and quartz–carbonate stage (IV). We present sulphur and lead isotope data on pyrite, and rhenium and osmium isotopes of molybdenite from the porphyry deposit and evaluate the timing and origin of ore formation. The δ³⁴S values of the pyrite range from ‐1.1‰ to −0.6‰, with an average of −0.875‰, suggesting origin from a mixture of magmatic/mantle sources and the basement rocks. The Pb isotope compositions of the pyrite show a range of 16.369 to 17.079 for ²⁰⁶Pb/²⁰⁴Pb, 15.201 to 15.355 for ²⁰⁷Pb/²⁰⁴Pb, and 36.696 to 37.380 for ²⁰⁸Pb/²⁰⁴Pb, indicating that the ore‐forming materials were derived from a mixture of lower crust (or basement rocks) and mantle. Rhenium contents in molybdenite samples from the main ore stage are between 74.73 to 254.43 ppm, with an average of 147.9 ppm, indicating a mixed crustal‐mantle source for the metal. Eight molybdenite separates yield model ages ranging from 124.17 to 130.80 Ma and a mean model age of 128.46 Ma. An isochron age of 126.7 ± 1.1 Ma (MSWD = 2.1, initial ¹⁸⁷Os = 0.0032 ± 0.0012 ppb) is computed, which reveals a close link between the Mo mineralization and the magmatism that generated the Mapeng batholith. The age is close to the zircon U–Pb age of ca. 130 Ma from the batholith reported in a recent study. The age is also consistent with the timing of mineralization in the Fuping ore cluster in the central NCC, as well as the peak time of lithosphere thinning and destruction of the NCC. We evaluate the spatio‐temporal distribution of the Mo deposits in the NCC and identify three important molybdenum provinces along the northern and southern margins of the craton formed during three distinct episodes: Middle to Late Triassic (240–220 Ma), Early Jurassic (190–175 Ma), and Late Jurassic to Early Cretaceous (150–125 Ma). The third period is considered to mark the most important metallogenic event, coinciding with the peak of lithosphere thinning and craton destruction in the NCC. Copyright © 2014 John Wiley & Sons, Ltd.