Major shifts in biogeographic regions of freshwater fishes as evidence of the Anthropocene epoch

Animals and plants worldwide are structured in global biogeographic regions, which were shaped by major geologic forces during Earth history. Recently, humans have changed the course of events by multiplying global pathways of introduction for nonindigenous species and propagating local species extirpations. Here, we report on how introductions and extirpations have changed the distributions of freshwater fishes worldwide and how it affected their natural biogeographic regions. We found major shifts in natural regions, with the emergence of an intercontinental region arising from the fusion of multiple faunas, which we named Pan-Anthropocenian Global North and East Asia (PAGNEA). The PAGNEA region is evocative of the Pangea supercontinent, as flows of introductions show that dispersal has become possible again across multiple continents, suggesting that human activities have superseded natural geological forces. Our results constitute evidence on the expected modification of biostratigraphic boundaries based on freshwater fish, which are abundant in the fossil record, thereby supporting the concept of the Anthropocene epoch.

Oklo, our scientific adventure with R. Naudet

Following several years of studies of the Francevillian basin, I was a specialist of the Gabonese geology. Since the discovery of the Oklo phenomenon, I worked with Roger Naudet, Jean-Paul Pfiffelmann and François Gauthier-Lafaye on this exceptional scientific case. This article relates the beginning of this adventure.

Characterization of the genesis and spatial variability in tectonic fractures in the Gaosong ore field (South China)

This study examines fractures in the Gaosong ore field to determine the main factors affecting the spatial variability in the fracture structure. The attributes of fractures, including the fracture orientation, intensity and intersection density, in the Wuzishan anticline and near the Lotus mountain fault in the Gaosong ore field in the GeJiu ore district were extracted by using a modified circular scanning line method. The fracture intensity and intersection density were analyzed based on the semivariance geostatistics function by using the volume of variation and the amount of relative variability. These parameters quantitatively describe the spatial variability in the fracture structure. The mean and standard variance of fracture intensity and intersection density in the ore field decrease with distance from the Lotus mountain fault, while the coefficient of variation increases. The spatial anisotropy is closely related to the axial direction of the Wuzishan anticline and the orientation of the Lotus mountain fault. The main factors affecting the spatial variability in the fault structure can be determined with the semivariance geostatistics function, and the results are useful for studying the geology of the mining area and can help to construct an accurate structural model to serve the needs of mine production.

Geological control of terrestrial background radiation in Garhwal Himalaya, India

Activity concentrations of 226Ra, 232Th, and 40K were measured in soil samples from several areas of Garhwal Himalaya, Northern India, by gamma-ray spectrometry. In this region, which extends around the Himalayan Main Central Thrust, a tectonic line that separates several geological provinces, background levels of natural terrestrial radiation were assessed. The maximum levels of radium, 285 Bq/kg and 136 Bq/kg, respectively, were found in the Budhakedar and Uttarkashi regions, exceeding the world average value of 35 Bq/kg. The mean radiation levels were found to be different between the areas, which reflects the geological diversity in the region. The overall absorbed dose rate owing to radionuclide presence in the Uttarkashi area ranged from 79 to 188 nGyh-1, with an average of 118 nGyh-1. That is more than UNSCEAR's world-populated weighted average value of 59 nGyh-1. The present investigation indicated that the absorbed dose rates are greater in Uttarkashi and Budhakedar than in other places. The multiple comparison analysis between geology and absorbed dose rate shows that the geology of Uttarkashi and Budhakedar are statistically similar. According to several hazard indices, terrestrial background radiation is not of radiological concern in the investigated region.

The effect of silvicultural systems on soil function depends on bedrock geology and altitude

Soil quality and function in forest environments are influenced by the interaction of soil-forming parameters and silvicultural systems. Hyrcanian forests were recently accepted as a UNESCO World Heritage Site, which extends across an area of approximately 1.8 million hectares and ascend to an elevation of 2800 m above sea level (m.a.s.l). In these woodlands, Oriental Beech (Fagus orientalis Lipsky) is the predominant tree species and could be observed at 700-1500 m.a.s.l., and occur on different parent rocks. Shelterwood and single-tree selection techniques have been the primary management methods for beech forests for the past forty years. Studies investigating the impacts of silvicultural systems have not yet been done on soil and forest floor features on different bedrock geology and altitudes. Therefore, in this study, we examined the influence of single-tree selection and shelterwood methods, 25 years after employing those methods, on soil quality and function compared to control areas (intact forests) in Hyrcanian beech stands. For this purpose, 15 forest floor (30 × 30 cm) and topsoil (0-10 cm depth) samples in each silvicultural systems (i.e., single-tree selection and shelterwood methods and control zones) × 4 regions (including Rasht, Nowshahr, Sari and Gorgan) × 4 altitude levels (with averages of 800, 1000, 1200 and 1400 m.a.s.l.) were considered. According to our findings, the investigated forest regions, forest floor and soil characteristics across various locations spots could be separated by principal component analysis output, and more than 85% of the variance was explained by the first and second axes. The structural equation model showed that the region, altitude and silvicultural systems had an effective role in the changes in soil biological activities by influencing the forest floor, and the soil physicochemical features. Based upon the network model, the C/N ratio, the sand content, the soil aggregate stability, the available K, the fulvic acid, and the Acarina density were found to be prominent factors with regard to soil function. In the control sites, increased soil organic material fractions, microbial/enzyme and biota activities were detected, particularly at the lower altitudes of the Nowshahr site, which had geological formations of dolomite and calcic layers. Taken together, it seems that the single-tree method, commonly referred to as the close-to-nature technique produces more suitable conditions for soil functioning compared to the shelterwood management approach. Silvicultural systems, bedrock geology and altitude can have major detrimental effects on soil function and fertility, over the long-term, impacts may increase with harvest intensity.

Linking health to geology-a new assessment and zoning model based on the frame of medical geology

With the growing concerns about the Earth's environment and human health, there has been a surge in research focused on the intersection of health and geology. This study quantitatively assesses the relationship between human health and geological factors using a new framework. The framework considers four key geological environment indicators related to health: soil, water, geological landform, and atmosphere. Results indicate that the atmospheric and water resource indicators in the study area were generally favorable, while the scores of geological landforms varied based on topography. The study also found that the selenium content in the soil greatly exceeded the local background value. Our research underscores the importance of geological factors on human health, establishes a new health-geological assessment model, and provides a scientific foundation for local spatial planning, water resource development, and land resource management. However, due to varying geological conditions worldwide, the framework and indicators for health geology may need to be adjusted accordingly.

Biased, Spasmodic, and Ridiculously Incomplete: Sequence Stratigraphy and the Emergence of a New Approach to Stratigraphic Complexity in Paleobiology, 1973-1995

This paper examines the emergence of a new approach to stratigraphic complexity, first in geology and then, following its creative appropriation, in paleobiology. The approach was associated with a set of models that together transformed stratigraphic geology in the decades following 1970. These included the influential models of depositional sequences developed by Peter Vail and others at Exxon. Transposed into paleobiology, they gave researchers new resources for studying the incompleteness of the fossil record and for removing biases imposed by the processes of sedimentary accumulation. In addition, they helped reconfigure the cultural landscape of paleobiology, consolidating a growing emphasis on fieldwork and eroding the barrier that had been erected in the 1970s between "paleontology" and "paleobiology." This paper traces these developments, paying special attention to the simulation models of stratigraphic paleobiologist Steven Holland. It also considers how the integration of sequence and event stratigraphy and paleobiology has begun to influence long-running discussions of incompleteness and bias in the fossil record.

Stability and probabilistic assessment of a large asymmetric bedding slope with a weak interlayer

The mining of open-pit mines has a significant influence on the surrounding ecological environment and safety. Once the slope of an open pit is unstable, it will bring extremely severe disasters to the surrounding environment. The south slope of the Fushun West open-pit mine is a large bedding slope with weak interlayer and asymmetrical features between the eastern and western regions. The different deformation mechanisms are worth further discussion, and the deformation mechanism of the eastern region needs to be systematically investigated. In this study, a comprehensive method combining a numerical model, the limit equilibrium method, and a probabilistic approach is used. The different deformation mechanisms of the slope are explored. Thereafter, the sensitivity of the geometric and mechanical parameters of the slope is determined. Finally, the probability of failure (PF) and slope stability are evaluated. The results illustrate that the difference in geological characteristics and constraint conditions is the fundamental reason for the difference in deformation features of the slope. The sensitivity analysis results show that the friction angle of the weak layer exhibits the most sensitivity and that the unit weight is irrelevant to the slope stability. Taking the friction angle of the weak layer as a variable, when the thickness of the basalt layer is less than 140 m and the height of the slope is more than 360 m, the PF of the slope reaches a high value. When the inclination of the weak layer exceeds 30°, the PF of the slope increases to a high value. When the Ru coefficient exceeds 0.2, the slope presents a high PF. The results illustrate that decision-makers should pay attention to the strength degradation, transitional mining of the basalt layer, and establish a good drainage system. The comprehensive method can also provide an effective way to assess the stability and PF of similar large-scale bedding slopes under complex geological conditions.

Bird clades with less complex appendicular skeletons tend to have higher species richness

Species richness is strikingly uneven across taxonomic groups at all hierarchical levels, but the reasons for this heterogeneity are poorly understood. It is well established that morphological diversity (disparity) is decoupled from taxonomic diversity, both between clades and across geological time. Morphological complexity has been much less studied, but there is theory linking complexity with differential diversity across groups. Here we devise an index of complexity from the differentiation of the fore and hind limb pairs for a sample of 983 species of extant birds. We test the null hypothesis that this index of morphological complexity is uncorrelated with clade diversity, revealing a significant and negative correlation between the species richness of clades and the mean morphological complexity of those clades. Further, we find that more complex clades tend to occupy a smaller number of dietary and habitat niches, and that this proxy for greater ecological specialisation correlates with lower species richness. Greater morphological complexity in the appendicular skeleton therefore appears to hinder the generation and maintenance of species diversity. This may result from entrenchment into morphologies and ecologies that are less capable of yielding further diversity.

Paramagnetic characterization of fossil mollusc shells at eastern part of the old Konya lake: its importance for EPR dating

Fossil mollusc shells are used for dating geological materials because they are well preserved throughout geological time. In this study, the radicals in the structure of fossil mollusc shells (Dreissena iconica, Valvata piscinalis, Bithynia tentaculate, Unio pictorum) collected from the Eastern Part of Old Konya Lake in Türkiye were investigated by EPR technique. For all fossil shells, microwave and temperature dependence of the signals were examined, and the signals suitable for dating are discussed. Characteristic features of intrinsic and impurity-related radicals were identified and the importance of paleontological evaluation of molluscs to get a reliable equivalent dose in EPR dating studies was emphasised.

Evaluating and prioritizing the healthcare waste disposal center locations using a hybrid multi-criteria decision-making method

Healthcare waste disposal center location (HCWDCL) impacts the environment and the health of living beings. Different and sometimes contradictory criteria in determining the appropriate site location for disposing of healthcare waste (HCW) complicate the decision-making process. This research presents a hybrid multi-criteria decision-making (MCDM) method, named PROMSIS, to determine the appropriate HCWDCL in a real case. The PROMSIS is the combination of two well-known MCDM methods, namely TOPSIS and PROMETHEE. Moreover, fuzzy theory is used to describe the uncertainties of the problem parameters. To provide a reliable decision on selecting the best HCWDCL, a comprehensive list of criteria is identified through a literature review and experts' opinions obtained from the case study. In total, 40 criteria are identified and classified into five major criteria, namely economic, environmental, social, technical, and geological. The weight of the considered criteria is determined by the Analytical Hierarchy Process (AHP) method. Then, the score of the alternative HCWDCLs in each considered criterion is obtained. Finally, the candidate locations for disposing of HCWs are ranked by the proposed fuzzy PROMSIS method. The results show that the most important criteria in ranking the alternatives in the studied case are economic, environmental, and social, respectively. Moreover, the sub-criteria of operating cost, transportation cost, and pollution are identified as the most important sub-criteria, respectively.

Investigation on large-scale 3D seepage characteristics of a pumped-storage power station: a case study in Zhejiang Province, China

Pumped-storage power stations (PSPSs) have higher requirements for anti-seepage compared with regular power stations. As a result, investigating the seepage distributions of PSPSs is particularly important. However, existing researches remain limited in assessing engineering needs such as ensuring the efficiency of a power station. Taking the Qingyuan PSPS as a typical case, this study aims to investigate the large-scale seepage field distribution while exploring the efficiency of the anti-seepage system. Considering the geological characteristics and structural location, a 3D finite element model is established. Based on the continuous medium model while combined with seepage control measures, the change in leakage while the anti-seepage system failed is further assessed. It is concluded that the operation status of anti-seepage measures will have a certain impact on the leakage volumes of each part. Using a comprehensive assessment, anti-seepage measures can effectively prevent seepage. When failure occurs on anti-seepage curtains, the leakage volume at the corresponding position will show an obvious growth. In summary, the findings of this study highlight the significance of avoiding excessive leakage caused by anti-seepage structure failure, the effective operation of anti-seepage measures must be ensured. The abovementioned results can provide scientific support for the seepage optimization design of PSPSs.

Ensemble models based on radial basis function network for landslide susceptibility mapping

Ensemble learning techniques have shown promise in improving the accuracy of landslide models by combining multiple models to achieve better predictive performance. In this study, several ensemble methods (Dagging, Bagging, and Decorate) and a radial basis function classifier (RBFC) were combined to predict landslide susceptibility in the Trung Khanh district of the Cao Bang Province, Vietnam. The ensemble models were developed using a geospatial database containing 45 historical landslides (1074 points) and thirteen influencing variables characterizing the topography, geology, land use/cover, and human activities of the study area. The performance of the models was evaluated based on the area under the receiver operating characteristic curve (AUC) and several other performance metrics, including positive predictive value (PPV), negative predictive value (NPV), sensitivity (SST), specificity (SPF), accuracy (ACC), and root mean square error (RMSE). The Bagging-RBFC model with PPV = 86%, NPV = 95%, SST = 95%, SPF = 87%, ACC = 91%, RMSE = 0.297, and AUC = 98% was found to be the most accurate model for the prediction of landslide susceptibility, followed by the Dagging-RBFC, Decorate-RBFC, and single RBFC models. The study demonstrates the efficacy of ensemble learning techniques in developing reliable landslide predictive models, which can ultimately save lives and reduce infrastructure damage in landslide-prone regions worldwide.

Colima volcano's archive of observations: The invention of a geological history from Johann Mortiz Rugendas to Paul Waitz

In 1936 the Austrian geologist Paul Waitz published a seminal bibliographical, historical essay on Colima volcano, Mexico. His article exemplifies well the paths by which geology became what Lorraine Daston has termed sciences of the archive, that is, the manner in which scientific disciplines became concerned with archival work. Waitz's historical description of studies of Colima volcano built a genealogy of observations, ultimately constructing a history of the volcano itself. By bringing attention not only to Waitz's discourse but also to his treatment of visual objects, such as pictorial and photographic landscapes, my article points out how long-term aesthetics, such as the picturesque and the sublime, functioned as tropes which enabled a standardized perception, essential to visualize a clear history of scientific observations, from the landscape paintings of the nineteenth-century artist Johann Moritz Rugendas to Waitz's own photographs.

Elevational surveys of Sulawesi herpetofauna 1: Gunung Galang, Gunung Dako Nature Reserve

The Indonesian island of Sulawesi has a unique geology and geography, which have produced an astoundingly diverse and endemic flora and fauna and a fascinating biogeographic history. Much biodiversity research has focused on the regional endemism in the island's Central Core and on its four peninsulas, but the biodiversity of the island's many upland regions is still poorly understood for most taxa, including amphibians and reptiles. Here, we report the first of several planned full-mountain checklists from a series of herpetological surveys of Sulawesi's mountains conducted by our team. In more than 3 weeks of work on Gunung Galang, a 2,254 m peak west of the city of Tolitoli, Sulawesi Tengah Province, on Sulawesi's Northern Peninsula, we recovered nearly fifty species of reptiles and amphibians, more than a dozen of which are either new to science or known but undescribed. The incompleteness of our sampling suggests that many more species remain to be discovered on and around this mountain.

Biogeography: The origin and spread of bee lineages

Where and when bees originated and how they dispersed and diversified across ancient continents has remained ambiguous. A new study that combines phylogenetics with fossil data reconstructs the origin and diversification of bees across geological time and space.

Understanding the exchange process between ground and surface water using mini drive point piezometer and mathematical models to identify suitable managed aquifer recharge sites

In the current study, subsurface characteristics within the complex formation of the Shilabati basin system of West Bengal, India, extending over an area of 3888 km2, have been estimated using a cost-effective piezometer and MIKE FEFLOW package based on a steady-state numerical model. Pore size and fine particle content of streambeds are affected by two opposing flow contraptions. Such opposite flow conditions are likely to affect the hydraulic conductivity of the streambed. However, analogies of the hydraulic conductivity (Kh) of streambeds for losing and gaining streams have not been well documented in the recent past. The Kh value from the piezometer has been highest at the Dakshin Pairachali site (6.765 m/day), with the stream gaining water from the discharge of the local aquifer. Analysis of the stream-aquifer interaction using the FEFLOW model has allowed us to understand the groundwater water head of the basin ranging from 160.33 to 0.32 m.a.s.l (meters above sea level). The present study also constitutes the first attempt for the identification of suitable sites for the implementation of managed aquifer recharge (MAR) technology in West Bengal, India, to manage extreme drought events. The suitable sites have been identified by means of three fuzzy multi-criteria decision analysis based on nine criteria: river discharge, moisture content, porosity, drainage type, rainfall, land use type, geology, aquifer material, and hydraulic conductivity. To design a radial collector well and infiltration gallery for the selected site in an anisotropic, homogeneous, unconfined, and semi-infinite aquifer near a fully penetrating stream, a pumping test has been conducted to optimize a safe yield of 12.096 MLD (megaliters per day).

Enhanced inner core fine-scale heterogeneity towards Earth's centre

Earth's inner core acquires texture as it solidifies within the fluid outer core. The size, shape and orientation of the mostly iron grains making up the texture record the growth of the inner core and may evolve over geologic time in response to geodynamical forces and torques1. Seismic waves from earthquakes can be used to image the texture, or fabric, of the inner core and gain insight into the history and evolution of Earth's core2-6. Here, we observe and model seismic energy backscattered from the fine-scale (less than 10 km) heterogeneities7 that constitute inner core fabric at larger scales. We use a novel dataset created from a global array of small-aperture seismic arrays-designed to detect tiny signals from underground nuclear explosions-to create a three-dimensional model of inner core fine-scale heterogeneity. Our model shows that inner core scattering is ubiquitous, existing across all sampled longitudes and latitudes, and that it substantially increases in strength 500-800 km beneath the inner core boundary. The enhanced scattering in the deeper inner core is compatible with an era of rapid growth following delayed nucleation.

Rift-induced disruption of cratonic keels drives kimberlite volcanism

Kimberlites are volatile-rich, occasionally diamond-bearing magmas that have erupted explosively at Earth's surface in the geologic past1-3. These enigmatic magmas, originating from depths exceeding 150 km in Earth's mantle1, occur in stable cratons and in pulses broadly synchronous with supercontinent cyclicity4. Whether their mobilization is driven by mantle plumes5 or by mechanical weakening of cratonic lithosphere4,6 remains unclear. Here we show that most kimberlites spanning the past billion years erupted about 30 million years (Myr) after continental breakup, suggesting an association with rifting processes. Our dynamical and analytical models show that physically steep lithosphere-asthenosphere boundaries (LABs) formed during rifting generate convective instabilities in the asthenosphere that slowly migrate many hundreds to thousands of kilometres inboard of rift zones. These instabilities endure many tens of millions of years after continental breakup and destabilize the basal tens of kilometres of the cratonic lithosphere, or keel. Displaced keel is replaced by a hot, upwelling mixture of asthenosphere and recycled volatile-rich keel in the return flow, causing decompressional partial melting. Our calculations show that this process can generate small-volume, low-degree, volatile-rich melts, closely matching the characteristics expected of kimberlites1-3. Together, these results provide a quantitative and mechanistic link between kimberlite episodicity and supercontinent cycles through progressive disruption of cratonic keels.

Climatic stability and geological history shape global centers of neo- and paleoendemism in seed plants

Assessing the distribution of geographically restricted and evolutionarily unique species and their underlying drivers is key to understanding biogeographical processes and critical for global conservation prioritization. Here, we quantified the geographic distribution and drivers of phylogenetic endemism for ~320,000 seed plants worldwide and identified centers and drivers of evolutionarily young (neoendemism) and evolutionarily old endemism (paleoendemism). Tropical and subtropical islands as well as tropical mountain regions displayed the world's highest phylogenetic endemism. Most tropical rainforest regions emerged as centers of paleoendemism, while most Mediterranean-climate regions showed high neoendemism. Centers where high neo- and paleoendemism coincide emerged on some oceanic and continental fragment islands, in Mediterranean-climate regions and parts of the Irano-Turanian floristic region. Global variation in phylogenetic endemism was well explained by a combination of past and present environmental factors (79.8 to 87.7% of variance explained) and most strongly related to environmental heterogeneity. Also, warm and wet climates, geographic isolation, and long-term climatic stability emerged as key drivers of phylogenetic endemism. Neo- and paleoendemism were jointly explained by climatic and geological history. Long-term climatic stability promoted the persistence of paleoendemics, while the isolation of oceanic islands and their unique geological histories promoted neoendemism. Mountainous regions promoted both neo- and paleoendemism, reflecting both diversification and persistence over time. Our study provides insights into the evolutionary underpinnings of biogeographical patterns in seed plants and identifies the areas on Earth with the highest evolutionary and biogeographical uniqueness-key information for setting global conservation priorities.

Gravel grabs: The rocky foundations of Indigenous geologic power in the Arctic

Infrastructure development cannot take place without gravel, which is scarce in the North American Arctic. Conditioning where development can occur, the commodity has become the target of Indigenous actors seeking to secure land and resource bases and their material futures, too. In Alaska, decades of litigation pitting Indigenous surface versus subsurface corporate landholders has contested gravel's legal location. In Canada, contrastingly, Inuvialuit land claims negotiators successfully secured access to granular resources. In both locales, legal processes have resulted in certain Indigenous actors' accumulation of geologic power. Rooted in the subterranean, this power enables them to transform the surface of the Earth. Contributing to research on geologic power and political geology and drawing on fieldwork and a review of court cases, policy documents and reports, this article critiques how gravel has become an Arctic resource lucrative to local communities rather than global markets and a key source of Indigenous political and economic agency. Going forward, struggles over Indigenous rights may concern securing ownership over not only the land base, but the land column.

Estimating geothermal and background radiation hotspots from primordial radionuclide concentrations in geology of South Africa

Naturally occurring radionuclides are the main generator of geothermal energy in the Earth's crust and mantle. The generated energy is consequently directly proportional to the concentrations of the three main naturally occurring radionuclides (uranium, thorium and potassium), which are primordial in origin. Concentrations of these naturally occurring radionuclides were extracted for all the different geological rock units in South Africa. The radionuclide concentrations were then mapped and integrated by using QGIS. The results were used to estimate and map the geothermal energy production rates for the rock units. The radionuclide concentrations in the rock units were also used to identify regions with high radiation background. These radiation hotspots were plotted and investigated. The estimated geothermal energy and background radiation hotspots were compared to measurements and projections of other studies and good corelations were found.

Distributed Thermal Response Multi-Source Modeling to Evaluate Heterogeneous Subsurface Properties

A thorough assessment of thermal properties in heterogeneous subsurface is necessary in design of low-temperature borehole heat exchangers (BHEs). A distributed thermal response test (DTRT), which combines distributed temperature sensing (DTS) with a conventional thermal response test (TRT), was conducted in a U-bend geothermal loop installed in an open borehole at the University of Illinois at Urbana-Champaign to estimate thermal properties by analyzing the thermal response of different geologic materials while applying a constant heat input rate. Fiber-optic cables in the DTRT were deployed both inside the U-bend geothermal loop and in the center of the borehole to improve the accuracy of calculated heat-loss rates and borehole temperature profile measurements. To assess the subsurface thermal conductivity during the heating phase of the DTRT, a single-source model and a multi-source model, both based on the infinite line source method, were developed using the borehole temperature data and temperatures inside and along the outside of the loop, separately. The two models returned similar thermal conductivity values. The multi-source modeling has the advantage of predicting the thermal conductivity of heterogeneous geologic materials from borehole temperature profiles during the DTRT heating phase. In addition, based on the distributed thermal conductivity measured in the borehole, estimates were made for both radial thermal impacts and the rate of heat loss in the BHE.

Prediction of elevated groundwater fluoride across India using multi-model approach: insights on the influence of geologic and environmental factors

Elevated fluoride in groundwater is a severe problem in India due to its extensive occurrence and detrimental health impacts on the large population that thrives on groundwater. Although fluoride is primarily a geogenic pollutant, existing model-based studies lack the amalgamation of the influence of geologic factors, specifically tectonics, for identifying groundwater fluoride distribution. This drawback encourages the present study to investigate the association of the tectonic framework with fluoride in a multi-model approach. We have applied three machine learning models (random forest, boosted regression tree, and logistic regression) to predict elevated groundwater fluoride based on fluoride measurements across India. The random forest model outperformed other models with an accuracy of 93%. Tectonics was found to be one of the most important predictors alongside "depth to water table." Two major areas of high risk identified were the northwest parts and the south-southeast cratonic peninsular region. The random forest model also performed significantly well over the validation dataset. We estimate that nearly 257 million people are exposed to elevated fluoride risk in India. We endeavor that the findings of our study would be an effective tool for identifying the areas at risk of elevated fluoride and also assist in undertaking effective groundwater management strategies.