Search NASA⌕ Search

SEARCH · Search NASA

Results for “Model Coupling”

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 739 records · Page 41

Characterization and performance of a field aligned ion cyclotron range of frequency antenna in Alcator C-Mod

Ion cyclotron range of frequency (ICRF) heating is expected to provide auxiliary heating for ITER and future fusion reactors where high Z metallic plasma facing components (PFCs) are being considered. Impurity contamination linked to ICRF antenna operation remains a major challenge particularly for devices with high Z metallic PFCs. Here, we report on an experimental investigation to test whether a field aligned (FA) antenna can reduce impurity contamination and impurity sources. We compare the modification of the scrape of layer (SOL) plasma potential of the FA antenna to a conventional, toroidally aligned (TA) antenna, in order to explore the underlying physics governing impurity contamination linked to ICRF heating. The FA antenna is a 4-strap ICRF antenna where the current straps and antenna enclosure sides are perpendicular to the total magnetic field while the Faraday screen rods are parallel to the total magnetic field. In principle, alignment with respect to the total magnetic field minimizes integrated E|| (electric field along a magnetic field line) via symmetry. A finite element method RF antenna model coupled to a cold plasma model verifies that the integrated E|| should be reduced for all antenna phases. Monopole phasing in particular is expected to have the lowest integrated E||. Consistent with expectations, we observed that the impurity contamination and impurity source at the FA antenna are reduced compared to the TA antenna. In both L and H-mode discharges, the radiated power is 20%–30% lower for a FA-antenna heated discharge than a discharge heated with the TA-antennas. However, inconsistent with expectations, we observe RF induced plasma potentials (via gas-puff imaging and emissive probes to be nearly identical for FA and TA antennas when operated in dipole phasing). Moreover, the highest levels of RF-induced plasma potentials are observed using monopole phasing with the FA antenna. Thus, while impurity contamination and sources are indeed reduced with the FA antenna configuration, the mechanism determining the SOL plasma potential in the presence of ICRF and its impact on impurity contamination and sources remains to be understood.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The effect of urban design parameters on the local microclimate

Two-dimensional steady-state simulations have been performed using the standard k-e turbulence model coupled with the heat transfer models available in the CFD software FLUENT 6.1, in order to examine the impact of radiation on the Urban Heat Island phenomenon. Specifically, the impact of radiation in three typical urban areas of Cyprus during the summer period is examined. The first geometry considered represents a typical suburban area and is termed as the reference geometry. The second geometry represents an area at the centre of a town with higher buildings and relatively narrower roads. The third geometry, on the other hand, describes a suburban area with wider roads and larger houses than the reference model. Computed values for air temperature in the urban street canyon have indicated that the increase in temperature associated with radiative heat transfer can be reduced by optimising the canyon geometry and, ultimately, help to mitigate the human thermal discomfort. The present study has also revealed that the selection of construction materials can be optimised to offer further reductions in the air temperature of the urban environment. It can be concluded that the combined effect of these remedies can lead to reductions in the energy consumption for building air-conditioning over the summer period.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Multilayer interface tracking model of pure tungsten oxidation

Here, we present a numerical model to predict oxide scale growth on tungsten surfaces under exposure to oxygen at high temperatures. The model captures the formation of four thermodynamically-compatible oxide sublayers, WO 2 , WO 2.72 , WO 2.9 , and WO 3 , on top of the metal substrate. Oxide layer growth is simulated by tracking the oxide/oxide and oxide/metal interfaces using a sharp-interface Stefan model coupled to diffusion kinetics. The model is parameterized using selected experimental measurements and electronic structure calculations of the diffusivities of all the oxide subphases involved. We simulate oxide growth at temperatures of 600°C and above, extracting the power law growth exponents in each case, which we find to deviate from classical parabolic growth in several cases. We conduct a comparison of the model predictions with an extensive experimental data set, with reasonable agreement at most temperatures. While many gaps in our understanding still exist, this work is a first attempt at embedding the thermodynamic and kinetic complexity of tungsten oxide growth into a comprehensive mesoscale kinetic model that attempts to capture the essential features of tungsten oxidation to fill existing knowledge gaps and guide and enhance future tungsten oxidation models.

36 MATERIALS SCIENCE↗

Characterizing W sources in the all-W wall, all-RF WEST tokamak environment

In this work, experimental data, together with interpretive modeling tools, are examined to study trends in the tungsten (W) source in the all-W environment of the WEST tokamak, both from the divertor and from the main chamber. In particular, a poloidal limiter protecting an ion cyclotron resonance heating (ICRH) antenna is used as proxy for main chamber sourcing. The key study is carried out by stepping up lower hybrid current drive (LHCD) power, as the only auxiliary power source. Limiter and divertor W sources exhibit a qualitatively similar proportionality to the total power crossing the separatrix, P SEP , although the main chamber source remains substantially lower than the divertor source, for the range of P SEP accessible in the experiments. Intepretive modeling of the limiter source is carried out with a particle-in-cell (PIC) sheath model coupled to a surface sputtering model. Oxygen is used as a proxy for all light impurity species allowing for characterization of the critical W erosion regions. To get a good quantitative match to the data, it is necessary to assume that the oxygen arrives at the surface mostly at high ionization stages (4+ and above). A separate simulation with SOLEDGE-EIRENE, constrained to measured upstream scrape-off-layer plasma profiles, gives oxygen fractional abundances that are compatible with the PIC simulation result. This is understood to arise from transport processes that dominate over recombination. Substituting the LHCD by ICRH, in an equivalent experiment, the local W source exhibits a 3× enhancement. This can be matched by the simulation, by assuming local RF electric field rectification, based on ~100 eV peak-to-peak, near-antennna electric field. This work has highlighted the particular importance of understanding the ion charge state balance of light impurities as these are most likely the dominant sputtering species in fusion devices with high-Z walls.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

On the implementation of the spherical collapse model for dark energy models

In this work we review the theory of the spherical collapse model and critically analyse the aspects of the numerical implementation of its fundamental equations. By extending a recent work by [1], we show how different aspects, such as the initial integration time, the definition of constant infinity and the criterion for the extrapolation method (how close the inverse of the overdensity has to be to zero at the collapse time) can lead to an erroneous estimation (a few per mill error which translates to a few percent in the mass function) of the key quantity in the spherical collapse model: the linear critical overdensity δ{sub c}, which plays a crucial role for the mass function of halos. We provide a better recipe to adopt in designing a code suitable to a generic smooth dark energy model and we compare our numerical results with analytic predictions for the EdS and the ΛCDM models. We further discuss the evolution of δ{sub c} for selected classes of dark energy models as a general test of the robustness of our implementation. We finally outline which modifications need to be taken into account to extend the code to more general classes of models, such as clustering dark energy models and non-minimally coupled models.

79 ASTRONOMY AND ASTROPHYSICS↗

Towards an Open (Data) Science Analytics-Hub for Reproducible Multi-Model Climate Analysis at Scale

Open Science is key to future scientific research and promotes a deep transformation in the whole scientific research process encouraging the adoption of transparent and collaborative scientific approaches aimed at knowledge sharing. Open Science is increasingly gaining attention in the current and future research agenda worldwide. To effectively address Open Science goals, besides Open Access to results and data, it is also paramount to provide tools or environments to support the whole research process, in particular the design, execution and sharing of transparent and reproducible experiments, including data provenance (or lineage) tracking. This work introduces the Climate Analytics-Hub, a new component on top of the Earth System Grid Federation (ESGF), which joins big data approaches and parallel computing paradigms to provide an Open Science environment for reproducible multi-model climate change data analytics experiments at scale. An operational implementation has been set up at the SuperComputing Centre of the Euro- Mediterranean Center on Climate Change, with the main goal of becoming a reference Open Science hub in the climate community regarding the multi-model analysis based on the Coupled Model Intercomparison Project (CMIP).

analytics-hub↗

Documenting Climate Models and Their Simulations

The results of climate models are of increasing and widespread importance. No longer is climate model output of sole interest to climate scientists and researchers in the climate change impacts and adaptation fields. Now nonspecialists such as government officials, policy makers, and the general public all have an increasing need to access climate model output and understand its implications. For this host of users, accurate and complete metadata (i.e., information about how and why the data were produced) is required to document the climate modeling results. We describe a pilot community initiative to collect and make available documentation of climate models and their simulations. In an initial application, a metadata repository is being established to provide information of this kind for a major internationally coordinated modeling activity known as CMIP5 (Coupled Model Intercomparison Project, Phase 5). We expected that for a wide range of stakeholders, this and similar community-managed metadata repositories will spur development of analysis tools that facilitate discovery and exploitation of Earth system simulations.

54 ENVIRONMENTAL SCIENCES↗

Modeling the Mechanisms of Linear and Nonlinear ENSO Responses to the Pacific Meridional Mode

Interactions between the Pacific meridional mode (PMM) and El Niño–Southern Oscillation (ENSO) are investigated using the National Center for Atmospheric Research (NCAR) Community Earth System Model (CESM) and an intermediate coupled model (ICM). The two models are configured so that the CESM simulates the PMM but not ENSO, and the ICM simulates ENSO but not the PMM, allowing for a clean separation between the PMM evolution and the subsequent ENSO response. An ensemble of CESM simulations is run with an imposed surface heat flux associated with the North Pacific Oscillation (NPO) generating a sea surface temperature (SST) and wind response representative of the PMM. The PMM wind is then applied as a forcing to the ICM to simulate the ENSO response. The positive (negative) ensemble-mean PMM wind forcing results in a warm (cold) ENSO event although the responses are not symmetric (warm ENSO events are larger in amplitude than cold ENSO events), and large variability between ensemble members suggests that any individual ENSO event is strongly influenced by natural variability contained within the CESM simulations. Sensitivity experiments show that 1) direct forcing of Kelvin waves by PMM winds dominates the ENSO response, 2) seasonality of PMM forcing and ENSO growth rates influences the resulting ENSO amplitude, 3) ocean dynamics within the ICM dominate the ENSO asymmetry, and 4) the nonlinear relationship between PMM wind anomalies and surface wind stress may enhance the La Niña response to negative PMM variations. Implications for ENSO variability are discussed.

54 ENVIRONMENTAL SCIENCES↗

Blocking Simulations in GFDL GCMs for CMIP5 and CMIP6

Abstract The frequency of atmospheric blocking has been largely underestimated by general circulation models (GCMs) participating in the Coupled Model Intercomparison Project (CMIP). Errors in the onset, persistence, barotropicity, geographical preference, seasonality, intensity, and moving speed of global blocking were diagnosed in 10 Geophysical Fluid Dynamics Laboratory (GFDL) GCMs for recent CMIP5 and CMIP6 using a detection approach that combines zonal eddies and the reversal of zonal winds. The blocking frequency, similar at 500 and 250 hPa, is underestimated by 50% in the Atlantic–Europe region during December–February but is overestimated by 60% in the Pacific–North America region during that season and by 70% in the southwest Pacific during July–August. These blocking biases at 500 hPa were investigated in the five CMIP6 models that showed improvements over the CMIP5 versions. The Atlantic–Europe underestimate corresponds to lower instantaneous blocking rates, lower persistent blocking rates, and higher persistent stationary ridge rates; the number of blocks with a duration of 4–5 days is only 40%–65% of that in observations. In contrast, the overestimate consists of excessive blocks with a duration longer than 12 days in the Pacific–North America and up to twice as many 4–6-day events in the southwest Pacific. Simulated December–February blocks up to 12 days in the Pacific–North America region tend to be stronger and to move more slowly than those in observations. Diagnostic sensitivity tests indicated that the zonal mean and zonal eddy components of the mean state play a key role, as replacing each with that of observations substantially reduced many of the outstanding biases in these GCMs.

54 ENVIRONMENTAL SCIENCES↗

Next Generation Climate Change Experiments Needed to Advance Knowledge and for Assessment of CMIP6

The Aspen Global Change Institute hosted a technical science workshop entitled, “Next generation climate change experiments needed to advance knowledge and for assessment of CMIP6,” on August 4-9, 2013 in Aspen, CO. Jerry Meehl (NCAR), Richard Moss (PNNL), and Karl Taylor (LLNL) served as co-chairs for the workshop which included the participation of 32 scientists representing most of the major climate modeling centers for a total of 160 participant days. In August 2013, AGCI gathered a high level meeting of representatives from major climate modeling centers around the world to assess achievements and lessons learned from the most recent generation of coordinated modeling experiments known as the Coupled Model Intercomparison Project – 5 (CMIP5) as well as to scope out the science questions and coordination structure desired for the next anticipated phase of modeling experiments called CMIP6. The workshop allowed for reflection on the coordination of the CMIP5 process as well as intercomparison of model results, such as were assessed in the most recent IPCC 5th Assessment Report, Working Group 1. For example, this slide from Masahiro Watanabe examines performance on a range of models capturing Atlantic Meridional Overturning Circulation (AMOC).

54 ENVIRONMENTAL SCIENCES↗

Quantifying the effect of nighttime interactions between roots and canopy physiology and their control of water and carbon cycling on feedbacks between soil moisture and terrestrial climatology under variable environmental conditions

The primary objective of this project is to characterize and quantify how the temporal variability of hydraulic redistribution (HR) and its physiological regulation in unmanaged and complex forests is affecting current water and carbon exchange and predict how future climate scenarios will affect these relationships and potentially feed back to the climate. Specifically, a detailed study of ecosystem water uptake and carbon exchange in relation to root functioning was proposed in order to quantify the mechanisms controlling temporal variability of soil moisture dynamic and HR in three active AmeriFlux sites, and to use published data of two other inactive AmeriFlux sites. Furthermore, data collected by our research group at the Duke Free Air CO2 enrichment (FACE) site was also being utilized to further improve our ability to forecast future environmental impacts of elevated CO2 concentration on soil moisture dynamic and its effect on carbon sequestration and terrestrial climatology. The overarching objective being to forecast, using a soil:plant:atmosphere model coupled with a biosphere:atmosphere model, the impact of root functioning on land surface climatology. By comparing unmanaged sites to plantations, we also proposed to determine the effect of land use change on terrestrial carbon sequestration and climatology through its effect on soil moisture dynamic and HR. Our simulations of HR by roots indicated that in some systems HR is an important mechanism that buffers soil water deficit, affects energy and carbon cycling; thus having significant implications for seasonal climate. HR maintained roots alive and below 70% loss of conductivity and our simulations also showed that the increased vapor pressure deficit at night under future conditions was sufficient to drive significant nighttime transpiration at all sites, which reduced HR. This predicted reduction in HR under future climate conditions played an important regulatory role in land atmosphere interactions by affecting whole ecosystem carbon and water balance. Under future climatic scenarios, HR was reduced thus affecting negatively plant water use and carbon assimilation. The discrepancy between the predicted and actual surface warming and atmospheric water vapor caused by the persistence of evapotranspiration during the dry season, increasing energy transfer in the form of latent heat. Under those simulations, we also evaluated how the hydraulic properties of soil and xylem limited the rate of carbon uptake, and carbon net ecosystem exchange. The multilayered hydraulically driven soil vegetation atmosphere carbon and water transfer model was designed to represent processes common to vascular plants, so that ecosystem atmosphere exchange could be captured by the same processes at different sites. Those models shown to be well suited for investigating the impact of drought on forest ecosystems because of its explicit treatment of water transport to leaves. This modeling work also confirmed that unmanaged, mixed hardwood site are more resilient to climatic variations than an adjacent pine plantation, but that future climatic conditions will reverse this trends.

54 ENVIRONMENTAL SCIENCES↗

Mid-Latitude Circulation and Extremes in a Changing Climate

Mid-latitude extreme weather events are responsible for a large part of climate-related damage. Yet large uncertainties remain in climate model projections of heat waves, droughts, and heavy rain/snow events on regional scales, limiting our ability to effectively use these projections for climate adaptation and mitigation. These uncertainties can be attributed to both the lack of spatial resolution in the models, and to the lack of a dynamical understanding of these extremes. The approach of this project is to relate the fine-scale features to the large scales in current climate simulations, seasonal re-forecasts, and climate change projections in a very wide range of models, including the atmospheric and coupled models of ECMWF over a range of horizontal resolutions (125 to 10 km), aqua-planet configuration of the Model for Prediction Across Scales and High Order Method Modeling Environments (resolutions ranging from 240 km – 7.5 km) with various physics suites, and selected CMIP5 model simulations. The large scale circulation will be quantified both on the basis of the well tested preferred circulation regime approach, and very recently developed measures, the finite amplitude Wave Activity (FAWA) and its spectrum. The fine scale structures related to extremes will be diagnosed following the latest approaches in the literature. The goal is to use the large scale measures as indicators of the probability of occurrence of the finer scale structures, and hence extreme events. These indicators will then be applied to the CMIP5 models and time-slice projections of a future climate.

54 ENVIRONMENTAL SCIENCES↗

Evaluation of a Potential Groundwater Tracer Test in the Ringold Upper Mud Aquifer at the 100-H Area of the Hanford Site

This document presents an evaluation of a potential groundwater tracer test in the uppermost RUM (Ringold Formation member of Wooded Island – upper mud unit) aquifer in the 100-H Area of the Hanford Site. The results from field electrical resistivity tomography (ERT) investigations, numerical flow and transport (F&T) modeling, coupled ERT-F&T model simulations, and laboratory resin experiments provide the technical basis for evaluating tracer test scenarios that minimize impacts to the 100-HX pump and treat (HX P&T) operations, and maximize hydrologic data acquisitions including the use of advanced surface geophysical techniques for monitoring tracer transport.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

CMIP6-based Multi-model Streamflow Projections over the Conterminous US

This dataset presents an ensemble of streamflow projections based on the hydroclimate projections dataset supporting the SECURE Water Act Section 9505 Assessment for the US Department of Energy (DOE) Water Power Technologies Office (WPTO). The six-member General Climate Model (GCM) ensemble from the Coupled Models Intercomparison Project phase 6 (CMIP6) downscaled using statistical (i.e., DBCCA) and dynamical (i.e., RegCM) and bias-corrected using two meteorological reference observations (Daymet & Livneh) are driven through two calibrated hydrologic models (VIC & PRMS) to simulate projected future hydrologic responses. This leads to the production of an ensemble of hydroclimate projections, including total runoff (surface runoff and baseflow) for 1980–2019 baseline and 2020–2059 near-term future periods under multiple emission scenarios at 1/24° (~4 km) spatial resolution across the CONUS. The total runoff projections are routed through the Routing Application for Parallel computatIon of Discharge (RAPID) routing model to produce an ensemble of streamflow projections for both periods across 2.7 million NHDPlusV2 stream reaches in the CONUS.

13 HYDRO ENERGY↗

CMIP6-based Multi-model Hydroclimate Projection over the Conterminous US

We present a suite of high-resolution downscaled hydro-climate projections over the conterminous United States (CONUS) based on a six-member General Climate Model (GCM) ensemble from the Coupled Models Intercomparison Project phase 6 (CMIP6). The CMIP6 GCMs are downscaled using two different downscaling approaches (statistical-based DBCCA & dynamical-based RegCM) based on two meteorological reference observations (Daymet & Livneh), and then fed to two calibrated hydrologic models (VIC & PRMS) to simulate projected future hydrologic responses. Each ensemble member covers 1980–2019 in the historic period and 2020–2059 in the near-term future period under the high-end (SSP585) emission scenario. Major variables such as daily maximum temperature (tmax), daily minimum temperature (tmin), daily total precipitation (prcp), daily average wind speed (wind), and daily total runoff (runoff) at 1/24° (~4 km) spatial resolution across the CONUS are provided through this data portal. This dataset is derived to support the SECURE Water Act Section 9505 Assessment for the US Department of Energy (DOE) Water Power Technologies Office (WPTO). Further details of this dataset can be referred to Kao et al. (2022) and Rastogi et al. (2022).

13 HYDRO ENERGY↗

Evaluating Uncertainties in Marine Biogeochemical Models: Benchmarking Aerosol Precursors

The effort to accurately estimate global radiative forcing has long been hampered by a degree of uncertainty in the tropospheric aerosol contribution. Reducing uncertainty in natural aerosol processes, the baseline of the aerosol budget, thus becomes a fundamental task. The appropriate representation of aerosols in the marine boundary layer (MBL) is essential to reduce uncertainty and provide reliable information on offsets to global warming. We developed an International Ocean Model Benchmarking package to assess marine biogeochemical process representations in Earth System Models (ESMs), and the package was employed to evaluate surface ocean concentrations and the sea–air fluxes of dimethylsulfide (DMS). Model performances were scored based on how well they captured the distribution and variability contained in high-quality observational datasets. Results show that model-data biases increased as DMS enters the MBL, but unfortunately over three-quarters of the models participating in the fifth Coupled Model Intercomparison Project (CMIP5) did not have a dynamic representation of DMS. When it is present, models tend to over-predict sea surface concentrations in the productive region of the eastern tropical Pacific by almost a factor of two, and the sea–air fluxes by a factor of three. Systematic model-data benchmarking as described here will help to identify such deficiencies and subsequently lead to improved subgrid-scale parameterizations and ESM development.

58 GEOSCIENCES↗

EVOLUTIONARY TRACKS OF THE CLIMATE OF EARTH-LIKE PLANETS AROUND DIFFERENT MASS STARS

The climatic evolution of the Earth depends strongly on the evolution of the insolation from the Sun and the amount of the greenhouse gasses, especially CO{sub 2} in the atmosphere. Here, we investigate the evolution of the climate of hypothetical Earths around stars whose masses are different from the solar mass with a luminosity evolution model of the stars, a mantle degassing model coupled with a parameterized convection model of the planetary interiors, and an energy balance climate model of the planetary surface. In the habitable zone (HZ), the climate of the planets is initially warm or hot, depending on the orbital semimajor axes. We found that, in the inner HZ, the climate of the planets becomes hotter with time owing to the increase in the luminosity of the central stars, while, in the outer HZ, it becomes colder and eventually globally ice-covered owing to the decrease in the CO{sub 2} degassing rate of the planets. The orbital condition for maintaining the warm climate similar to the present Earth becomes very limited, and more interestingly, the planet orbiting in the outer HZ becomes globally ice-covered after a certain critical age (∼3 Gyr for the hypothetical Earth with standard parameters), irrespective of the mass of the central star. This is because the critical age depends on the evolution of the planets and planetary factors, rather than on the stellar mass. The habitability of the Earth-like planet is shown to be limited with age even though it is orbiting within the HZ.

79 ASTRONOMY AND ASTROPHYSICS↗

Modulation of springtime surface sensible heating over the Tibetan Plateau on the interannual variability of East Asian dust cycle

Previous observational evidence and numerical simulations have revealed that the surface sensible heating in spring (March–April–May, MAM) over the Tibetan Plateau (TPSH) can affect the Asian regional hydrological cycle, surface energy balance, and climate through altering atmospheric heat source of the Tibetan Plateau (TP). This study aims to investigate the impacts of MAM TPSH on the interannual variability of East Asian dust cycle by using CAM4-BAM (version 4 of the Community Atmosphere Model coupled to a bulk aerosol model), MERRA-2 (version 2 of the Modern-Era Retrospective analysis for Research and Applications) surface dust concentration, and TPSH measurements. Our simulations show that the surface dust concentrations over the East Asian (EA) dust source region and over the northwestern Pacific (NP) in MAM are significantly positively correlated with TPSH, with regionally averaged correlation coefficients of 0.49 for EA and 0.44 for NP. Similar positive correlations are also shown between the MAM TPSH measurements averaged over the 73 observation sites and the surface dust concentration from MERRA-2. Simulation-based comparisons between strongest and weakest TPSH years reveal that, the MAM surface dust concentration in the strongest TPSH years increases with relative differences of 13.1 % over EA and 36.9 % over NP. These corresponding differences are found in MERRA-2 with 22.9 % and 13.3 % over EA and NP, respectively. Further simulated results show that the processes of whole dust cycles (e.g., dust loading, emission, and transport, as well as dust deposition) are also significantly enhanced during the strongest TPSH years over EA and NP. Through enhancing the TP heat source, stronger TPSH in MAM generates an anticyclonic anomaly in middle and upper troposphere over the TP and over the downstream Pacific region, respectively. These atmospheric circulation anomalies induced by the increased TPSH result in increasing the westerly winds over both EA and NP, which in turn increases dust emissions over the dust source, and dust transport over these two regions, as well as the regional dust cycles. These results suggest that addressing the East Asian dust changes in the future climates requires considering not only increasing greenhouse gas emissions but also the variations of the TP's heat source under global warming.

54 ENVIRONMENTAL SCIENCES↗