Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “terrestrial models”

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 217 records · Page 12

Isolating the Impacts of Soil Moisture Initialization on Land Carbon Fluxes at the Sub-Seasonal to Seasonal Scale

Soil moisture, through its impact on transpiration, plays a significant role in controlling the productivity of terrestrial vegetation. Initial soil moisture conditions in a seasonal forecast may therefore affect the forecast of land carbon uptake. Here we investigate the relative impacts of soil moisture initialization and carbon reservoir initialization on forecasts of land carbon fluxes at the sub-seasonal to seasonal (S2S) scale. The bias-corrected, retrospective meteorology of NASA Global Modeling and Assimilation Office (GMAO)’s S2S ensemble forecast was used to force the stand-alone Catchment-CN model and thereby estimate terrestrial carbon responses out to nine lead months. Our results show that soil moisture initialization is a major contributor (approximately 44%) to the high global carbon uptake forecast skill seen during the first three lead months. The carbon reservoir initialization explains roughly another half of the monthly carbon forecast skill during this period and becomes relatively more important at longer leads (while the overall forecast skill decreases after the 3 lead month), suggesting a slower but longer-lasting influence of carbon reservoir initialization on carbon fluxes. Our results highlight the significance of a good soil moisture initialization for improved forecasts of carbon fluxes at leads of several months, further support for the usefulness of assimilating the satellite based soil moisture information into terrestrial biosphere models and of a short term carbon forecast to understand current events with often a lag in the availability of flux estimates.

Land carbon↗

Assimilation of GRACE Terrestrial Water Storage Data into a Land Surface Model

The NASA Gravity Recovery and Climate Experiment (GRACE) system of satellites provides observations of large-scale, monthly terrestrial water storage (TWS) changes. In. this presentation we describe a land data assimilation system that ingests GRACE observations and show that the assimilation improves estimates of water storage and fluxes, as evaluated against independent measurements. The ensemble-based land data assimilation system uses a Kalman smoother approach along with the NASA Catchment Land Surface Model (CLSM). We assimilated GRACE-derived TWS anomalies for each of the four major sub-basins of the Mississippi into the Catchment Land Surface Model (CLSM). Compared with the open-loop (no assimilation) CLSM simulation, assimilation estimates of groundwater variability exhibited enhanced skill with respect to measured groundwater. Assimilation also significantly increased the correlation between simulated TWS and gauged river flow for all four sub-basins and for the Mississippi River basin itself. In addition, model performance was evaluated for watersheds smaller than the scale of GRACE observations, in the majority of cases, GRACE assimilation led to increased correlation between TWS estimates and gauged river flow, indicating that data assimilation has considerable potential to downscale GRACE data for hydrological applications. We will also describe how the output from the GRACE land data assimilation system is now being prepared for use in the North American Drought Monitor.

Reichle, Rolf H.↗

The Martian dust cycle: A proposed model

Despite more than a decade of study of martian dust storms, many of their characteristics and associated processes remain enigmatic, including the mechanisms for dust raising, modes of settling, and the nature of dust deposits. However, observations of Mars dust, considerations of terrestrial analogs, theoretical models, and laboratory simulations permit the formulation of a Martian Dust Cycle Model, which consists of three main processes: (1) suspension threshold, (2) transportation, and (3) deposition; two associated processes are also included: (4) dust removal and (5) the addition of new dust to the cycle. Although definitions vary, dust includes particles less than 4 to approx. 60 microns in diameter, which by terrestrial usage includes silt, loess, clay, and aerosolic dust particles. The dust cycle model is explained.

Greeley, Ronald↗

Atmospheric Constraints on Landing Site Selection

The Martian atmosphere is a significant part of the environment that the Mars Exploration Rovers (MER) will encounter. As such, it imposes important constraints on where the rovers can and cannot land. Unfortunately, as there are no meteorological instruments on the rovers, there is little atmospheric science that can be accomplished, and no scientific preference for landing sites. The atmosphere constrains landing site selection in two main areas, the entry descent and landing (EDL) process and the survivability of the rovers on the surface. EDL is influenced by the density profile and boundary layer winds (up to altitudes of 5 to 10 km). Surface survivability involves atmospheric dust, temperatures and winds. During EDL, the atmosphere is used to slow the lander down, both ballistically and on the parachute. This limits the maximum elevation of the landing site to -1.3 km below the MOLA reference aeroid. The landers need to encounter a sufficiently dense atmosphere to be able to stop, and the deeper the landing site, the more column integrated atmosphere the lander can pass through before reaching the surface. The current limit was determined both by a desire to be able to reach the hematite region and by a set of atmosphere models we developed for EDL simulations. These are based on Thermal Emission Spectrometer (TES) atmospheric profile measurements, Ames Mars General Circulation Model (MGCM) results, and the 1-D Ames GCM radiative/convective model by J. Murphy. The latter is used for the near surface diurnal cycle. The current version of our model encompasses representative latitude bands, but we intend to make specific models for the final candidate landing sites to insure that they fall within the general envelope. The second constraint imposed on potential landing sites through the EDL process is the near surface wind. The wind in the lower approximately 5 km determines the horizontal velocity that the landers have when they land. Due to the mechanics of the landing process, the total velocity (including both the horizontal and vertical components) determines whether or not the landers are successful. Unfortunately, the landing system has no easy way to nullify any horizontal velocity imparted by the wind, so the landing sites selected need to have as little wind as possible. In addition to the mean wind velocity, the landing system is sensitive to vertical wind shear in the lowest kilometer or so. Wind shear can deflect the retro rockets (RADs) from their nominal vertical orientation producing unwanted horizontal spacecraft velocities. Both mean velocity and wind shear are dominated by the the local topography and other surface properties (in particular albedo and thermal inertia which control the surface temperature). This is seen even in simplified 2-D mesoscale models. The effects in a fully 3-D model are expected to he even more topographically dependent. In particular there is potential for wind channeling in canyons and other terrain features. Boundary layer winds and wind shear are currently being modeled based on terrestrial data and boundary layer scaling laws modified for Martian conditions. We hope to supplement this with mesoscale model results (from several sources) once the number of landing sites is reduced to a manageable number.

Kass, David M.↗

Application of PLIF to Investigate the Hypersonic Wake of a LOFTID-Relevant Model at Mach 10

An experimental investigation of hypersonic wake flows using the planar laser-induced fluorescence (PLIF) measurement technique is summarized in this report. Off-body measurements, primarily flow visualization and velocity, were obtained for the wake of a blunt body model relevant for terrestrial re-entry and Martian entry applications within a Mach 10 hypersonic flow. The design of the model was based on the forebody aeroshell used for the LOFTID (Low-Earth Orbit Flight Test of an Inflatable Decelerator) flight vehicle. This document serves to document the experimental conditions and provide an initial analysis of select test conditions.

PLIF↗

Validation Assessment for the Soil Moisture Active Passive (SMAP) Level 4 Carbon (L4_C) Data Product Version 5

The post-launch Cal/Val phase of the SMAP mission is guided by two primary objectives for each science product team: 1) to calibrate, verify, and improve the performance of the science algorithms, and 2) validate accuracies of the science data products as specified in the SMAP Level-1 mission science requirements. Algorithm science and product maintenance activities during the SMAP extended mission phase have also involved periodic algorithm calibration and product refinements to maintain or enhance product consistency and performance as well as science utility. This report provides an assessment of the latest (Version 5) SMAP Level 4 Carbon (L4_C) product. The L4_C Version 5 (v5) global record now spans more than six years (March 2015 – present) of SMAP operations and has benefited from five major reprocessing updates to the operational product. These reprocessing events and L4_C product release updates have incorporated various algorithm refinements and calibration adjustments to account for similar refinements to the upstream GEOS land model assimilation system, SMAP brightness temperatures, and MODIS vegetation inputs used for L4_C processing. The SMAP L4_C algorithms utilize a terrestrial carbon flux model informed by daily surface and root zone soil moisture information contributed from the SMAP Level 4 Soil Moisture (L4_SM) product along with optical remote sensing-based (e.g. MODIS-based) land cover and canopy fractional photosynthetic active radiation (fPAR), and other ancillary biophysical data. The carbon flux model estimates global daily net ecosystem CO2 exchange (NEE) and the component carbon fluxes, namely, vegetation gross primary production (GPP) and soil heterotrophic respiration (Rh). Other L4_C product elements include surface (~0-5 cm depth) soil organic carbon (SOC) stocks and associated environmental constraints to these processes, including soil moisture-related controls on GPP and ecosystem respiration (Kimball et al. 2014, Jones et al. 2017). The L4_C product addresses SMAP carbon cycle science objectives by: 1) providing a direct link between terrestrial carbon fluxes and underlying freeze/thaw and soil moisture-related constraints to these processes, 2) documenting primary connections between terrestrial water, energy and carbon cycles, and 3) improving understanding of terrestrial carbon sink activity. The SMAP L4_C algorithms and operational product are mature and at a CEOS Validation Stage 4 level (Jackson et al. 2012) based on extensive validation of the multi-year record against a diverse array of independent benchmarks, well characterized global performance, and systematic refinements gained from five major reprocessing events. There are no Level-1 mission science requirements for the L4_C product; however, self-imposed requirements have been established focusing on NEE as the primary product field for validation, and on demonstrating L4_C accuracy and success in meeting product science requirements (Jackson et al. 2012). The other L4_C product fields also have strong utility for carbon science applications (e.g., Liu et al. 2019, Endsley et al. 2020); however, analysis of these other fields is considered secondary relative to primary validation activities focusing on NEE. The L4_C targeted accuracy requirements are to meet or exceed a mean unbiased root-mean-square error (ubRMSE, or standard deviation of the error) for NEE of 1.6 g C m-2 d-1 and 30 g C m-2 yr-1, emphasizing northern (≥45°N) boreal and arctic ecosystems; this accuracy is similar to that of tower eddy covariance measurement-based observations (Baldocchi 2008). Methods used for the latest v5 L4_C product performance and validation assessment have been established from the SMAP Cal/Val plan and previous studies (Jackson et al. 2012, Jones et al. 2017) and include: 1) consistency evaluations of the product fields against earlier product releases (version 4 or earlier); 2) comparisons of daily carbon flux estimates with independent tower eddy 4covariance measurement-based daily carbon (CO2) flux observations from core tower validation sites (CVS); and 3) consistency checks against other global carbon products, including soil carbon inventory records, global GPP records derived from tower observation upscaling methods, and satellite-based observations of canopy solar induced chlorophyll fluorescence (SIF) as a surrogate for GPP. Metrics used to evaluate relative agreement between L4_C product fields and observational benchmarks include correlation (r-value), RMSE differences, bias and model sensitivity diagnostics. Following these validation criteria, the present report provides a validation assessment of the latest L4_C product release (v5). Detailed descriptions of the L4_C algorithm and additional global product accuracy and performance results are given elsewhere (Jones et al. 2017, Endsley et al. 2020). The v5 L4_C product replaces earlier product versions and continues to show: (i) accuracy and performance levels meeting or exceeding SMAP L4_C science requirements; (ii) improvement over the previous product version (version 4); and (iii) suitability for a diversity of science applications. Example L4_C applications from the recent literature include clarifying environmental trends and controls on the northern terrestrial carbon sink (Liu et al. 2019), diagnosing drought-related impacts on ecosystem productivity (Li et al. 2020), and regional monitoring of cropland conditions for projecting annual yields (Wurster et al. 2020). 2EXPECTED L4_C ALGORITHM AND PRODUCTPERFORMANCE The L4_C algorithm performance, including variance and uncertainty estimates of model outputs, was determined during the mission pre-launch phase through spatially explicit model sensitivity studies using available model inputs similar to those currently being used for operational production and evaluating the resulting model simulations over the observed range of northern (≥45 °N) and global conditions (Kimball et al. 2012, Entekhabi et al. 2014). The L4_C algorithm options were also evaluated during the mission prelaunch phase, including deriving canopy fPAR from lower order NDVI (Normalized Difference Vegetation Index) inputs in lieu of using MODIS (MOD15) fPAR; and including an explicit model representation of boreal fire disturbance recovery impacts. These results indicated that the L4_C accuracy requirements (i.e., NEE ubRMSE ≤ 30 g C m-2 yr-1or ≤ 1.6 g C m-2 d-1) could be met from the baseline algorithms over more than 82% and 89% of global and northern vegetated land areas, respectively (Yi et al. 2013, Kimball et al. 2014). The global L4_C algorithm error budget for NEE derived during the mission prelaunch phase indicated that the estimated NEE ubRMSE uncertainty is proportional to GPP and is therefore larger in higher biomass productivity areas, including forests and croplands (Kimball et al. 2014). Likewise, NEE ubRMSE uncertainty is expected to be lower in less-productive areas, including grasslands and shrublands. Expected model NEE ubRMSE levels were also generally within targeted accuracy levels for characteristically less-productive boreal and Arctic biomes, even though relative model error as a proportion of total productivity (NEE RMSE / GPP) may be large in these areas. The estimated NEE uncertainty was lower than expected in some warmer tropical high biomass productivity areas (e.g. Amazon rainforest) because of reduced low temperature and moisture constraints to the L4_C respiration calculations so that the bulk of model uncertainty is contributed by GPP in these areas. Model NEE uncertainty in the African Congo was estimated to

SMAP↗

Re-evaluation of IIH as the Ideal Terrestrial Analog for Sans: Is There a Better Model to Consider?

While astronauts are returning from long duration spaceflight with multiple ocular signs that mimic those seen in terrestrial patients with elevated intracranial pressure (ICP), evidence has yet to prove a clinically significant increase in ICP during space.1 Preliminary research evidence may even suggest that ICP decreases in microgravity. Idiopathic intracranial hypertension (IIH) has long been considered the ideal terrestrial analogue to Spaceflight Associated Neuro-ocular Syndrome (SANS).1 However, there are several critical features of SANS that do not complement any reported case of IIH on Earth. These findings mandate a closer look at the accuracy of IIH as a terrestrial SANS analog.

Batliwala, Shehzad↗

The magnetic field of Mercury

The paper examines the magnetic field observations and their analyses relating to the determination of the Mercury magnetic field. Methods of analyzing data included: (1) comparison of bow shock and magnetopause relative positions at Mercury to the earth, (2) direct spherical harmonic analysis, (3) magnetosphere modeling by an image dipole, and (4) scaling of a mathematical model for the terrestrial magnetosphere. Dipole moments were determined using partial quadrupole and octupole terms to improve the least-square fit of models to observations; analyses by method (2) yield a convergent series of dipole moments values considered to best represent the intrinsic planetary field. Finally, it is suggested that the origin of the magnetic field of Mercury cannot be uniquely determined, but the sources of convective energy may be radiogenic decay and heat release, gravitational settling, and differentiation of processional torques.

Ness, N. F.↗

The Invasive Species Forecasting System

The Invasive Species Forecasting System (ISFS) provides computational support for the generic work processes found in many regional-scale ecosystem modeling applications. Decision support tools built using ISFS allow a user to load point occurrence field sample data for a plant species of interest and quickly generate habitat suitability maps for geographic regions of management concern, such as a national park, monument, forest, or refuge. This type of decision product helps resource managers plan invasive species protection, monitoring, and control strategies for the lands they manage. Until now, scientists and resource managers have lacked the data-assembly and computing capabilities to produce these maps quickly and cost efficiently. ISFS focuses on regional-scale habitat suitability modeling for invasive terrestrial plants. ISFS s component architecture emphasizes simplicity and adaptability. Its core services can be easily adapted to produce model-based decision support tools tailored to particular parks, monuments, forests, refuges, and related management units. ISFS can be used to build standalone run-time tools that require no connection to the Internet, as well as fully Internet-based decision support applications. ISFS provides the core data structures, operating system interfaces, network interfaces, and inter-component constraints comprising the canonical workflow for habitat suitability modeling. The predictors, analysis methods, and geographic extents involved in any particular model run are elements of the user space and arbitrarily configurable by the user. ISFS provides small, lightweight, readily hardened core components of general utility. These components can be adapted to unanticipated uses, are tailorable, and require at most a loosely coupled, nonproprietary connection to the Web. Users can invoke capabilities from a command line; programmers can integrate ISFS's core components into more complex systems and services. Taken together, these features enable a degree of decentralization and distributed ownership that have helped other types of scientific information services succeed in recent years.

Schnase, John↗

Using Orbiting Carbon Observatory-2 (OCO-2) column CO2 retrievals to rapidly detect and estimate biospheric surface carbon flux anomalies

The global carbon cycle is experiencing continued perturbations via increases in atmospheric carbon concentrations, which are partly reduced by terrestrial biosphere and ocean carbon uptake. Greenhouse gas satellites have been shown to be useful in retrieving atmospheric carbon concentrations and observing surface and atmospheric CO2 seasonal-to-interannual variations. However, limited attention has been placed on using satellite column CO2 retrievals to evaluate surface CO2 fluxes from the terrestrial biosphere without advanced inversion models at low latency. Such applications could be useful to monitor, in near real time, biosphere carbon fluxes during climatic anomalies like drought, heatwaves, and floods, before more complex terrestrial biosphere model outputs and/or advanced inversion modelling estimates become available. Here, we explore the ability of Orbiting Carbon Observatory-2 (OCO-2) column-averaged dry air CO2 (XCO2) retrievals to directly detect and estimate terrestrial biosphere CO2 flux anomalies using a simple mass-balance approach. An initial global analysis of surface–atmospheric CO2 coupling and transport conditions reveals that the western US, among a handful of other regions, is a feasible candidate for using XCO2 for detecting terrestrial biosphere CO2 flux anomalies. Using the CarbonTracker model reanalysis as a test bed, we first demonstrate that a well-established mass-balance approach can estimate monthly surface CO2 flux anomalies from XCO2 enhancements in the western United States. The method is optimal when the study domain is spatially extensive enough to account for atmospheric mixing and has favorable advection conditions with contributions primarily from one background region. We find that errors in individual soundings reduce the ability of OCO-2 XCO2 to estimate more frequent, smaller surface CO2 flux anomalies. However, we find that OCO-2 XCO2 can often detect and estimate large surface flux anomalies that leave an imprint on the atmospheric CO2 concentration anomalies beyond the retrieval error/uncertainty associated with the observations. OCO-2 can thus be useful for low-latency monitoring of the monthly timing and magnitude of extreme regional terrestrial biosphere carbon anomalies.

Andrew F. Feldman↗

Terrestrial spreading centers under Venus conditions - Evaluation of a crustal spreading model for Western Aphrodite Terra

The model of Reid and Jackson (1981) for terrestrial spreading centers is applied to Venus conditions. On the basis of spreading rate, mantle temperature, and surface temperature, the model predicts both isostatic topography and crustal thickness. The model and Pioneer Venus altimetry and gravity data are used to test the hypothesis of Head and Crumpler (1987) that Western Aphrodite Terra is the location of crustal spreading on Venus. It is concluded that a spreading center model for Ovda Regio in Western Aphrodite Terra could account for the observed topography and line-of-sight gravity anomalies found in the Pioneer data.

Sotin, C.↗

Similar explosive eruptions of lunar and terrestrial volcanoes

A mathematical model of one-dimensional, steady duct flow of a mixture of gas and small solid particles (rock), has been applied to the earth and the moon under geometrically and dynamically similar conditions to simulate explosive eruption conditions. Numerical results for equilibrium two-phase flows for lunar and terrestrial explosive eruptions under similar conditions indicate that: (1) a lunar explosive vent is much larger than the corresponding terrestrial vent, (2) the exit velocity from the lunar explosive flow may be higher than the lunar escape velocity, but the exit velocity of terrestrial explosive flow is much less than that of the lunar case (this result supports the hypothesis that Australian tektites came from the moon as a stream of a mixture of rock and gas of extremely high speed), and (3) the thermal effects on the lunar explosive flows are much larger than those of the terrestrial case.

Pai, S. I.↗

The geologic history of quartz-normative and olivine-normative basalts in the vicinity of Hadley Rille (Apollo 15)

The geologic history of the quartz normative (QNB) and olivine normative (ONB) basalt types at Hadley Rille are discussed. A model for the geology of the mare basalts was constructed from a combination of field observations, sample chemistry, sample petrology and personal bias from terrestrial experience. The model proposes that the QNBs are the only mare lava type that is present as outcrop in the area traversed by the astronauts during the Apollo 15 mission. The returned QNB samples formed during a single eruptive phase of the Hadley Rille lava tube system. The ONB lavas are an exotic component transported to the site by a cratering event, or the ONBs are samples excavated from older are bedrock that was partly covered by the QNB lavas.

Grove, T. L.↗

Origin and evolution of planetary atmospheres

This report concerns several research tasks related to the origin and evolution of planetary atmospheres and the large-scale distribution of volatile elements in the Solar System. These tasks and their present status are as follows: (1) we have conducted an analysis of the volatility and condensation behavior of compounds of iron, aluminum, and phosphorus in the atmosphere of Venus in response to publish interpretations of the Soviet Venera probe XRF experiment data, to investigate the chemistry of volcanic gases, injection of volatiles by cometary and asteroidal impactors, and reactions in the troposphere; (2) we have completed and are now writing up our research on condensation-accretion modeling of the terrestrial planets; (3) we have laid the groundwork for a detailed study of the effects of water transport in the solar nebula on the bulk composition, oxidation state, and volatile content of preplanetary solids; (4) we have completed an extensive laboratory study of cryovolcanic materials in the outer solar system; (5) we have begun to study the impact erosion and shock alteration of the atmosphere of Mars resulting from cometary and asteroidal bombardment; and (6) we have developed a new Monte Carlo model of the cometary and asteroidal bombardment flux on the terrestrial planets, including all relevant chemical and physical processes associated with atmospheric entry and impact, to assess both the hazards posed by this bombardment to life on Earth and the degree of cross-correlation between the various phenomena (NO(x) production, explosive yield, crater production, iridium signature, etc.) that characterize this bombardment. The purpose of these investigations has been to contribute to the developing understanding of both the dynamics of long-term planetary atmosphere evolution and the short-term stability of planetary surface environments.

Lewis, John S.↗

The magnetic field of Mercury

Data from Mariner 10 observations of Mercury indicate that there exists an intrinsic magnetic field of the planet, sufficiently strong at present to deflect the solar wind flow around the planet and to form a detached bow shock wave in the super Alfvenic solar wind. Four methods used to analyze the magnetic field data and derive quantitative values for the description of the planetary field include (1) comparison of bow shock and magnetopause relative positions at Mercury to those at Earth; (2) direct spherical harmonic analysis of the data; (3) modeling of the magnetosphere by an image dipole and infinite 2-D current sheet in addition to the planetary field; and (4) scaling of a mathematical model for the terrestrial magnetosphere. The results obtained yield dipole moments ranging from 2.4 to 5.1x10 to the 22d power, with the lower values associated cw cm with certain models using partial quadrupole and octupole terms to improve the least squares fitting of models to observations.

Ness, N. F.↗

Global changes in biogeochemical cycles in response to human activities

The main objective of our research was to characterize biogeochemical cycles at continental and global scales in both terrestrial and aquatic ecosystems. This characterization applied to both natural ecosystems and those disturbed by human activity. The primary elements of interest were carbon and nitrogen and the analysis sought to quantify standing stocks and dynamic cycling processes. The translocation of major nutrients from the terrestrial landscape to the atmosphere (via trace gases) and to fluvial systems (via leaching, erosional losses, and point source pollution) were of particular importance to this study. Our aim was to develop the first generation of Earth System Models. Our research was organized around the construction and testing of component biogeochemical models which treated terrestrial ecosystem processes, aquatic nutrient transport through drainage basins, and trace gas exchanges at the continental and global scale. A suite of three complementary models were defined within this construct. The models were organized to operate at a 1/2 degree latitude by longitude level of spatial resolution and to execute at a monthly time step. This discretization afforded us the opportunity to understand the dynamics of the biosphere down to subregional scales, while simultaneously placing these dynamics into a global context.

Moore, Berrien, III↗

The AE-8 trapped electron model environment

The machine sensible version of the AE-8 electron model environment was completed in December 1983. It has been sent to users on the model environment distribution list and is made available to new users by the National Space Science Data Center (NSSDC). AE-8 is the last in a series of terrestrial trapped radiation models that includes eight proton and eight electron versions. With the exception of AE-8, all these models were documented in formal reports as well as being available in a machine sensible form. The purpose of this report is to complete the documentation, finally, for AE-8 so that users can understand its construction and see the comparison of the model with the new data used, as well as with the AE-4 model.

Vette, James I.↗

Feeding Ten Billion People Is Possible Within Four Terrestrial Planetary Boundaries

Global agriculture puts heavy pressure on planetary boundaries, posing the challenge to achieve future food security without compromising Earth system resilience. On the basis of process-detailed, spatially explicit representation of four interlinked planetary boundaries (biosphere integrity, land-system change, freshwater use, nitrogen flows) and agricultural systems in an internally consistent model framework, we here show that almost half of current global food production depends on planetary boundary transgressions. Hotspot regions, mainly in Asia, even face simultaneous transgression of multiple underlying local boundaries. If these boundaries were strictly respected, the present food system could provide a balanced diet (2,355 kcal per capita per day) for 3.4 billion people only. However, as we also demonstrate, transformation towards more sustainable production and consumption patterns could support 10.2 billion people within the planetary boundaries analysed. Key prerequisites are spatially redistributed cropland, improved water–nutrient management, food waste reduction and dietary changes. Adoption of the Sustainable Development Goals by all nations in 2015 is the first ever commitment to a world development path that safeguards the stability of the Earth system as a prerequisite for meeting universal human standards1. The longstanding challenge of achieving food security through sustainable agriculture is particularly acute in this context as world agriculture is a leading cause for the current transgressions of multiple planetary boundaries (PBs) globally and regionally2–5. The PB framework is a comprehensive scientific attempt to synoptically define our planet’s biogeophysical limits to anthropogenic interference. It suggests bounds to nine interacting processes that together delineate a Holocene-like Earth system state. The Holocene is chosen as the reference state as it is the only period known to provide a safe operating space for a world population of several billion people, and according to a precautionary principle, the PBs are set in sufficient distance from processes that may critically undermine Earth system resilience and global sustainability. A challenging question, thus, is whether human development goals such as food security can be met while maintaining multiple PBs along with their subglobal manifestations. Further PB transgressions could jeopardize the chances of providing sufficient food for a world population projected to be wealthier and reach >9 billion by 2050. This conundrum portrays a tradeoff between Earth’s biophysical carrying capacity and humankind’s rising food demand, calling in response for radical rethinking of food production and consumption patterns6–9. Yield gap closures, avoidance of excessive input use, shifts towards less resource-demanding diets, food waste reductions and efficient international trade are crucial options for sustainably increasing the food supply10–15. For example, enhancing water-use efficiency on irrigated and rain-fed farms can triple or quadruple crop yields in low-performing systems, suggesting possible global gains of >20% (ref. 16). Even higher gains appear feasible through globally optimized configurations of the land-use pattern17, and cutting food losses by half could generate food for another billion people18. Thus, collective large-scale implementation of such options could sustain food for a further growing world population19. Yet achieving this within a safe operating space as defined by PBs requires not only a halt to but actually a reversal of existing PB transgressions. Previous studies suggest that such a reconciliation might be possible, but these were based on aggregate representations of PBs (not accounting for the spatial patterns of limits, transgressions and interactions) or considered only one boundary in isolation17,20–23. Here, we systematically quantify to what extent current food production depends on local to global transgressions of the PBs for biosphere integrity, land-system change, freshwater use and nitrogen (N) flows, along with the potential of a range of solutions to avoid these transgressions and still increase food supply (Table 1). To this end, we configured an internally consistent process-based model of the terrestrial biosphere including agriculture (LPJmL) with multiple spatially distributed PBs and their interactions. LPJmL is among the longest-established and best-evaluated biosphere models, showing robust performance regarding simulation of, for example, carbon, water and crop yield dynamics (Supplementary Figs. 1 and 2 and Supplementary Table 1; see ref. 24 for a comprehensive benchmarking and Supplementary Methods for more detail on model evaluations). In principle following established definitions4, we refine the computation of some PBs with respect to their regional patterns and interactions (Methods), providing globally gridded precautionary limits to human interference with the Earth system at a level of great detail. In particular, we account for the evidence that many PBs need to be represented spatially explicitly4 to cover their

Gerten, Dieter↗