Alkali-calcic and alkaline post-orogenic (PO) granite magmatism: petrologic constraints and geodynamic settings

The end of an orogenic Wilson cycle corresponds to amalgamation of terranes into a Pangaea and is marked by widespread magmatism dominated by granitoids. The post-collision event starts with magmatic processes still influenced by subducted crustal materials. The dominantly calc-alkaline suites show a shift from normal to high-K to very high-K associations. Source regions are composed of depleted and later enriched orogenic subcontinental lithospheric mantle, affected by dehydration melting and generating more and more K- and LILE-rich magmas. In the vicinity of intra-crustal magma chambers, anatexis by incongruent melting of hydrous minerals may generate peraluminous granitoids bearing mafic enclaves. The post-collision event ends with emplacement of bimodal post-orogenic (PO) suites along transcurrent fault zones. Two suites are defined, (i) the alkali-calcic monzonite–monzogranite–syenogranite–alkali feldspar granite association characterised by [biotite+plagioclase] fractionation and moderate [LILE+HFSE] enrichments and (ii) the alkaline monzonite–syenite–alkali feldspar granite association characterised by [amphibole+alkali feldspar] fractionation and displaying two evolutionary trends, one peralkaline with sodic mafic mineralogy and higher enrichments in HFSE than in LILE, and the other aluminous biotite-bearing marked by HFSE depletion relative to LILE due to accessory mineral precipitation. Alkali-calcic and alkaline suites differ essentially in the amounts of water present within intra-crustal magma chambers, promoting crystallisation of various mineral assemblages. The ultimate enriched and not depleted mantle source is identical for the two PO suites. The more primitive LILE and HFSE-rich source rapidly replaces the older orogenic mantle source during lithosphere delamination and becomes progressively the thermal boundary layer of the new lithosphere. Present rock compositions are a mixture of major mantle contribution and various crustal components carried by F-rich aqueous fluids circulating within convective cells created around magma chambers. In favourable areas, PO suites pre-date a new orogenic Wilson cycle.     

A semi-empirical wind set-up forecasting model for Lake Champlain

The precision of Lake Champlain's water level estimation is a key component in the flood forecasting process for the Richelieu River. Hydrological models do not typically take into consideration the effects of the wind on the water level (also known as the wind set-up). The objective of this study is to create an empirical wind set-up forecast model for Lake Champlain during high wind events. The proposed model uses wind speed and direction across the Lake, as well as wind gusts as inputs. The model is calibrated to a subset of observations and evaluated on an independent sample, considering four wind speed bins. It is tested and compared to a variant of the Zuider Zee equation on 20 wind set-up events that occurred between 2017 and 2019 using hindcast data from five different numerical weather prediction systems (GDPS, RDPS, HRDPS, NOAA and ECMWF). A quantile mapping-based forecast calibration scheme is implemented for each of the forecast products to correct their biases. Results show that events are successfully predicted by the proposed model at least 72 h in advance. These results are better than the other comparative models found in the literature and tested herein. Overall, significant improvements are obtained by including wind speed and wind gusts from different weather stations.     

Oxygen isotope disequilibrium during serpentinite dehydration

Talc + olivine in metaperidotites result from the serpentinite breakdown due to increasing temperature in the Bergell contact aureole. Jack‐straw olivine textures are present in close proximity to the serpentine breakdown reaction. As the intrusion is approached, the number of olivine crystals increases while the size of the crystals decreases; this feature documents increased overstepping with increased heating rates. Talc veinlets are observed in the outer parts of the talc–olivine zone and are interpreted to be pathways of fluid produced during devolatilization of serpentine. None of the talc–olivine oxygen isotope pairs analysed are in isotopic equilibrium with respect to the peak contact temperature. This implies that escaping fluids cannot be in equilibrium with both phases. Hence, the fluid produced by serpentine reaction does not directly reflect the protolith composition, and attention must be given to the reaction mechanism before interpreting fluid isotope composition.     

The Multiscale TROPIcal CatchmentS critical zone observatory M-TROPICS dataset II: Land use,hydrology and sediment production monitoring in Houay Pano,northern Lao PDR

Mountain regions of the humid tropics are characterized by steep slopes and heavy rains. These regions are thus prone to both high surface runoff and soil erosion. In Southeast Asia, uplands are also subject to rapid land-use change, predominantly as a result of increased population pressure and market forces. Since 1998, the Houay Pano site, located in northern Lao PDR (19.85°N 102.17°E) within the Mekong basin, aims at assessing the long-term impact of the conversion of traditional slash-and-burn cultivation systems to commercial perennial monocultures such as teak tree plantations, on the catchment hydrological response and sediment yield. The instrumented site monitors hydro-meteorological and soil loss parameters at both microplot (1 m²) and small catchment (0.6 km²) scales. The monitored catchment is part of the network of critical zone observatories named Multiscale TROPIcal CatchmentS (M-TROPICS). The data shared by M-TROPICS in Houay Pano are (1) rainfall, (2) air temperature, air relative humidity, wind speed, and global radiation, (3) catchment land use, (4) stream water level, suspended particulate matter, bed particulate matter and stones, (5) soil surface features, and (6) soil surface runoff and soil detachment. The dataset has already been used to interpret suspended particulate matter and bed particulate matter sources and dynamics, to assess the impact of land-use change on catchment hydrology, soil erosion, and sediment yields, to understand bacteria fate and weed seed transport across the catchment, and to build catchment-scale models focused on hydrology and water quality issues. The dataset may be further used to, for example, assess the role of headwater catchments in large tropical river basin hydrology, support the interpretation of new variables measured in the catchment (e.g., contaminants other than faecal bacteria), and assess the relative impacts of both climate and land-use change on the catchment.     

Contribution of magmatism,partial melting buffering and localized crustal thinning on the late Variscan thermal structure of the Agly massif (French Pyrenees)

Low pressure-high temperature (LPHT) metamorphism, with geothermal gradients in the order of 50–100°C/km, is a common feature of the late evolution of collisional orogens. These abnormal thermal conditions may be the results of complex interactions between magmatism, metamorphism and deformation. The Agly massif, in the French Pyrenees, preserves the metamorphic footprints of the late Variscan thermal structure of an almost continuous section from the upper and middle continental crust. The upper crust is characterized by a very high geothermal gradient of ~55°C/km, evolving from greenschist to amphibolite facies, while the middle crust, exposed in a gneissic core, exhibits granulite facies conditions with a near isothermal geothermal gradient (<8°C/km) between 740 and 790°C. The abnormal and discontinuous crustal geothermal gradient, dated at c. 305 Ma on syn-granulitic monazite by LA-ICP-MS, is interpreted to be the result of magmatic intrusions at different structural levels in the crust: the Ansignan charnockite (c. 305 Ma) in the deepest part of the gneissic core, the Tournefort granodiorite (c. 308 Ma) at the interface between the gneissic core and the upper crust and the Saint-Arnac granite (c. 304 Ma) in the upper section of the massif. The heat input from these magmas combined with the thermal buffering effect of the biotite dehydration-melting reaction resulted in the near isothermal geothermal gradient in the gneissic core (melt-enhanced geotherm). The higher geothermal gradient (>50°C/km) in the upper crust is only due to conduction between the hot middle crust and the Earth's surface. The estimated maximum finite pressure range suggests that ~10 to 12 km of crust are exposed in the Agly massif while the present-day thickness does not exceed 5–6 km. This pressure/depth gap is consistent with the presence of several normal mylonitic shear zones that could have contributed to the subtraction of ~5 km of the rock pile. Monazite U–Th–Pb ages carried out on monazite overgrowths from a highly mylonitized sample suggest that this vertical thinning of the massif occurred at c. 296–300 Ma. This later Variscan extension might have slightly perturbed the 305 Ma geothermal gradient, resulting in an apparent higher conductive geothermal gradient in the upper crust. Although the Agly massif has been affected by Cretaceous extension and Eocene Alpine compression, we suggest that most of the present-day thickness of the column rock was acquired by the end of the Palaeozoic.     

Spatially distributed long‐term hydrologic simulation using a continuous SCS CN method‐based hybrid hydrologic model

A continuous Soil Conservation Service (SCS) curve number (CN) method that considers time‐varied SCS CN values was developed based on the original SCS CN method with a revised soil moisture accounting approach to estimate run‐off depth for long‐term discontinuous storm events. The method was applied to spatially distributed long‐term hydrologic simulation of rainfall‐run‐off flow with an underlying assumption for its spatial variability using a geographic information systems‐based spatially distributed Clark's unit hydrograph method (Distributed‐Clark; hybrid hydrologic model), which is a simple few parameter run‐off routing method for input of spatiotemporally varied run‐off depth, incorporating conditional unit hydrograph adoption for different run‐off precipitation depth‐based direct run‐off flow convolution. Case studies of spatially distributed long‐term (total of 6 years) hydrologic simulation for four river basins using daily NEXRAD quantitative precipitation estimations demonstrate overall performances of Nash–Sutcliffe efficiency (E_NS) 0.62, coefficient of determination (R²) 0.64, and percent bias 0.33% in direct run‐off and E_NS 0.71, R² 0.72, and percent bias 0.15% in total streamflow for model result comparison against observed streamflow. These results show better fit (improvement in E_NS of 42.0% and R² of 33.3% for total streamflow) than the same model using spatially averaged gauged rainfall. Incorporation of logic for conditional initial abstraction in a continuous SCS CN method, which can accommodate initial run‐off loss amounts based on previous rainfall, slightly enhances model simulation performance; both E_NS and R² increased by 1.4% for total streamflow in a 4‐year calibration period. A continuous SCS CN method‐based hybrid hydrologic model presented in this study is, therefore, potentially significant to improved implementation of long‐term hydrologic applications for spatially distributed rainfall‐run‐off generation and routing, as a relatively simple hydrologic modelling approach for the use of more reliable gridded types of quantitative precipitation estimations.     

The fresnel interferometric imager

The Fresnel Interferometric Imager has been proposed to the European Space Agency (ESA) Cosmic Vision plan as a class L mission. This mission addresses several themes of the CV Plan: Exoplanet study, Matter in extreme conditions, and The Universe taking shape. This paper is an abridged version of the original ESA proposal. We have removed most of the technical and financial issues, to concentrate on the instrumental design and astrophysical missions. The instrument proposed is an ultra-lightweight telescope, featuring a novel optical concept based on diffraction focussing. It yields high dynamic range images, while releasing constraints on positioning and manufacturing of the main optical elements. This concept should open the way to very large apertures in space. In this two spacecraft formation-flying instrument, one spacecraft holds the focussing element: the Fresnel interferometric array; the other spacecraft holds the field optics, focal instrumentation, and detectors. The Fresnel array proposed here is a 3.6 ×3.6 m square opaque foil punched with 10⁵ to 10⁶ void “subapertures”. Focusing is achieved with no other optical element: the shape and positioning of the subapertures (holes in the foil) is responsible for beam combining by diffraction, and 5% to 10% of the total incident light ends up into a sharp focus. The consequence of this high number of subapertures is high dynamic range images. In addition, as it uses only a combination of vacuum and opaque material, this focussing method is potentially efficient over a very broad wavelength domain. The focal length of such diffractive focussing devices is wavelength dependent. However, this can be corrected. We have tested optically the efficiency of the chromatism correction on artificial sources (500 < λ < 750 nm): the images are diffraction limited, and the dynamic range measured on an artificial double source reaches 6.2 10^− 6. We have also validated numerical simulation algorithms for larger Fresnel interferometric arrays. These simulations yield a dynamic range (rejection factor) close to 10^− 8 for arrays such as the 3.6 m one we propose. A dynamic range of 10^− 8 allows detection of objects at contrasts as high as than 10^− 9 in most of the field. The astrophysical applications cover many objects in the IR, visible an UV domains. Examples are presented, taking advantage of the high angular resolution and dynamic range capabilities of this concept.     

Weathering rates from top to bottom in a carbonate environment

This study investigates U-series, Sr isotopes, major and trace elements in a chalk aquifer system located in Eastern France. Soil and rock samples were collected along depth profiles down to 45 m in four localities as an attempt to investigate the weathering processes in the soil, the unsaturated zone and the saturated zone of the aquifer. Interstitial water was extracted from soils and rocks by a centrifugation technique. U-series offer a powerful tool to calculate weathering rates because the relative mobility of the U- and Th-isotopes can be precisely measured and it does not require the determination of a reference state as in other approaches. As expected, the data show very large mobile element depletion in the soil with large ²³⁰Th excess relative to ²³⁸U, while the rocks show more limited but not insignificant mobile element depletion. The U-series data have been used to constrain weathering rates based on a 1-D reactive transport model. Weathering rates in the near surface are about 10–100 times faster than at depth. However, when integrated over the depth of the cores, including the unsaturated and the saturated zones, this underground weathering represents more than 30% of the total weathering flux, assuming congruent dissolution of carbonates. The (²³⁴U/²³⁸U) ratios in interstitial water are consistent with solid samples showing ²³⁴U depletion near the surface and an excess ²³⁴U at depth. A leaching experiment performed on chalk shows that the excess ²³⁴U in natural waters percolating through carbonate rocks results both from preferential ²³⁴U leaching and direct recoil in the interstitial water. A new approach was used to derive the recoil ejection factor based on BET measurements and the fractal dimension of chalk surface. Consideration of preferential leaching and recoil allows a more accurate modeling of weathering rates.     

Litho- and biofacies analysis of postglacial marine mud using CT-Scanning

Marine silty clay deposited during the Late-Wisconsinian postglacial marine transgression of eastern Québec (Goldthwait Sea) is ubiquitous in the sedimentary column of intertidal zones of the St-Lawrence Estuary. This mud is very compact and limits the penetration of organisms composing the modern Macoma balthica community. In order to describe the characteristics of intertidal sediments containing Goldthwait Sea mud, axial tomography (CT-Scan) is used. CT-Scan is a non-destructive method that can be used to describe sediment characteristics (grain size, mineralogy, primary and secondary sedimentary structures, fabric, shape and roundness, bedding contact), and to obtain high resolution, 3D representations of structures within sediment cores. Based on differences in the densities of analysed materials, the different lithologies, lithofacies, and organisms within the core can be discriminated, and a quantification of the volume occupied by the different components of the material can be made. Here, CT-Scan images provide information on the distribution, orientation and interweaving of thanatocœnosis shell beds that alternate with massive or faintly laminated postglacial marine mud beds, as well as on ichnofacies characteristics. In addition, we show 3D images of bioturbation structures within the recent sediment layer, which is distinguished from the underlying Goldthwait Sea mud. When coupled with conventional sedimentary (grain size statistics) and radiochronological (¹⁴C) analyses, these data provide information which is valuable for identifying depositional processes within sedimentary environments.