154 research outputs found
An Assessment of Data Transfer Performance for Large-Scale Climate Data Analysis and Recommendations for the Data Infrastructure for CMIP6
We document the data transfer workflow, data transfer performance, and other
aspects of staging approximately 56 terabytes of climate model output data from
the distributed Coupled Model Intercomparison Project (CMIP5) archive to the
National Energy Research Supercomputing Center (NERSC) at the Lawrence Berkeley
National Laboratory required for tracking and characterizing extratropical
storms, a phenomena of importance in the mid-latitudes. We present this
analysis to illustrate the current challenges in assembling multi-model data
sets at major computing facilities for large-scale studies of CMIP5 data.
Because of the larger archive size of the upcoming CMIP6 phase of model
intercomparison, we expect such data transfers to become of increasing
importance, and perhaps of routine necessity. We find that data transfer rates
using the ESGF are often slower than what is typically available to US
residences and that there is significant room for improvement in the data
transfer capabilities of the ESGF portal and data centers both in terms of
workflow mechanics and in data transfer performance. We believe performance
improvements of at least an order of magnitude are within technical reach using
current best practices, as illustrated by the performance we achieved in
transferring the complete raw data set between two high performance computing
facilities. To achieve these performance improvements, we recommend: that
current best practices (such as the Science DMZ model) be applied to the data
servers and networks at ESGF data centers; that sufficient financial and human
resources be devoted at the ESGF data centers for systems and network
engineering tasks to support high performance data movement; and that
performance metrics for data transfer between ESGF data centers and major
computing facilities used for climate data analysis be established, regularly
tested, and published
Quantile-based bias correction and uncertainty quantification of extreme event attribution statements
Extreme event attribution characterizes how anthropogenic climate change may
have influenced the probability and magnitude of selected individual extreme
weather and climate events. Attribution statements often involve quantification
of the fraction of attributable risk (FAR) or the risk ratio (RR) and
associated confidence intervals. Many such analyses use climate model output to
characterize extreme event behavior with and without anthropogenic influence.
However, such climate models may have biases in their representation of extreme
events. To account for discrepancies in the probabilities of extreme events
between observational datasets and model datasets, we demonstrate an
appropriate rescaling of the model output based on the quantiles of the
datasets to estimate an adjusted risk ratio. Our methodology accounts for
various components of uncertainty in estimation of the risk ratio. In
particular, we present an approach to construct a one-sided confidence interval
on the lower bound of the risk ratio when the estimated risk ratio is infinity.
We demonstrate the methodology using the summer 2011 central US heatwave and
output from the Community Earth System Model. In this example, we find that the
lower bound of the risk ratio is relatively insensitive to the magnitude and
probability of the actual event.Comment: 28 pages, 4 figures, 3 table
Quantifying statistical uncertainty in the attribution of human influence on severe weather
Event attribution in the context of climate change seeks to understand the
role of anthropogenic greenhouse gas emissions on extreme weather events,
either specific events or classes of events. A common approach to event
attribution uses climate model output under factual (real-world) and
counterfactual (world that might have been without anthropogenic greenhouse gas
emissions) scenarios to estimate the probabilities of the event of interest
under the two scenarios. Event attribution is then quantified by the ratio of
the two probabilities. While this approach has been applied many times in the
last 15 years, the statistical techniques used to estimate the risk ratio based
on climate model ensembles have not drawn on the full set of methods available
in the statistical literature and have in some cases used and interpreted the
bootstrap method in non-standard ways. We present a precise frequentist
statistical framework for quantifying the effect of sampling uncertainty on
estimation of the risk ratio, propose the use of statistical methods that are
new to event attribution, and evaluate a variety of methods using statistical
simulations. We conclude that existing statistical methods not yet in use for
event attribution have several advantages over the widely-used bootstrap,
including better statistical performance in repeated samples and robustness to
small estimated probabilities. Software for using the methods is available
through the climextRemes package available for R or Python. While we focus on
frequentist statistical methods, Bayesian methods are likely to be particularly
useful when considering sources of uncertainty beyond sampling uncertainty.Comment: 41 pages, 11 figures, 1 tabl
The effect of geographic sampling on evaluation of extreme precipitation in high resolution climate models
Traditional approaches for comparing global climate models and observational
data products typically fail to account for the geographic location of the
underlying weather station data. For modern high-resolution models, this is an
oversight since there are likely grid cells where the physical output of a
climate model is compared with a statistically interpolated quantity instead of
actual measurements of the climate system. In this paper, we quantify the
impact of geographic sampling on the relative performance of high resolution
climate models' representation of precipitation extremes in Boreal winter (DJF)
over the contiguous United States (CONUS), comparing model output from five
early submissions to the HighResMIP subproject of the CMIP6 experiment. We find
that properly accounting for the geographic sampling of weather stations can
significantly change the assessment of model performance. Across the models
considered, failing to account for sampling impacts the different metrics
(extreme bias, spatial pattern correlation, and spatial variability) in
different ways (both increasing and decreasing). We argue that the geographic
sampling of weather stations should be accounted for in order to yield a more
straightforward and appropriate comparison between models and observational
data sets, particularly for high resolution models. While we focus on the CONUS
in this paper, our results have important implications for other global land
regions where the sampling problem is more severe
Quantifying the effect of interannual ocean variability on the attribution of extreme climate events to human influence
In recent years, the climate change research community has become highly
interested in describing the anthropogenic influence on extreme weather events,
commonly termed "event attribution." Limitations in the observational record
and in computational resources motivate the use of uncoupled,
atmosphere/land-only climate models with prescribed ocean conditions run over a
short period, leading up to and including an event of interest. In this
approach, large ensembles of high-resolution simulations can be generated under
factual observed conditions and counterfactual conditions that might have been
observed in the absence of human interference; these can be used to estimate
the change in probability of the given event due to anthropogenic influence.
However, using a prescribed ocean state ignores the possibility that estimates
of attributable risk might be a function of the ocean state. Thus, the
uncertainty in attributable risk is likely underestimated, implying an
over-confidence in anthropogenic influence.
In this work, we estimate the year-to-year variability in calculations of the
anthropogenic contribution to extreme weather based on large ensembles of
atmospheric model simulations. Our results both quantify the magnitude of
year-to-year variability and categorize the degree to which conclusions of
attributable risk are qualitatively affected. The methodology is illustrated by
exploring extreme temperature and precipitation events for the northwest coast
of South America and northern-central Siberia; we also provides results for
regions around the globe. While it remains preferable to perform a full
multi-year analysis, the results presented here can serve as an indication of
where and when attribution researchers should be concerned about the use of
atmosphere-only simulations
Explaining the unexplainable: leveraging extremal dependence to characterize the 2021 Pacific Northwest heatwave
In late June, 2021, a devastating heatwave affected the US Pacific Northwest
and western Canada, breaking numerous all-time temperature records by large
margins and directly causing hundreds of fatalities. The observed 2021 daily
maximum temperature across much of the U.S. Pacific Northwest exceeded upper
bound estimates obtained from single-station temperature records even after
accounting for anthropogenic climate change, meaning that the event could not
have been predicted under standard univariate extreme value analysis
assumptions. In this work, we utilize a flexible spatial extremes model that
considers all stations across the Pacific Northwest domain and accounts for the
fact that many stations simultaneously experience extreme temperatures. Our
analysis incorporates the effects of anthropogenic forcing and natural climate
variability in order to better characterize time-varying changes in the
distribution of daily temperature extremes. We show that greenhouse gas
forcing, drought conditions and large-scale atmospheric modes of variability
all have significant impact on summertime maximum temperatures in this region.
Our model represents a significant improvement over corresponding
single-station analysis, and our posterior medians of the upper bounds are able
to anticipate more than 96% of the observed 2021 high station temperatures
after properly accounting for extremal dependence.Comment: 19 pages, 4 figures and 2 table
Evaluation of NASA\u27s MERRA Precipitation Product in Reproducing the Observed Trend and Distribution of Extreme Precipitation Events in the United States
This study evaluates the performance of NASA’s Modern-Era Retrospective Analysis for Research and Applications (MERRA) precipitation product in reproducing the trend and distribution of extreme precipitation events. Utilizing the extreme value theory, time-invariant and time-variant extreme value distributions are developed to model the trends and changes in the patterns of extreme precipitation events over the contiguous United States during 1979–2010. The Climate Prediction Center (CPC) U.S. Unified gridded observation data are used as the observational dataset. The CPC analysis shows that the eastern and western parts of the United States are experiencing positive and negative trends in annual maxima, respectively. The continental-scale patterns of change found in MERRA seem to reasonably mirror the observed patterns of change found in CPC. This is not previously expected, given the difficulty in constraining precipitation in reanalysis products. MERRA tends to overestimate the frequency at which the 99th percentile of precipitation is exceeded because this threshold tends to be lower in MERRA, making it easier to be exceeded. This feature is dominant during the summer months. MERRA tends to reproduce spatial patterns of the scale and location parameters of the generalized extreme value and generalized Pareto distributions. However, MERRA underestimates these parameters, particularly over the Gulf Coast states, leading to lower magnitudes in extreme precipitation events. Two issues in MERRA are identified: 1) MERRA shows a spurious negative trend in Nebraska and Kansas, which is most likely related to the changes in the satellite observing system over time that has apparently affected the water cycle in the central United States, and 2) the patterns of positive trend over the Gulf Coast states and along the East Coast seem to be correlated with the tropical cyclones in these regions. The analysis of the trends in the seasonal precipitation extremes indicates that the hurricane and winter seasons are contributing the most to these trend patterns in the southeastern United States. In addition, the increasing annual trend simulated by MERRA in the Gulf Coast region is due to an incorrect trend in winter precipitation extremes
Detecting Extreme Temperature Events Using Gaussian Mixture Models
Extreme temperature events have traditionally been detected assuming a
unimodal distribution of temperature data. We found that surface temperature
data can be described more accurately with a multimodal rather than a unimodal
distribution. Here, we applied Gaussian Mixture Models (GMM) to daily
near-surface maximum air temperature data from the historical and future
Coupled Model Intercomparison Project Phase 6 (CMIP6) simulations for 46 land
regions defined by the Intergovernmental Panel on Climate Change (IPCC). Using
the multimodal distribution, we found that temperature extremes, defined based
on daily data in the warmest mode of the GMM distributions, are getting more
frequent in all regions. Globally, a 10-year extreme temperature event relative
to 1985-2014 conditions will occur 13.6 times more frequently in the future
under 3.0{\deg}C of Global Warming Levels (GWL). The frequency increase can be
even higher in tropical regions, such that 10-year extreme temperature events
will occur almost twice a week. Additionally, we analysed the change in future
temperature distributions under different GWL and found that the hot
temperatures are increasing faster than cold temperatures in low latitudes,
while the cold temperatures are increasing faster than the hot temperatures in
high latitudes. The smallest changes in temperature distribution can be found
in tropical regions, where the annual temperature range is small. Our method
captures the differences in geographical regions and shows that the frequency
of extreme events will be even higher than reported in previous studies.Comment: 32 pages, 10 figure
Detected changes in precipitation extremes at their native scales derived from in situ measurements
The gridding of daily accumulated precipitation -- especially extremes --
from ground-based station observations is problematic due to the fractal nature
of precipitation, and therefore estimates of long period return values and
their changes based on such gridded daily data sets are generally
underestimated. In this paper, we characterize high-resolution changes in
observed extreme precipitation from 1950 to 2017 for the contiguous United
States (CONUS) based on in situ measurements only. Our analysis utilizes
spatial statistical methods that allow us to derive gridded estimates that do
not smooth extreme daily measurements and are consistent with statistics from
the original station data while increasing the resulting signal to noise ratio.
Furthermore, we use a robust statistical technique to identify significant
pointwise changes in the climatology of extreme precipitation while carefully
controlling the rate of false positives. We present and discuss seasonal
changes in the statistics of extreme precipitation: the largest and most
spatially-coherent pointwise changes are in fall (SON), with approximately 33%
of CONUS exhibiting significant changes (in an absolute sense). Other seasons
display very few meaningful pointwise changes (in either a relative or absolute
sense), illustrating the difficulty in detecting pointwise changes in extreme
precipitation based on in situ measurements. While our main result involves
seasonal changes, we also present and discuss annual changes in the statistics
of extreme precipitation. In this paper we only seek to detect changes over
time and leave attribution of the underlying causes of these changes for future
work
Recommended from our members
Western North Pacific tropical cyclone model tracks in present and future climates
Western North Pacific tropical cyclone (TC) model tracks are analyzed in two large multimodel ensembles, spanning a large variety of models and multiple future climate scenarios. Two methodologies are used to synthesize the properties of TC tracks in this large data set: cluster analysis and mass moment ellipses. First, the models' TC tracks are compared to observed TC tracks' characteristics, and a subset of the models is chosen for analysis, based on the tracks' similarity to observations and sample size. Potential changes in track types in a warming climate are identified by comparing the kernel smoothed probability distributions of various track variables in historical and future scenarios using a Kolmogorov-Smirnov significance test. Two track changes are identified. The first is a statistically significant increase in the north-south expansion, which can also be viewed as a poleward shift, as TC tracks are prevented from expanding equatorward due to the weak Coriolis force near the equator. The second change is an eastward shift in the storm tracks that occur near the central Pacific in one of the multimodel ensembles, indicating a possible increase in the occurrence of storms near Hawaii in a warming climate. The dependence of the results on which model and future scenario are considered emphasizes the necessity of including multiple models and scenarios when considering future changes in TC characteristics
- …