Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 1,243 records · Page 69
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Effective radiative forcing and adjustments in CMIP6 models

Abstract. The effective radiative forcing, which includes the instantaneous forcing plus adjustments from the atmosphere and surface, has emerged as the key metric of evaluating human and natural influence on the climate. We evaluate effective radiative forcing and adjustments in 17 contemporary climate models that are participating in the Coupled Model Intercomparison Project (CMIP6) and have contributed to the Radiative Forcing Model Intercomparison Project (RFMIP). Present-day (2014) global-mean anthropogenic forcing relative to pre-industrial (1850) levels from climate models stands at 2.00 (±0.23) W m−2, comprised of 1.81 (±0.09) W m−2 from CO2, 1.08 (± 0.21) W m−2 from other well-mixed greenhouse gases, −1.01 (± 0.23) W m−2 from aerosols and −0.09 (±0.13) W m−2 from land use change. Quoted uncertainties are 1 standard deviation across model best estimates, and 90 % confidence in the reported forcings, due to internal variability, is typically within 0.1 W m−2. The majority of the remaining 0.21 W m−2 is likely to be from ozone. In most cases, the largest contributors to the spread in effective radiative forcing (ERF) is from the instantaneous radiative forcing (IRF) and from cloud responses, particularly aerosol–cloud interactions to aerosol forcing. As determined in previous studies, cancellation of tropospheric and surface adjustments means that the stratospherically adjusted radiative forcing is approximately equal to ERF for greenhouse gas forcing but not for aerosols, and consequentially, not for the anthropogenic total. The spread of aerosol forcing ranges from −0.63 to −1.37 W m−2, exhibiting a less negative mean and narrower range compared to 10 CMIP5 models. The spread in 4×CO2 forcing has also narrowed in CMIP6 compared to 13 CMIP5 models. Aerosol forcing is uncorrelated with climate sensitivity. Therefore, there is no evidence to suggest that the increasing spread in climate sensitivity in CMIP6 models, particularly related to high-sensitivity models, is a consequence of a stronger negative present-day aerosol forcing and little evidence that modelling groups are systematically tuning climate sensitivity or aerosol forcing to recreate observed historical warming.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Continental United States climate projections based on thermodynamic modification of historical weather

Abstract Regional climate models can be used to examine how past weather events might unfold under different climate conditions by simulating analogue versions of those events with modified thermodynamic conditions (i.e., warming signals). Here, we apply this approach by dynamically downscaling a 40-year sequence of past weather from 1980–2019 driven by atmospheric re-analysis, and then repeating this 40-year sequence a total of 8 times using a range of time-evolving thermodynamic warming signals that follow 4 80-year future warming trajectories from 2020–2099. Warming signals follow two emission scenarios (SSP585 and SSP245) and are derived from two groups of global climate models based on whether they exhibit relatively high or low climate sensitivity. The resulting dataset, which contains 25 hourly and over 200 3-hourly variables at 12 km spatial resolution, can be used to examine a plausible range of future climate conditions in direct reference to previously observed weather and enables a systematic exploration of the ways in which thermodynamic change influences the characteristics of historical extreme events.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

Are atmospheric updrafts a key to unlocking climate forcing and sensitivity?

Both climate forcing and climate sensitivity persist as stubborn uncertainties limiting the extent to which climate models can provide actionable scientific scenarios for climate change. A key, explicit control on cloud–aerosol interactions, the largest uncertainty in climate forcing, is the vertical velocity of cloud-scale updrafts. Model-based studies of climate sensitivity indicate that convective entrainment, which is closely related to updraft speeds, is an important control on climate sensitivity. Updraft vertical velocities also drive many physical processes essential to numerical weather prediction. Vertical velocities and their role in atmospheric physical processes have been given very limited attention in models for climate and numerical weather prediction. The relevant physical scales range down to tens of meters and are thus frequently sub-grid and require parameterization. Many state-of-science convection parameterizations provide mass fluxes without specifying vertical velocities, and parameterizations that do provide vertical velocities have been subject to limited evaluation against what have until recently been scant observations. Atmospheric observations imply that the distribution of vertical velocities depends on the areas over which the vertical velocities are averaged. Distributions of vertical velocities in climate models may capture this behavior, but it has not been accounted for when parameterizing cloud and precipitation processes in current models. New observations of convective vertical velocities offer a potentially promising path toward developing process-level cloud models and parameterizations for climate and numerical weather prediction. Taking account of the scale dependence of resolved vertical velocities offers a path to matching cloud-scale physical processes and their driving dynamics more realistically, with a prospect of reduced uncertainty in both climate forcing and sensitivity.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

Are Atmospheric Updrafts a Key to Unlocking Climate Forcing and Sensitivity?

Both climate forcing and climate sensitivity persist as stubborn uncertainties limiting the extent to which climate models can provide actionable scientific scenarios for climate change. A key, explicit control on cloud-aerosol interactions, the largest uncertainty in climate forcing, is the vertical velocity of cloud-scale updrafts. Model-based studies of climate sensitivity indicate that convective entrainment, which is closely related to updraft speeds, is an important control on climate sensitivity. Updraft vertical velocities also drive many physical processes essential to numerical weather prediction. Vertical velocities and their role in atmospheric physical processes have been given very limited attention in models for climate and numerical weather prediction. The relevant physical scales range down to tens of meters and are thus frequently sub-grid and require parameterization. Many state-of-science convection parameterizations provide mass fluxes without specifying vertical velocities, and parameterizations which do provide vertical velocities have been subject to limited evaluation against what have until recently been scant observations. Atmospheric observations imply that the distribution of vertical velocities depends on the areas over which the vertical velocities are averaged. Distributions of vertical velocities in climate models may capture this behavior, but it has not been accounted for when parameterizing cloud and precipitation processes in current models. New observations of convective vertical velocities offer a potentially promising path toward developing process-level cloud models and parameterizations for climate and numerical weather prediction. Taking account of scale-dependence of resolved vertical velocities offers a path to matching cloud-scale physical processes and their driving dynamics more realistically, with a prospect of reduced uncertainty in both climate forcing and sensitivity.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2003

Observed and Modeled Trends in Southern Ocean Sea Ice

Conceptual models and global climate model (GCM) simulations have both indicated the likelihood of an enhanced sensitivity to climate change in the polar regions, derived from the positive feedbacks brought about by the polar abundance of snow and ice surfaces. Some models further indicate that the changes in the polar regions can have a significant impact globally. For instance, 37% of the temperature sensitivity to a doubling of atmospheric CO2 in simulations with the GCM of the Goddard Institute for Space Studies (GISS) is attributable exclusively to inclusion of sea ice variations in the model calculations. Both sea ice thickness and sea ice extent decrease markedly in the doubled CO, case, thereby allowing the ice feedbacks to occur. Stand-alone sea ice models have shown Southern Ocean hemispherically averaged winter ice-edge retreats of 1.4 deg latitude for each 1 K increase in atmospheric temperatures. Observations, however, show a much more varied Southern Ocean ice cover, both spatially and temporally, than many of the modeled expectations. In fact, the satellite passive-microwave record of Southern Ocean sea ice since late 1978 has revealed overall increases rather than decreases in ice extents, with ice extent trends on the order of 11,000 sq km/year. When broken down spatially, the positive trends are strongest in the Ross Sea, while the trends are negative in the Bellingshausen/Amundsen Seas. Greater spatial detail can be obtained by examining trends in the length of the sea ice season, and those trends show a coherent picture of shortening sea ice seasons throughout almost the entire Bellingshausen and Amundsen Seas to the west of the Antarctic Peninsula and in the far western Weddell Sea immediately to the east of the Peninsula, with lengthening sea ice seasons around much of the rest of the continent. This pattern corresponds well with the spatial pattern of temperature trends, as the Peninsula region is the one region in the Antarctic with a strong record of temperature increases. Still, although the patterns of the temperature and ice changes match fairly well, there is a substantial ways to go before these patterns are understood (and can be modeled) in the full context of global change.

Parkinson, Claire L.↗
NASA NTRS2018

The Present and Future of Secondary Organic Aerosol Direct Forcing on Climate

Secondary organic aerosols (SOA), a subset of organic aerosols that are chemically produced in the atmosphere, are included in climate modeling calculations using very simple parameterizations. Estimates on their shortwave forcing on climate span almost two orders of magnitude, being potentially comparable to sulfate direct forcing. In the longwave, a neglected part of the spectrum when it comes to SOA, the direct SOA forcing could exceed that of sulfate and black carbon, although in absolute values it is much weaker than the shortwave forcing. Critical for these estimates is the vertical distribution of the climate active agents, pointing to SOA temperature-dependent volatility. Over the last few years, research also revealed the highly oxidized character of organic aerosol and its chemical aging in the atmosphere that partially leads to the formation of brown carbon, an absorbing form of organic aerosol. This review summarizes critical advances in the understanding of SOA behavior and properties relevant to direct climate forcing and puts them in perspective with regard to primary organic aerosol and brown carbon. These findings also demonstrate an emerging dynamic picture of organic aerosol that has not yet been integrated in climate modeling. The challenges for the coming years in order to reduce uncertainties in the direct organic aerosol climate impact are discussed. High priority for future model development should be given to the dynamic link between “white” and “brown” organic aerosol and between primary and secondary organic aerosol. The SOA temperature-dependent volatility parameterizations and wavelength-dependent refractive index should be also included.

Kostas Tsigaridis↗
DOE OSTI2014

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↗
DOE OSTI2025

Using multiple high-resolution datasets to benchmark the energy exascale earth system model (E3SM) for renewable resource assessment

The United States is accelerating its shift toward a renewable energy system. However, renewable resources, which harness energy from the Earth system, are susceptible to both present-day climate variability and future climate change. For example, variations in regional climate can alter renewable energy production patterns and site viability. The use of high-resolution climate model projections can therefore facilitate and may be critical to long-term planning of renewable energy investments. However, climate models must first be validated for renewable resource assessment. This research employs multiple high-spatiotemporal-resolution datasets to assess the capability of the Department of Energy’s (DOE) Energy Exascale Earth System Model version 2 North American Regionally Refined Model (E3SMv2-NARRM) for predicting multi-year climatological values of solar and wind energy capacity factors in the continental U.S., with a focus on regional and seasonal variability. Present-day E3SMv2-NARRM simulations are compared with reported utility-scale production data obtained from the Energy Information Administration (EIA). In addition, E3SMv2-NARRM data are evaluated against non-climate benchmark models from the National Renewable Energy Laboratory, including the Wind Integration National Dataset Toolkit and the National Solar Radiation Database (NSRDB), as well as three wind energy datasets from PLUSWIND. Our analysis indicates that solar capacity factors from E3SM closely match those from the NSRDB dataset. However, both datasets tend to overestimate values by 10% in comparison to EIA data. Furthermore, biases in wind capacity factors within E3SM are notably pronounced in the West Coast regions, where the seasonal cycle diverges from EIA data.

Energy forecasting, Capacity factor, Renewable ene↗
DOE OSTI2013

Simulating the temperature and precipitation signal in an Alpine ice core

Accumulation and δ 18 O data from Alpine ice cores provide information on past temperature and precipitation. However, their correlation with seasonal or annual mean temperature and precipitation at nearby sites is often low. This is partly due to the irregular sampling of the atmosphere by the ice core (i.e. ice cores almost only record precipitation events and not dry periods) and the possible incongruity between annual layers and calendar years. Using daily meteorological data from a nearby station and reanalyses, we replicate the ice core from the Grenzgletscher (Switzerland, 4200 m a.s.l.) on a sample-by-sample basis by calculating precipitation-weighted temperature (PWT) over short intervals. Over the last 15 yr of the ice core record, accumulation and δ 18 O variations can be well reproduced on a sub-seasonal scale. This allows a wiggle-matching approach for defining quasi-annual layers, resulting in high correlations between measured quasi-annual δ 18 O and PWT. Further back in time, the agreement deteriorates. Nevertheless, we find significant correlations over the entire length of the record (1938–1993) of ice core δ 18 O with PWT, but not with annual mean temperature. This is due to the low correlations between PWT and annual mean temperature, a characteristic which in ERA-Interim reanalysis is also found for many other continental mid-to-high-latitude regions. The fact that meteorologically very different years can lead to similar combinations of PWT and accumulation poses limitations to the use of δ 18 O from Alpine ice cores for temperature reconstructions. Rather than for reconstructing annual mean temperature, δ 18 O from Alpine ice cores should be used to reconstruct PWT over quasi-annual periods. This variable is reproducible in reanalysis or climate model data and could thus be assimilated into conventional climate models.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2010

Quantifying Uncertainty in Projections of Stratospheric Ozone Over the 21st Century

Future stratospheric ozone concentrations will be determined both by changes in the concentration of ozone depleting substances (ODSs) and by changes in stratospheric and tropospheric climate, including those caused by changes in anthropogenic greenhouse gases (GHGs). Since future economic development pathways and resultant emissions of GHGs are uncertain, anthropogenic climate change could be a significant source of uncertainty for future projections of stratospheric ozone. In this pilot study, using an ensemble of opportunity of chemistry-climate model (CCM) simulations, the contribution of scenario uncertainty from different plausible emissions pathways for 10 ODSs and GHGs to future ozone projections is quantified relative to the contribution from model uncertainty and internal variability of the chemistry-climate system. For both the global, annual mean ozone concentration and for ozone in specific geographical regions, differences between CCMs are the dominant source of uncertainty for the first two-thirds of the 21 st century, up-to and after the time when ozone concentrations 15 return to 1980 values. In the last third of the 21st century, dependent upon the set of greenhouse gas scenarios used, scenario uncertainty can be the dominant contributor. This result suggests that investment in chemistry-climate modelling is likely to continue to refine projections of stratospheric ozone and estimates of the return of stratospheric ozone concentrations to pre-1980 levels.

Charlton-Perez, A. J.↗
NASA NTRS1990

Model test of CCN-cloud albedo climate forcing

Cloud condensation nuclei (CCN) influence cloud albedo through their effect on the cloud droplet size distribution. A number of studies have evaluated the climatic impact of the CCN-cloud albedo feedback, but all have assumed that cloud distributions, cloud thicknesses, and cloud liquid water contents would remain constant as the climate adjusted. This assumption has been tested using the Livermore version of the National Center for Atmospheric Research Community Climate Model. The results indicate that there are no significant compensating changes in cloud properties that would counteract the 1.7 percent global albedo increase resulting from a fourfold increase in marine CCN concentration. Furthermore, when ocean surface temperatures are decreased 4 C in a manner broadly consistent with the enhanced cloud albedos, an increase in cloud fraction of 3.5 percent and a reduction in cloud altitude are predicted, leading to a positive feedback from clouds that would imply a climate impact roughly double that calculated from cloud droplet size distribution change alone.

Ghan, S. J.↗
NASA NTRS2000

Model Study of the Sources of Upper Tropospheric Ozone Over The Atlantic, with Emphasis on Interhemispheric Asymmetries

Measurements of ozone and aerosols from ships transecting the Atlantic Ocean reveal considerable interhemispheric asymmetries in the vertical structure of the ozone distribution. There is much more upper tropospheric ozone to the South of the Equator than in the northern hemisphere. These profiles suggest that the air in the southern hemispheric tropics and subtropics has a source in the stratosphere, while that in the northern hemisphere originates in the boundary layer. A model experiment designed to investigate this problem will be discussed. They are based on a new general circulation model, based on a Finite-Volume dynamical core and the physical parametrizations from the Community Climate Model at the National Center for Atmospheric Research (NCAR). The model simulates a climatological state in reasonable agreement with observations, comparable to other contemporary climate models, but has a much more accurate transport capability than most other present-day models. The experiment involves using a suite of tracers, designed to highlight air from different source regions. Some of these (water vapor and ozone) represent real atmospheric quantities, while others are markers which indicate whether the air originates in the stratosphere, the boundary layer, or has traversed regions of convection (for example). The experiments were run for periods of up to one month, initialized from snapshots from a long climate integration of the model. A sample of results will be presented, chosen to illustrate the different sources of air in the upper troposphere and the physical mechanisms by which these air parcels arrived there.

Pawson, Steven↗
DOE OSTI2018

Computational Benefit of GPU Optimization for the Atmospheric Chemistry Modeling

Abstract Global chemistry‐climate models are computationally burdened as the chemical mechanisms become more complex and realistic. Optimization for graphics processing units (GPU) may make longer global simulation with regional detail possible, but limited study has been done to explore the potential benefit for the atmospheric chemistry modeling. Hence, in this study, the second‐order Rosenbrock solver of the chemistry module of CAM4‐Chem is ported to the GPU to gauge potential speed‐up. We find that on the CPU, the fastest performance is achieved using the Intel compiler with a block interleaved memory layout. Different combinations of compiler and memory layout lead to ~11.02× difference in the computational time. In contrast, the GPU version performs the best when using a combination of fully interleaved memory layout with block size equal to the warp size, CUDA streams for independent kernels, and constant memory. Moreover, the most efficient data transfer between CPU and GPU is gained by allocating the memory contiguously during the data initialization on the GPU. Compared to one CPU core, the speed‐up of using one GPU alone reaches a factor of ~11.7× for the computation alone and ~3.82× when the data transfer between CPU and GPU is considered. Using one GPU alone is also generally faster than the multithreaded implementation for 16 CPU cores in a compute node and the single‐source solution (OpenACC). The best performance is achieved by the implementation of the hybrid CPU/GPU version, but rescheduling the workload among the CPU cores is required before the practical CAM4‐Chem simulation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

CAUSES: On the Role of Surface Energy Budget Errors to the Warm Surface Air Temperature Error Over the Central United States

Many weather forecasting and climate models simulate a warm surface air temperature (T2m) bias over mid-latitude continents during the summertime, especially over the Great Plains. We present here one of a series of papers from a multi-model intercomparison project (CAUSES: Cloud Above the United States and Errors at the Surface), which aims to evaluate the role of cloud, radiation, and precipitation biases in contributing to T2m bias using a short-term hindcast approach with observations mainly from the Atmospheric Radiation Measurement (ARM) Southern Great Plains (SGP) site during the period of April to August 2011. The present study examines the contribution of surface energy budget errors to the bias. All participating models simulate higher net shortwave and longwave radiative fluxes at the surface but there is no consistency on signs of biases in latent and sensible heat fluxes over the Central U.S. and ARM SGP. Nevertheless, biases in net shortwave and downward longwave fluxes, as well as surface evaporative fraction (EF) are the main contributors to T2m bias. Radiation biases are largely affected by cloud simulations, while EF is affected by soil moisture modulated by seasonal accumulated precipitation and evaporation. An approximate equation is derived to further quantify the magnitudes of radiation and EF contributions to T2m bias. Our analysis suggests that radiation errors are always an important source of T2m error for long-term climate runs with EF errors either of equal or lesser importance. However, for the short-term hindcasts, EF errors are more important provided a model has a substantial EF bias.

Ma, H. -Y.↗
NASA NTRS2013

Radiative Forcing in the ACCMIP Historical and Future Climate Simulations

A primary goal of the Atmospheric Chemistry and Climate Model IntercomparisonProject (ACCMIP) was to characterize the short-lived drivers of preindustrial to 2100climate change in the current generation of climate models. Here we evaluate historicaland 5 future radiative forcing in the 10 ACCMIP models that included aerosols, 8 of whichalso participated in the Coupled Model Intercomparison Project phase 5 (CMIP5).The models generally reproduce present-day climatological total aerosol opticaldepth (AOD) relatively well. components to this total, however, and most appear to underestimate AOD over East10 Asia. The models generally capture 1980-2000 AOD trends fairly well, though theyunderpredict AOD increases over the YellowEastern Sea. They appear to strongly underestimate absorbing AOD, especially in East Asia, South and Southeast Asia, SouthAmerica and Southern Hemisphere Africa.We examined both the conventional direct radiative forcing at the tropopause (RF) and the forcing including rapid adjustments (adjusted forcing AF, including direct andindirect effects). The models calculated all aerosol all-sky 1850 to 2000 global meanannual average RF ranges from 0.06 to 0.49 W m(sup -2), with a mean of 0.26 W m(sup -2) and a median of 0.27 W m(sup -2. Adjusting for missing aerosol components in some modelsbrings the range to 0.12 to 0.62W m(sup -2), with a mean of 0.39W m(sup -2). Screen20ing the models based on their ability to capture spatial patterns and magnitudes ofAOD and AOD trends yields a quality-controlled mean of 0.42W m(sup -2) and range of0.33 to 0.50 W m(sup -2) (accounting for missing components). The CMIP5 subset of ACCMIPmodels spans 0.06 to 0.49W m(sup -2), suggesting some CMIP5 simulations likelyhave too little aerosol RF. A substantial, but not well quantified, contribution to histori25cal aerosol RF may come from climate feedbacks (35 to 58). The mean aerosol AF during this period is 1.12W m(sup -2) (median value 1.16W m(sup -2), range 0.72 to1.44W m(sup -2), indicating that adjustments to aerosols, which include cloud, water vaporand temperature, lead to stronger forcing than the aerosol direct RF.

Climate Models↗
DOE OSTI2024

Quantifying and Modeling the Impact of Phase State on the Ice Nucleation Abilities of 2-Methyltetrols as a Key Component of Secondary Organic Aerosol Derived from Isoprene Epoxydiols

Organic aerosols (OAs) may serve as ice-nucleating particles (INPs), impacting the formation and properties of cirrus clouds when their phase state and viscosity are in the semisolid to glassy range. However, there is a lack of direct parameterization between aerosol viscosity and their ice nucleation capabilities. In this study, we experimentally measured the ice nucleation rate of 2-methyltetrols (2-MT) aerosols, a key component of isoprene-epoxydiol-derived secondary organic aerosols (IEPOX-SOA), at different viscosities. These results demonstrate that the phase state has a significant impact on the ice nucleation abilities of OA under typical cirrus cloud conditions, with the ice nucleation rate increasing by 2 to 3 orders of magnitude when the phase state changes from liquid to semisolid. An innovative parametric model based on classical nucleation theory was developed to directly quantify the impact of viscosity on the heterogeneous nucleation rate. This model accurately represents our laboratory measurement and can be implemented into climate models due to its simple, equation-based form. Based on data collected from the ACRIDICON-CHUVA field campaign, our model predicts that the INP concentration from IEPOX-SOA can reach the magnitude of 1 to tens per liter in the cirrus cloud region impacted by the Amazon rainforest, consistent with recent field observations and estimations. This novel parameterization framework can also be applied in regional and global climate models to further improve representations of cirrus cloud formation and associated climate impacts.

2-methyltetrol↗
NASA NTRS2022

Automatic variable selection in ecological niche modeling: A case study using Cassin’s Sparrow (Peucaea cassinii)

MERRA/Max provides a feature selection approach to dimensionality reduction that enables direct use of global climate model outputs in ecological niche modeling. The system accomplishes this reduction through a Monte Carlo optimization in which many independent MaxEnt runs, operating on a species occurrence file and a small set of randomly selected variables in a large collection of variables, converge on an estimate of the top contributing predictors in the larger collection. These top predictors can be viewed as potential candidates in the variable selection step of the ecological niche modeling process. MERRA/Max’s Monte Carlo algorithm operates on files stored in the underlying filesystem, making it scalable to large data sets. Its software components can run as parallel processes in a high-performance cloud computing environment to yield near real-time performance. In tests using Cassin’s Sparrow (Peucaea cassinii) as the target species, MERRA/Max selected a set of predictors from Worldclim’s Bioclim collection of 19 environmental variables that have been shown to be important determinants of the species’ bioclimatic niche. It also selected biologically and ecologically plausible predictors from a more diverse set of 86 environmental variables derived from NASA’s Modern-Era Retrospective Analysis for Research and Applications Version 2 (MERRA-2) reanalysis, an output product of the Goddard Earth Observing System Version 5 (GEOS-5) modeling system. We believe these results point to a technological approach that could expand the use global climate model outputs in ecological niche modeling, foster exploratory experimentation with otherwise difficult-to-use climate data sets, streamline the modeling process, and, eventually, enable automated bioclimatic modeling as a practical, readily accessible, low-cost, commercial cloud service.

John L. Schnase↗
DOE OSTI2018

Spatially Dependent Multiple Testing Under Model Misspecification, With Application to Detection of Anthropogenic Influence on Extreme Climate Events

The Weather Risk Attribution Forecast (WRAF) is a forecasting tool that uses output from global climate models to make simultaneous attribution statements about whether and how greenhouse gas emissions have contributed to extreme weather across the globe. However, in conducting a large number of simultaneous hypothesis tests, the WRAF is prone to identifying false “discoveries.” A common technique for addressing this multiple testing problem is to adjust the procedure in a way that controls the proportion of true null hypotheses that are incorrectly rejected, or the false discovery rate (FDR). Unfortunately, generic FDR procedures suffer from low power when the hypotheses are dependent, and techniques designed to account for dependence are sensitive to misspecification of the underlying statistical model. In this paper, we develop a Bayesian decision-theoretical approach for dependent multiple testing and a nonparametric hierarchical statistical model that flexibly controls false discovery and is robust to model misspecification. We illustrate the robustness of our procedure to model error with a simulation study, using a framework that accounts for generic spatial dependence and allows the practitioner to flexibly specify the decision criteria. Finally, we apply our procedure to several seasonal forecasts and discuss implementation for the WRAF workflow. Lastly, supplementary materials for this article, including a standardized description of the materials available for reproducing the work, are available as an online supplement.

54 ENVIRONMENTAL SCIENCES↗