Search

Abstract Subterranean ecosystems (comprising terrestrial, semi-aquatic, and aquatic components) are increasingly threatened by human activities; however, the current network of surface-protected areas is inadequate to safeguard subterranean biodiversity. Establishing protected areas for subterranean ecosystems is challenging. First, there are technical obstacles in mapping three-dimensional ecosystems with uncertain boundaries. Second, the rarity and endemism of subterranean organisms, combined with a scarcity of taxonomists, delays the accumulation of essential biodiversity knowledge. Third, establishing agreements to preserve subterranean ecosystems requires collaboration among multiple actors with often competing interests. This perspective addresses the challenges of preserving subterranean biodiversity through protected areas. Even in the face of uncertainties, we suggest it is both timely and critical to assess general criteria for subterranean biodiversity protection and implement them based on precautionary principles. To this end, we examine the current status of European protected areas and discuss solutions to improve their coverage of subterranean ecosystems.

This study investigates the long-term trends (1992–2022) of nitrogen and phosphorus loadings exported by the River Po to the Adriatic Sea, to better analyse how changes in hydrology are affecting the timing and magnitude of river nutrient loads. We used 30 years of monitoring data in order to (a) identify the main temporal patterns and their interactions at a decadal, annual and seasonal scale, (b) estimate precipitation effects on load formation and evaluate whether and to which extent the hydrological regime affects nutrient export across the years and (c) analyse the nutrient export regime at a monthly scale and the main transport dynamic of N and P chemical species (hydrological vs. biogeochemical control). The long-term analysis shows a general decrease of both P and N loadings, but the trends are different between the elements and their chemical species, as well as undergoing different seasonal variations. We found a statistically significant relationships between precipitation and loads, which demonstrates that precipitation patterns drive the exported load at the intra- and interannual time scales considered in this study. Precipitation-induced load trends trigger seasonal changes in nutrient deliveries to the sea, peaking in spring and autumn. The nitrogen decrease is mainly concentrated in the summer dry period, while total phosphorus diminishes mainly in spring and autumn. This mismatch of N and P results in variable molar N:P ratios within the year. The effects of extreme drought and flood events, along with the progressive decrease in the snowmelt contribution to water fluxes, are expected to exacerbate the variability in the N and P loadings, which in turn is expected to perturbate the biodiversity, food webs and trophic state of the Northern Adriatic Sea.

Bundera sinkhole, located in north-western Australia, is the only known continental anchialine system in the Southern Hemisphere. Anchialine environments are characterised by stratified water columns with complex physicochemical profiles spanning hypoxic and anoxic regions, often displaying high levels of endemism. Research on these systems has focused on eukaryotic fauna, however interest in the microbial diversity of these environments is growing, enabled by next-generation DNA sequencing. Here we report detailed analyses of the microbial communities across a depth profile within Bundera sinkhole (from 2 to 28 m), involving parallel physicochemical measurements, cell population counts and 16S rRNA amplicon analyses. We observed clear shifts in microbial cell counts, community diversity, structure and membership across the depth profile, reflecting changing levels of light, organic and inorganic energy sources as well as shifts in pH and salinity. While Proteobacteria were the most abundant phylum found, there was a high degree of taxonomic novelty within these microbial communities, with 13,028 unique amplicon sequence variants (ASVs) identified, belonging to 67 identifiable bacterial and archaeal phyla. Of these ~4,600, more than one third of the total, were unclassified below family level. A small number of ASVs were highly abundant at select depths, all of which were part of the set not classified below family level. The 2 m and 6 m samples had in common two highly abundant ASVs, belonging to the Ectothiorhodospiraceae and Thiotrichaceae families, while the 8 m community contained a single predominant ASV belonging to family Thioglobaceae. At lower depths a different Ectothiorhodospiraceae ASV comprised up to 68% relative abundance, peaking at 26 and 28 m. Canonical correspondence analyses indicated that community structure was strongly influenced by differences in key physicochemical parameters, particularly salinity, dissolved organic and inorganic carbon, phosphate and sulphate concentrations. This work highlights the potential for anchialine systems to house considerable microbial novelty, potentially driven by adaptations to the specific physicochemical makeup of their local environment. As only a small number of anchialine systems have been examined via microbial community studies to date, this work is particularly valuable, contributing new insight regarding the microbial residents of these important and sensitive environments.

Abstract Groundwaters host vital resources playing a key role in the near future. Subterranean fauna and microbes are crucial in regulating organic cycles in environments characterized by low energy and scarce carbon availability. However, our knowledge about the functioning of groundwater ecosystems is limited, despite being increasingly exposed to anthropic impacts and climate change-related processes. In this work we apply novel biochemical and genetic techniques to investigate the ecological dynamics of an Australian calcrete under two contrasting rainfall periods (LR—low rainfall and HR—high rainfall). Our results indicate that the microbial gut community of copepods and amphipods experienced a shift in taxonomic diversity and predicted organic functional metabolic pathways during HR. The HR regime triggered a cascade effect driven by microbes (OM processors) and exploited by copepods and amphipods (primary and secondary consumers), which was finally transferred to the aquatic beetles (top predators). Our findings highlight that rainfall triggers ecological shifts towards more deterministic dynamics, revealing a complex web of interactions in seemingly simple environmental settings. Here we show how a combined isotopic-molecular approach can untangle the mechanisms shaping a calcrete community. This design will help manage and preserve one of the most vital but underrated ecosystems worldwide.

Groundwater environments interact with and support subterranean biota as well as superficial aquatic and terrestrial ecosystems. However, knowledge of subterranean energy flows remains incomplete. Cross-boundary investigations are needed to better understand the trophic structures of groundwater ecosystems and their reliance on carbon inputs from aboveground. In this study we used carbon and nitrogen stable isotope analyses combined with radiocarbon fingerprints to characterise organic flows in groundwater ecosystems. We coupled these data with DNA metabarcoding of the gut contents of stygofauna to further elucidate organic matter (OM) sources and shifts in diet preferences. Samples were collected from the arid zone Sturt Meadows calcrete aquifer under low rainfall (LR) and high rainfall (HR) conditions. Bayesian modelling of Δ14C, δ13C, and δ15N data indicated that primary consumers (copepods) incorporated mainly particulate organic carbon (POC) under LR but during HR shifted to root derived material (either exudates or direct root grazing). By contrast, diets of secondary consumers (amphipods) were dominated by root material under both LR and HR. Our DNA metabarcoding-based results indicate that amphipods relied primarily on root inputs from perennial trees (likely Eucalyptus and Callitris) during the dry season (LR). Under HR, diets of both amphipods and copepods also included organic material derived from a broad range of more shallow rooted shrubs, and ephemeral herbs and grasses. Our findings illustrate the complexity of functional linkages between groundwater biota and surface terrestrial ecosystems in environments where aboveground productivity, diversity and OM flux to groundwater are intimately linked to often episodic rainfall.

Abstract Food web dynamics are vital in shaping the functional ecology of ecosystems. However, trophic ecology is still in its infancy in groundwater ecosystems due to the cryptic nature of these environments. To unravel trophic interactions between subterranean biota, we applied an interdisciplinary Bayesian mixing model design (multi‐factor BMM) based on the integration of faunal C and N bulk tissue stable isotope data (δ13C and δ15N) with radiocarbon data (Δ14C), and prior information from metagenomic analyses. We further compared outcomes from multi‐factor BMM with a conventional isotope double proxy mixing model (SIA BMM), triple proxy (δ13C, δ15N, and Δ14C, multi‐proxy BMM), and double proxy combined with DNA prior information (SIA + DNA BMM) designs. Three species of subterranean beetles (Paroster macrosturtensis, Paroster mesosturtensis, and Paroster microsturtensis) and their main prey items Chiltoniidae amphipods (AM1: Scutachiltonia axfordi and AM2: Yilgarniella sturtensis), cyclopoids and harpacticoids from a calcrete in Western Australia were targeted. Diet estimations from stable isotope only models (SIA BMM) indicated homogeneous patterns with modest preferences for amphipods as prey items. Multi‐proxy BMM suggested increased—and species‐specific—predatory pressures on amphipods coupled with high rates of scavenging/predation on sister species. SIA + DNA BMM showed marked preferences for amphipods AM1 and AM2, and reduced interspecific scavenging/predation on Paroster species. Multi‐factorial BMM revealed the most precise estimations (lower overall SD and very marginal beetles' interspecific interactions), indicating consistent preferences for amphipods AM1 in all the beetles' diets. Incorporation of genetic priors allowed crucial refining of the feeding preferences, while integration of more expensive radiocarbon data as a third proxy (when combined with genetic data) produced more precise outcomes but close dietary reconstruction to that from SIA + DNA BMM. Further multidisciplinary modeling from other groundwater environments will help elucidate the potential behind these designs and bring light to the feeding ecology of one the most vital ecosystems worldwide.

Thermal tolerance limits in animals are often thought to be related to temperature and thermal variation in their environment. Recently, there has been a focus on studying upper thermal limits due to the likelihood for climate change to expose more animals to higher temperatures and potentially extinction. Organisms living in underground environments experience reduced temperatures and thermal variation in comparison to species living in surface habitats, but how these impact their thermal tolerance limits are unclear. In this study, we compare the thermal critical maximum (CTmax) of two subterranean diving beetles (Dytiscidae) to that of three related surface-dwelling species. Our results show that subterranean species have a lower CTmax (38.3-39.0°C) than surface species (42.0-44.5°C). The CTmax of subterranean species is ∼10°C higher than the highest temperature recorded within the aquifer. Groundwater temperature varied between 18.4°C and 28.8°C, and changes with time, depth and distance across the aquifer. Seasonal temperature fluctuations were 0.5°C at a single point, with the maximum heating rate being ∼1000x lower (0.008°C/hour) than that recorded in surface habitats (7.98°C/hour). For surface species, CTmax was 7-10°C higher than the maximum temperature in their habitats, with daily fluctuations from ∼1°C to 16°C and extremes of 6.9°C and 34.9°C. These findings suggest that subterranean dytiscid beetles are unlikely to reach their CTmax with a predicted warming of 1.3-5.1°C in the region by 2090. However, the impacts of long-term elevated temperatures on fitness, different life stages and other species in the beetle's trophic food web are unknown.

Freshwater ecosystems play a key role in shaping the global carbon cycle and maintaining the ecological balance that sustains biodiversity worldwide. Surficial water bodies are often interconnected with groundwater, forming a physical continuum, and their interaction has been reported as a crucial driver for organic matter (OM) inputs in groundwater systems. However, despite the growing concerns related to increasing anthropogenic pressure and effects of global change to groundwater environments, our understanding of the dynamics regulating subterranean carbon flows is still sparse. We traced carbon composition and transformations in an arid zone calcrete aquifer using a novel multidisciplinary approach that combined isotopic analyses of dissolved organic carbon (DOC) and inorganic carbon (DIC) (δ13CDOC, δ13CDIC, 14CDOC and 14CDIC) with fluorescence spectroscopy (Chromophoric Dissolved OM (CDOM) characterisation) and metabarcoding analyses (taxonomic and functional genomics on bacterial 16S rRNA). To compare dynamics linked to potential aquifer recharge processes, water samples were collected from two boreholes under contrasting rainfall: low rainfall ((LR), dry season) and high rainfall ((HR), wet season). Our isotopic results indicate limited changes and dominance of modern terrestrial carbon in the upper part (northeast) of the bore field, but correlation between HR and increased old and 13C-enriched DOC in the lower area (southwest). CDOM results show a shift from terrestrially to microbially derived compounds after rainfall in the same lower field bore, which was also sampled for microbial genetics. Functional genomic results showed increased genes coding for degradative pathways-dominated by those related to aromatic compound metabolisms-during HR. Our results indicate that rainfall leads to different responses in different parts of the bore field, with an increase in old carbon sources and microbial processing in the lower part of the field. We hypothesise that this may be due to increasing salinity, either due to mobilisation of Cl- from the soil, or infiltration from the downstream salt lake during HR. This study is the first to use a multi-technique assessment using stable and radioactive isotopes together with functional genomics to probe the principal organic biogeochemical pathways regulating an arid zone calcrete system. Further investigations involving extensive sampling from diverse groundwater ecosystems will allow better understanding of the microbiological pathways sustaining the ecological functioning of subterranean biota.

Hypersaline ecosystems—aquatic environments where concentration of salt exceeds 35 g L−1—host microbial communities that are highly specialised to cope with these extreme conditions. However, our knowledge on the taxonomic diversity and functional metabolisms characterising microbial communities in the water columns of hypersaline ecosystems is still limited, and this may compromise the future preservation of these unique environments. DNA metabarcoding provides a reliable and affordable tool to investigate environmental dynamics of aquatic ecosystems, and its use in brine can be highly informative. Here, we make use of bacterial 16S metabarcoding techniques combined with hydrochemical analyses to investigate the microbial patterns (diversity and functions) from five hypersaline lakes located at Rottnest Island (WA). Our results indicate lake-driven microbial aquatic assemblages that are characterised by taxonomically and functionally moderately to extremely halophilic groups, with TDS (total dissolved solids) and alkalinity amongst the most influential parameters driving the community patterns. Overall, our findings suggest that DNA metabarcoding allows rapid but reliable ecological assessment of the hypersaline aquatic microbial communities at Rottnest Island. Further studies involving different hypersaline lakes across multiple seasons will help elucidate the full extent of the potential of this tool in brine.

Subterranean ecosystems host highly adapted aquatic invertebrate biota which play a key role in sustaining groundwater ecological functioning and hydrological dynamics. However, functional biodiversity studies in groundwater environments, the main source of unfrozen freshwater on Earth, are scarce, probably due to the cryptic nature of the systems. To address this, we investigate groundwater trophic ecology via stable isotope analysis, employing δ13C and δ15N in bulk tissues, and amino acids. Specimens were collected from a shallow calcrete aquifer in the arid Yilgarn region of Western Australia: a well-known hot-spot for stygofaunal biodiversity. Sampling campaigns were carried out during dry (low rainfall: LR) and the wet (high rainfall: HR) periods. δ13C values indicate that most of the stygofauna shifted towards more 13C-depleted carbon sources under HR, suggesting a preference for fresher organic matter. Conversion of δ15N values in glutamic acid and phenylalanine to a trophic index showed broadly stable trophic levels with organisms clustering as low-level secondary consumers. However, mixing models indicate that HR conditions trigger changes in dietary preferences, with increasing predation of amphipods by beetle larvae. Overall, stygofauna showed a tendency towards opportunistic and omnivorous habits-typical of an ecologically tolerant community-shaped by bottom-up controls linked with changes in carbon flows. This study provides baseline biochemical and ecological data for stygofaunal trophic interactions in calcretes. Further studies on the carbon inputs and taxa-specific physiology will help refine the interpretation of the energy flows shaping biodiversity in groundwaters. This will aid understanding of groundwater ecosystem functioning and allow modelling of the impact of future climate change factors such as aridification.

Abstract Ammonia is a crucial component in the production of fertilizers and various nitrogen-based compounds. Now, the lithium-mediated nitrogen reduction reaction (Li-NRR) has emerged as a promising approach for ammonia synthesis at ambient conditions. The proton shuttle plays a critical role in the proton transfer process during Li-NRR. However, the structure-activity relationship and design principles for effective proton shuttles have not yet been established in practical Li-NRR systems. Here, we propose a general procedure for verifying a true proton shuttle and established design principles for effective proton shuttles. We systematically evaluate several classes of proton shuttles in a continuous-flow reactor with hydrogen oxidation at the anode. Among the tested proton shuttles, phenol exhibits the highest Faradaic efficiency of 72 ± 3% towards ammonia, surpassing that of ethanol, which has been commonly used so far. Experimental investigations including operando isotope-labelled mass spectrometry proved the proton-shuttling capability of phenol. Further mass transport modeling sheds light on the mechanism.