Cenozoic metallogeny of Greece and potential for precious,critical and rare metals exploration

The Cenozoic metallogeny in Greece includes numerous major and minor hydrothermal mineral deposits, associated with the closure of the Western Tethyan Ocean and the collision with the Eurasian continental plate in the Aegean Sea, which started in the Cretaceous and is still ongoing. Mineral deposits formed in four main periods: Oligocene (33–25 Ma), early Miocene (22–19 Ma), middle to late Miocene (14–7 Ma), and Pliocene-Pleistocene (3–1.5 Ma). These metallogenic periods occurred in response to slab-rollback and migration of post-collisional calc-alkaline to shoshonitic magmatism in a back-arc extensional regime from the Rhodopes through the Cyclades, and to arc-related magmatism along the active south Aegean volcanic arc. Invasion of asthenospheric melts into the lower crust occurred due to slab retreat, and were responsible for partial melting of metasomatized lithosphere and lower crustal cumulates. These geodynamic events took place during the collapse of the Hellenic orogen along large detachment faults, which exhumed extensive metamorphic core complexes in mainly two regions, the Rhodopes and the Cyclades. The detachment faults and supra-detachment basins controlled magma emplacement, fluid circulation, and mineralization.The most significant mineralization styles comprise porphyry, epithermal, carbonate-replacement, reduced intrusion-related gold, intrusion-related Mo-W and polymetallic veins. Porphyry and epithermal deposits are commonly associated with extensive hydrothermal alteration halos, whereas in other cases alteration is of restricted development and mainly structurally controlled. Porphyry deposits include Cu-Au-, Cu-Mo-Au-Re, Mo-Re, and Mo-W variants. Epithermal deposits include mostly high- and intermediate-sulfidation (HS and IS) types hosted in volcanic rocks, although sedimentary and metamorphic rock hosted mineralized veins, breccias, and disseminations are also present. The main metal associations are Cu-Au-Ag-Te and Pb-Zn-Au-Ag-Te in HS and IS epithermal deposits, respectively. Major carbonate-replacement deposits in the Kassandra and Lavrion mining districts are rich in Au and Ag, and together with reduced intrusion-related gold systems played a critical role in ancient economies. Finally hundreds of polymetallic veins hosted by metamorphic rocks in the Rhodopes and Cyclades significantly add to the metal endowment of Greece.     

Combined zircon and cassiterite U–Pb dating of the Piaotang granite-related tungsten–tin deposit,southern Jiangxi tungsten district,China

The large low-grade Piaotang W–Sn deposit in the southern Jiangxi tungsten district of the eastern Nanling Range, South China, is related to a hidden granite pluton of Jurassic age. The magmatic-hydrothermal system displays a zonation from an inner greisen zone to quartz veins and to peripheral veinlets/stringers (Five-floor zonation model). Most mineralization is in quartz veins with wolframite > cassiterite. The hidden granite pluton in underground exposures comprises three intrusive units, i.e. biotite granite, two-mica granite and muscovite granite. The latter unit is spatially associated with the W–Sn deposit.Combined LA-MC-ICP-MS U–Pb dating of igneous zircon and LA-ICP-MS U–Pb dating of hydrothermal cassiterite are used to constrain the timing of granitic magmatism and hydrothermal mineralization. Zircon from the three granite units has a weighted average ²⁰⁶Pb/²³⁸U age of 159.8 ± 0.3 Ma (2 σ, MSWD = 0.3). The cathodoluminescence (CL) textures indicate that some of the cassiterite crystals from the wolframite-cassiterite quartz vein system have growth zonations, i.e. zone I in the core and zone II in the rim. Dating on cassiterite (zone II) yields a weighted average ²⁰⁶Pb/²³⁸U age of 159.5 ± 1.5 Ma (2 σ, MSWD = 0.4), i.e. the magmatic and hydrothermal systems are synchronous. This confirms the classical model of granite-related tin–tungsten mineralization, and is against the view of a broader time gap of >6 Myr between granite magmatism and W–Sn mineralization which has been previously proposed for the southern Jiangxi tungsten district. The elevated trace element concentrations of Zr, U, Nb, Ta, W and Ti suggest that cassiterite (zone II) formed in a high-temperature quartz vein system related to the Piaotang granite pluton.     

Dating the giant Zhuxi W–Cu deposit (Taqian–Fuchun Ore Belt) in South China using molybdenite Re–Os and muscovite Ar–Ar system

The recently discovered Zhuxi W–Cu ore deposit is located within the Taqian–Fuchun Ore Belt in the southeastern edge of the Yangtze Block, South China. Its inferred tungsten resources, based on new exploration data, are more than 280 Mt by 2016. At least three paragenetic stages of skarn formation and ore deposition have been recognized: prograde skarn stage; retrograde stage; and hydrothermal sulfide stage. Secondly, greisenization, marmorization and hornfels formation are also observed. Scheelite and chalcopyrite are the dominant metal minerals in the Zhuxi deposit and their formation was associated with the emplacement of granite stocks and porphyry dykes intruded into the surrounding Carboniferous carbonate sediments (Huanglong and Chuanshan formations) and the Neoproterozoic slate and phyllites. The scheelite was mostly precipitated during the retrograde stage, whereas the chalcopyrite was widely precipitated during the hydrothermal sulfide stage. A muscovite ⁴⁰Ar/³⁹Ar plateau age of about 150 Ma is interpreted as the time of tungsten mineralization and molybdenite Re–Os model ages ranging from 145.9 ± 2.0 Ma to 148.7 ± 2.2 Ma (for the subsequent hydrothermal sulfide stage of activity) as the time of the copper mineralization. Our new molybdenite Re–Os and muscovite ⁴⁰Ar/³⁹Ar dating results, along with previous zircon U–Pb age data, indicate that the hydrothermal activity from the retrograde stage to the last hydrothermal sulfide stage lasted up to 5 Myr, from 150.6 ± 1.5 to 145.9 ± 1 Ma, and is approximately coeval or slightly later than the emplacement of the associated granite porphyry and biotite granite. The new ages reported here confirm that the Zhuxi tungsten deposit represents one of the Mesozoic magmatic–hydrothermal mineralization events that took place in South China in a setting of lithospheric extension during the Late Jurassic (160–150 Ma). It is suggested that mantle material played a role in producing the Zhuxi W–Cu mineralization and associated magmatism.     

The Cenomanian–Turonian Boundary Event (CTBE) on the southern slope of the Subalpine Basin (SE France) and its bearing on a probable tectonic pulse on a larger scale

The Cenomanian–Turonian Boundary Event (CTBE) event is not associated with a transgression on the southern margin of the Subalpine Basin, but with a steady shallowing-up trend beginning in the lower half of the δ¹³C positive shift. The SW–NE Rouaine Fault had a complex role, first in isolating a black shale basin to the west and a large, deep submarine plateau devoid of black shale to the east, then by a strike-slip movement that induced a forced progradation to the north of the southern platform in the eastern compartment. This compressive tectonic reactivation of the southern margin began around the deposition of the local equivalent of the Plenus bed of boreal basins, as shown by correlation supported by both isotope and palaeontological data. Other local data are pieced together to suggest that the whole of SE France underwent a short-lived transpressive tectonic pulse around the Cenomanian–Turonian boundary, probably connected with the early compressive movement of Africa vs. Europe. On a larger scale, other published data suggest that this pulse could be a global one. It is coeval with renewed thrust loading, volcanism and transgression in the North-American Western Interior, local emergences during the event along the eastern Atlantic margin, suggesting a slight tendency to inversion of the margin, and a tilting to the east of the North-Africa plate that could explain the large transgression recorded from Morocco to Tunisia on the Saharan Craton.New isotope and palaeontological (coiling ratio of Muricohedbergella delrioensis) data from SE France suggest that two coolings of suprabasinal importance occurred just before and during the build-up of the d¹³C shift, including the boreal “Plenus Marls“, especially its middle limestone bed and its SE France equivalent.Regarding the extinction of the genus Thalmaninella and Rotalipora and during the event, neither anoxia nor climate changes can fully explain the palaeontological crisis, given that Rotalipora cushmani crosses the first phase of anoxia without harm, as well as the two coolings, not only in SE France but on a large scale, as shown by the correlation of the published data. This extinction needs alternative explanations as we challenge both anoxia and climate as major causes.     

Zircon U–Pb ages of granitic and mafic dikes associated with granitoids and metamorphic rocks of the Ryoke Belt in the eastern Kii Peninsula,southwest Japan

Zircon U–Pb dates for felsic and intermediate to mafic dikes intruding into the Ryoke granitoids and metamorphic rocks at selected outcrops in the Takamiyama area of the eastern Kii Peninsula, southwest Japan, were determined along with their geology and petrography to reveal the history of Cretaceous magmatism. At each outcrop, the felsic and intermediate to mafic dikes exhibit specific structures that are indicative of magma intermingling and have coeval intrusion ages of ca. 81–77 Ma. Our zircon U–Pb data complement previously published data, suggesting that the mafic magmatism continued intermittently from 83 to 76 Ma in the Takamiyama area and that magmatism migrated eastward within the Ryoke Belt. A comparison of intrusion ages between a dike and a host Ryoke granitoid at one outcrop indicates that the host rock experienced ductile deformation at ~88 to ~83 Ma. Judging from the small number of zircons and the concordant date distributions, we didn't recognize the evidence suggesting the partial melting of the host rocks, as Nakajima et al. (Journal of the Geological Society of Japan, 2021, 127, 69–78) reported.     

Spatial rainfall model using a pattern classifier for estimating missing daily rainfall data

Missing data in daily rainfall records are very common in water engineering practice. However, they must be replaced by proper estimates to be reliably used in hydrologic models. Presented herein is an effort to develop a new spatial daily rainfall model that is specifically intended to fill in gaps in a daily rainfall dataset. The proposed model is different from a convectional daily rainfall generation scheme in that it takes advantage of concurrent measurements at the nearby sites to increase the accuracy of estimation. The model is based on a two-step approach to handle the occurrence and the amount of daily rainfalls separately. This study tested four neural network classifiers for a rainfall occurrence processor, and two regression techniques for a rainfall amount processor. The test results revealed that a probabilistic neural network approach is preferred for determining the occurrence of daily rainfalls, and a stepwise regression with a log-transformation is recommended for estimating daily rainfall amounts.     

Effects of antecedent moisture and macroporosity on infiltration and water flow in frozen soil

Infiltration into frozen soil plays an important role in soil freeze–thaw and snowmelt-driven hydrological processes. To better understand the complex thermal energy and water transport mechanisms involved, the influence of antecedent moisture content and macroporosity on infiltration into frozen soil was investigated. Ponded infiltration experiments on frozen macroporous and non-macroporous soil columns revealed that dry macroporous soil produced infiltration rates reaching 10³ to 10⁴ mm day⁻¹, two to three orders of magnitude larger than dry non-macroporous soil. Results suggest that rapid infiltration and drainage were a result of preferential flow through initially air-filled macropores. Using recorded flow rates and measured macropore characteristics, calculations indicated that a combination of both saturated flow and unsaturated film flow likely occurred within macropores. Under wet conditions, regardless of the presence of macropores, infiltration was restricted by the slow thawing rate of pore ice, producing infiltration rates of 2.8 to 5.0 mm day⁻¹. Reduced preferential flow under wet conditions was attributed to a combination of soil swelling, due to smectite-rich clay (that reduced macropore volume), and pore ice blockage within macropores. In comparison, dry soil column experiments demonstrated that macropores provided conduits for water and thermal energy to bypass the frozen matrix during infiltration, reducing thaw rates compared with non-macroporous soils. Overall, results showed the dominant control of antecedent moisture content on the initiation, timing, and magnitude of infiltration and flow in frozen macroporous soils, as well as the important role of macropore connectivity. The study provides an important data set that can aid the development of hydrological models that consider the interacting effects of soil freeze–thaw and preferential flow on snowmelt partitioning in cold regions.     

Characterizing hydrological processes within kettle holes using stable water isotopes in the Uckermark of northern Brandenburg,Germany

Understanding the hydrologic connectivity between kettle holes and shallow groundwater, particularly in reaction to the highly variable local meteorological conditions, is of paramount importance for tracing water in a hydro(geo)logically complex landscape and thus for integrated water resource management. This article is aimed at identifying the dominant hydrological processes affecting the kettle holes' water balance and their interactions with the shallow groundwater domain in the Uckermark region, located in the north-east of Germany. For this reason, based on the stable isotopes of oxygen (δ¹⁸O ) and hydrogen (δ²H ), an isotopic mass balance model was employed to compute the evaporative loss of water from the kettle holes from February to August 2017. Results demonstrated that shallow groundwater inflow may play the pivotal role in the processes taking part in the hydrology of the kettle holes in the Uckermark region. Based on the calculated evaporation/inflow (E/I) ratios, most of the kettle holes (86.7%) were ascertained to have a partially open, flow-through-dominated system. Moreover, we identified an inverse correlation between E/I ratios and the altitudes of the kettle holes. The same holds for electrical conductivity (EC) and the altitudes of the kettle holes. In accordance with the findings obtained from this study, a conceptual model explaining the interaction between the shallow groundwater and the kettle holes of Uckermark was developed. The model exhibited that across the highest altitudes, the recharge kettle holes are dominant, where a lower ratio of E/I and a lower EC was detected. By contrast, the lowest topographical depressions represent the discharge kettle holes, where a higher ratio of E/I and EC could be identified. The kettle holes existing in between were categorized as flow-through kettle holes through which the recharge takes place from one side and discharge from the other side.     

Response of a scale-model pile group for a jacket foundation of an offshore wind turbine in liquefiable ground during shaking table tests

To investigate the seismic response of a pile group during liquefaction, shaking table tests on a 1/25 scale model of a 2 × 2 pile group were conducted, which were pilot tests of a test project of a scale-model offshore wind turbine with jacket foundation. A large laminar shear box was utilized as the soil container to prepare a liquefiable sandy ground specimen. The pile group model comprising four slender aluminum piles with their pile heads connected by a rigid frame was designed with similitude considerations focusing on soil–pile interaction. The input motions were 2-Hz sinusoids with various acceleration amplitudes. The excess pore water pressure generation indicated that the upper half of the ground specimen reached initial liquefaction under the 50-gal-amplitude excitation, whereas in the 75-gal-amplitude test, almost entire ground was liquefied. Accelerations in soil, on the movable frames composing the laminar boundary of the shear box, and along the pile showed limited difference at the same elevation before liquefaction. After liquefaction, the soil and the movable-frame accelerations that represented the ground response considerably reduced, whereas both the movable frames and the piles exhibited high-frequency jitters other than 2-Hz sinusoid, and meantime, remarkable phase difference between the responses of the pile group and the ground was observed, all probably due to the substantial degradation of liquefied soil. Axial strains along the pile implied its double-curvature bending behavior, and the accordingly calculated moment declined significantly after liquefaction. These observations demonstrated the interaction between soil and piles during liquefaction.