Search NASA⌕ Search

SEARCH · Search NASA

Results for “CLIMATE 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 487 records · Page 27

Reconciling divergent trends and millennial variations in Holocene temperatures

Cooling during most of the past two millennia has been widely recognized1,2 and has been inferred to be the dominant global temperature trend of the past 11,700 years (the Holocene epoch)3. However, long-term cooling has been difficult to reconcile with global forcing4, and climate models consistently simulate long-term warming4. The divergence between simulations and reconstructions emerges primarily for northern mid-latitudes, for which pronounced cooling has been inferred from marine and coastal records using multiple approaches3. Here we show that temperatures reconstructed from sub-fossil pollen from 642 sites across North America and Europe closely match simulations, and that long-term warming, not cooling, defined the Holocene until around 2,000 years ago. The reconstructions indicate that evidence of long-term cooling was limited to North Atlantic records. Early Holocene temperatures on the continents were more than two degrees Celsius below those of the past two millennia, consistent with the simulated effects of remnant ice sheets in the climate model Community Climate System Model 3 (CCSM3)5. CCSM3 simulates increases in ‘growing degree days’—a measure of the accumulated warmth above five degrees Celsius per year—of more than 300 kelvin days over the Holocene, consistent with inferences from the pollen data. It also simulates a decrease in mean summer temperatures of more than two degrees Celsius, which correlates with reconstructed marine trends and highlights the potential importance of the different subseasonal sensitivities of the records. Despite the differing trends, pollen- and marine-based reconstructions are correlated at millennial-to-centennial scales, probably in response to ice-sheet and meltwater dynamics, and to stochastic dynamics similar to the temperature variations produced by CCSM3. Although our results depend on a single source of palaeoclimatic data (pollen) and a single climate-model simulation, they reinforce the notion that climate models can adequately simulate climates for periods other than the present-day. They also demonstrate that amplified warming in recent decades increased temperatures above the mean of any century during the past 11,000 years.

Science & Technology - Other Topics↗

WATER TRAPPING ON TIDALLY LOCKED TERRESTRIAL PLANETS REQUIRES SPECIAL CONDITIONS

Surface liquid water is essential for standard planetary habitability. Calculations of atmospheric circulation on tidally locked planets around M stars suggest that this peculiar orbital configuration lends itself to the trapping of large amounts of water in kilometers-thick ice on the night side, potentially removing all liquid water from the day side where photosynthesis is possible. We study this problem using a global climate model including coupled atmosphere, ocean, land, and sea ice components as well as a continental ice sheet model driven by the climate model output. For a waterworld, we find that surface winds transport sea ice toward the day side and the ocean carries heat toward the night side. As a result, nightside sea ice remains O(10 m) thick and nightside water trapping is insignificant. If a planet has large continents on its night side, they can grow ice sheets O(1000 m) thick if the geothermal heat flux is similar to Earth's or smaller. Planets with a water complement similar to Earth's would therefore experience a large decrease in sea level when plate tectonics drives their continents onto the night side, but would not experience complete dayside dessiccation. Only planets with a geothermal heat flux lower than Earth's, much of their surface covered by continents, and a surface water reservoir O(10%) of Earth's would be susceptible to complete water trapping.

79 ASTRONOMY AND ASTROPHYSICS↗

Characterizing Sierra Nevada Snowpack Using Variable-Resolution CESM

The location, timing, and intermittency of precipitation in California make the state integrally reliant on winter-season snowpack accumulation to maintain its economic and agricultural livelihood. Of particular concern is that winter-season snowpack has shown a net decline across the western United States over the past 50 years, resulting in major uncertainty in water-resource management heading into the next century. Cutting-edge tools are available to help navigate and preemptively plan for these uncertainties. This paper uses a next-generation modeling technique-variable-resolution global climate modeling within the Community Earth System Model (VR-CESM)-at horizontal resolutions of 0.125° (14 km) and 0.25° (28 km). VR-CESM provides the means to include dynamically large-scale atmosphere-ocean drivers, to limit model bias, and to provide more accurate representations of regional topography while doing so in a more computationally efficient manner than can be achieved with conventional general circulation models. This paper validates VR-CESM at climatological and seasonal time scales for Sierra Nevada snowpack metrics by comparing them with the "Daymet," "Cal-Adapt," NARR, NCEP, and North American Land Data Assimilation System (NLDAS) reanalysis datasets, the MODIS remote sensing dataset, the SNOTEL observational dataset, a standard-practice global climate model (CESM), and a regional climate model (WRF Model) dataset. Overall, given California's complex terrain and intermittent precipitation and that both of the VR-CESM simulations were only constrained by prescribed sea surface temperatures and data on sea ice extent, a 0.68 centered Pearson product-moment correlation, a negative mean SWE bias of < 7 mm, an interquartile range well within the values exhibited in the reanalysis datasets, and a mean December-February extent of snow cover that is within 7% of the expected MODIS value together make apparent the efficacy of the VR-CESM framework.

54 ENVIRONMENTAL SCIENCES↗

Fast domain-aware neural network emulation of a planetary boundary layer parameterization in a numerical weather forecast model

Abstract. Parameterizations for physical processes in weather and climate models are computationally expensive. We use model output from the Weather Research Forecast (WRF) climate model to train deep neural networks (DNNs) and evaluate whether trained DNNs can provide an accurate alternative to the physics-based parameterizations. Specifically, we develop an emulator using DNNs for a planetary boundary layer (PBL) parameterization in the WRF model. PBL parameterizations are used in atmospheric models to represent the diurnal variation in the formation and collapse of the atmospheric boundary layer – the lowest part of the atmosphere. The dynamics and turbulence, as well as velocity, temperature, and humidity profiles within the boundary layer are all critical for determining many of the physical processes in the atmosphere. PBL parameterizations are used to represent these processes that are usually unresolved in a typical climate model that operates at horizontal spatial scales in the tens of kilometers. We demonstrate that a domain-aware DNN, which takes account of underlying domain structure (e.g., nonlocal mixing), can successfully simulate the vertical profiles within the boundary layer of velocities, temperature, and water vapor over the entire diurnal cycle. Results also show that a single trained DNN from one location can produce predictions of wind speed, temperature, and water vapor profiles over nearby locations with similar terrain conditions with correlations higher than 0.9 when compared with the WRF simulations used as the training dataset.

54 ENVIRONMENTAL SCIENCES↗

Quantification of marine aerosol subgrid variability and its correlation with clouds based on high‐resolution regional modeling

Abstract One limitation of most global climate models (GCMs) is that with the horizontal resolutions they typically employ, they cannot resolve the subgrid variability (SGV) of clouds and aerosols, adding extra uncertainties to the aerosol radiative forcing estimation. To inform the development of an aerosol subgrid variability parameterization, here we analyze the aerosol SGV over the southern Pacific Ocean simulated by the high‐resolution Weather Research and Forecasting model coupled to Chemistry. We find that within a typical GCM grid, the aerosol mass subgrid standard deviation is 15% of the grid‐box mean mass near the surface on a 1 month mean basis. The fraction can increase to 50% in the free troposphere. The relationships between the sea‐salt mass concentration, meteorological variables, and sea‐salt emission rate are investigated in both the clear and cloudy portion. Under clear‐sky conditions, marine aerosol subgrid standard deviation is highly correlated with the standard deviations of vertical velocity, cloud water mixing ratio, and sea‐salt emission rates near the surface. It is also strongly connected to the grid box mean aerosol in the free troposphere (between 2 km and 4 km). In the cloudy area, interstitial sea‐salt aerosol mass concentrations are smaller, but higher correlation is found between the subgrid standard deviations of aerosol mass and vertical velocity. Additionally, we find that decreasing the model grid resolution can reduce the marine aerosol SGV but strengthen the correlations between the aerosol SGV and the total water mixing ratio (sum of water vapor, cloud liquid, and cloud ice mixing ratios).

Lin, Guangxing↗

Assessing Clouds Using Satellite Observations Through Three Generations of Global Atmosphere Models

Abstract Clouds are parameterized in climate models using quantities on the model grid‐scale to approximate the cloud cover and impact on radiation. Because of the complexity of processes involved with clouds, these parameterizations are one of the key challenges in climate modeling. Differences in parameterizations of clouds are among the main contributors to the spread in climate sensitivity across models. In this work, the clouds in three generations of an atmosphere model lineage are evaluated against satellite observations. Satellite simulators are used within the model to provide an appropriate comparison with individual satellite products. In some respects, especially the top‐of‐atmosphere cloud radiative effect, the models show generational improvements. The most recent generation, represented by two distinct branches of development, exhibits some regional regressions in the cloud representation; in particular the southern ocean shows a positive bias in cloud cover. The two branches of model development show how choices during model development, both structural and parametric, lead to different cloud climatologies. Several evaluation strategies are used to quantify the spatial errors in terms of the large‐scale circulation and the cloud structure. The Earth mover's distance is proposed as a useful error metric for the passive satellite data products that provide cloud‐top pressure‐optical depth histograms. The cloud errors identified here may contribute to the high climate sensitivity in the Community Earth System Model, version 2 and in the Energy Exascale Earth System Model, version 1.

54 ENVIRONMENTAL SCIENCES↗

Evaluating and Improving Convective Parameterization for GCMs Using ARM Observations. Final report

This technical report summarizes the achievements during the project period funded by the Atmospheric System Research (ASR) program. The general goal of the project is to improve the representation of atmospheric deep convection in global climate models. By using DOE Atmospheric Radiation Measurement (ARM) program observations and a hierarchy of models from cloud-resolving to global climate models, we evaluated and improved many aspects of convection parameterization schemes, such as the trigger function for convection and closure condition to determine the amount of convection. We also incorporated the stochasticity of convection into a convection scheme. These improvements, when implemented into the National Center for Atmospheric Research (NCAR) Community Earth System Model (CESM) and the DOE Exascale Energy Earth System Model (E3SM), led to improved simulation of precipitation characteristics, including precipitation intensity, Intertropical Convergence Zone, and Madden-Julian oscillation. The project also improved the representation of microphysical processes in convective clouds. It enables the interaction of cloud microphysics with aerosols in weather and climate models.

54 ENVIRONMENTAL SCIENCES↗

Improved representation of black carbon mixing structures suggests stronger direct radiative heating

Black carbon significantly influences the Earth system because of its strong solar radiation absorption. However, its direct radiative effect remains poorly understood in current climate models, partly because current climate models oversimplify the diverse structures formed when black carbon mixes with other atmospheric components. Here we show that incorporating more realistic, multi-mixing-structure representations of black carbon increases the direct radiative effect. We find that aged black carbon particles, with thicker coatings and higher embedded fractions, enhance the direct radiative effect more efficiently. Using machine learning alongside the Community Earth System Model, we show that the direct radiative effect at the top of the atmosphere in regions with heavy black carbon pollution is 31.6% greater when multi-mixing structures are considered. These findings highlight the importance of modeling complex mixing structures of particle-resolved black carbon to accurately capture their warming impacts on global atmosphere, particularly in highly polluted regions.

54 ENVIRONMENTAL SCIENCES↗

Observations and Projections of Heat Waves in South America

Although Heat Waves (HWs) are expected to increase due to global warming, they are a regional phenomenon that demands for local analyses. In this paper, we assess four HW metrics (HW duration, HW frequency, HW amplitude, and number of HWs per season) as well as the share of extremely warm days (TX95, according to the 95th percentile) in South America (SA). Our analysis included observations as well as simulations from global and regional models. In particular, Regional Climate Models (RCMs) from the Coordinated Regional Climate Downscaling Experiment (CORDEX), and Global Climate Models (GCMs) from the Coupled Model Intercomparison Project Phase 5 (CMIP5) were used to project both TX95 estimates and HW metrics according to two representative concentration pathways (RCP4.5 and RCP8.5). We found that in recent decades the share of extremely warm days has at least doubled over the period December–January–February (DJF) in northern SA; less significant increases have been observed in southern SA. We also found that by midcentury, under the RCP4.5 scenario, extremely warm DJF days (as well as the number of HWs per season) are expected to increase by 5–10 times at locations close to the Equator and in the Atacama Desert. Increases are expected to be less pronounced in southern SA. Projections under the RCP8.5 scenario are more striking, particularly in tropical areas where half or more of the days could be extremely warm by midcentury.

54 ENVIRONMENTAL SCIENCES↗

Bias Correction and Statistical Downscaling of Solar Radiation Using NA-CORDEX and the NSRDB

The current state-of-art for estimating long-term PV production uses long-term estimates of solar radiation variables, such as global horizontal irradiance (GHI), from previous years. This data is used in models such as the System Advisor Model (SAM) or PYSyst to predict annual production for a PV plant. This information is then used to estimate the production over the next 20 years (a typical plant lifetime) under the assumption that the variability over the current period is representative of the future. As the PV industry moves to extend plant lifetimes to 50 years the current assumptions of representativeness of weather may not be appropriate. This is especially true as our climate changes rapidly. To assess long-term PV production, future projections for solar radiation based on projected carbon emissions are readily available in regional and global climate models. However, climate model projections contain inherent biases that may need to be corrected for accurate analysis of future projections of climate variables. Several studies have analyzed projections of solar radiation for future years, however the accuracy of the model output compared to current and historic data has not been widely studied. Chen (2021) showed that available climate models do not accurately represent solar radiation in some cases, over-projecting GHI at the surface while under-projecting its obstructions, such as clouds and aerosols. This works aims to (1) increase understanding of the accuracy of solar radiation currently available in global and regional climate models and (2) implement bias correction through linear models based on reanalysis data compared to observed solar radiation. The latter aim will be conducted using available observed solar radiation data and modeled data from several regional climate models (RCMs). The bias correction method will be applied to projections of solar radiation resulting in a more accurate representation of the future of solar production.

climate data↗

Interactive Correlation Analysis and Visualization of Climate Data

The relationship between our ability to analyze and extract insights from visualization of climate model output and the capability of the available resources to make those visualizations has reached a crisis point. The large volume of data currently produced by climate models is overwhelming the current, decades-old visualization workflow. The traditional methods for visualizing climate output also have not kept pace with changes in the types of grids used, the number of variables involved, and the number of different simulations performed with a climate model or the feature-richness of high-resolution simulations. This project has developed new and faster methods for visualization in order to get the most knowledge out of the new generation of high-resolution climate models. While traditional climate images will continue to be useful, there is need for new approaches to visualization and analysis of climate data if we are to gain all the insights available in ultra-large data sets produced by high-resolution model output and ensemble integrations of climate models such as those produced for the Coupled Model Intercomparison Project. Towards that end, we have developed new visualization techniques for performing correlation analysis. We have also introduced highly scalable, parallel rendering methods for visualizing large-scale 3D data. This project was done jointly with climate scientists and visualization researchers at Argonne National Laboratory and NCAR.

54 ENVIRONMENTAL SCIENCES↗

Fusing time-varying mosquito data and continuous mosquito population dynamics models

Climate change is arguably one of the most pressing issues affecting the world today and requires the fusion of disparate data streams to accurately model its impacts. Mosquito populations respond to temperature and precipitation in a nonlinear way, making predicting climate impacts on mosquito-borne diseases an ongoing challenge. Data-driven approaches for accurately modeling mosquito populations are needed for predicting mosquito-borne disease risk under climate change scenarios. Many current models for disease transmission are continuous and autonomous, while mosquito data is discrete and varies both within and between seasons. This study uses an optimization framework to fit a non-autonomous logistic model with periodic net growth rate and carrying capacity parameters for 15 years of daily mosquito time-series data from the Greater Toronto Area of Canada. The resulting parameters accurately capture the inter-annual and intra-seasonal variability of mosquito populations within a single geographic region, and a variance-based sensitivity analysis highlights the influence each parameter has on the peak magnitude and timing of the mosquito season. This method can easily extend to other geographic regions and be integrated into a larger disease transmission model. This method addresses the ongoing challenges of data and model fusion by serving as a link between discrete time-series data and continuous differential equations for mosquito-borne epidemiology models.

97 MATHEMATICS AND COMPUTING↗

Important ice processes are missed by the Community Earth System Model (CESM) in Southern Ocean mixed–phase clouds: Bridging SOCRATES observations to model developments

Global climate models (GCMs) are challenged by difficulties in simulating cloud phase and cloud radiative effect over the Southern Ocean (SO). Some of the new generation GCMs predict too much liquid and too little ice in mixed-phase clouds. This misrepresentation of cloud phase in GCMs results in weaker negative cloud feedback over the SO and a higher climate sensitivity. Based on a model comparison with observational data obtained during the Southern Ocean Cloud Radiation and Aerosol Transport Experimental Study (SOCRATES), this study addresses a key uncertainty in the Community Earth System Model version 2 (CESM2) related to cloud phase, namely ice formation in pristine remote SO clouds. It is found that sea spray organic aerosols (SSOA) are the most important type of ice nucleating particles (INPs) over the SO with concentrations one order of magnitude higher than those of dust INPs based on measurements and CESM2 simulations. Secondary ice production (SIP) which includes riming splintering, rain droplet shattering, and ice-ice collisional fragmentation as implemented in CESM2 is the dominant ice production process in moderately cold clouds with cloud temperatures greater than –20°C. SIP enhances the in-cloud ice number concentrations (Ni) by 1-3 orders of magnitude and predicts more mixed-phase (with percentage occurrence increased from 15% to 21%), in better agreement with the observations. Finally, this study highlights the importance of accurately representing cloud phase over the pristine remote SO by considering the ice nucleation of SSOA and SIP processes, which are currently missing in most GCM cloud microphysics parameterizations.

54 ENVIRONMENTAL SCIENCES↗

A unified parameterization of clouds and turbulence using CLUBB and subcolumns in the Community Atmosphere Model

Most global climate models parameterize separate cloud types using separate parameterizations. This approach has several disadvantages, including obscure interactions between parameterizations and inaccurate triggering of cumulus parameterizations. Alternatively, a unified cloud parameterization uses one equation set to represent all cloud types. Such cloud types include stratiform liquid and ice cloud, shallow convective cloud, and deep convective cloud. Vital to the success of a unified parameterization is a general interface between clouds and microphysics. One such interface involves drawing Monte Carlo samples of subgrid variability of temperature, water vapor, cloud liquid, and cloud ice, and feeding the sample points into a microphysics scheme.This study evaluates a unified cloud parameterization and a Monte Carlo microphysics interface that has been implemented in the Community Atmosphere Model (CAM) version 5.3. Results describing the mean climate and tropical variability from global simulations are presented. The new model shows a degradation in precipitation skill but improvements in short-wave cloud forcing, liquid water path, long-wave cloud forcing, precipitable water, and tropical wave simulation. Also presented are estimations of computational expense and investigation of sensitivity to number of subcolumns.

58 GEOSCIENCES↗

Pacific decadal oscillation remotely forced by the equatorial Pacific and the Atlantic Oceans

Abstract The Pacific Decadal Oscillation (PDO), the leading mode of Pacific decadal sea surface temperature variability, arises mainly from combinations of regional air-sea interaction within the North Pacific Ocean and remote forcing, such as from the tropical Pacific and the Atlantic. Because of such a combination of mechanisms, a question remains as to how much PDO variability originates from these regions. To better understand PDO variability, the equatorial Pacific and the Atlantic impacts on the PDO are examined using several 3-dimensional partial ocean data assimilation experiments conducted with two global climate models: the CESM1.0 and MIROC3.2m. In these partial assimilation experiments, the climate models are constrained by observed temperature and salinity anomalies, one solely in the Atlantic basin and the other solely in the equatorial Pacific basin, but are allowed to evolve freely in other regions. These experiments demonstrate that, in addition to the tropical Pacific’s role in driving PDO variability, the Atlantic can affect PDO variability by modulating the tropical Pacific climate through two proposed processes. One is the equatorial pathway, in which tropical Atlantic sea surface temperature (SST) variability causes an El Niño-like SST response in the equatorial Pacific through the reorganization of the global Walker circulation. The other is the north tropical pathway, where low-frequency SST variability associated with the Atlantic Multidecadal Oscillation induces a Matsuno-Gill type atmospheric response in the tropical Atlantic-Pacific sectors north of the equator. These results provide a quantitative assessment suggesting that 12–29% of PDO variance originates from the Atlantic Ocean and 40–44% from the tropical Pacific. The remaining 27–48% of the variance is inferred to arise from other processes such as regional ocean-atmosphere interactions in the North Pacific and possibly teleconnections from the Indian Ocean.

Meteorology & Atmospheric Sciences↗

Multifidelity Monte Carlo estimation for efficient uncertainty quantification in climate-related modeling

Abstract. Uncertainties in an output of interest that depends on the solution of a complex system (e.g., of partial differential equations with random inputs) are often, if not nearly ubiquitously, determined in practice using Monte Carlo (MC) estimation. While simple to implement, MC estimation fails to provide reliable information about statistical quantities (such as the expected value of the output of interest) in application settings such as climate modeling, for which obtaining a single realization of the output of interest is a costly endeavor. Specifically, the dilemma encountered is that many samples of the output of interest have to be collected in order to obtain an MC estimator that has sufficient accuracy – so many, in fact, that the available computational budget is not large enough to effect the number of samples needed. To circumvent this dilemma, we consider using multifidelity Monte Carlo (MFMC) estimation which leverages the use of less costly and less accurate surrogate models (such as coarser grids, reduced-order models, simplified physics, and/or interpolants) to achieve, for the same computational budget, higher accuracy compared to that obtained by an MC estimator – or, looking at it another way, an MFMC estimator obtains the same accuracy as the MC estimator at lower computational cost. The key to the efficacy of MFMC estimation is the fact that most of the required computational budget is loaded onto the less costly surrogate models so that very few samples are taken of the more expensive model of interest. We first provide a more detailed discussion about the need to consider an alternative to MC estimation for uncertainty quantification. Subsequently, we present a review, in an abstract setting, of the MFMC approach along with its application to three climate-related benchmark problems as a proof-of-concept exercise.

54 ENVIRONMENTAL SCIENCES↗

The effect of temperature on photosystem II efficiency across plant functional types and climate

Modeling terrestrial gross primary productivity (GPP) is central to predicting the global carbon cycle. Much interest has been focused on the environmentally induced dynamics of photosystem energy partitioning and how improvements in the description of such dynamics assist the prediction of light reactions of photosynthesis and therefore GPP. The maximum quantum yield of photosystem II (Φ PSIImax ) is a key parameter of the light reactions that influence the electron transport rate needed for supporting the biochemical reactions of photosynthesis. Φ PSIImax is generally treated as a constant in biochemical photosynthetic models even though a constant Φ PSIImax is expected only for non-stressed plants. We synthesized reported Φ PSIImax values from pulse-amplitude-modulated fluorometry measurements in response to variable temperatures across the globe. We found that Φ PSIImax is strongly affected by prevailing temperature regimes with declined values in both hot and cold conditions. To understand the spatiotemporal variability in Φ PSIImax , we analyzed the temperature effect on Φ PSIImax across plant functional type (PFT) and habitat climatology. The analysis showed that temperature's impact on Φ PSIImax is shaped more by climate than by PFT for plants with broad latitudinal distributions or in regions with extreme temperature variability. There is a trade-off between the temperature range within which Φ PSIImax remains maximal and the overall rate of decline of Φ PSIImax outside the temperature range such that species cannot be simultaneously tolerant and resilient to extreme temperatures. Our study points to a quantitative approach for improving electron transport and photosynthetic productivity modeling under changing climates at regional and global scales.

54 ENVIRONMENTAL SCIENCES↗

Conceptual Model of Climate Change Impacts at LANL

Goal 9 of the LANL FY15 Site Sustainability Plan (LANL 2014a) addresses Climate Change Adaptation. As part of Goal 9, the plan reviews many of the individual programs the Laboratory has initiated over the past 20 years to address climate change impacts to LANL (e.g. Wildland Fire Management Plan, Forest Management Plan, etc.). However, at that time, LANL did not yet have a comprehensive approach to climate change adaptation. To fill this gap, the FY15 Work Plan for the LANL Long Term Strategy for Environmental Stewardship and Sustainability (LANL 2015) included a goal of (1) establishing a comprehensive conceptual model of climate change impacts at LANL and (2) establishing specific climate change indices to measure climate change and impacts at Los Alamos. Establishing a conceptual model of climate change impacts will demonstrate that the Laboratory is addressing climate change impacts in a comprehensive manner. This paper fulfills the requirement of goal 1. The establishment of specific indices of climate change at Los Alamos (goal 2), will improve our ability to determine climate change vulnerabilities and assess risk. Future work will include prioritizing risks, evaluating options/technologies/costs, and where appropriate, taking actions. To develop a comprehensive conceptual model of climate change impacts, we selected the framework provided in the National Oceanic and Atmospheric Administration (NOAA) Climate Resilience Toolkit (http://toolkit.climate.gov/).

54 ENVIRONMENTAL SCIENCES↗