Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Planetary Science”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Quasi-stationary Planetary Scale Waves in Modern Climate Models

We evaluate climate model simulations of geopotential height at the largest space-time scales, considering the first Fourier components (in longitude) of long-term climatological means. Quasi-stationary planetary scale waves emerge from this analysis. Variations of these waves are important components of climate change, for example associated with regional drought. Systematic examination of their long-term climatology provides context for understanding their evolution in time. We compare four reanalysis datasets with “historical runs” from the latest version of the Coupled Model Intercomparison Project (CMIP) and find general agreement, but we also identify model errors outside the range of both observational uncertainty and the uncertainty implied by interannual variability. Errors in wave phase have possible implications for climate model projections of regional climate and of shifting weather patterns. One model in the set we examine, originally created nearly 20 years before the others, exhibits noticeably greater errors.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Constraining Mornings and Evenings on Distant Worlds: A new Semianalytical Approach and Prospects with Transmission Spectroscopy

The technique of transmission spectroscopy—the variation of a planetary radius with wavelength due to opacity sources in the planet’s terminator region—has been one of the most successful in the characterization of exoplanet atmospheres to date, providing key insights into the composition and structure of these distant worlds. A common assumption made when using this technique, however, is that the variations are the same in the entire terminator region. In reality, the morning and evening terminators might have distinct temperature, pressure, and thus compositional profiles due to the inherent 3D nature of the planet, which would, in turn, give rise to different spectra on each side of it. Constraining those might be fundamental for our understanding of not only the weather patterns in these distant worlds but also the planetary formation signatures that might only be possible to extract once these features are well understood. Motivated by this physical picture, in this work, we perform a detailed study of the observational prospects of detecting this effect. We present an open-source semianalytical framework with which this information can be extracted directly from transit light curves and perform a detailed study of the prospects of detecting the effect with current missions, such as TESS, and upcoming ones, such as JWST. Our results show that these missions show great promise for the detection of this effect. Transmission spectroscopy studies with JWST, in particular, could provide spectra of each of the limbs, allowing us to convey 3D information previously accessible only via phase curves.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

Air quality impact of the Northern California Camp Fire of November 2018

The Northern California Camp Fire that took place in November 2018 was one of the most damaging environmental events in California history. Here, we analyze ground-based station observations of airborne particulate matter that has a diameter <2.5 µm (PM 2.5 ) across Northern California and conduct numerical simulations of the Camp Fire using the Weather Research and Forecasting model online coupled with chemistry (WRF-Chem). Simulations are evaluated against ground-based observations of PM 2.5 , black carbon, and meteorology, as well as satellite measurements, such as Tropospheric Monitoring Instrument (TROPOMI) aerosol layer height and aerosol index. The Camp Fire led to an increase in Bay Area PM 2.5 to over 50 µg m –3 for nearly 2 weeks, with localized peaks exceeding 300 µg m –3 . Using the Visible Infrared Imaging Radiometer Suite (VIIRS) high-resolution fire detection products, the simulations reproduce the magnitude and evolution of surface PM 2.5 concentrations, especially downwind of the wildfire. The overall spatial patterns of simulated aerosol plumes and their heights are comparable with the latest satellite products from TROPOMI. WRF-Chem sensitivity simulations are carried out to analyze uncertainties that arise from fire emissions, meteorological conditions, feedback of aerosol radiative effects on meteorology, and various physical parameterizations, including the planetary boundary layer model and the plume rise model. Downwind PM 2.5 concentrations are sensitive to both flaming and smoldering emissions over the fire, so the uncertainty in the satellite-derived fire emission products can directly affect the air pollution simulations downwind. Our analysis also shows the importance of land surface and boundary layer parameterization in the fire simulation, which can result in large variations in magnitude and trend of surface PM 2.5 . Inclusion of aerosol radiative feedback moderately improves PM 2.5 simulations, especially over the most polluted days. Results of this study can assist in the development of data assimilation systems as well as air quality forecasting of health exposures and economic impact studies.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Aquaplanets as a Framework for Examination of Aerosol Effects

Although fundamental to the planetary radiative balance, aerosol impacts are highly uncertain in climate simulations because of the uneven distribution of aerosol sources and the complex interactions with radiation and clouds that are difficult to represent in climate models. This study proposes that aquaplanet configurations represent an idealized framework to investigate aerosol effects. As a simple demonstration, a series of aquaplanet simulations with the Community Atmosphere Model version 6 shows that the spatial distribution of aerosol emissions changes the aerosol effective radiative forcing even with unchanged total emissions. Some statistical properties of the simulations are presented to show that relatively short model integrations yield robust results. Much of the aerosol effect is shown to arise from aerosol–cloud interactions, especially through rapid adjustments associated with the aerosol lifetime effect that alter the cloud optical thickness.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Origin and effect of film sub-stoichiometry on ultraviolet, ns-laser damage resistance of hafnia single layers

Understanding the origin of laser damage-prone precursors in high index materials such as hafnia holds the key to the development of laser damage-resistant multilayer dielectric coated optics for high power and energy laser systems. In this study, we investigate the source of sub-stoichiometry, a potent laser damage precursor, in hafnia films produced by an ion beam sputtering (IBS) deposition method and the effect of such defects on the film performance upon ns ultraviolet (UV) laser (8 ns, 355 nm) exposure. Chemical analysis of data obtained via Rutherford backscattering spectroscopy (RBS) suggests that hafnia films deposited at two different planetary locations from the same deposition run exhibit anisotropic and location-dependent stoichiometries. While the oxygen-to-hafnium ratio is at the stoichiometric value of 2 for the hafnia film at the edge location, the ratio is significantly deviated and is 1.7 for that deposited at the planetary center. The sub-stoichiometric hafnia films display a much lower 1-on-1 damage onset at 1.6 ± 0.2 J/cm 2 compared to 2.3 ± 0.2 J/cm 2 in a stoichiometric film. The low damage performance films also have an over three times higher damage density at fluences above initiation. Coupled with Monte Carlo simulations, we reveal that sub-stoichiometry is primarily attributed to preferential removal of oxygen during film deposition by the bombardment of energetic reflected argon neutrals. The resulting oxygen deficiencies create the sub-bandgap states which facilitate the strong laser energy coupling and reduce the resistance to laser-induced damage in the hafnia single layer films.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Tropical Cyclone Precipitation Response to Surface Warming in Aquaplanet Simulations With Uniform Thermal Forcing

While many modeling studies have attempted to estimate how tropical cyclone (TC) precipitation is impacted by climate change, the multitude of analysis techniques and methodologies have resulted in varying conclusions. Simplified models may be able to help overcome this problem. Radiative-convective equilibrium (RCE) model simulations have been used in various configurations to study fundamental aspects of Earth's climate. While many RCE modeling studies have focused on TC genesis, intensification, and size, limited work has been done using RCE to study TC precipitation. Here, in this study, the response of TC precipitation to sea surface temperature (SST) change is analyzed in global Community Atmosphere Model (CAM) aquaplanet simulations run with Radiative-Convective Equilibrium Model Intercomparison Project protocols, with the addition of planetary rotation. We expect that the insight gained about how TC precipitation responds to SST warming will help predict how TCs in the real world respond to climate change. In the CAM RCE simulations, the warmer SST simulations have less TCs on average, but the TCs tend to be larger in outer size and more intense. As simulation SST increases, more extreme precipitation rates occur within TCs, and more of the TC precipitation comes from these extreme rates. For extreme (99th percentile) TC precipitation, SST, and TC intensity increases dominate the 8.6% per K increase, while TC outer size changes have little impact. For accumulated TC precipitation, SST, and TC intensity contributions are still the majority, but TC outer size changes also contribute to the 6.6% per K increase.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

High-resolution analysis of clay minerals and amorphous materials in martian analog environments

Mudrocks and mud-sized sediments (i.e., silt to clay) dominate the surface of Earth and Mars. These fine-grained sediments preserve a rich history of sedimentary processes from source to sink and shed light on ancient climates. However, both the physical and chemical nature of these materials make them difficult to fully characterize with traditional laboratory techniques. Here, we explore a cross-disciplinary and high-resolution approach using synchrotron radiation for X-ray diffraction, pair distribution function analysis, and submicrometer-scale X-ray fluorescence, combined with transmission electron microscopy, to better understand the nanostructure and composition of mud-sized sediments from a glacio-fluvial watershed in southwest Iceland. Our results demonstrate that sediments in the cold and wet climate of Iceland are more altered than previously thought, as evidenced by the identification of kaolinite and mixed-layer kaolinite-smectite. Additionally, sediments are enriched in amorphous materials and nanocrystalline phases, as determined from grain morphologies and compositions consistent with allophane, hisingerite, ferrihydrite, and halloysite. These alteration products are present as intimate mixtures that vary across depositional sites, demonstrating the dynamic nature of the secondary assemblage from source to sink. This work has implications for Mars, where, for example, basalt-sourced sedimentary rocks from Gale crater are abundant in clay minerals and amorphous materials. Finally, this work underpins the importance of using high-resolution techniques, a coordinated methodology, and developing innovative approaches for future planetary sample return missions (e.g., Mars sample return).

36 MATERIALS SCIENCE↗
DOE OSTI2023

Fulfilling global climate pledges can lead to major increase in forest land on Tibetan Plateau

The UN Climate Change Conference in Glasgow spawned the enhancement and updating of many nations’ climate pledges. Previous research has investigated the effects of these pledges on limiting planetary warming, but their spatially explicit effects on land use/cover are unknown. Here, we linked the Glasgow pledges and the spatially explicit responses of the Tibetan Plateau’s land systems. We found that while fulfilling global climate pledges may not significantly affect the global shares of forestland, grassland/pasture, shrubland, and cropland, it needs a 9.4% increase in the forest area of the Tibetan Plateau. This need is an area 11.4 times the increase of the plateau’s forest in the 2010s, or greater than the size of Belgium. The new forest comes mainly from the medium-density grassland in the Yangtze River basin, calling for more proactive environmental management for the headwaters area of this longest river in Asia.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Interfacial Reactivity and Speciation Emerging from Na-Montmorillonite Interactions with Water and Formic Acid at 200 °C: Insights from Reactive Molecular Dynamics Simulations, Infrared Spectroscopy, and X-ray Scattering Measurements

Reactive organic fluid - mineral interactions at elevated temperatures contribute to the evolution of planetary matter. One of the less studied but important transformations in this regard involves the reactions of formic acid with naturally occurring clays such as sodium montmorillonite. To advance a mechanistic understanding of these interactions, we use ReaxFF reactive molecular dynamics simulations in conjunction with infrared (IR) spectroscopy and X-ray scattering experiments to investigate the speciation behavior of water-formic acid mixtures on sodium montmorillonite interfaces at 473 K and 1 atm. Using a newly developed reactive forcefield, we show that the experimental IR spectra of unreacted and reacted mixture can be accurately reproduced by ReaxFF/MD. We further benchmark the simulation predictions of sodium carbonate and bicarbonate formation in the clay interlayers using Small and Wide-Angle X-ray Scattering measurements. Subsequently, leveraging the benchmarked forcefield, we interrogate the pathway of speciation reactions with emphasis on carbonate, formate, and hydroxide groups elucidating the energetics, transition states, intermediates, and preferred products. Further, we also delineate the differences in reactivities and catalytic effects of clay edges, facets, and interlayers owing to their local chemical environments, which have far reaching consequences in their speciation behavior. The experimental and simulation approaches described in this study and the transferable forcefields can be applied translationally to advance the science of clay-fluid interactions for several applications including subsurface fluid storage and recovery and clay-pollutant dynamics.

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

Turbulent properties from Doppler lidar

Permanent and mobile facilities deployed by the ARM program cover a wide range of climatic regions, with a suite of active remote-sensing of cloud and aerosol profiles complemented by comprehensive in-situ measurements of the chemical, physical and optical properties of aerosol measured at the surface. To be able to fully exploit these datasets requires an understanding of the transport from the surface and within the atmosphere, especially within the planetary boundary layer. This requires a method of identifying when the surface and atmosphere are dynamically-linked. Determining the state of the boundary layer, and defining a well-mixed boundary-layer, is therefore suitably described in terms of mixing; here, we will assume that mixing essentially means turbulent mixing. This necessitates accurate measurements of the winds and turbulent mixing in the atmosphere. Doppler lidar is an excellent tool for observing turbulent properties, and ARM has deployed commercially available Doppler lidars at many of the fixed and mobile facilities. Doppler lidars have been shown to provide accurate retrievals of wind and turbulent properties within the boundary layer , however, robust characterization and calculation of the Doppler lidar measurement uncertainties is required before deriving higher level data products.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

A Shallow-Water Model for Convective Self-Aggregation

Randomly distributed convective storms can self-aggregate in the absence of large-scale forcings. Here we present a 1D shallow-water model to study the convective self-aggregation. This model simulates the dynamics of the planetary boundary layer and represents convection as a triggered process. Once triggered, convection lasts for finite time and occupies finite length. Here we show that the model can successfully simulate self-aggregation, and that the results are robust to a wide range of parameter values. In the simulations, convection excites gravity waves. The gravity waves then form a standing wave pattern, separating the domain into convectively active and inactive regions. We analyze the available potential energy (APE) budget and show that convection generates APE, providing energy for self-aggregation. By performing dimensional analysis, we develop a scaling theory for the size of convective aggregation, which is set by the gravity wave speed, damping time scale, and number density of convective storms. This paper provides a simple modeling framework to further study convective self-aggregation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

A Verification Suite of Test Cases for the Barotropic Solver of Ocean Models

Abstract The development of any atmosphere or ocean model warrants a suite of test cases (TCs) to verify its spatial and temporal discretizations, order of accuracy, stability, reproducibility, portability, scalability, etc. In this paper, we present a suite of shallow water TCs designed to verify the barotropic solver of atmosphere and ocean models. These include the non‐dispersive coastal Kelvin wave; the dispersive inertia‐gravity wave; the dispersive planetary and topographic Rossby waves; the barotropic tide; and a non‐linear manufactured solution. These TCs check the implementation of the linear pressure gradient term; the linear constant or variable‐coefficient Coriolis and bathymetry terms; and the non‐linear advection terms. Simulation results are presented for a variety of time‐stepping methods as well as two spatial discretizations: a mimetic finite volume method based on the TRiSK scheme, and a high‐order discontinuous Galerkin spectral element method. The experimental procedure for conducting these numerical experiments is detailed. It underscores several key considerations that vary depending on the chosen spatial discretization method. Finally, convergence studies of every TC are conducted with refinement in both space and time, only in space, and only in time. The convergence slopes match the expected theoretical predictions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Observed Sensitivity of Low-Cloud Radiative Effects to Meteorological Perturbations over the Global Oceans

Understanding how marine low clouds and their radiative effects respond to changing meteorological conditions is crucial to constrain low-cloud feedbacks to greenhouse warming and internal climate variability. In this paper, we use observations to quantify the low-cloud radiative response to meteorological perturbations over the global oceans to shed light on physical processes governing low-cloud and planetary radiation budget variability in different climate regimes. We assess the independent effect of perturbations in sea surface temperature, estimated inversion strength, horizontal surface temperature advection, 700-hPa relative humidity, 700-hPa vertical velocity, and near-surface wind speed. Stronger inversions and stronger cold advection greatly enhance low-level cloudiness and planetary albedo in eastern ocean stratocumulus and midlatitude regimes. Warming of the sea surface drives pronounced reductions of eastern ocean stratocumulus cloud amount and optical depth, and hence reflectivity, but has a weaker and more variable impact on low clouds in the tropics and middle latitudes. By reducing entrainment drying, higher free-tropospheric relative humidity enhances low-level cloudiness. At low latitudes, where cold advection destabilizes the boundary layer, stronger winds enhance low-level cloudiness; by contrast, wind speed variations have weak influence at midlatitudes where warm advection frequently stabilizes the marine boundary layer, thus inhibiting vertical mixing. These observational constraints provide a framework for understanding and evaluating marine low-cloud feedbacks and their simulation by models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Global simulations of monoterpene-derived peroxy radical fates and the distributions of highly oxygenated organic molecules (HOMs) and accretion products

Abstract. We evaluate monoterpene-derived peroxy radical (MT-RO2) unimolecular autoxidation and self- and cross-reactions with other RO2 species in the GEOS-Chem global chemical transport model. The formation of associated highly oxygenated organic molecules (HOMs) and accretion products are tracked in competition with other bimolecular reactions. Autoxidation is the dominant fate up to 6–8 km for first-generation MT-RO2, which can undergo unimolecular H shifts. Reaction with NO can be a more common fate for H-shift rate constants < 0.1 s−1 or at altitudes higher than 8 km due to the imposed Arrhenius temperature dependence of unimolecular H shifts. For MT-derived HOM-RO2, generated by multistep autoxidation of first-generation MT-RO2, reaction with other RO2 species is predicted to be the major fate throughout most of the boreal and tropical forest regions, whereas reaction with NO dominates in the temperate and subtropical forests of the Northern Hemisphere. The newly added reactions result in an approximate 4 % global average decrease in HO2 and RO2, mainly due to faster self-/cross-reactions of MT-RO2, but the impact upon HO2, OH, and NOx abundances is only important in the planetary boundary layer (PBL) over portions of tropical forests. Predicted HOM concentrations in MT-rich regions and seasons can exceed total organic aerosol predicted by the standard version of the GEOS-Chem model depending on the parameters used. Comparisons to observations reveal that large uncertainties remain for key reaction parameters and processes, especially with respect to the photochemical lifetime and volatility of HOMs as well as the rates and branching of associated RO2 accretion products. Further observations and laboratory studies related to MT-RO2-derived HOMs and gas-phase RO2 accretion product formation kinetics – especially their atmospheric fate, such as gas–particle partitioning, multiphase chemistry, and net secondary organic aerosol formation – are needed.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Recreating the California New Year's Flood Event of 1997 in a Regionally Refined Earth System Model

Abstract The 1997 New Year's flood event was the most costly in California's history. This compound extreme event was driven by a category 5 atmospheric river that led to widespread snowmelt. Extreme precipitation, snowmelt, and saturated soils produced heavy runoff causing widespread inundation in the Sacramento Valley. This study recreates the 1997 flood using the Regionally Refined Mesh capabilities of the Energy Exascale Earth System Model (RRM‐E3SM) under prescribed ocean conditions. Understanding the processes causing extreme events informs practical efforts to anticipate and prepare for such events in the future, and also provides a rich context to evaluate model skill in representing extremes. Three California‐focused RRM grids, with horizontal resolution refinement of 14 km down to 3.5 km, and six forecast lead times, 28 December 1996 at 00Z through 30 December 1996 at 12Z, are assessed for their ability to recreate the 1997 flood. Planetary to synoptic scale atmospheric circulations and integrated vapor transport are weakly influenced by horizontal resolution refinement over California. Topography and mesoscale circulations, such as the Sierra barrier jet, are better represented at finer horizontal resolutions resulting in better estimates of storm total precipitation and storm duration snowpack changes. Traditional time‐series and causal analysis frameworks are used to examine runoff sensitivities state‐wide and above major reservoirs. These frameworks show that horizontal resolution plays a more prominent role in shaping reservoir inflows, namely the magnitude and time‐series shape, than forecast lead time, 2‐to‐4 days prior to the 1997 flood onset.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

TOI-1634 b: An Ultra-short-period Keystone Planet Sitting inside the M-dwarf Radius Valley

Studies of close-in planets orbiting M dwarfs have suggested that the M-dwarf radius valley may be well explained by distinct formation timescales between enveloped terrestrials and rocky planets that form at late times in a gas-depleted environment. This scenario is at odds with the picture that close-in rocky planets form with a primordial gaseous envelope that is subsequently stripped away by some thermally driven mass-loss process. These two physical scenarios make unique predictions of the rocky/enveloped transition’s dependence on orbital separation such that studying the compositions of planets within the M-dwarf radius valley may be able to establish the dominant physics. Here, we present the discovery of one such keystone planet: the ultra-short-period planet TOI-1634 b (P = 0.989 days, F=121F{sub ⊕}, r{sub p}=1.790{sub -0.081}{sup +0.080} R {sub ⊕}) orbiting a nearby M2 dwarf (K {sub s} = 8.7, R {sub s} = 0.450 R {sub ⊙}, M {sub s} = 0.502 M {sub ⊙}) and whose size and orbital period sit within the M-dwarf radius valley. We confirm the TESS-discovered planet candidate using extensive ground-based follow-up campaigns, including a set of 32 precise radial velocity measurements from HARPS-N. We measure a planetary mass of 4.91{sub -0.70}{sup +0.68} M {sub ⊕}, which makes TOI-1634 b inconsistent with an Earth-like composition at 5.9σ and thus requires either an extended gaseous envelope, a large volatile-rich layer, or a rocky composition that is not dominated by iron and silicates to explain its mass and radius. The discovery that the bulk composition of TOI-1634 b is inconsistent with that of Earth supports the gas-depleted formation mechanism to explain the emergence of the radius valley around M dwarfs with M{sub s}≲0.5 M {sub ⊙}.

36 MATERIALS SCIENCE↗
DOE OSTI2024

A Study of Slow-mode Shocks in the Near-Earth Magnetotail with MMS Observations and Hybrid Simulations

The structure of the magnetic reconnection boundary, particularly the presence of slow-mode shocks in the near-Earth magnetotail was studied by using magnetospheric multiscale (MMS) observations and 2.5D hybrid simulations. A total of 51 crossings of MMS from 2017 to 2021 were analyzed. We found that the detection percentage of slow-mode shocks in the near-Earth magnetotail is 41%–55%. Previous studies have only reported one slow-mode shock event in the near-Earth magnetotail and a slow-mode shock detection percentage of 10% or lower in the mid-to-distant magnetotail. It was observed that if the high-energy beam region data is removed from the slow-mode shock downstream observations then the detection of slow-mode shocks reduces, implying that the kinetic effects play an important role in the detection of slow-mode shocks. For the crossings where the interface was not identified as a slow-mode shock, it was found that the turbulence in those crossings can change the mass flux values and disrupt the detection of slow-mode shock. However, the macroscopic slow-mode shock-like structure stably exists around the magnetic reconnection interface, as most of the conditions for slow-mode shocks were satisfied. This result suggests that slow-mode shocks are a general feature of magnetic reconnection geometry. We find that the lack of detection of slow-mode shocks in previous observations and simulations can be explained by taking into account the kinetic structure of slow-mode shocks and the presence of turbulence.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Short-term impacts of 2017 western North American wildfires on meteorology, the atmosphere’s energy budget, and premature mortality

Abstract Western North American fires have been increasing in magnitude and severity over the last few decades. The complex coupling of fires with the atmospheric energy budget and meteorology creates short-term feedbacks on regional weather altering the amount of pollution to which Americans are exposed. Using a combination of model simulations and observations, this study shows that the severe fires in the summer of 2017 increased atmospheric aerosol concentrations leading to a cooling of the air at the surface, reductions in sensible heat fluxes, and a lowering of the planetary boundary layer height over land. This combination of lower-boundary layer height and increased aerosol pollution from the fires reduces air quality. We estimate that from start of August to end of October 2017, ∼400 premature deaths occurred within the western US as a result of short-term exposure to elevated PM 2.5 from fire smoke. As North America confronts a warming climate with more fires the short-term climate and pollution impacts of increased fire activity should be assessed within policy aimed to minimize impacts of climate change on society.

Bernstein, Diana N. (ORCID:0000000308279336)↗