Parameterization of the Ångström–Prescott formula based on machine learning benefit estimation of reference crop evapotranspiration with missing solar radiation data

Accurately estimated reference evapotranspiration (ET₀) is essential to regional water management. The FAO recommends coupling the Penman–Monteith (P-M) model with the Ångström–Prescott (A-P) formula as the standard method for ET₀ estimation with missing R_s measurements. However, its application is usually restricted by the two fundamental coefficients (a and b) of the A-P formula. This paper proposes a new method for estimating ET₀ with missing R_s by combining machine learning with physical-based P-M models (PM-ET₀). The benchmark values of the A-P coefficients were first determined at the daily, monthly, and yearly scales, and further evaluated in R_s and ET₀ estimates at 80 national R_s measuring stations. Then, three empirical models and four machine-learning methods were evaluated in estimating the A-P coefficients. Machine learning methods were also used to estimate ET₀ (ML-ET₀) to compare with the PM-ET₀. Finally, the optimal estimation method was used to estimate the A-P coefficients for the 839 regular weather stations for ET₀ estimation without R_s measurement for China. The results demonstrated a descending trend for coefficient a from northwest to southeast China, with larger values in cold seasons. However, coefficient b showed the opposite distribution as the coefficient a. The FAO has recommended a larger a but a smaller b for southeast China, which produced the region's largest R_s and ET₀ estimation errors. Additionally, the A-P coefficients calibrated at the daily scale obtained the best estimation accuracy for both R_s and ET₀, and slightly outperformed the monthly and yearly coefficients without significant difference in most cases. The machine learning methods outperformed the empirical methods for estimating the A-P coefficients, especially for the sites with extreme values. Further, ML-ET₀ outperformed the PM-ET₀ with yearly A-P coefficients but underperformed those with daily and monthly ones. This study indicates an exciting potential for combining machine learning with physical models for estimating ET₀. However, we found that using the A-P coefficients with finer time scales is unnecessary to deal with the missing R_s measurements.     

Effect of seasonal freeze–thaw process on spatial and temporal distribution of soil water and its infiltration to recharge groundwater

Clarifying the distribution and dynamics of soil moisture during the freeze–thaw process is crucial for surface ecology and is an objective requirement to investigate the mechanism of changes during the groundwater recharge process in a freeze–thaw zone. Based on the monitoring data of soil moisture and temperature in the Changbai Mountain area, the freeze–thaw process is classified into four periods. This study investigates the hydrothermal migration processes during different periods. The simultaneous heat and water model is used to simulate and analyse the infiltration of soil moisture into groundwater under five precipitation insurance rates. The results are as follows: (1) The smaller the soil depth, the stronger is the correlation between soil temperature and air temperature during the freeze–thaw process. (2) The redistribution of soil moisture before and after freeze–thaw is significantly affected by the soil texture, and soil permeability affects the recharge of soil moisture from the upper region to the lower region during the thawing period. (3) Groundwater receives vertical infiltration recharge mainly during non-freezing and is supplied by freezing and snowmelt recharge during the stable thawing period. The percentage of soil water infiltration during the stable thawing period in the total annual infiltration increases gradually with the precipitation insurance rate.     

‘I hate sand … it gets everywhere’—Phanerozoic sedimentary recycling from NW Africa

Over 10,000 published detrital zircon ages have been reprocessed (applying a +10% normal and reverse concordance range) and analysed to understand the evolution of the detrital zircon record of north-west Africa during the Phanerozoic. Using dissimilarity and clustering analysis, shifts in detrital zircon populations allow interpretation of the evolution of source regions and source to sink systems throughout the Phanerozoic within the West Gondwana superfan. Previous thermochronology and field studies conducted across north-west Africa indicate significant and sustained shifts in source regions in Meso-Cenozoic times which are not recorded in the detrital zircon geochronology record. This discrepancy is most notable for Mesozoic to modern source to sink studies focused on the evolution of the Atlasic rift and opening of the Atlantic and Tethyan Oceans to the west and north respectively. Our results indicate a high degree of similarity between samples from Cambrian times onwards due to successive phases of sediment recycling. This highlights the need to integrate detrital zircon analysis with other techniques to provide confident reconstruction of sediment routing systems across Morocco. This systematic review also reveals the ubiquitous occurrence of Mesoproterozoic zircons within Moroccan sediment. No basement of this age is known from north-west Africa—often described as the ‘Mesoproterozoic Gap’, which was thought to be a diagnostic feature of sediment derived from the West African Craton. However, zircons of this age form 7% of all analysed zircons and are present in sediments from at least 700 Ma. The presence of this population is interpreted as strongly diagnostic of provenance from either the Amazonian Craton or the Eastern Gondwana Orogen within Central Africa. Their presence in the Moroccan detrital record from the Neoproterozoic onwards raises questions about the position of the West African Craton in the Proterozoic, and for the spatial extent of Mesoproterozoic orogeny within north Africa.     

Mesozoic–Cenozoic multistage tectonic evolution of the Pamir: Detrital fission-track constraints from the Tajik Basin

Knowledge of the tectonic history of the Pamir contributes to our understanding of both the evolution of collisional orogenic belts as well as factors controlling Central Asian aridification. It is, however, not easy to decipher the Mesozoic–Cenozoic tectonics of the Pamir due to extensive Neogene deformation in an orogen that remains largely understudied. This study reports detrital apatite and zircon fission-track (FT) ages from both the eastern Tajik Basin sedimentary rocks and Pamir modern river sands. These FT data, supported by vitrinite reflectance and zircon and apatite U–Pb double dating, suggest that the majority of the FT ages are unreset and record exhumation stages of the Pamir, which has served as the source terrane of the Tajik Basin since the Cretaceous. Furthermore, we combine the new data with a compilation of published detrital apatite and zircon FT data from both the Tajik Basin sedimentary rocks and Pamir modern river sands, to explore the Mesozoic–Cenozoic tectonic history of Pamir. Deconvolved FT Peak Ages document two major Mesozoic exhumation events associated with the Late Triassic–Early Jurassic Cimmerian orogeny that reflects accretion of the Pamir terranes, as well as the Early–early Late Cretaceous deformation associated with the northward subduction of the Neo-Tethys Ocean beneath Pamir. The compiled data also show significant Late Eocene–Neogene exhumation associated with the ongoing formation of the Pamir, which peaks at ca. 36, 25, 14 and 7 Ma.     

新疆喀喇昆仑大红柳滩一带锂矿光谱特征 及其找矿指示意义

近年来,新疆喀喇昆仑大红柳滩地区锂铍矿取得重大找矿发现,显示良好的成矿潜力。受其高寒深切割环境的制约,矿区外围的地质找矿工作推进不易。前人利用高分、多光谱遥感数据及高光谱技术在该区花岗伟晶岩型锂铍矿找矿工作中发挥出积极作用,发现了505、507、俘虏沟南1号、俘虏沟南2号等多处大型锂铍矿产地。笔者通过对区内花岗伟晶岩型锂矿床中岩、矿石开展地面光谱测试,系统对比不同岩性赋矿地层、岩体、含矿伟晶岩、不含矿伟晶岩的光谱特征。通过光谱解算剖析其矿物组合,分析显示含矿伟晶岩中含锂矿物以锂辉石、锂云母为主,少量锂绿泥石。含矿伟晶岩具高反射率,发育强Al-OH吸收峰,受其接触交代蚀变影响,发育Mg-OH吸收峰,二者共同构成区内含矿伟晶岩的标志性遥感异常组合。采用国产高光谱遥感卫星数据,利用匹配滤波方法提取上述异常组合,圈定找矿有利区7处,为该区下一步锂铍找矿工作提供了依据。     

Metamorphic stages in mountain belts during a Wilson cycle: A case study in the central Sanandaj–Sirjan zone (Zagros Mountains,Iran)

Polymetamorphic units are important constituents of continent–continent collisional orogens, and rift metamorphic assemblages are often overprinted by subsequent metamorphism during subduction and collision. This study reports the metamorphic conditions and evolution of the Dorud–Azna metamorphic units in the central part of the Sanandaj–Sirjan zone (SSZ), Iran. Here, new geothermobarometry results are integrated with ⁴⁰Ar/³⁹Ar mineral and Th–U–Pb monazite and thorite ages to provide new insight of polyphase metamorphism in the two different basement units of the SSZ, the lower Galeh-Doz orthogneiss and higher Amphibolite-Metagabbro units. In the Amphibolite-Metagabbro unit, staurolite micaschist underwent a prograde P–T evolution from 640 ± 20 °C/6.2 ± 0.8 kbar in garnet cores (M1) to 680 ± 20 °C/7.2 ± 1.0 kbar in garnet rims (M2). Three Th–U–Pb monazite ages of 306 ± 5 Ma, 322 ± 28 Ma and 336 ± 39 Ma from the garnet-micaschists testify the Carboniferous age of M1 metamorphism. In the same unit, the metagabbro records P–T conditions of 4.0 ± 0.8 kbar and 580 ± 50 °C in the (magmatic) amphibole core (Late Carboniferous intrusion) to 7.5 ± 0.7 kbar and 700 ± 20 °C in the amphibole rim indicating a prograde P–T path during subsequent burial (M1). New ⁴⁰Ar/³⁹Ar dating of white mica from the staurolite micaschist yielded a staircase pattern ranging from 36 ± 12 Ma to 170 ± 2 Ma. This implies polymetamorphism with a minimum Late Jurassic cooling age through the Ar retention temperature of ca. 425 ± 25 °C after M2 metamorphism and a Paleogene low-grade metamorphic overprint (M3), while ⁴⁰Ar/³⁹Ar white mica dating of garnet micaschist yielded a plateau age of 137.84 ± 0.65 Ma. We therefore interpret the amphibolite-grade metamorphism M2 to have predated 170 Ma and is likely between 180 and 200 Ma. Furthermore, it is overprinted at about 36 Ma under retrogressive low-grade M3 metamorphism (at temperatures of ~350–240 °C) during final shortening and exhumation. In the underlying Galeh-Doz unit, the Panafrican granitic orthogneiss intruded at P–T conditions of 3.2 ± 4 kbar and 700 ± 20 °C, then it was metamorphosed and deformed at 600 ± 50 °C and 2.0 ± 0.8 kbar (metamorphic stage M1) prior to Late Carboniferous intrusion of mafic dikes. ⁴⁰Ar/³⁹Ar dating of amphibole from the Galeh-Doz orthogneiss gave plateau-like steps between 260 and 270 Ma, representing the age of cooling through ca. 500 °C after the M1 metamorphic event. Interestingly, the results of this study demonstrate polyphase metamorphic histories in both the Galeh-Doz orthogneiss and Amphibolite-Metagabbro units at different P–T conditions and final thick-skinned Paleogene emplacement of these units over the underlying low-grade metamorphic June Complex. Our findings suggest that both units are affected by high-T/low-P Late Carboniferous orogenic metamorphism along with the bimodal magmatism, as result of rifting. We propose that the Early Jurassic amphibolite-grade M2 metamorphism of the SSZ is correlated with the initial subduction of the Neotethyan Ocean. Eventually, the investigated units reflect various stages of a Wilson cycle, from rifting to initiation of the subduction in final plate collision.     

Geochronology and geochemistry of Cretaceous–Eocene granites,Tengchong Block (SW China): Petrogenesis and implications for Mesozoic-Cenozoic tectonic evolution of Eastern Tethys

The Early Cretaceous–Early Eocene granitoids in the Tengchong Block record the evolutionary history of the Mesozoic-Cenozoic tectono-magmatic evolution of Eastern Tethys. (a) The Early Cretaceous granitoids with relatively low (⁸⁷Sr/⁸⁶Sr)_i ratios of 0.7090–0.7169 and ε_Nd(t) values of ?9.8 to ?7.8 display metaluminous, calc-alkaline dominated by I-type granite affinity and hybrid mantle–crust geochemical signatures. They may have been derived from melting of the subducted Meso-Tethyan Bangong-Nujiang oceanic crust with terrigenous sediments in an arc-continent collisional setting. (b) The Late Cretaceous–Paleocene granitoids with relatively high (⁸⁷Sr/⁸⁶Sr)_i ratios of 0.7109–0.7627, and ε_Nd(t) values of ?12.1 to ?7.9 exhibit metaluminous to peraluminous, calc-alkaline dominated by S-type granite affinity and hybrid Lower–Upper crust geochemical signatures, which may be originated from partial melting of the Meso-Proterozoic continental crust in the collision setting between the Tengchong Block and Baoshan Block. (c) The Early Eocene granitoids have metaluminous, calc-alkaline I-type and S-type granites dual affinity, with relatively high (⁸⁷Sr/⁸⁶Sr)_i ratios of 0.711–0.736, ε_Nd(t) values of ?9.4 to ?4.7, showing crust-mantle mixing geochemical signatures. They may have been originated from partial melting of the late Meso-Proterozoic upper crustal components mixed with some upper mantle material during the ascent process of mantle magma caused by the subduction of the Neo-Tethyan Putao–Myitkyian oceanic crust, and collision between the Western Burma Block and the Tengchong Block. It is these multi-stage subductions and collisions that caused the spatial and temporal distribution of the granitic rocks in the Tengchong Block.     

Crust–mantle interaction recorded by Early Cretaceous volcanic rocks within the southern margin of the North China Craton

This study presents new geochronological and geochemical data for Early Cretaceous volcanic rocks in the southern margin of the North China Craton (NCC), to discuss the crust–mantle interaction. The studied rocks include pyroxene andesites from Daying Formation, hornblende andesites and andesites from Jiudian Formation, and rhyolites from a hitherto unnamed Formation. These rocks formed in Early Cretaceous (138–120 Ma), with enrichment in light rare earth elements (REE), depletion in heavy REE and arc-like trace elements characteristics. Pyroxene andesites show low SiO₂ contents and enriched Sr–Nd–Pb–Hf isotopic compositions, with orthopyroxene phenocryst and Paleoproterozoic (2320–1829 Ma) inherited zircons, suggesting that they originated from lithospheric mantle after metasomatism with NCC lower crustal materials. Hornblende andesites have low SiO₂ contents and high Mg# (Mg# = 100 Mg/(Mg + Fe²⁺)) values, indicating a lithospheric-mantle origin. Considering the distinct whole-rock Sr isotopic compositions we divide them into two groups. Among them, the low (⁸⁷Sr/⁸⁶Sr)_i andesites possess amount inherited Neoarchean to Neoproterozoic (2548–845 Ma) zircons, indicating the origin of lithospheric mantle with addition of Yangtze Craton (YC) and NCC materials. In comparison, the high (⁸⁷Sr/⁸⁶Sr)_i andesites, with abundant Neoarchean–Paleozoic inherited zircons (3499–261 Ma), are formed by partial melting of lithospheric mantle with incorporation of NCC supracrustal rocks and YC materials. Rhyolites have extremely high SiO₂ (77.63–82.52 wt.%) and low total Fe₂O₃, Cr, Ni contents and Mg# values, combined with ancient inherited zircon ages (2316 and 2251 Ma), suggesting an origin of NCC lower continental crust. Considering the presence of resorption texture of quartz phenocryst, we propose a petrogenetic model of ‘crystal mushes’ for rhyolites prior to their eruption. These constraints record the intense crust–mantle interaction in the southern margin of the NCC. Given the regional data and spatial distribution of Early Cretaceous rocks within NCC, we believe that the formation of these rocks is related to the contemporaneous far-field effect of the Paleo-Pacific Plate.