A document by Shiheng Duan. Click on the document to view its contents.

We analyze the large-scale statistically meaningful patterns (LSMPs), also called large-scale meteorological patterns, that precede extreme precipitation (PEx) events over Northern California (NorCal). We find LSMPs by applying k-means clustering to the two leading principal components of daily 500hPa geopotential height anomalies persisting two days before the onset. A statistical significance test based on the Monte Carlo simulations suggests the existence of a minimum of four statistically distinguished LSMP clusters. The four LSMP clusters are characterized as the NW continental negative height anomaly, the Eastward positive “PNA”, the Westward negative “PNA”, and the Prominent Alaskan ridge. These four clusters, shown in multiple atmospheric and oceanic variables, evolve very differently and have distant links to the Arctic and tropical Pacific regions. Using binary forecast skill measures and a new copula-based framework for predicting PEx events, we show that the LSMP indices are useful predictors of NorCal PEx events, with the moisture-based variables being the best predictors of PEx events at least six days before the onset, and the lower atmospheric variables being better than their upper atmospheric counterparts any day in advance.

The 1997 New Year’s flood event was the most costly in California’s history. This compound extreme event was driven by a category 5 atmospheric river that led to widespread snowmelt. Extreme precipitation, snowmelt, and saturated soils produced heavy runoff causing widespread inundation in the Sacramento Valley. This study recreates the 1997 flood using the Regionally Refined Mesh capabilities of the Energy Exascale Earth System Model (RRM-E3SM) under prescribed ocean conditions. Understanding the processes causing extreme events inform practical efforts to anticipate and prepare for such events in the future, and also provides a rich context to evaluate model skill in representing extremes. Three California-focused RRM grids, with horizontal resolution refinement of 14km down to 3.5km, and six forecast lead times, 28 December 1996 at 00Z through 30 December 1996 at 12Z, are assessed for their ability to recreate the 1997 flood. Planetary to synoptic scale atmospheric circulations and integrated vapor transport are weakly influenced by horizontal resolution refinement over California. Topography and mesoscale circulations, such as the Sierra barrier jet, are prominently influenced by horizontal resolution. The finest resolution RRM-E3SM simulation best represents storm total precipitation and storm duration snowpack changes. Traditional time-series and causal analysis frameworks are used to examine runoff sensitivities state-wide and above major reservoirs. These frameworks show that horizontal resolution plays a more prominent role in shaping reservoir inflows, namely the magnitude and time-series shape, than forecast lead time, 2-to-4 days prior to the 1997 flood onset.

The water cycle is an important component of the earth system and it plays a key role in many facets of society, including energy production, agriculture, and human health and safety. In this study, the Energy Exascale Earth System Model version 1 (E3SMv1) is run with low-resolution (roughly 110 km) and high-resolution (roughly 25 km) configurations — as established by the High Resolution Model Intercomparison Project protocol — to evaluate the atmospheric and terrestrial water budgets over the conterminous United States (CONUS) at the large watershed scale. The water cycle slows down in the HR experiment relative to the LR, with decreasing fluxes of precipitation, evapotranspiration, atmospheric moisture convergence, and runoff. The reductions in these terms exacerbate biases for some watersheds, while reducing them in others. For example, precipitation biases are exacerbated at HR over the Eastern and Central CONUS watersheds, while precipitation biases are reduced at HR over the Western CONUS watersheds. The most pronounced changes to the water cycle come from reductions in precipitation and evapotranspiration, the latter of which results from decreases in evaporative fraction. While the HR simulation is warmer than the LR, moisture convergence decreases despite the increased atmospheric water vapor, suggesting circulation biases are an important factor. Additional exploratory metrics show improvements to water cycle extremes (both in precipitation and streamflow), fractional contributions of different storm types to total precipitation, and mountain snowpack.

The Atmospheric River (AR) Tracking Method Intercomparison Project (ARTMIP) is a community effort to systematically assess how the uncertainties from AR detectors (ARDTs) impact our scientific understanding of ARs. This study describes the ARTMIP Tier 2 experimental design and initial results using the Coupled Model Intercomparison Project (CMIP) Phases 5 and 6 multi-model ensembles. We show that AR statistics from a given ARDT in CMIP5/6 historical simulations compare remarkably well with the MERRA-2 reanalysis. In CMIP5/6 future simulations, most ARDTs project a global increase in AR frequency, counts, and sizes, especially along the western coastlines of the Pacific and Atlantic oceans. We find that the choice of ARDT is the dominant contributor to the uncertainty in projected AR frequency when compared with model choice. These results imply that new projects investigating future changes in ARs should explicitly consider ARDT uncertainty as a core part of the experimental design.

In 2021, Nepal underwent its most severe fire season, resulting in a fire rate 10 times greater than the historical average in many areas of the country with record-high air pollution levels. Leading the fire outbreaks in March of 2021, the country experienced an extreme precipitation deficit and drought in the post-monsoon season. Current community forest management practices and resultant forest growth may have exacerbated the conflagration, but an analysis using observational, reanalysis, and climate model ensemble data indicates that climate variability and climate change induced severe drought conditions that resulted in the anomalous fire season. While warning of the likely re-occurrence of extremely active fire seasons in Nepal through the end of the 21st century, this research also proposes a statistical model for sub-seasonal prediction that could help mitigate the projected effects of the drought-fire paradigm.

Atmospheric rivers (ARs) are long and narrow filaments of vapor transport responsible for most poleward moisture transport outside of the tropics. Many AR detection algorithms have been developed to automatically identify ARs in climate data. The diversity of these algorithms has introduced appreciable uncertainties in quantitative measures of AR properties and thereby impedes the construction of a unified and internally consistent climatology of ARs. This paper compares eight global AR detection algorithms from the perspective of AR life cycles following the propagation of ARs from origin to termination in the MERRA2 reanalysis over the period 1980-2017. Uncertainties related to lifecycle characteristics, including number, lifetime, intensity, and frequency distribution are discussed. Notably, the number of AR events per year in the Northern Hemisphere can vary by a factor of 5 with different algorithms. Although all algorithms show that the maximum origin (termination) frequency locates over the northwestern (northeastern) Pacific, significant disagreements occur in regional distribution. Spreads are large in AR lifetime and intensity. The number of landfalling AR events produced by the algorithms can vary from 16 to 78 events per cool season, i.e. by almost a factor of five, although the agreement improves for stronger ARs. By examining the AR’s connection with the Madden-Julian Oscillation and El Niño Southern Oscillation, we find that the overall responses of ARs (such as changes in AR frequency, origin, and landfalling activity) to low-frequency climate variabilities are consistent among algorithms.

Increasing wildfire and declining snowpacks in mountain regions threaten water availability. We combine satellite-based fire detection with snow seasonality classifications to examine fire activity in California’s seasonal and ephemeral snow areas. We find a nearly tenfold increase in fire activity during 2020 and 2021 compared to 2001-2019 as measured by satellite data. Accumulation season snow albedo declined 17-77% in two burned sites as measured by in-situ data relative to un-burned conditions, with greater declines associated with increased soil burn severity. By enhancing snowpack susceptibility to melt, decreased snow albedo drove mid-winter melt during a multi-week midwinter dry spell in 2022. Despite similar meteorological conditions in 2013 and 2022, which we link to persistent high pressure weather regimes, minimal melt occurred in 2013. Post-fire differences are confirmed with satellite measurements. Our findings suggest larger areas of California’s snowpack will be increasingly impacted by the compounding effects of dry spells and wildfire.