Implications of climate model selection for projections of decision-relevant metrics: A case study of chill hours in California
Not Available
SEARCH · Engineering Papers
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.
Not Available
Understanding terrestrial ecosystems and their response to anthropogenic climate change requires quantification of land-atmosphere carbon exchange. However, top-down and bottom-up estimates of large-scale land-atmosphere fluxes, including the northern extratropical growing season net flux (GSNF), show significant discrepancies. We developed a data-driven metric for the GSNF using atmospheric carbon dioxide concentration observations collected during the High-Performance Instrumented Airborne Platform for Environmental Research Pole-to-Pole Observations and Atmospheric Tomography Mission flight campaigns. This aircraft-derived metric is bias-corrected using three independent atmospheric inversion systems. We estimate the northern extratropical GSNF to be 5.7 ± 0.3 Pg C and use it to evaluate net biosphere productivity from the Coupled Model Intercomparison Project phase 5 and 6 (CMIP5 and CMIP6) models. While the model-to-model spread in the GSNF has decreased in the CMIP6 models relative to that of the CMIP5 models, there is still disagreement on the magnitude and timing of seasonal carbon uptake with most models underestimating the GSNF and overestimating the length of the growing season relative to the observations. We also use an emergent constraint approach to estimate annual northern extratropical gross primary productivity to be 56 ± 17 Pg C, heterotrophic respiration to be 25 ± 13 Pg C, and net primary productivity to be 28 ± 12 Pg C. The flux inferred from these aircraft observations provides an additional constraint on large-scale gross fluxes in prognostic Earth system models that may ultimately improve our ability to accurately predict carbon-climate feedbacks.
Evaluation methodologies in visualization have mostly focused on how well the tools and techniques cater to the analytical needs of the user. While this is important in determining the effectiveness of the tools and advancing the state-of-the-art in visualization research, a key area that has mostly been overlooked is how well established visualization theories and principles are instantiated in practice. This is especially relevant when domain experts, and not visualization researchers, design visualizations for analysis of their data or for broader dissemination of scientific knowledge. There is very little research on exploring the synergistic capabilities of cross-domain collaboration between domain experts and visualization researchers. To fill this gap, in this paper we describe the results of an exploratory study of climate data visualizations conducted in tight collaboration with a pool of climate scientists. The study analyzes a large set of static climate data visualizations for identifying their shortcomings in terms of visualization design. The outcome of the study is a classification scheme that categorizes the design problems in the form of a descriptive taxonomy. The taxonomy is a first attempt for systematically categorizing the types, causes, and consequences of design problems in visualizations created by domain experts. We demonstrate the use of the taxonomy for a number of purposes, such as, improving the existing climate data visualizations, reflecting on the impact of the problems for enabling domain experts in designing better visualizations, and also learning about the gaps and opportunities for future visualization research. We demonstrate the applicability of our taxonomy through a number of examples and discuss the lessons learnt and implications of our findings.
To assess deep convective parameterizations in a variety of GCMs and examine the fast-time-scale convective transition, a set of statistics characterizing the pickup of precipitation as a function of column water vapor (CWV), PDFs and joint PDFs of CWV and precipitation, and the dependence of the moisture–precipitation relation on tropospheric temperature is evaluated using the hourly output of two versions of the GFDL Atmospheric Model, version 4 (AM4), NCAR CAM5 and superparameterized CAM (SPCAM). The 6-hourly output from the MJO Task Force (MJOTF)/GEWEX Atmospheric System Study (GASS) project is also analyzed. Contrasting statistics produced from individual models that primarily differ in representations of moist convection suggest that convective transition statistics can substantially distinguish differences in convective representation and its interaction with the large-scale flow, while models that differ only in spatial–temporal resolution, microphysics, or ocean–atmosphere coupling result in similar statistics. Most of the models simulate some version of the observed sharp increase in precipitation as CWV exceeds a critical value, as well as that convective onset occurs at higher CWV but at lower column RH as temperature increases. While some models quantitatively capture these observed features and associated probability distributions, considerable intermodel spread and departures from observations in various aspects of the precipitation–CWV relationship are noted. For instance, in many of the models, the transition from the low-CWV, nonprecipitating regime to the moist regime for CWV around and above critical is less abrupt than in observations. Additionally, some models overproduce drizzle at low CWV, and some require CWV higher than observed for strong precipitation. For many of the models, it is particularly challenging to simulate the probability distributions of CWV at high temperature.
Herein the relationship between African dust and the climatology of tropical cyclones (TCs) in the North Atlantic is explored using the Community Atmosphere Model at a global horizontal resolution of 28 km. A simulation in which the aerosol model is modified to significantly reduce the amount of airborne dust is compared to a standard simulation. The simulation with reduced dust increases TC frequency globally, with the largest increase occurring in the North Atlantic. The increase in TC activity in the North Atlantic is consistent with an environment that is more conducive for the genesis and intensification of storms. TCs are more frequent (27%) and on average significantly longer lived (13%) in the low dust configuration but only slightly stronger (3%). This results in a 57% increase in accumulated cyclone energy per hurricane season on average. This work has implications for projections of future climate and resulting changes in TC activity.
Not Available
We introduce the Clouds Above the United States and Errors at the Surface (CAUSES) project with its aim of better understanding the physical processes leading to warm screen-temperature biases over the American Midwest in many numerical models. In this first of four companion papers, 11 different models, from 9 institutes, perform a series of 5-day hindcasts, each initialised from reanalyses. After describing the common experimental protocol and detailing each model configuration, a gridded temperature data set is derived from observations and used to show that all the models have a warm bias over parts of the Midwest. Additionally, a strong diurnal cycle in the screen-temperature bias is found in most models. In some models the bias is largest around midday, while in others it is largest during the night. At the Department of Energy Atmospheric Radiation Measurement Southern Great Plains (SGP) site, the model biases are shown to extend several kilometers into the atmosphere. Finally, to provide context for the companion papers, in which observations from the SGP site are used to evaluate the different processes contributing to errors there, it is shown that there are numerous locations across the Midwest where the diurnal cycle of the error is highly correlated with the diurnal cycle of the error at SGP. This suggests that conclusions drawn from detailed evaluation of models using instruments located at SGP will be representative of errors that are prevalent over a larger spatial scale.
Abstract Atmospheric carbon dioxide (CO 2 ) concentrations have increased as a direct result of human activity and are at their highest level over the last 2 million years, with profound impacts on the Earth system. However, the magnitude and future dynamics of land and ocean carbon sinks are not well understood; therefore, the amount of anthropogenic fossil fuel emissions that remain in the atmosphere (the airborne fraction) is poorly constrained. This work aims to quantify the sources and controls of atmospheric CO 2 , the fate of anthropogenic CO 2 over time, and the likelihood of a trend in the airborne fraction. We use Hector v3.0, a coupled simple climate and carbon cycle model with the novel ability to explicitly track carbon as it flows through the Earth system. We use key model parameters in a Monte Carlo analysis of 15 000 model runs from 1750 to 2300. Results are filtered for physical realism against historical observations and CMIP6 projection data, and we calculate the relative importance of parameters controlling how much anthropogenic carbon ends up in the atmosphere. Modeled airborne fraction was roughly 52%, consistent with observational studies. The overwhelming majority of model runs exhibited a negative trend in the airborne fraction from 1960–2020, implying that current-day land and ocean sinks are proportionally taking up more carbon than the atmosphere. However, the percentage of atmospheric CO 2 derived from anthropogenic origins can be much higher because of Earth system feedbacks. We find it peaks at over 90% between 2010–2050. Moreover, when looking at the destination of anthropogenic fossil fuel emissions, only a quarter ends up in the atmosphere while more than half of emissions are taken up by the land sink on centennial timescales. This study evaluates the likelihood of airborne fraction trends and provides insights into the dynamics of anthropogenic CO 2 in the Earth system.
We report that earth system models are complex and represent a large number of processes, resulting in a persistent spread across climate projections for a given future scenario. Owing to different model performances against observations and the lack of independence among models, there is now evidence that giving equal weight to each available model projection is suboptimal. This Perspective discusses newly developed tools that facilitate a more rapid and comprehensive evaluation of model simulations with observations, process-based emergent constraints that are a promising way to focus evaluation on the observations most relevant to climate projections, and advanced methods for model weighting. Lastly, these approaches are needed to distill the most credible information on regional climate changes, impacts, and risks for stakeholders and policy-makers.
Model projections of permafrost thaw during the next century diverge widely. Here we used ecosys to examine how climate change will affect permafrost thaw in a polygonal tundra at Barrow AK. The model was tested against diurnal and seasonal variation in energy exchange, soil heat flux, soil temperature ( T s ), and active layer depth (ALD) measured during 2014 and 2015, and interannual variation in ALD measured from 1991 to 2015. During RCP 8.5 climate change from 2015 to 2085, increases in T a and precipitation ( P ) to 6.2 °C and 27% above current values, and in atmospheric CO 2 concentrations ( C a ) to 763 μmol mol -1 , altered energy exchange by increasing leaf area index of dominant sedge relative to that of moss. Increased Ca and sedge leaf area index imposed greater stomatal control of transpiration and reduced soil heat fluxes, slowing soil warming, limiting increases in evapotranspiration, and thereby causing gradual soil wetting. Consequently, increases in surface T s and ALD of 2.4–4.7 °C and 21–24 cm above current values were modeled after 70 years. ALD increase was slowed if model boundary conditions were altered to improve landscape drainage. These rates were smaller than those of earlier modeling studies, some of which did not account for changes in vegetation, but are closer to those derived from current studies of warming impacts in the region. Therefore, accounting for climate change effects on vegetation density and composition, and consequent effects on surface energy budgets, will cause slower increases in T s and ALD to be modeled during climate change simulations.
Abstract The present study aims to identify the precipitation bias associated with the interactions among fast physical processes in the Community Atmospheric Model, version 5 (CAM5), during the abrupt onset of the South China Sea (SCS) summer monsoon, a key precursor of the overall East Asia summer monsoon (EASM). The multiyear hindcast approach is utilized to obtain the well-constrained synoptic-scale horizontal circulation each year during the onset period from the years 1998 to 2012. In the pre-onset period, the ocean precipitation over the SCS is insufficiently suppressed in CAM5 hindcasts and thus weaker land–ocean precipitation contrasts. This is associated with the weaker and shallower convection simulated over the surrounding land, producing weaker local circulation within the SCS basin. In the post-onset period, rainfall of the organized convection over the Philippine coastal ocean is underestimated in the hindcasts, with overestimated upper-level heating. These biases are further elaborated as the underrepresentation of the convection diurnal cycle and coastal convection systems, as well as the issue of precipitation sensitivity to environmental moisture during the SCS onset period. The biases identified in hindcasts are consistent with the general bias of the EASM in the climate simulation of CAM5. The current results highlight that the appropriate representation of land–ocean–convection interactions over coastal areas can potentially improve the simulation of seasonal transition over the monsoon regions.
A modal aerosol module (MAM) has been developed for the Community Atmosphere Model version 5 (CAM5), the atmospheric component of the Community Earth System Model version 1 (CESM1). MAM is capable of simulating the aerosol size distribution and both internal and external mixing between aerosol components, treating numerous complicated aerosol processes and aerosol physical, chemical and optical properties in a physically based manner. Two MAM versions were developed: a more complete version with seven-lognormal modes (MAM7), and a three-lognormal mode version (MAM3) for the purpose of long-term (decades to centuries) simulations. Major approximations in MAM3 include assuming immediate mixing of primary organic matter (POM) and black carbon (BC) with other aerosol components, merging of the MAM7 fine dust and fine sea salt modes into the accumulation mode, merging of the MAM7 coarse dust and coarse sea salt modes into the single coarse mode, and neglecting the explicit treatment of ammonia and ammonium cycles. Simulated sulfate and secondary organic aerosol (SOA) mass concentrations are remarkably similar between MAM3 and MAM7 as most ({approx}90%) of these aerosol species are in the accumulation mode. Differences of POM and BC concentrations between MAM3 and MAM7 are also small (mostly within 10%) because of the assumed hygroscopic nature of POM, so that freshly emitted POM and BC are wet-removed before mixing internally with soluble aerosol species. Sensitivity tests with the POM assumed to be hydrophobic and with slower aging process increase the POM and BC concentrations, especially at high latitudes (by several times). The mineral dust global burden differs by 10% and sea salt burden by 30-40% between MAM3 and MAM7 mainly due to the different size ranges for dust and sea salt modes and different standard deviations of log-normal size distribution for sea salt modes between MAM3 and MAM7. The model is able to qualitatively capture the observed geographical and temporal variations of aerosol mass and number concentrations, size distributions, and aerosol optical properties. However, there are noticeable biases, e.g., simulated sulfate and mineral dust concentrations at surface over the oceans are too low. Simulated BC concentrations are significant low in the Arctic. There is a low bias in modeled aerosol optical depth on the global scale, especially in the developing counties. There biases in aerosol simulations clearly indicate the need for improvements of aerosol processes (e.g., emission fluxes of anthropogenic aerosols and precursor gases in developing countries, boundary layer nucleation) and properties (e.g., primary aerosol emission size, POM hygroscopicity). In addition the critical role of cloud properties (e.g., liquid water content, cloud fraction) responsible for the wet scavenging of aerosol is highlighted.
Two types of Arctic mixed-phase clouds observed during the ISDAC and M-PACE field campaigns are simulated using a 3-dimensional cloud-resolving model (CRM) with size-resolved cloud microphysics. The modeled cloud properties agree reasonably well with aircraft measurements and surface-based retrievals. Cloud properties such as the probability density function (PDF) of vertical velocity (w), cloud liquid and ice, the regime of ice growth at the expense of liquid water (i.e., Wegener-Bergeron-Findeisen (WBF) process), and the inherent relationships among cloud properties/processes in the mixed-phase layers are examined to gain insights for improving the representation of the mixed-phase processes in General Circulation Models (GCMs). We find that, the WBF process only occurs in about 50% of the mixed-phase regime with the vast majority occurring in the downdrafts. In updrafts both liquid and ice grow simultaneously. But in GCMs, it is not necessary to treat the WBF process at the subgrid scale. Our CRM results produce a w distribution well represented by a Gaussian normal function, validating, at least for arctic clouds, the subgrid treatment used in GCMs. Our CRM results also support the assumption frequently used in GCMs that mixed phase clouds maintain water vapor very near liquid saturation. A Gamma function with a fixed variance does not accurately represent the subgrid variability of cloud liquid. The PDFs of cloud liquid and cloud ice can be fitted with Gamma functions, and a normal function can be used for total water, but the variance should not be fixed. The relationship between the ice depositional growth rate and cloud ice strongly depends on the capacitance of ice particles. The assumption for the capacitance of ice particles (e.g., 1.0 for spheres) used in GCMs could lead to a large deviation in ice depositional growth. At large sales, the maximum overlap assumption looks appropriate.
In this DOE project the improvements to parameterization of marine primary organic matter (POM) emissions, hygroscopic properties of marine POM, marine isoprene derived secondary organic aerosol (SOA) emissions, surfactant effects, new cloud droplet activation parameterization have been implemented into Community Atmosphere Model (CAM 5.0), with a seven mode aerosol module from the Pacific Northwest National Laboratory (PNNL)'s Modal Aerosol Model (MAM7). The effects of marine aerosols derived from sea spray and ocean emitted biogenic volatile organic compounds (BVOCs) on microphysical properties of clouds were explored by conducting 10 year CAM5.0-MAM7 model simulations at a grid resolution 1.9° by 2.5° with 30 vertical layers. Model-predicted relationship between ocean physical and biological systems and the abundance of CCN in remote marine atmosphere was compared to data from the A-Train satellites (MODIS, CALIPSO, AMSR-E). Model simulations show that on average, primary and secondary organic aerosol emissions from the ocean can yield up to 20% increase in Cloud Condensation Nuclei (CCN) at 0.2% Supersaturation, and up to 5% increases in droplet number concentration of global maritime shallow clouds. Marine organics were treated as internally or externally mixed with sea salt. Changes associated with cloud properties reduced (absolute value) the model-predicted short wave cloud forcing from -1.35 Wm -2 to -0.25 Wm -2 . By using different emission scenarios, and droplet activation parameterizations, this study suggests that addition of marine primary aerosols and biologically generated reactive gases makes an important difference in radiative forcing assessments. All baseline and sensitivity simulations for 2001 and 2050 using global-through-urban WRF/Chem (GU-WRF) were completed. The main objective of these simulations was to evaluate the capability of GU-WRF for an accurate representation of the global atmosphere by exploring the most accurate configuration of physics options in GWRF for global scale modeling in 2001 at a horizontal grid resolution of 1° x 1°. GU-WRF model output was evaluated using observational datasets from a variety of sources including surface based observations (NCDC and BSRN), model reanalysis (NCEP/ NCAR Reanalysis and CMAP), and remotely-sensed data (TRMM) to evaluate the ability of GU-WRF to simulate atmospheric variables at the surface as well as aloft. Explicit treatment of nanoparticles produced from new particle formation in GU-WRF/Chem-MADRID was achieved by expanding particle size sections from 8 to 12 to cover particles with the size range of 1.16 nm to 11.6m. Simulations with two different nucleation parameterizations were conducted for August 2002 over a global domain at a 4º by 5º horizontal resolution. The results are evaluated against field measurement data from the 2002 Aerosol Nucleation and Real Time Characterization Experiment (ANARChE) in Atlanta, Georgia, as well as satellite and reanalysis data. We have also explored the relationship between clean marine aerosol optical properties and ocean surface wind speed using remotely sensed data from the Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) on board the CALIPSO satellite and the Advanced Microwave Scanning Radiometer (AMSR-E) on board the AQUA satellite. Detailed data analyses were carried out over 15 regions selected to be representative of different areas of the global ocean for the time period from June 2006 to April 2011. We show that for very low (less than 4 m s -1 ) and very high (more than 12 m s -1 ) wind speed conditions the mean CALIPSO-derived aerosol optical depth (AOD) has little dependency on the surface wind speed. For an intermediate (between 4 and 12 m s -1 ) marine AOD was linearly correlated with the surface wind speed values, with a slope of 0.0062 s m -1 . Results of our study suggest that considerable improvements to both optical properties of marine aerosols and their production mechanisms can be achieved by discriminating clean marine aerosols (or sea salt particles) from all other types of aerosols present over the ocean.
Models predicting ecosystem carbon dioxide (CO 2 ) exchange under future climate change rely on relatively few real-world tests of their assumptions and outputs. Here, we demonstrate a rapid and cost-effective method to estimate CO 2 exchange from intact vegetation patches under varying atmospheric CO 2 concentrations.We find that net ecosystem CO 2 uptake (NEE) in a boreal forest rose linearly by 4.7 0.2% of the current ambient rate for every 10 ppm CO 2 increase, with no detectable influence of foliar biomass, season, or nitrogen (N) fertilization. The lack of any clear short-term NEE response to fertilization in such an N-limited system is inconsistent with the instantaneous downregulation of photosynthesis formalized in many global models. Incorporating an alternative mechanism with consider-able empirical support diversion of excess carbon to storage compounds into an existing earth system model brings the model output into closer agreement with our field measurements. A global simulation incorporating this modified model reduces a long-standing mismatch between the modeled and observed seasonal amplitude of atmospheric CO 2 . Wider application of this chamber approach would provide critical data needed to further improve modeled projections of biosphere atmosphere CO 2 exchange in a changing climate.
Abstract. The Scenario Model Intercomparison Project (ScenarioMIP) defines and coordinates the main set of future climate projections, based on concentration-driven simulations, within the Coupled Model Intercomparison Project phase 6 (CMIP6). This paper presents a range of its outcomes by synthesizing results from the participating global coupled Earth system models. We limit our scope to the analysis of strictly geophysical outcomes: mainly global averages and spatial patterns of change for surface air temperature and precipitation. We also compare CMIP6 projections to CMIP5 results, especially for those scenarios that were designed to provide continuity across the CMIP phases, at the same time highlighting important differences in forcing composition, as well as in results. The range of future temperature and precipitation changes by the end of the century (2081–2100) encompassing the Tier 1 experiments based on the Shared Socioeconomic Pathway (SSP) scenarios (SSP1-2.6, SSP2-4.5, SSP3-7.0 and SSP5-8.5) and SSP1-1.9 spans a larger range of outcomes compared to CMIP5, due to higher warming (by close to 1.5 ∘C) reached at the upper end of the 5 %–95 % envelope of the highest scenario (SSP5-8.5). This is due to both the wider range of radiative forcing that the new scenarios cover and the higher climate sensitivities in some of the new models compared to their CMIP5 predecessors. Spatial patterns of change for temperature and precipitation averaged over models and scenarios have familiar features, and an analysis of their variations confirms model structural differences to be the dominant source of uncertainty. Models also differ with respect to the size and evolution of internal variability as measured by individual models' initial condition ensemble spreads, according to a set of initial condition ensemble simulations available under SSP3-7.0. These experiments suggest a tendency for internal variability to decrease along the course of the century in this scenario, a result that will benefit from further analysis over a larger set of models. Benefits of mitigation, all else being equal in terms of societal drivers, appear clearly when comparing scenarios developed under the same SSP but to which different degrees of mitigation have been applied. It is also found that a mild overshoot in temperature of a few decades around mid-century, as represented in SSP5-3.4OS, does not affect the end outcome of temperature and precipitation changes by 2100, which return to the same levels as those reached by the gradually increasing SSP4-3.4 (not erasing the possibility, however, that other aspects of the system may not be as easily reversible). Central estimates of the time at which the ensemble means of the different scenarios reach a given warming level might be biased by the inclusion of models that have shown faster warming in the historical period than the observed. Those estimates show all scenarios reaching 1.5 ∘C of warming compared to the 1850–1900 baseline in the second half of the current decade, with the time span between slow and fast warming covering between 20 and 27 years from present. The warming level of 2 ∘C of warming is reached as early as 2039 by the ensemble mean under SSP5-8.5 but as late as the mid-2060s under SSP1-2.6. The highest warming level considered (5 ∘C) is reached by the ensemble mean only under SSP5-8.5 and not until the mid-2090s.
Several independent measurements of warm-season soil moisture and surface atmospheric variables recorded at the ARM Southern Great Plains (SGP) research facility are used to estimate the terrestrial component of land-atmosphere coupling (LAC) strength and its regional uncertainty. The observations reveal substantial variation in coupling strength, as estimated from three soil moisture measurements at a single site, as well as across six other sites having varied soil and land cover types. The observational estimates then serve as references for evaluating SGP terrestrial coupling strength in the Community Atmospheric Model coupled to the Community Land Model. These coupled model components are operated in both a free-running mode and in a controlled configuration, where the atmospheric and land states are reinitialized daily, so that they do not drift very far from observations. Although the controlled simulation deviates less from the observed surface climate than its free-running counterpart, the terrestrial LAC in both configurations is much stronger and displays less spatial variability than the SGP observational estimates. Preliminary investigation of vegetation leaf area index (LAI) substituted for soil moisture suggests that the overly strong coupling between model soil moisture and surface atmospheric variables is associated with too much evaporation from bare ground and too little from the vegetation cover. Lastly, these results imply that model surface characteristics such as LAI, as well as the physical parameterizations involved in the coupling of the land and atmospheric components, are likely to be important sources of the problematical LAC behaviors.
Abstract Understanding internal variability of the climate system is critical when isolating internal and anthropogenically forced signals. Here, we investigate the modes of Atlantic Meridional Overturning Circulation ( AMOC ) variability using perturbation experiments with the Institut Pierre‐Simon Laplace's (IPSL) coupled model and compare them to Coupled Model Intercomparison Project Phase 6 (CMIP6) pre‐industrial control simulations. We identify two characteristic modes of variability—decadal‐to‐multidecadal ( DMD var ) and centennial ( CEN var ). The former is driven largely by temperature anomalies in the subpolar North Atlantic, while the latter is driven by salinity in the western subpolar North Atlantic. The amplitude of each mode scales linearly with the mean AMOC strength in the IPSL experiments. The DMD var amplitude correlates well with the AMOC mean strength across CMIP6 models, while the CEN var mode does not. These findings suggest that the strength of DMD var depends robustly on the North Atlantic mean state, while the CEN var mode may be model‐dependent.