Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reservoir operation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 685 records · Page 38
All sourcesNASA NTRSDOE OSTI
NASA NTRS2003

Zeolite Formation and Weathering Processes Within the Martian Regolith: An Antarctic Analog

As more information is obtained about the nature of the surface compositions and processes operating on Mars, it is clear that significant erosional and depositional features are present on the surface. Apparent aqueous or other fluid activity on Mars has produced many of the erosional and outflow features observed. Evidence of aqueous activity on Mars has been reported by earlier studies. Gooding and colleagues championed the cause of pre-terrestrial aqueous alteration processes recorded in Martian meteorites. Oxygen isotope studies on Martian meteorites by Karlsson et al. and Romenek et al. gave evidence for two separate water reservoirs on Mars. The oxygen isotopic compositions of the host silicate minerals was different from the oxygen isotopic composition of the secondary alteration products within the SNC meteorites. This implied that the oxygen associated with fluids which produced the secondary alteration was from volatiles which were possibly added to the planetary inventory after formation of the primary silicates from which the SNC s were formed. The source of the oxygen may have been from a cometary or volatile-rich veneer added to the planet in its first 600 million years.

Gibson, E. K.↗
DOE OSTI2022

Impact of Organic and Volatile Compounds in Produced Water from Unconventional Reservoirs on Direct Contact Membrane Distillation Permeate Quality

The expansion of oil and gas extraction from unconventional reservoirs has led to an increase in the amount of produced water that has to be managed by this industry. Direct contact membrane distillation (DCMD) is a promising technology for treatment of produced water to enable water recovery and reduce the environmental footprint of this industry. The feasibility of DCMD for the treatment of highly saline produced water from the Permian Basin in TX with commercially available polyethylene and polytetrafluoroethylene membranes was evaluated in this study. An increase in water recovery by a DCMD system operated in the batch (concentrating) mode led to an increase in permeate conductivity. Partial removal of organic compounds from the produced water by biodegradation, chemical oxidation, and/or activated carbon adsorption could not resolve deterioration in permeate quality, and none of the organics observed in the permeate contributed to its conductivity. The observed increase in permeate conductivity was attributed to the passage of ammonia vapor from the feed side followed by protonation on the permeate side. Furthermore, this study revealed that permeate conductivity may not always be a reliable indicator of membrane wetting and underscores the importance of understanding the interactions between specific solutes and membrane materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

High Temperature Ternary Chloride Molten Salt Pump

A ternary chloride (MgCl2-NaCl-KCl) molten salt pump has been developed to operate up to 720°C. The design is the result of a thermal, structural and material study of critical components. The vertical overhung pump incorporates first of a kind salt-wetted hydrodynamic bearings. Several bearing material pairings were tested which include: colmonoy grades, stellite grades, and a novel NiWC-based alloy (i.e. HybrimetTM NiWC3b) used as a High Velocity Oxygen Fuel (HVOF) coating. The coating withstands the corrosion effects of the molten salt, and is also used to protect a 316L stainless steel pump reservoir. The pump bowl assembly components are made of Inconel 625. A Nitrogen flushed carbon ring seal is used to keep salt vapors from interacting with external ambient air. Preliminary testing shows a promising pump design for use in the Concentrating Solar Power Gen 3 systems and sets the ground for further development of the HybrimetTM NiWC3b as a competitive alternative to available superalloys.

Handy-Cardenas, Lewis (ORCID:0009000744006972)↗
NASA NTRS1979

A fundamental study of liquid spray evaporation at a heated surface

The influence of feedwater temperature and the influence of pulsing nozzle flow on the heat transfer characteristics in the spray evaporative cooling process were investigated. A commercial water circulation heater was installed in the feedwater line between the supply reservoir and the nozzle. A pass through the heater provided water at temperatures up to 50 C at the nozzle; higher temperatures were obtained by multiple heater passes. A thermocouple installed in the water line just upstream of the nozzle provided a measure of the feedwater temperature. Initial tests were run with belljar exhaust orifice number two which provided a pressure in a range about 7mm Hg abs. The feedwater temperature was heated to a temperature of 45 C. Results indicate that for a given system operating pressure, there is an upper limit on the feedwater temperature above which the spray will flash to vapor before hitting the surface to be cooled.

Wierum, F. A.↗
DOE OSTI2020

Coupled Time-Lapse Full-Waveform Inversion for Subsurface Flow Problems Using Intrusive Automatic Differentiation

We describe a novel framework for estimating subsurface properties, such as rock permeability and porosity, from time-lapse observed seismic data by coupling full-waveform inversion (FWI), subsurface flow processes, and rock physics models. For the inverse modeling, we handle the back propagation of gradients by an intrusive automatic differentiation strategy that offers three levels of user control: (1) At the wave physics level, we adopted the discrete adjoint method in order to use our existing high-performance FWI code; (2) at the rock physics level, we used built-in automatic differentiation operators from the TensorFlow backend; (3) at the flow physics level, we implemented customized partial differential equation (PDE) operators for the multiphase flow equations. The three-level coupled inversion strategy strikes a good balance between computational efficiency and programming efforts, and when the gradients are chained together, it constitutes a coupled inverse system. Our numerical experiments demonstrate that the three-level coupled inverse problem is superior in terms of accuracy to a traditional decoupled inversion strategy. Additionally, our method is able to simultaneously invert for parameters in empirical relationships such as the rock physics models. Our proposed inverted model can be used for reservoir performance prediction and reservoir management/optimization purposes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Utah FORGE 1-2409: Zonal Isolation Solution for Geothermal Wells - Workshop Presentation

This is a presentation on the Zonal Isolation Solution for Geothermal Wells project by PetroQuip Energy Services, presented by VP of operations Robert Coon. The project's objective was to design and develop a multi-stage system for zonally isolating fluids inside and outside of geothermal well casings. This presentation was featured in the Utah FORGE R&D Annual Workshop on September 7, 2023. The workshop provided a valuable opportunity to explore the progress made in each of the 17 Research and Development projects funded under Solicitation 2020-1 which aim to enhance our understanding of the crucial factors influencing the development of Enhanced Geothermal Systems (EGS) reservoirs and resources.

15 GEOTHERMAL ENERGY↗
DOE OSTI2022

Computational Tools and Workflows for Quantitative Risk Assessment and Decision Support for Geologic Carbon Storage Sites: Progress and Insights from the U.S. DOE’s National Risk Assessment Partnership

The 2005 Intergovernmental Panel on Climate Change (IPCC) Special Report on CCS raised the profile of CO2 capture and storage (CCS) as an important technology for reducing greenhouse gas (GHG) emissions. CCS is now recognized as a key component of most climate change mitigation scenarios. Since publication of that report the international research, development, and deployment (RD&D) community has advanced key technical aspects, clarified regulatory requirements, explored value chain and infrastructure solutions, and developed incentive paradigms to enable and promote large-scale deployment of CCS. These efforts have included research to better characterize geologic storage resources, to improve injection performance and storage efficiency, to assess and manage subsurface environmental risks, and to advance monitoring technologies to assure system conformance. These efforts have helped to build confidence in the viability of geologic carbon storage (GCS), but stakeholder concerns about long-term risks and liability associated with GCS remain a hurdle to broad acceptance and large-scale deployment of CCS. Since 2010, the U.S. DOE’s National Risk Assessment Partnership (NRAP) – a research collaboration between five contributing national laboratories – has worked to establish and demonstrate methods and tools to quantify and manage the subsurface environmental risks associated with GCS, amidst uncertainty. This work supports the Office of Fossil Energy and Carbon Management Carbon Transport and Storage Program’s goal of advancing safe and secure commercial-scale GCS deployment. To address the technical challenge of simulating the physical response of the GCS site to large-scale CO2 injection, NRAP has adopted an approach that relies on coupling computationally efficient reduced-order and/or data-driven proxy models of important system components (i.e., storage reservoir, sealing caprock, leakage pathways, intermediate formations, overlying groundwater aquifers, and the atmosphere) in integrated assessment framework. That integrated model of the physical system is complemented with fit-for purpose functionality to support site characterization and risk-related decisions. The recently released NRAP Phase II toolset includes the Open-Source Integrated Assessment Model (NRAP-Open-IAM) for evaluation of trends in leakage risk and potential impact, tools to support monitoring design optimization (Designs for Risk Evaluation and Management – DREAM v3.0 and Passive Seismic Monitoring Tool - PSMT), and tools for state of stress evaluation (State-of-Stress Analysis Tool - SOSAT) and forecasting induced seismicity risk. The NRAP team has also released a pair of reports describing conceptual workflows to incorporate physics-based, quantitative risk assessment into many of the design, planning, operation, and closure decisions for GCS projects. An online catalogue highlights published studies where these tools and methods are demonstrated. In this presentation, the utility of these products to assess risks and address key stakeholder questions will be highlighted through examples, and related insights about the safety and security of geologic carbon storage in qualified storage sites will be discussed. The prospect of rapid, large-scale deployment of GCS technology to aggressively reduce anthropogenic CO2 emissions requires careful consideration of interference between multiple commercial-scale storage projects within a basin. Going forward, NRAP is expanding and adapting site-scale risk quantification tools and methods to enable assessment of risks and inform management decisions for basin-scale deployment. Increasingly, this work will leverage next-generation approaches for surrogate modelling, fast prediction, and advanced visualization enabled by machine learning and artificial intelligence to promote virtual learning, scenario evaluation, and augment risk-based decision making.

quantitative risk assessment, geologic carbon stor↗
NASA NTRS2000

Integrated Global Observation Strategy - Ozone and Atmospheric Chemistry Project

The "Long Term Continuity of Stratospheric Ozone Measurements and Atmospheric Chemistry" project was one of six established by the Committee on Earth Observing Satellites (CEOS) in response to the Integrated Global Observing Strategy (IGOS) initiative. IGOS links satellite and ground based systems for global environmental observations. The strategy of this project is to develop a consensus of user requirements including the scientific (SPARC, IGAC, WCRP) and the applications community (WMO, UNEP) and to develop a long-term international plan for ozone and atmospheric chemistry measurements. The major components of the observing system include operational and research (meeting certain criteria) satellite platforms planned by the space faring nations which are integrated with a well supported and sustained ground, aircraft, and balloon measurements program for directed observations as well satellite validation. Highly integrated and continuous measurements of ozone, validation, and reanalysis efforts are essential to meet the international scientific and applications goals. In order to understand ozone trends, climate change, and air quality, it is essential to conduct long term measurements of certain other atmospheric species. These species include key source, radical, and reservoir constituents.

Hilsenrath, Ernest↗
DOE OSTI2026

A Demonstration System for Geological Thermal Energy Storage of Concentrating Solar Thermal

Energy storage is increasingly necessary as variable energy technologies are deployed. Seasonal energy storage can shift energy generation from the summer to the winter, but these technologies must have extremely large energy capacities and low costs. Geological Thermal Energy Storage (GeoTES) is proposed as a solution for long-term energy storage [1]. Excess thermal energy can be stored in permeable reservoirs such as aquifers and depleted hydrocarbon reservoirs for several months. The energy capacity cost of GeoTES is very low which makes it suitable for both daily- and seasonal storage of Concentrating Solar Thermal (CST) energy, thus enabling CST to provide value to electricity markets and thermal energy off-takers. A CST-GeoTES demonstration system has been funded by the U.S. Department of Energy, Solar Energy Technology Office. In this article, we will describe this demonstration system and progress that has been made in its development. The demonstration system will comprise a 2 MWth parabolic trough with an 8m aperture developed by Gossamer Space Frames, seven wells, and a 100 kWe power cycle. The demonstration system will be deployed in Kern County, California by Premier Resource Management. A techno-economic model for CST-GeoTES systems has also been developed [3] and is applied to the demonstration system and its planned future expansion. The model integrates the output of specialist models of each subsystem, which enables the performance and cost of both the subsurface and surface systems to be captured. Off-design models enable the performance to be evaluated at each hour of the year, before being aggregated to evaluate the economic potential of CST-GeoTES. Initial analysis indicates that CST-GeoTES can provide long duration energy storage capabilities with low marginal costs of energy capacity - leading to a low value of Levelized Cost of Storage (LCOS) compared to alternative technologies, see Figure 2. In this article, we investigate the cost and performance of the specific demonstration site being developed by PRM and explore a range of operational profiles that can deliver different value streams, such as daily and seasonal storage, capacity, and resiliancy.

14 SOLAR ENERGY↗
DOE OSTI2024

Advancing spatiotemporal forecasts of CO 2 plume migration using deep learning networks with transfer learning and interpretation analysis

Accurate and timely forecasts of CO 2 plume distribution throughout the injection and post-injection phases are crucial for detecting plume migration, assessing leakage risks, and supporting operational decisions in geologic carbon storage (GCS). Current convolutional neural network-based approaches primarily focus on spatial information and overlook temporal dependencies in plume distributions, thus limiting their ability to capture dynamic movement effects and provide accurate predictions of plume migration. In this work, we propose two deep learning models, Auto-Encoder (AE)-LSTM and Encoder-Decoder (ED)-ConvLSTM, each uniquely designed to capture both spatial and temporal features. We apply the proposed methods to forecast the dynamic distribution of CO 2 plumes based on 108 reservoir simulations over a 30-year injection and a 30-year post-injection period. The results indicate that the ED-ConvLSTM model outperforms the AE-LSTM model in accurately predicting the spatiotemporal dynamics of CO 2 plume migration, achieving R 2 values above 0.99. To provide a deeper understanding of these model predictions, we employ a gradient-based explanation method on the trained models. This approach provides insights into the influence of input variables on plume migration forecasts and uncovers the underlying prediction mechanisms of the proposed models. Furthermore, we introduce a transfer learning technique, enabling fast and accurate plume migration forecasting in the post-injection phase by leveraging the trained model during the injection phase. This reduces the necessity for extensive data collection or re-training. In conclusion, the methods proposed in our work enhances the performance and interpretability of CO 2 plume migration forecasts, thereby facilitating informed decision-making throughout the entire lifecycle of GCS applications.

58 GEOSCIENCES↗
NASA NTRS2022

Mapping Firn Saturation over Greenland using NASA’s Soil Moisture Active Passive Satellite

Mapping the spatial extent of recently identified englacial hydrological features (i.e., ice slabs and perennial firn aquifers) formed by meters-thick water-saturated firn layers over the percolation facies of the Greenland Ice Sheet using L-band microwave radiometry has recently been demonstrated. However, these initial maps are binary, and do not provide a parameter to estimate the spatial variability in the thickness and volumetric fraction of meltwater stored within the firn pore space. Here, we exploit enhanced-resolution vertical-polarization L-band brightness temperature (T^(B)_(V)) imagery (2015–2019) generated using observations collected over Greenland by NASA’s Soil Moisture Active Passive (SMAP) satellite and a simple two-layer L-band brightness temperature model. We map water-saturated firn layers via a “firn saturation” parameter, and interpret our results together with ice slab and perennial firn aquifer spatial extents, estimates of snow accumulation simulated via the Regional Atmospheric Climate Model (RACMOp2.3), and airborne radar surveys collected via NASA’s Operation IceBridge (OIB) campaigns. We find that variable firn saturation parameter values are mapped in lower snow accumulation ice slab areas in western, northern, and northeastern Greenland, where firn is colder and water-saturated firn layers seasonally refreeze as solid-ice. Higher firn saturation parameter values are mapped in higher snow accumulation perennial firn aquifer areas in southeastern, southern, and northwestern Greenland, where firn is near the melting point, and meters-thick water-saturated firn layers exist. Our results have implications for identifying expansive englacial reservoirs that store significant volumes of meltwater in locations that are vulnerable to meltwater-induced hydrofracturing and accelerated outlet glacier flow year-round.

Greenland ice sheet↗
DOE OSTI2020

Analysis of Microseismicity and Reactivated Fault Size to Assess the Potential for Felt Events by CO 2 Injection in the Illinois Basin

The results of monitoring of carbon dioxide (CO 2 ) injection at the Illinois Basin—Decatur Project (IBDP) and the companion Illinois Industrial Carbon Capture and Sequestration Sources (IL-ICCS) project—have shown that reservoir response to fluid pressure changes can vary significantly at different injection locations within the same reservoir. Predrill reservoir characterization is important to identify potentially seismogenic faults. However, interpretations of newly reprocessed 3D seismic reflection data illustrate the challenges related to their identification in a region dominated by faulting with small vertical offsets. Faults interpreted in the 3D seismic volume range from ~300 to 1200 m wide and are in the same size range as faults that could have been the source of historical events up to Mw 2.7 in central Illinois. The array of monitoring sensors that was installed for the IBDP continues to collect data, as injection operates in IL-ICCS, the second injection well. CO 2 injection rates for the IL-ICCS well are on average 1.7 times the rates injected in the IBDP well, but a significantly reduced rate of induced seismicity is observed. This article presents results of passive seismic monitoring for the duration of the project to date, integrating active and passive seismic data to develop a new interpretation of the subsurface structure at the Decatur site that explicitly identifies pathways for fluid flow into the basement leading to induced seismicity, and provides a geological explanation for the sharp reduction of induced seismicity during injection at higher rates into the second well. The use of seismic moment to estimate the length of seismogenic slip planes in the local subsurface suggests that faults large enough to produce felt seismicity are unlikely to be present at or near the Decatur site.

Geochemistry & Geophysics↗
DOE OSTI2020

Optimizing Cornell’s future geothermal district heating performance through systems engineering and simulation

Cornell University’s intention to lower its carbon footprint has motivated this engineering evaluation of using low-temperature geothermal energy to supply heat to the campus district energy system. Optimal selection and operation of heat pumps can significantly improve system performance. To support this analysis, we model the quantitative relationships between heat pump configurations and source flow, supply and distribution temperatures, facility heating design options, output heat generation, and carbon abatement outcomes. A systems approach is used for analysis, troubleshooting, and improvement of the model. The result of this effort is a dynamic tool that appropriately responds to hourly thermal demand and communicates energy response for techno-economic analysis. Simulations indicate that a single well-pair in the local Basement Contact Zone reservoir can satisfy almost 68% of annual thermal demand at nearly 20 MWth average capacity while remaining financially competitive compared to conventional heating with levelized cost of heating (LCOH) as low as $4.55/MMBTU ($15.53/MWh th ). A direct-use scheme with heat pump augmentation provided more carbon abatement and 50 to 150% more annual heat supply than one without augmentation. Exploration of three and four well-pair scenarios prove the capability to achieve 50 MW th baseload heating capacity with competitive LCOH and significant carbon reductions.

15 GEOTHERMAL ENERGY↗
NASA NTRS1987

Infrared spectroscopic measurements of halogenated source gases in the stratosphere with the ATMOS instrument

The volume mixing ratios of the six most important halogenated source species (CH3Cl, CF2Cl2, CFCl3, CHF2Cl, CCl4, and CF4) have been retrieved over the 10- to 30-km altitude range from the analysis of 0.01/cm resolution infrared solar occultation spectra recorded near 30 deg N and 47 deg S latitudes with the ATMOS (Atmospheric Trace Molecule Spectroscopy) instrument, operating from on board Spacelab 3 (April-May 1985). The results reported here, although in satisfactory agreement with recent in situ values obtained from air sampling techniques, are limited in accuracy by the limited absorption representative of most of the species and by uncertainties in the spectroscopic parameters currently available for these gases. They demonstrate, however, the power of the IR remote sensing approach for evaluating on a global scale the total chlorine budget of the atmosphere, and they provide an independent set of simultaneous data acquired for the important source and reservoir halogenated molecular species in the upper atmosphere.

Zander, R.↗
DOE OSTI2023

Autonomous quantum error correction and fault-tolerant quantum computation with squeezed cat qubits

We propose an autonomous quantum error correction scheme using squeezed cat (SC) code against excitation loss in continuous-variable systems. Through reservoir engineering, we show that a structured dissipation can stabilize a two-component SC while autonomously correcting the errors. The implementation of such dissipation only requires low-order nonlinear couplings among three bosonic modes or between a bosonic mode and a qutrit. While our proposed scheme is device independent, it is readily implementable with current experimental platforms such as superconducting circuits and trapped-ion systems. Compared to the stabilized cat, the stabilized SC has a much lower dominant error rate and a significantly enhanced noise bias. Furthermore, the bias-preserving operations for the SC have much lower error rates. In combination, the stabilized SC leads to substantially better logical performance when concatenating with an outer discrete-variable code. The surface-SC scheme achieves more than one order of magnitude increase in the threshold ratio between the loss rate κ 1 and the engineered dissipation rate κ 2 . Under a practical noise ratio κ 1 /κ 2 = 10 -3 , the repetition-SC scheme can reach a 10 -15 logical error rate even with a small mean excitation number of 4, which already suffices for practically useful quantum algorithms.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1986

Apparatus for detecting at least one predetermined condition and providing an informational signal in response thereto in a medication infusion system

A medication infusion system provides redundant safety and includes condition detecting and informational alarm signal generating apparatus for indicating if (1) a fluid leak occurs in different portions of the system; (2) a programmable input from a patient or physician would result in exceeding a safe dosage limit; (3) the reservoir containing medication has been filled; (4) the intended medication pumping does not correlate with the pumping actually effected; (5) battery voltage is low; (6) the medication reserve is low; and (7) the system has been switched off. The apparatus may provide subcutaneous electrical, thermal, or audible stimulation to the patient and also provides a signal which a physician may monitor. The stimulation may be coded to separately identify each above-listed deviation in nominal system performance. In addition, the number of medication requests are correlated with actual medication dispensing to assure proper operation. An identification scheme is provided which matches the patient with his or her corresponding medication.

Fischell, Robert E.↗
DOE OSTI2020

Advanced supercritical water-based process concepts for treatment and beneficial reuse of brine in oil/gas production

Produced water generated from oil/gas reservoirs is a cost burden for oil/gas producers, with over 21 billion barrels of the waste generated in 2019. Average disposal costs range from 4.00–8.00·bbl -1 , with up to 20·bbl -1 spent on transportation alone. Dissolved solids content is a contributing factor which limits economical treatment options for this waste stream. This work proposes a novel technique that can handle high salinity waste, employing favorable properties past the critical point of water. In this study, process simulation of two supercritical water desalination (SCWD) scenarios was completed. Zero liquid discharge (ZLD) and brine concentration operating scenarios were compared, weighing the associated economics and benefits for each case. The results were shown to be economically feasible for brines with a high dissolved solids content, ranging from 3.49 to 17.28·bbl -1 in an expanded sensitivity analysis.

42 ENGINEERING↗
DOE OSTI2022

Self-consistent dispersal puts tight constraints on the spatiotemporal organization of species-rich metacommunities

Dispersal can be critical to the maintenance of ecosystems as it allows local communities to be recolonized after extinction. However, it remains unclear whether the extinction-mitigating effect of dispersal persists when the number of competing species is large. Based on a spatially explicit mathematical description of metacommunities, we show that when many species coexist, each species operates near its extinction threshold, barely surviving due to dispersal. This has general consequences for spatiotemporal abundance patterns. For short-range dispersal, species organize into fractal spatiotemporal extinction patterns characteristic of a directed percolation phase transition. As species approach their extinction threshold, biodiversity is very sensitive to perturbation, suggesting that dispersal within a metacommunity puts tight constraints on the robustness and evolution of species-rich metacommunities. Biodiversity is often attributed to a dynamic equilibrium between the immigration and extinction of species. This equilibrium forms a common basis for studying ecosystem assembly from a static reservoir of migrants—the mainland. Yet, natural ecosystems often consist of many coupled communities (i.e., metacommunities), and migration occurs between these communities. The pool of migrants then depends on what is sustained in the ecosystem, which, in turn, depends on the dynamic migrant pool. This chicken-and-egg problem of survival and dispersal is poorly understood in communities of many competing species, except for the neutral case—the “unified neutral theory of biodiversity.” Employing spatiotemporal simulations and mean-field analyses, we show that self-consistent dispersal puts rather tight constraints on the dynamic migration–extinction equilibrium. When the number of species is large, species are pushed to the edge of their global extinction, even when competition is weak. As a consequence, the overall diversity is highly sensitive to perturbations in demographic parameters, including growth and dispersal rates. When dispersal is short range, the resulting spatiotemporal abundance patterns follow broad scale-free distributions that correspond to a directed percolation phase transition. The qualitative agreement of our results for short-range and long-range dispersal suggests that this self-organization process is a general property of species-rich metacommunities. Our study shows that self-sustaining metacommunities are highly sensitive to environmental change and provides insights into how biodiversity can be rescued and maintained.

54 ENVIRONMENTAL SCIENCES↗