Geosciences

The Kuznetsk Basin (Kuzbass) is one of the largest coal basins in Siberia and a reference area for the ancient Angaraland continent. The proximity of the Kuzbass and Siberian Platform caused their biotic similarities in the Late Palaeozoic. However, due to biota endemism, the Kuzbass Upper Palaeozoic does not correlate directly with the International Chronostratigraphic Chart (ICC). This paper discusses radioisotopic (CA-ID-TIMS) dating of zircons from a volcanic tuff located in the Starokuznetsk Formation (Fm). This level matches the interval of the Balakhonka/Kolchugino (B/K) floral change in Kuzbass, i.e., the gradual replacement of cordaitoid-dominated wet forests (Balakhonka flora) with more arid fern–pteridosperm–cordaitoid assemblages (Kolchugino flora). New age (276.9 ± 0.4 Ma) directly correlates the Starokuznetsk Fm with the Upper Kungurian of the ICC. We compared the Kuzbass data with data of the Western Verkhoyanie, where Middle Permian ammonoids (Sverdrupites assemblage) occur in strata recording the B/K floral change. The available (ICC) and new datings indicate the lag between the B/K floral change in low (Kuzbass) and high (Verkhoyanie) latitudes of Angaraland. The B/K floral change in the Kuzbass began in the early Late Kungurian and was completed by the end of this age. In contrast, the B/K floral change in Verkhoyanie began at the end of the Late Kungurian and was completed in the Late Wordian. The delay in the floral changes at different latitudes of Angaraland suggests that existing interregional correlations need further improvement. Full article

On 3 November 2023, a moment magnitude (M_W) 5.7 (Local Magnitude, M_L6.4) earthquake struck the western region of Nepal, one of the most powerful seismic events since 1505 in the region. Even though the earthquake was of moderate magnitude, it caused significant damage to several masonry buildings and caused slope failures in some regions. The field reconnaissance carried out on 6–9 November by the study team, following the earthquake, conducted the first-hand preliminary damage assessment in the three most affected districts—Jajarkot; West Rukum; and Salyan. This study covers the observed typical structural failures and geotechnical case studies from the field study. To have a robust background understanding, this paper examines the seismotectonic setting and regional seismic activity in the region. The observations of earthquake damage suggest that most of the affected buildings were made of stone or brick masonry without seismic consideration, while most of the reinforced concrete (RC) buildings remained intact. Case histories of damaged buildings, the patterns, and the failure mechanisms are discussed briefly in this paper. Significant damage to Khalanga Durbar, a historical monument in Jajarkot, was also observed. Medium- to large-scale landslides and rockfalls were recorded along the highway. The motorable bridge in the Bheri River suffered from broken bolts, rotational movement at the expansion joint, and damage to the stoppers. The damage observations suggest that, despite the existence of building codes, their non-implementation could have contributed to the heavy impact in the region. This study highlights that the local population faces a potential threat of subsequent disasters arising from earthquakes and earthquake-induced landslides. This underscores the necessity for proactive measures in preparedness for future disasters. Full article

Shallow Geothermal Energy (SGE) exchanges heat with the ground. In continuous, long-term operation, the initial temperature field can be disturbed, and subsurface thermal changes can be developed. In this paper, the thermal impact of a SGE system under a Mediterranean climate is handled. Temperature monitoring was conducted on 15 investigation boreholes equipped with a total of 92 thermal sensors placed at specific depths. Investigation boreholes were drilled 1–2 m from SGE system borehole heat exchangers installed in a university building. The analysis handles a one-year monitoring period of SGE system operation. Temperature depth profiles, reaching up to 140 m depth, were registered with a 10 min time step, resulting in a large amount of data. Ground thermal conductivity was estimated experimentally and semi-empirically, allowing us to obtain, using a numerical model, the initial undisturbed ground temperature profiles and compare them with the monitored values. Climate data were recorded by the university meteorological station. Globally, the measured and computed data were coherent, and a non-negligible impact of the SGE system operation in the first year was observed. The building orientation as well as the nearby departments had significant impacts on the shallow ground temperature. Maximum ground temperature changes observed at depths higher than 10–20 m, ranging from 2 to 3 °C as observed in different boreholes, indicate that the system is operating efficiently. Full article

The occurrence of sinkhole phenomena in Italy is a prevalent and very uncertain class of geological hazards that pose a significant threat to human infrastructure and individuals. These events are characterized by their unpredictability and the challenges associated with their accurate forecasting. Both natural and anthropic factors influence the occurrence of these events; therefore, accurate identification of the above factors is critical for effective proactive and predictive efforts. The work presented in this paper refers to a collapse that occurred in a volcanic region in northern Latium (central Italy) on 31 January 2023. The area has been monitored using drones since the early stages of the sinkhole’s formation and has continued to date. Then, the collapse and the neighboring area were examined via geophysical and geochemical investigations to identify potential underlying factors. Geophysical and geochemical data were combined to provide a preliminary hypothesis on the collapse’s genesis. The obtained data indicate that the structural collapse can be attributable to the fluctuation in groundwater levels as well as the development of instabilities along its banks, leading to a growth in its dimensions. Full article

Scaly clays are structurally complex clay formations found throughout the world. Their typical fissured structure, the low shear strength and the high swelling potential often make them unsuitable for earthworks in road and railway infrastructure. This research has attempted to extend the possibilities of using this geomaterial in this field after appropriate lime treatment. A laboratory test programme was carried out to evaluate the response of the treated geomaterial to typical loads acting on road infrastructures. Unconfined and confined compression tests as well as cyclic triaxial tests, in undrained conditions, were carried out to investigate the static and dynamic mechanical behaviour. The results show that lime treatment induces significant improvements in the geomechanical properties and limits the swelling behaviour upon saturation of the geomaterial. Dynamic tests showed that, after only 28 days of curing, the treated scaly clay became insensitive to the damaging cyclic loading caused by vehicular traffic. The collected results show that the scaly clay can be properly used as a subgrade and embankment layer in road and railway construction with limited economic and environmental costs, after accurate treatment with lime. These results are significant for researchers and practitioners to increase sustainability in the construction of linear infrastructures involving excavations in scaly clays and to avoid landfill, which in some cases represented the only option. Full article

Studies highlight researchers’ concerns about how science should be taught today. It is recognised that teachers have difficulty involving and motivating students to learn about different complex topics, such as geology. Schools must promote skills development to develop citizens who can be active and informed in society. One way of undertaking this is to use active methodologies such as educational games, in which students play an essential role in developing activities. Games encourage changes in conceptions regarding the relevance of this scientific topic that is often undervalued by students. Games have gained space in recent years in several disciplines, and it is essential that this tool is thought out and planned within a consistent pedagogical proposal. This educational resource is used to increase motivation for learning, as well as enhance and strengthen the effects of learning. An intervention plan can be framed within game-based teaching. Teachers have been underrepresented in the game-based learning literature, with more emphasis on games’ effects. However, the pedagogical issue of games has been particularly relevant in recent decades. The current investigation used a survey given to geology teachers (n = 112) from public and private middle and secondary schools in Portugal. Its purpose was to assess teachers’ perceptions regarding game-based teaching and its potential to promote active learning. Our sample ages ranged from 24 to 64 years (average of 48.9 years old); 81.3% were women and 18.8% were men. The analysis of the results seems to confirm that although they do not always use games to promote learning in geology, most teachers still recognise their potential to motivate, enhance, and reinforce the learning of geological content, with digital games being the preferred option. They emphasise the importance of teacher training in this area and the inclusion of game applications in school textbooks to approach different geology-related themes. Our results seem to indicate some lack of consistency in teachers’ opinions on the impact of games on student learning. Full article

Here, we characterize the statistical behaviour of the Mt. Vesuvius seismicity using distinct available catalogues. Our analysis confirms that for this area, the GR distribution exhibited two scaling regimes of the b-value, not commonly observed for the standard frequency-magnitude distribution of earthquakes. By assuming a physical cause, we tested four different hypotheses for the source of the break in the scaling: finite size effect, depth variations in the b-value, radial dependence in the b-value, and different b-values for swarm and non-swarm events. None of the above reasons are able to explain the observation. Thus, we investigated the possibility of some pitfalls in magnitude estimation. Based on our analysis, we suggest there is a bias in the duration magnitude the catalogues are based on. This is due to the arbitrary extrapolation to smaller magnitudes of a linear regression derived for earthquakes with $m\ge 3.0$. When a suitable correction is applied to the estimated magnitude, the GR distribution assumes the usual shape, with a b-value closer to that usually observed in volcanic areas. Finally, the analysis of the time variation of some statistical parameters reveals that the state of the volcano appears to be stationary over the entire analysed period, possibly with only a slight decrease in the b-value, indicating a small reduction in differential stress. Full article

As a transboundary river with rich and unique wetland types, the Irtysh River faces various challenges and threats from human activities and climate change, which affect area, type, and function of wetland. To accurately obtain information on the spatial and temporal distribution of wetlands in this basin, this study compares and evaluates the consistency and accuracy of a total of eleven remote sensing (RS) based land use/land cover (LULC), and wetland products. The information extraction effect of each RS product was examined through methods such as wetland area and type description, thematic map comparison, and similarity coefficient and Kappa coefficient calculations, which can reflect the wetland distribution characteristics and differences among the RS products in the Irtysh River Basin. The results show that although there is a consensus among the products in the major wetland distribution areas, there are still obvious deviations in detail depiction due to differences in factors such as data sources and methods. The products of Global 30 m Wetland Fine Classification Data (GWL_FCS30) and Global 30 m Land Cover Data (GLC_FCS30-2020) released by the Institute of Space and Astronautical Information Innovation (ISAI) of the Chinese Academy of Sciences (CAS) have a clear advantage in extracting spatial morphology features of wetlands due to the use of multi-source data, while the Esri Global 10 m Land Cover Data (ESRI_Global-LULC_10m) and products such as the global 10 m land cover data (FROM_GLC10_2017) from Tsinghua University have higher classification consistency. Moreover, data resolution, classification scheme design, and validation methods are key factors affecting the quality of wetland information extraction in the Irtysh River Basin. In practical terms, the findings of this study hold significant implications for informed decision-making in wetland conservation and management within the Irtysh River Basin. By advancing wetland monitoring technologies and addressing critical considerations in information extraction, this research effectively bridges the gap between remote sensing technology and practical applications, offering valuable insights for regional wetland protection efforts. Full article

Research was conducted on plutonic rocks, previously referred to as the Hai Van granitic rocks, distributed in Da Nang and Quang Nam provinces in Central Vietnam. The granitic rocks in Da Nang Province have low magnetic susceptibilities and have geochemical signatures typical of S-type. Additionally, a negative Eu anomaly suggests that the source rock is an organic matter-bearing sedimentary rock. The granitic rocks were likely formed during the collision between the Indochina and South China blocks. In contrast, plutonic rocks in Quang Nam Province have high magnetic susceptibilities and have geochemical signatures of I-type. No Eu anomaly was observed, and they are adakitic rocks in nature. Based on these findings, the plutonic rocks in Quang Nam Province are distinctly different from the Hai Van granitic rocks in Da Nang Province, but they are Van Canh plutonic rocks. The Hai Van granitic rocks in Da Nang Province and the Van Canh plutonic rocks in Quang Nam Province are located in the Truong Son Fold Belt. The Van Canh plutonic rocks are located farther away from the Song Ma Suture Zone than the Hai Van granitic rocks. The Van Canh plutonic rocks were generated due to the subduction of the hot Song Ma Ocean beneath the Indochina Block. The Hai Van granitic rocks are understood to have been generated in a compressional field where the Song Ma Ocean was pushing against the Indochina Block; however, the Van Canh plutonic rocks are supposed to have been generated in an extensional field, like in a back-arc-like environment generated by the subduction of the Song Ma Ocean beneath the Indochina Block. Full article

This paper presents an initial study of teachers’ perspectives and practices on students’ assessment methods regarding the subject of Geology-Geography taught in junior secondary education in Greece. While the application of descriptive assessment is recorded in other European countries, the main focus of the research is the current situation in Greece, as well as the willingness of science teachers to adopt an alternative form of evaluation and the definition of the respective circumstances. The methodology that was followed in this research included questionnaires that were sent to science teachers. Their answers were processed with the IBM SPSS 23 statistical software. The results reveal that most teaching staff surveyed positively acknowledge the need for an effective appraisal system. Their views and opinions on how student assessment contributes to the overall educational procedure were also recorded during the survey. Statistical analysis of the raw data shows that the views and opinions of teaching staff were not differentiated due to their level of experience or their level of education but only due to their field of specialization. The analysis also points to a range of additional factors that impact teaching staff’s preferences and perceptions of the effectiveness of different student appraisal and evaluation approaches. Such factors create skepticism among teachers towards alternative assessment methods, such as descriptive assessment. While this is preliminary research, it raises important issues regarding effective evaluation methods that would promote students’ development. Full article

New, detailed geological/structural mapping and field-based structural analysis were carried out to investigate the deformation pattern of well-preserved high-pressure rocks of the Blueschist Unit exposed in SE Syros (Cyclades, Greece). Geological mapping revealed the occurrence of extensive alternations between different rock groups, as well as interfingering patterns in map-scale that are possibly the result of folding. The earlier ductile deformation phase recognized in the mapped area is associated with the development of a penetrative foliation, which was formed at eclogite/blueschist-facies conditions under peak metamorphism. The subsequent main deformation phase occurred under blueschist facies conditions synchronous with the early stages of exhumation of the high-pressure rocks. This phase is mainly associated with the formation of WNW-trending folds and a pervasive axial planar foliation linked with ESE-directed shearing. The main deformation ceased under blueschist-facies conditions, and exhumation of the rocks to greenschist-facies conditions took place under very weak and localized deformation. Greenschist retrogression observed in the southwestern part of the mapped area seems to be controlled by fluids, rather than by intense deformation and formation of major syn-greenschist shear zones. Full article

In this paper, the Hellenic orogenic belt’s main geological structure and architecture of deformation are presented in an attempt to achive a better interpretation of its geotectonic evolution during Alpine orogeny. This study was based not only on recent research that I and my collaborators conducted on the deformational history of the Hellenides but also on more modern views published by other colleagues concerning the Alpine geotectonic reconstruction of the Hellenides. The structural evolution started during the Permo–Triassic time with the continental breaking of the supercontinent Pangea and the birth of the Neotethyan ocean realm. Bimodal magmatism and A-type granitoid intrusions accompanied the initial stages of continental rifting, followed by Triassic–Jurassic multiphase shallow- and deep-water sediment deposition on both formed continental margins. These margins were the Apulian margin, containing Pelagonia in the western part of the Neotethyan Ocean, and the European margin, containing continental parts of the Serbo-Macedonian and Rhodope massifs in the eastern part of the Neotethyan ocean. Deformation and metamorphism are recorded in six main deformational stages from the Early–Middle Jurassic to the present day, beginning with Early–Middle Jurassic Neotethyan intra-oceanic subduction and ensimatic island arc magmatism, as well as the formation of a suprasubduction oceanic lithosphere. Compression, nappe stacking, calc-alkaline magmatism, and high-pressure metamorphic events related to subduction processes alternated successively over time with extension, orogenic collapse, medium- to high-temperature metamorphism, adakitic and calc-alkaline magmatism, and partial migmatization related to the uplift and exhumation of deep crustal levels as tectonic windows or metamorphic core complexes. A S- to SW-ward migration of dynamic peer compression vs. extension is recognized during the Tertiary Alpine orogenic stages in the Hellenides. It is suggested that all ophiolite belts in the Hellenides originated from a single source, and this was the Neotethyan Meliata/Maliac-Axios/Vardar ocean basin, parts of which obducted during the Mid–Late Jurassic on both continental margins, Apulian (containing Pelagonia) and European (containing units of the Serbo-Macedonian/Rhodope nappe stack), W-SW-ward and E-NE-ward, respectively. In this case, the ophiolite nappes should be considered far-traveled nappes on the continental parts of the Hellenides associated with the deposition of Middle–Late Jurassic ophiolitic mélanges in basins at the front of the adjacent ophiolite thrust sheets. The upper limit of the ophiolite emplacement are the Mid–Upper Jurassic time(Callovian–Oxfordian), as shown by the deposition of the Kimmeridgian–Tithonian Upper Jurassic sedimentary carbonate series on the top of the obducted ophiolite nappes. The lowermost Rhodope Pangaion unit is regarded as a continuation of the marginal part of the Apulian Plate (External Hellenides) which was underthrust during the Paleocene–Eocene time below the unified Sidironero–Kerdylia unit and the Pelagonian nappe, following the Paleocene–Eocene subduction and closure of a small ocean basin in the west of Pelagonia (the Pindos–Cyclades ocean basin). It preceded the Late Cretaceous subduction of the Axios/Vardar ocean remnants below the European continental margin and the final closure of the Axios/Vardar ocean during the Paleocene–Eocene time, which was associated with the overthrusting of the European origins Vertiskos–Kimi nappe on the Sidironero–Kerdylia nappe and, subsequently, the final collision of the European margin and the Pelagonian fragment. Subsequently, during a synorogenic Oligocene–Miocene extension associated with compression and new subduction processes at the more external orogenic parts, the Olympos–Ossa widow and the Cyclades, together with the lower-most Rhodope Pangaion unit, were exhumed as metamorphic core complexes. Full article

In this paper we present textural and compositional data for columbite group minerals (CGMs) and associated Nb-Ta-Sn oxides from lithium-beryllium-tantalum pegmatite deposits of the Kolmozero–Voronja belt, NW Russia, with the aim of deciphering these characteristics for minerals from deposits with different mineral signatures and lithium ore grade. Minerals from four deposits, including two of world-class (Kolmozero and Polmostundra), are examined. The main controlling factors for CGM compositional ranges are the diversity and rate of magmatic fractionation, hydrothermal overprint and mineral paragenesis, following the specific geochemical signature of the different pegmatite deposits. CGM from Kolmozero include several mineral species (columbite-(Fe), columbite-(Mn), tantalite-(Fe), and tantalite-(Mn)), showing large compositional variations, mainly controlled by Nb-Ta fractionation (Ta/(Ta + Nb) = 0.16–0.70; Mn/(Mn + Fe) = 0.45–0.63). Textural patterns are various (oscillatory, homogeneous, and patchy); spongy domains and overgrowing Ta-rich rims are also observed somewhere. This indicates the involvement of numerous magmatic and hydrothermal processes. The Polmostundra CGMs are represented by columbite-(Fe) with Ta/(Ta + Nb), ranging from 0.05 to 0.39; some crystals are homogenous, and others present normal, oscillatory, mottled and reverse-zoning patterns. The Okhmylk CGMs are irregular normal, patchy and homogeneous columbite-(Fe) and columbite-(Mn), with Ta/(Ta + Nb) = 0.09–0.24 and Mn/(Mn + Fe) = 0.29–0.92, indicating the suppressed magmatic fractionation and iron drop due to precipitation of Fe minerals. Columbite-(Fe) and columbite-(Mn) from the Be-Ta Shongui deposit are less evolved, with Ta/(Ta + Nb) = 0.07–0.23 and Mn/(Mn + Fe) = 0.31–0.55. The minerals are characterized by progressive normal, oscillatory, homogeneous and irregular reverse patterns. Associated pyrochlore minerals occur both as early magmatic (Kolmozero) and late hydrothermal (Polmostundra, Okhmylk). Cassiterite is found only in the Okhmylk dykes, and is apparently of hydrothermal origin. CGM from Li pegmatites have impurities of Ti (0.01–0.05 apfu) and W (up to 0.02 apfu), whereas CGM from Be pegmatites contains elevated Ti (up to 0.09 apfu). The mineral system analysis presented here is relevant for exploration. Full article

A pilot study employing a template-matching approach on the detection of repeating earthquakes (or repeaters) in a small volume of the Irpinia fault (Southern Italy) is presented here. A catalog of repeaters nucleating over about a decade with local magnitudes ranging between $M_L$ 0.4 and $M_L$ 1.7 is generated. The events nucleate at depths of about 11 km over a fault patch with a size of about 600 m. Assuming a constant stress drop for all repeaters, the yearly coseismic slip rate of the patch is estimated in the range of 1–4.8 mm with a mean value of about 3.4 mm, which is close to the fault slip rate estimated by geodetic approaches. Given the mechanical properties of the medium in the source region of the repeaters revealed by previous tomographic studies, and in particular the high Vp/Vs ratio (close to 1.95), it is postulated that these events are fluid triggered, consistently with the mechanism proposed for the generation of most of the microseismicity in the area. The encouraging outcomes of this study suggest to focus future measurement campaigns on this area and open perspectives to extend this approach to other segments of the Irpinia fault to reveal small-scale frictional properties in an area capable to generate M7 earthquakes. Full article

The zircon from the pegmatite of the Adui granitic massif displays the unique the rare earth element (REE) distribution spectrum with the tetrad effect in REE fractionation. The tetrad effect often occurs in granitoid rocks, but it is rarely encountered in minerals, e.g., zircon. Fluid saturated with volatiles, water and trace elements is a factor responsible for the tetrad effect in the zircon. The detailed isotopic-geochemical study of the zircon has revealed several zones differing in internal structure (in the back-scattered electron (BSE) image), composition and REE distribution. The zones indicate changes in the crystallization environment provoked by the evolution of the pegmatite-forming melt. They occur as the gradually growing changes in composition from the unaltered zones that are light-colored in BSE to the altered zones that are dark-colored in BSE. The unaltered zones are consistent in composition and geochemical features with magmatic zircons. The high content of trace (U, Th, REE) and volatile elements (F, Cl) in water suggests its crystallization from the fluid-saturated magmatic melt. The altered zircon zones occur as recrystallized zones with high content of non-formula elements (Y, Ca, Sr, Nb, P), a non-differentiated REE distribution spectrum and an absent Ce anomaly. These features are consistent with those of hydrothermal-metasomatic zircon. Full article

A destructive Mw = 6.2 earthquake struck NW Croatia on 29 December 2020, ranking among the most hazardous events in the country in the last century. This study assesses the impact of the earthquake on the permanent geodetic stations of the Croatian national positioning system. Using a novel approach, we combined positioning station data with satellite radar data for the first time in a geodynamic study across Croatia. The analysis involved collecting daily solutions data from 11 nearby network stations as well as multi-temporal satellite images before and after the earthquake. The results reveal a significant co-seismic shift, with the Sisak reference station experiencing the largest displacement of 5 cm to the southeast. Stations up to 100 km from the epicenter exhibited horizontal shifts between 1 and 2.5 cm and vertical shifts between 2 and 4 cm. Satellite images confirmed the co-seismic effects on the permanent stations. The study establishes a strong correlation between geodetic and geological findings, suggesting that vertical displacements result from superficial deposit consolidation in the river plains due to the earthquake shaking, while positional displacement reflects tectonic block movements along the activated strike-slip fault system. The late 2020 earthquake may have further relaxed the accumulated strain on the Mt. Medvednica thrust fault system, which had ruptured in early 2020 with a lower magnitude than expected. Full article

The endorheic Konya Basin is a vast aggradational plain in Central Anatolia, Türkiye. It occupies a significant portion of Konya Province, covering approximately 50,000 km². The basin is subjected to intense groundwater withdrawal and extensive agricultural activities with excessive irrigation. These activities have led to human-induced hazards, such as sinkholes and regional land subsidence. Although sinkhole occurrence mainly occurs in the Karapınar area, land subsidence is primarily observed in the central sector of Konya city, with 2 million inhabitants, as well as in various parts of the basin. This study focuses on determining the extent and rate of land subsidence throughout the basin, understanding sinkhole formation, and unraveling their relationship with anthropogenic activities. For this purpose, Interferometric Synthetic Aperture Radar (InSAR) analysis of Sentinel-1 data from 2014 to 2022 was conducted to identify and assess land subsidence. We also used the land cover data and groundwater-level information to better understand the spatial and temporal patterns of land subsidence and sinkhole occurrence. Additionally, the land cover data were used to resolve spatial–temporal variations in the cultivated area and urbanization, which are the main factors governing groundwater exploitation in the region. Our study identified widespread subsidence zones with rates as high as 90 mm/y. Groundwater overexploitation to sustain extensive agricultural operations is the main cause of the high rate of land subsidence. Additionally, it was discovered that the number of sinkholes has substantially increased due to anthropogenic influences, currently amounting to as many as 660. Full article

Restoration of open-pit mines may utilize waste rock for landscape reconstruction, which can include the construction of backfill aquifers. Weathering and contaminant transport may be different in backfill aquifers compared to the surrounding aquifer because of newly available mineral surfaces and transportable nano- to micro-scale particles generated during mining. Waste rock from the Cordero Rojo open-pit coal mine in the Powder River Basin was exposed to benchtop leachate experiments for 20 weeks at temperatures of 5 °C and 20 °C. Collected leachate was analyzed for Eh, pH, specific conductance, alkalinity, and cation and anion concentrations as unfiltered and 0.45-μm and 0.2-μm filtered concentrations. During the experiment, leachate Eh and pH substantially varied during the first 55 days, which corresponds to a period of high specific conductance (>1000 µS/cm) and alkalinity (>200 mg/L). Correspondingly, anion and cation concentrations were the largest during this early weathering stage, and the filter fractions indicated multiple forms of transported elements. After this early weathering stage, column leachate evolved towards a weathering equilibrium of neutral, oxidizing, and low solute conditions indicated by positive Eh values, pH near 7, and specific conductance <500 μS/cm. This evolution was reflected in the decline and stabilization or non-detection of metal(loid) concentrations reflective of a shift to primarily bulk aluminosilicate weathering when coal- and salt-associated elements, such as arsenic, cadmium, and selenium, were not detected or at minimal concentrations. Over the course of the experiment, the solute trend of certain elements indicated particular weathering processes—cadmium and nanoparticle transport, selenium and salt dissolution, and arsenic and pyrite oxidation. The mining of overburden formations and use of the waste rock for backfill aquifers as part of landscape reconstruction will create newly available mineral surfaces and nanoparticles that will weather to produce solute concentrations not typically found in groundwater associated with the original overburden. Full article

Groundwater resource management requires understanding the groundwater basin’s hydrogeology and would be improved with the development of a three-dimensional hydrogeologic framework model (HFM). A wide range of methods and software exist to quantify the extent, structure, and properties of geologic systems. However, most geologic software is proprietary and cost-prohibitive for use in developing countries. GemPy is a Python-based, open-source (no-cost) tool for generating three-dimensional geological models. This study uses available data and GemPy to develop the Kobo Valley Hydrogeologic Framework Model (KV-HFM), a three-dimensional HFM for Kobo Valley in northern Ethiopia, which is part of the East African Rift System. The KV-HFM is a conceptual model that comprises the hydrostratigraphy, structural features, and hydraulic properties of the Kobo Valley groundwater system. The limited data described the extent and altitude of the hydrostratigraphic units using the GemPy implicit potential–field interpolation. The KV-HFM showed the existence of an east-to-west, structural-based groundwater divide composed of volcanic rock and clay. This divide splits the catchment into two groundwater systems with limited interconnected flow. This study illustrates the use of open-source software for developing an HFM using sparse, existing geologic data. Full article