Climate Variability of Atmospheric Rivers and Droughts over the West Coast of the United States from 2006 to 2019

: Water resources are crucial to the livelihood and sustainability of the general public across the western United States. This study covers the timespan of both the third driest drought in Californian history between 2012 and 2015 as well as the extreme atmospheric river year in 2016–2017. The evaluation of vertical moisture proﬁles using Constellation Observing System for Meteorology, Ionosphere, and Climate (COSMIC) Radio Occultation (RO) data, National Centers for Environmental Prediction (NCEP)/National Center for Atmospheric Research (NCAR) Reanalysis of 500 hPa geopotential heights, 1000–500 hPa thickness, Optimum Interpolation (OI) Sea Surface Temperature (SST), NOAA/NDBC buoy data, and NASA, MEaSUREs, Gridded Sea Surface Height Anomalies (SSHA) were performed. The daily COSMIC time evolution from 2006 through 2015 showed a ﬂat to slightly upward trend of both temperature and water vapor proﬁles through the entirety of the western US drought. Subsequently, a signiﬁcant increase of temperatures and water vapor were recorded in early 2016 before the extreme Atmospheric River (AR) season of 2016–2017. The quantitative analyses suggest that warmer SST and higher SSHA lead to an increase of heat ﬂuxes from the ocean into the troposphere, which forces thickness changes and thus the position of troughs in the geopotential height ﬁeld changes afterwards, consequently pushing the trough eastward over the Paciﬁc Northwest and potentially leading to an active AR year in the western US. It appears that regional COSMIC RO moisture proﬁles, seasonal SST, and SLH anomalies may serve as a precursor for seasonal or sub-seasonal precipitation outlook along the western US. we managers verify Paciﬁc wave-pattern A future research study with COSMIC-2 and monthly climatology of SST and SLH data promises future developments for seasonal forecasting. The warranted study could include the relationship between several climate indices (i.e., El Niño Southern Oscillation (ENSO), Arctic Oscillation (AO), and Paciﬁc Decadal Oscillation (PDO)) and COSMIC mixing ratios. The potential linkage with climatological indices could help improve the forecast of enhanced AR years and droughts nearly three months ahead of a given season. In summary, we recommend continued research using COSMIC-2 data within the Weather Research and Forecasting (WRF) model to provide a lasting beneﬁt and better understanding of AR forecasting.


Introduction
Atmospheric Rivers (ARs) are responsible for over 90% of the poleward water vapor transport, while only covering about 10% of the global longitudinal area (Neiman et al. 2008b; Rutz et al. 2014) [1,2]. Roughly 60-74% and 33-40% of droughts in the Pacific Northwest and California, respectively, have ended due to the arrival of AR storms (Dettinger 2013) [3]. Atmospheric River events are an essential source of moisture and precipitation that are observed along the west coast and are a key factor in the prevalence of droughts and floods. Nevertheless, water resources are critical issues that California is facing today and will continue to face amidst the threat of global climate change. In California, three-fourths of the population and their water demands are primarily in the central and southern parts of the state, while two-thirds of the precipitation and runoff occurs in the northern one-third of the state (Dettinger 2011) [4]. Furthermore, global environmental changes have potentially increased flooding threats along the western United States after multi-year droughts (Dettinger 2011) [4]. Considering the ability to end droughts and the geographic disparity of California's resources, understanding the strength and climatology of ARs and their impact on the west coast is essential for water management, reservoir operations, and the mitigation of flood risks.
California from 2006 to 2019 in order to discern if the strength of ARs are increasing compared to previous years and also investigate the variability of droughts. Comparisons will then be made to a 13-year record of COSMIC data along the eastern Pacific coast. The methods used to provide these comparisons will be described in Section 2, followed by data and results in Section 3, and our concluding results and recommendations for future work can be seen in Section 4.

COSMIC Vertical Profiles
University Corporation for Atmospheric Research (UCAR's) Constellation Observing System for Meteorology, Ionosphere, and Climate (FORMOSAT-3/COSMIC) mission was launched in 2006 as a joint research project by Taiwan and the United States [15]. COSMIC consists of a constellation of six small polar orbiting satellites that are distributed equally around the Earth, providing users with extensive data coverage over data-sparse locations including oceans and inhabited land. Furthermore, Ho et al. (2019a) [16] states that, through 2016, COSMIC has offered users approximately 1000 soundings a day and one sounding for each 5 • × 5 • latitudinal box. Radio occultation (RO) employs active satellites that use Global Navigation Satellite System Radio Occultation (GNSS-RO) measurements to accurately vertically profile the Earth's atmosphere (Ho et al. 2019a) [16].
Each satellite consists of three instruments, which include the following: a GPS radio occultation receiver that receives radio waves from a GPS satellite, an ionospheric photometer that measures the strength of electromagnetic waves throughout the Earth's ionosphere, and a tri-band beacon (e.g., Very High Frequency, Ultra High Frequency, and L-Band) which is a device that provides the satellite with a bearing by locating the radio waves Constellation Observing System for Meteorology Ionosphere, and Climate (COSMIC) 2020; UCAR Community Programs 2020) [17,18]. COSMIC satellites intercept radio signals from other GNSS satellites in orbit and measures their bending and signal delay (UCAR Community Programs 2020) [18]. After refracting the radio waves, COSMIC then relays the data down to the COSMIC Data Analysis and Archive Center (CDAAC) with vertical profiles of the following parameters: temperature, pressure, water vapor, bending angles, and refractivity (UCAR Community Programs 2020) [18].
The target area that was analyzed in this study, as shown in Figure 1, consists of the eastern Pacific Ocean and the western mid-latitudes from the Washington/Canadian border through California ( (1)) values from the lowest level to 5 km above ground level.
The lowest level to 5 km profiles were evaluated in this study (denoted as 0-5 km) to depict how the moisture and temperature profiles changed between the surface and the midatmosphere. In order to improve interpretation in the study, each COSMIC plot included the daily value (marked as red dots) and running means of weekly, monthly, seasonal, and annual trends. COSMIC data after March 2019 was excluded in this study because only one of the six satellites remained online, consequently leading to a plethora of missing data and an abundance of inaccuracies in the dataset. COSMIC data was made available by the Data Analysis and Archive Center Constellation Observing System for Meteorology Ionosphere and Climate (University Corporation for Atmospheric Research) [19].

Regional Climate Characteristics
Global reanalysis from the NCEP/NCAR Reanalysis Project, Kalnay et al. (1996) [20], were used in this study to plot and analyze 500 hPa geopotential heights and 1000-500 hPa thickness. Each variable studied will portray autumn (SON) and winter (DJF) seasonal anomalies for the purpose of comparing seasonal and sub-seasonal variations in the dataset. In order to develop a 30-year anomaly dataset, the researched seasonal anomalies were calculated by taking the mean of each seasonal group (i.e., SON and DJF) from January 1989 through March 2019. Quantitative tables were developed to average out each of the variables in their domains to provide future statistical analysis. Additionally, the original target area was then divided into four sub-regions consisting of a 10 • × 10 • latitudinal box to display the spatial distributions of each analyzed variable (See Figure 1) [21] were employed to evaluate the sea surface height variations over the target area. MEaSUREs provides a high resolution global (0.17 • × 0.17 • ) analysis dataset with a 5-day temporal resolution of sea surface height and sea surface height anomalies. The NASA JPL MEaSUREs sea level heights are constructed by the kriging method from two simultaneous sets of altimetric satellites (MEaSUREs 2015) [21]. The first set of simultaneous satellites include the TOPEX/Poseidon, Jason-1, Jason-2, and Jason-3, while the second set of satellites include the ERS-1, ERS-2 Envisat, AltiKa, Cryosat-2, and Sentinel 3A. Furthermore, in terms of accuracy of sea level heights, MEaSUREs (2015) [21] shows a reliable solution of approximately −1 to 0 cm error term from 1993 to 1999 and a positive 0 to 1 cm error through 2015.
From an in-situ observations perspective, NOAA National Data Buoy Center (NDBC) offers quality measurements across broad marine environments. NOAA NDBC (2020) [22] states that, since the early 1980's, buoys have been deployed to spatially cover the western Atlantic to the northeastern Pacific. Furthermore, NOAA NDBC (2020) [22] has installed 50 Coastal-Marine Automated Network (C-MAN) buoys and over 100 meteorological buoys that provide information on barometric pressure, wind direction, wind speed, wind gust, and air temperature.
The NOAA 1/4 • daily Optimum Interpolation (OI) Sea Surface Temperature (SST) V2 data (Zlotnicki 2019) [23] was used to evaluate the annual and seasonal differential of sea surface temperatures throughout the domain. OISST data is constructed by combining observations from different platforms (satellites, ships, buoys, and Argo floats) on a regular global grid. A spatially complete SST map is produced by interpolating to fill in gaps. The methodology includes bias adjustment of satellite and ship observations (referenced to buoys) to compensate for platform differences and sensor biases (Zlotnicki 2019) [23].
The Climate Prediction Center (2020) [24] ENSO amplitude archive was used to monitor El Niño and La Niña conditions. ENSO amplitudes are found by taking a 3-month running mean of SST anomalies between 5 • N and 5 • S to 120 • W and 170 • W (Climate Prediction Center 2020) [24]. In addition to the amplitudes, El Niño/La Niña is broken down into categories to portray the severity of the event. Therefore, a weak ENSO amplitude has 0.5 to 0.9 SST anomaly, moderate 1 to 1.4, strong 1.5 to 1.9, and very strong when SST is greater than 2 (Null 2020) [25].
The National Weather Service (NWS) (Advanced Hydrologic Prediction Center 2020) [26] rain gauges were used to evaluate the daily precipitation accumulation around the San Francisco Bay Area (e.g., region 4 in Figure 1). These rain gauges are well maintained and accurate; however, inaccuracies can occur due to the following: freezing precipitation, windy conditions, and obstructions (Advanced Hydrologic Prediction Center 2020) [26].

COSMIC Moisture and Temperature Observations and Annual Features
Time evolution of COSMIC moisture profiles were conducted in association with the seasonal and sub-seasonal variability of mixing ratio and temperature characteristics of the target region. Additionally, heat maps of moisture were added to portray any seasonal trends or patterns that were observed during the study.
The first COSMIC variable investigated in this study was mixing ratio. When visualizing the monthly averages, the seasonality of the mixing ratio is clearly observed, with the highest peak found in the summer months. The summer month peaks are due to an increase of solar radiation to the surface, consequently forcing an increase of latent heat release into the atmosphere. As shown in Figure 2, the historical time series depicts a slight decreasing trend~−0.19 g/kg (~−0.047 g/kg/year) of vertical profile mixing ratio from July 2006 through early 2013. This slight decrease of mixing ratio between 2011 through 2014 is analogous to the extreme western United States drought, which greatly reduced the moisture field upstream of the resilient ridge in the Pacific Northwest. Between late 2015 and 2016, the mixing ratio trend line showed the steepest climb of~0.32 g/kg (~0.162 g/kg/year). Afterwards, starting in October 2016 through early 2017, the annual mixing ratio trend line peaked and remained constant at~2.8 g/kg. In 2018, the trend line appeared to decrease again by~0.2 g/kg and flatlined around~2.6 g/kg. The observed decrease of mixing ratio in 2018 was directly attributable to a persistent ridge hovering along the entire west coast, forcing storm tracks northward.
The magnitude of mixing ratio was processed via heatmaps in order to further analyze the trends and patterns of mixing ratios in the target region. The analysis includes both Figure 3a normalization of mixing ratios and Figure 3b mixing ratio anomalies. Both heatmaps were limited to October through March in order to focus on how the mixing ratio changed during the climatological west coast rain season. As shown in Figure 3a The key finding observed in the heatmaps (red dashed line) was a possible wavelike oscillation. Future investigation will be needed to determine if the observed wave maintains the same wave period and if the wave can be directly influenced or linked to any climatological oscillations. This future research could provide a possible summary of when the western United States may see the next drought and/or extreme AR rain season.  The second COSMIC variable that was investigated was the 0-5 km (i.e., lowest level to 5 km above ground level (AGL)) temperatures ( Figure 4). Similar to patterns seen in the moisture field, temperatures in the domain decreased by~−2.10°C (~−0.52 • C/year) from June 2006 through late 2013. In order to compare the average 5 km profile with in-situ surface observations, eight sets of buoy data from NDBC were analyzed ( Figure 5). Buoy data showed the warmest surface temperatures in 2014 through 2015, however, the COSMIC temperature profile observed the warmest temperatures in 2016. According to Phillips and Leslie (2017) [27], 2016 global average temperatures witnessed the warmest year in the last 137 years of recorded data. Furthermore, they suggested that a strong El Niño in early 2016 directly increased the global temperature by approximately 0.5°C from the 1981-2010 mean (Phillips and Leslie 2017) [27]. The enhanced global warming in 2016 helped induce an increase of latent heat release and consequently affected western United States precipitation in the 2016-2017 rain season. After the above average rain season in 2016-2017, both COSMIC and buoy temperatures showed a decrease by approximately 1 • C. Potentially, the cooler temperatures observed by COSMIC in late 2017 forced a decrease of latent heat into the atmosphere. Therefore, with less moisture and a persistent 500 hPa ridge along the west coast, California saw a tremendous reduction of moisture leading to a below-average rain season in 2017-2018. Temperatures within the target area rebounded in early 2018, consequently evaporating more water into the atmosphere, leading to another above average rain season in 2018-2019.

Regional Analysis
This section explores the reanalysis data for a few important large scale meteorological factors. Analysis will include seasonal anomaly maps from 2006 through 2019 of 500 hPa geopotential heights with 1000-500 hPa thickness, sea level heights (SLH), and sea surface temperatures (SSTs). Additionally, for each variable analyzed, a quantitative table was constructed to organize the data and to highlight potential trends and patterns.

500 hPa Geopotential Heights and 1000-500 hPa Thickness
The results of 500 hPa geopotential heights with 1000-500 hPa thickness [20] overlaid are shown in Figure 6a [29] suggested that the downstream trough that brought cold outbreaks in the Midwest, which was seen during the winter months in 2013 ( Figure 6), was deepened by the Pacific ridge. In the autumn of 2013, weak negative height and thickness anomalies were found in California due to the remnants of a late summer monsoon; however, drought conditions returned to the desert southwest in DJF and persisted into the spring months. In 2015, geopotential height anomalies dropped back to 2006-2011 climatological average values. An eastern push of the 500 hPa trough was observed in the autumn of 2016 west of Washington, with anomalies bottoming out at −6.66 decameter (dm) in region 1, −3.98 dm in region 2, and −1.97 dm in region 3, and remaining near normal (−0.09 dm) in region 4. (Table 1a). Thickness anomalies in SON 2016 saw a −2.78 dm anomaly in region 1, −1.75 dm in region 2, −1.10 dm in region 3, and −0.10 dm in region 4 (Table 1c). This eastward track of the upper-level longwave trough was crucial in the development of an above average rain season in 2016. In DJF 2016, low geopotential height and thickness anomalies persisted along the northwestern United States with −3.78 dm in region 1, −5.67 dm in region 2, −2.03 dm in region 3, and −2.00 dm in region 4 (Table 1b). In terms of thickness during DJF 2016, region 1 observed a −4.07 dm anomaly, while this was −4.09 dm in region 2, −0.59 dm in region 3, and −0.69 dm in region 4 (Table 1d). In 2017, SON observed near normal height anomalies (Table 1a); however, there was the return of a resilient ridge with positive anomalies in DJF around positive 5 to 7 dm. The resilience of the winter ridge in 2017-2018 inevitably pushed the storm track further into Canada, hence dramatically decreasing the western United States' precipitation. Subsequently, in 2018, near-zero anomalies were found in autumn; however, a negative 2-4 dm anomaly developed. The wintertime trough helped drop the storm track further south to the California-Oregon border in late winter and spring, which yielded multiple atmospheric river events and an above average rain season in 2018-2019.

Sea Level Heights
The second variable analyzed within the target area was NASA MEaSUREs [21] sea level heights (Figure 7a,b and Table 2). In this section, the target area was broken into four sub-regions to quantify how the sea level heights changed and to explain which location may be critical in enhancing western US precipitation.  (Table 2b). To further analyze the sub-regional changes in sea level heights per year, the average seasonal percentage change was calculated ( Table 3).
As shown in Table 3, a positive percentage change rate between 2012 and 2015 coincides with the beginning of the 3-year long west coast drought. In 2013, sea level heights observed the biggest percentage change of 8650%, and the ocean surface warmed. This result is consistent with the regional average buoy and COSMIC temperatures rise during the west coast drought. In the period 2014-2015, both SON (0.111 m) and DJF (0.061 m) observed the highest sea level heights, with an annual average of 0.086 m. In 2016, SON average sea level height anomalies measured 0.072 m; however, in contrast to high heights in 2015, 2016 DJF decreased substantially to 0.001 m. The tremendous decrease of sea level heights in DJF was directly related to an extreme atmospheric river rain season due to an increase in the total amount of storms. Heights continued to fall by another 40% through the end of 2017. After a dry and warm 2018, another bump of 78% in sea level heights can be seen, which positively influenced another above average precipitation year in 2018-2019.
To verify the patterns and trends within the sea level height data, each of the four subdomains were compared using an annual time series (Figure 8). A noteworthy point is that the sea level heights annual trend line shows a similar possible oscillation to the COSMIC mixing ratio (Figure 3). This observation was crucial for showing a potential link between sea level heights and COSMIC mixing ratio due to the thermal expansion/compression of the ocean (e.g., Widlansky et al. 2020) [30]. In order to explain the variations of SLH

Sea Surface Temperatures
The third and final variable, analyzed in Figure 9a,b, is sea surface temperatures (SSTs) from the NOAA 1/4 • daily Optimum Interpolation (OI) Sea Surface Temperature (SST) Version 2 data. Sea surface temperatures between 2006 and 2012 observed negative anomalies between 0 and −1 • C. In the period 2013-2014, the large-scale pattern in both autumn and winter shifted from negative anomalies and became positive (Figure 9). The eastern Pacific detected a push of warmer water poleward; however, it persisted away from the coast in regions 1 and 3. In 2014, SST anomalies for both autumn and winter warmed between positive 1-2 • C and persisted through the DJF of 2016 (Table 4. El Niño typically leads to above average precipitation in California; however, the west coast observed near normal precipitation amounts in 2014 and 2015. According to the Oceanic Niño Index (Table 5) [31,32] suggested that one possible reason for this was due to the internal atmospheric variability. The peak of SST occurred in 2015, which was verified with the highest annual average amplitude of 2.5 in the ENSO Oscillation [24]. In 2016, a return of cooler SSTs (−1 to −1.5 • C) in regions 1 and 2 and negative (−0.7 to −0.3) ENSO amplitudes were found. Furthermore, the Oceanic Niño Index (Table 5) [24], 2016-2017 saw a weak La Niña (SST anomalies −0.5 to −0.9), and a strong La Niña (SST anomalies −1.5 to −1.9) was observed in 2017-2018. The drop off of SST prior to the extreme 2016-2017 rain season can be explained by a more atmospheric river events and the occurrence of storms, which consequently led to more latent heat release into the atmosphere from the ocean surface. Therefore, the extra heat from the ocean surface that was transferred into the atmosphere would directly lead to an increase of moisture into the atmosphere, which is consistent with the COSMIC vertical profiles shown earlier (Figure 2). In 2017, the eastern Pacific observed another push of warm positive anomalies poleward; however, the mid-level ridge ( Figure 6) did allow the transport of moisture to the west. On the other hand, in 2018, similar warm SST anomalies were found within the target region, but the mid-level trough tapped into the subtropical moisture, consequently producing another active rain season in 2018-2019.

Statistics-Linear Regression
Linear regression was performed to identify the statistical significance of each variable that may have influenced the western United States precipitation. Analysis in this section compares the annual average of the four subdomains for sea surface temperature v. sea level heights (Figure 10a) and the NASA MEaSUREs sea level heights v. COSMIC mixing ratio (Figure 10b). Sea surface temperatures and sea level heights have a strong correlation coefficient of R 2 = 0.50. The direct influence of SST and SLH are related to expansion and/or contraction as ocean temperature increases and/or decreases. Sea level height and mixing ratio had the highest correlation coefficient of R 2 = 0.65. This finding suggests that, as sea level heights increase, COSMIC's mixing ratio increases, verifying that there is a direct impact to the overall moisture field forced by more latent heat released into the atmosphere.

Atmospheric Rivers and Precipitation
The final set of analysis incorporates the Hecht and Ralph (2020) [33][34][35] classification of the total number of AR storms sorted by the location where they made landfall since 2016 and comparing it to the reanalysis data ( Table 6). As explained using the geopotential height and thickness fields, in 2016 the eastward push of the 500 hPa trough to the Pacific Northwest ( Figure 6), directly affected the total number of AR storms and their landfalling location.    (Figure 6). In addition to a persistent ridge, the target area observed a decrease in both sea surface temperatures ( Figure 9) and sea level heights (Figure 7). In 2018-2019, another eastward push of the 500 hPa geopotential heights and thicker and warmer SSTs allowed for another above average rain season, with a total of 41 AR events occurring. The 2018-2019 rain season saw similar results to the extreme AR year of 2016-2017 in terms of predominant landfalling locations in northern California and Oregon. This can be explained by the dip southward of the mid-level trough, which allowed ARs to track further south ( Figure 6).  To further verify precipitation and AR events, a comparison between the total number of ARs with the total annual precipitation for five northern and central California locations from 1 October through 31 March was analyzed (NCEI 2020) [36]. In Figure 11, each of the five locations observed substantially wide ranges of total annual precipitation due to a variety of orographic effects. San Jose, which is located in a rain shadow of the Santa Cruz Mountains, received about 10.7 inches on average over the 13-year period (Figure 11a). Next, Santa Cruz, a beach city southwest of the Santa Cruz Mountains that observes some orographic enhancement, received approximately 24.1 inches annually (Figure 11b). Similar to Santa Cruz's orographic effects, Monterey, a beach city along the California central coast, collected around 13.9 inches (Figure 11c). San Francisco International Airport (KSFO), located along the San Francisco Peninsula, obtained roughly 15.4 inches per year (Figure 11d). Finally, Big Sur, a mountainous part of the California central coast, acquired about 33.2 inches of precipitation every year (Figure 11e). A key finding revealed from annual precipitation data was that, when the annual running mean was calculated, a similar wave-like pattern could be observed, which is similar to the COSMIC mixing ratio heatmap ( Figure 3) and the northeastern Pacific sea level heights (Figure 8).

Discussion and Summary
Annual precipitation, regional reanalysis of geopotential height, 1000-500 hPa thickness, sea surface height anomalies, sea surface temperature, and COSMIC moisture and temperature profiles were evaluated in this study from 2006 to 2019. The study included the third worst drought, from 2012-2014, and compared and contrasted the dynamical changes to the extreme AR wet seasons of 2016-2017 in California. This study was designed to evaluate the COSMIC data and how it can be used to improve our current understanding of AR occurrence along the west coast. COSMIC has provided a quantification of high-resolution vertical profiles in terms of pressure, vapor pressure, and temperature throughout the mid-latitudes and tropics. This quantification assists and builds on the scientific knowledge concerning how likely it is that the regional changing climate has directly affected the moisture and temperature fields. Understanding the prior thirteen years Results suggested that increasing global temperatures tend to lead to an increase of heat fluxes into the atmosphere, consequently increasing moisture within the target area. Therefore, moisture from an increase of latent heat release by warmer SSTs, as suggested by Chen and Leung (2020) and Dai et al. (1998) [8,9] and the total frequency and strength of AR could potentially increase, as was seen on the west coast in 2016-2017. In 2018, mixing ratios dipped bỹ 0.2 g/kg and stayed constant till the end of the studied period in March of 2019. Another interesting finding using the COSMIC mixing ratio heatmaps from October to March was that they revealed a wave-like pattern. In order to show how increasing moisture changes AR frequency and strength, the total number of ARs and the total annual precipitation for five California cities were compared to the COSMIC mixing ratio data ( Figure 11). It can be explained by the winter trough during the active wet seasons. In above average rain seasons, the trough transports moisture into the atmosphere from warmer SSTs and less moisture and convection from cooler SSTs during droughts. A noteworthy feature from the annual precipitation analysis showed another similar wave pattern to the COSMIC mixing ratios (Figure 3).
The analysis of climatological regional scale factors found that the most efficient and useful variable relating to AR frequency was sea level height anomalies. Furthermore, sea level height anomalies had the highest correlation coefficient of R 2 = 0.65 with the COSMIC mixing ratio. Consistent with meteorological theory, sea level heights rose due to the enhanced warming during the drought, consequently increasing evaporation and providing more vapor to the atmosphere. The link between sea surface temperature and mixing ratio has been asserted for years. This study found that 50% of the variation of sea surface temperature (R 2 = 0.50) was attributed to sea level heights. Because sea surface temperatures lag behind the seasonal air temperature due to the high specific heat capacity of water, it appears that having a warmer summer and autumn can allow more moisture to be lofted into the atmosphere for the following rain season. Sea surface temperature can directly impact the location of mid-level troughs and ridges, as observed within the geopotential height and thickness fields. In terms of the drought, cooler sea surface temperatures prevailed, which allowed a persistent 500 hPa ridge to advect over the western United States and lead to a below average rain season. On the other hand, during the strong AR year in 2016-2017, the west observed warmer sea surface temperatures that pushed the mid-level trough eastward toward the Pacific Northwest, allowing for the surface low to propagate over the west coast, leading to an above-average rain season.
The findings from this study provide new information to the scientific community about which atmospheric factors are beneficial to droughts and ARs. In this experiment, we acknowledge that 13 years of COSMIC satellite data may not be sufficient to analyze the climate variability of droughts and ARs. In order to extend the studied timespan, future experiments should include the new COSMIC-2 satellite, which was launched in June 2019. Potentially in the future, this research should help forecasters and water managers develop seasonal forecasts. Therefore, continued work moving the domain to another location to verify the northeastern Pacific sinusoidal wave-pattern is recommended. A future research study with COSMIC-2 and monthly climatology of SST and SLH data promises future developments for seasonal forecasting. The warranted study could include the relationship between several climate indices (i.e., El Niño Southern Oscillation (ENSO), Arctic Oscillation (AO), and Pacific Decadal Oscillation (PDO)) and COSMIC mixing ratios. The potential linkage with climatological indices could help improve the forecast of enhanced AR years and droughts nearly three months ahead of a given season. In summary, we recommend continued research using COSMIC-2 data within the Weather Research and Forecasting (WRF) model to provide a lasting benefit and better understanding of AR forecasting.