Asian Long-Range Transport in Relation to Atmospheric Rivers in Northern California

: The study investigates the e ﬀ ect of aerosol long-range transport on precipitation over Northern California during atmospheric river (AR) events in the 2017 cold season (January–April). ARs in 2017 were one of the strongest to date, and the intense precipitation associated with the ARs resulted in ﬂooding, destruction of property, and contamination of water supplies. The Aerosol Optical Depth (AOD) from Moderate Resolution Imaging Spectroradiometer (MODIS) data shows Asian dust traveling across the Northern Paciﬁc Ocean along with AR events. Aerosol measurements in California, provided by the Interagency Monitoring of Protected Visual Environments (IMPROVE), show that more Asian dust tends to be observed over the coast, while non-Asian / localized dust is observed inland. A mixture of Asian and localized dust is observed over the mountains, although higher amounts of both are observed in the spring (March–April). Back trajectory analysis conﬁrms that Asian aerosols are transported along the air parcels, and each AR event has its own transport pattern in terms of horizontal advection and vertical lifting. Correlation between precipitation and aerosols is low. This suggests that aerosols contribute little to the decrease of local precipitation during the 2017 AR events.


Introduction
In April 2017, California lifted its State of Emergency that had been declared in January 2014 in response to a five-year drought [1]. The 2017 water year alone (October 2016-April 2017) saw significant amounts of rain, with the majority of station readings in California reading above 130% of normal levels. For example, Los Angeles International Airport reached 133% of its normal rain total, compared to the previous year at 55% of its normal level [2]. While the rain alleviated drought conditions within the state, the intense precipitation resulted in a flooding of cities and freeways, mudslides and mountain trail erosions, destruction of property, and contamination of water supplies; all of which threaten public safety.
The motivation for the study is to improve precipitation forecasts and prepare for future extreme events by understanding the dynamics and characteristics of atmospheric phenomena that contribute to US west coast precipitation amounts. To start, most of California's rainfall arrives during the cool months (October-March), when low level jets draw heat and moisture from the tropics, near Hawaii, northeastward to the midlatitudes [3]. This strip of heat and moisture to the North American West Coast is commonly known as the Pineapple Express, or an Atmospheric River (AR) [4]. ARs are usually characterized by a length of around 2000 km and a width up to 1000 km. The air column typically has at least 2 cm of integrated water vapor (IWV), or greater than 250 kg m −1 s −1 of integrated water vapor transport (IVT), within the first 2.5 km from the surface, and a center of strong winds at 2 km [5,6].
Despite the discrepancies in dust aerosol effects in various studies, they all agree that aerosol transport is the strongest in the spring, when Asian dust and pollution emissions are the highest [8,10,[16][17][18]. However, little to no literature has described the extent and effects of Asian aerosol transport on the US West Coast during the cool months (October-March) in terms of atmospheric processes (precipitation, development, etc.). That being said, Ault et al. (2011) [16] detected mineral dust from Asia in orographic precipitation over California during the CalWater Early Start campaign, which examined two extratropical cyclone systems between 22 February and 11 March 2009. The enhancement of horizontal IVT and IWV concentrations that resulted from the two systems were attributed to AR conditions. Still, the question of how the long-range transport of Asian aerosols affect precipitation amounts, system location, and AR frequency, particularly for the 2017 events, remains. The study will compare the properties of six 2017 AR events to a ten year (2006)(2007)(2008)(2009)(2010)(2011)(2012)(2013)(2014)(2015)(2016) climatology, in order to (1) understand the effect of Asian long-range aerosol transport on different locations over Northern California during ARs, (2) examine how the long-range transport of aerosols affect precipitation intensity, location, and the frequency of ARs over Northern California, and (3) assess how aerosols are correlated with precipitation associated with AR events. It is hypothesized that there is a positive correlation between aerosols and precipitation, in that aerosols will have contributed to the precipitation enhancement of the 2017 AR events.

AR Events and Synoptic Conditions in 2017
For the study, the months of January to April were taken into consideration, with a focus on the months of January and February, as the rainfall was particularly heavy during this time. The highest amounts of observed rainfall occurred at the coasts and mountains, as seen in Figure 1a [19], with amounts exceeding 20 inches (508 mm). Inland measurements were up to 6-10" (152-254 mm). Rainfall amounts were higher overall in northern California, and decreased towards the south as well as further inland on the other side of the Sierras. This trend is also observed during the March and April months. In comparison to climatological rainfall norms collected by NOAA since 2005, the majority of California experienced a 110% to 600+% increase during January and February 2017 (Figure 1b), with the highest percentage changes in the mountains. For March and April 2017, rainfall was at or slightly below normal amounts for central California, while southern California amounts were below 50%. Most of the damages, as a result of flooding and erosion, occurred during January and February. Overall, California experienced at least 150% of its normal rainfall amounts for 2017. For this reason, case studies were taken from these two months.

January 2017
The first series of AR events occurred from 01/07-11, bringing as much as 20" (508 mm) of rainfall to areas such as Strawberry Valley, located at the western base of the Sierras. On 01/08, the Uvas Creek in Santa Clara County flooded roads and properties up to three feet (914.40 mm), due to

January 2017
The first series of AR events occurred from 01/07-11, bringing as much as 20" (508 mm) of rainfall to areas such as Strawberry Valley, located at the western base of the Sierras. On 01/08, the Uvas Creek in Santa Clara County flooded roads and properties up to three feet (914.40 mm), due to the Uvas Reservoir overflowing and spilling into the creek. The reservoir typically receives a monthly average of 4.88" (124.95 mm) of rainfall during January, but experienced a max 48 h rainfall total of 7.05" (179.07 mm) in 2017. Other incidents that occurred during this event were road flooding in Morgan Hill, and 48 h rain gauge amounts of 8.74" (222 mm) at the Guadalupe Watershed, 9.49" (241.05 mm) at Uvas Canyon County Park, and 6.73" (170.94 mm) at the Lexington Reservoir. January rainfall amounts at these areas are generally 3.47", 3.63", and 3.97", (88.14 mm, 92.02 mm, 100.84 mm) respectively.
The second series of AR events occurred between 01/17 and 01/19 where southern California received up to 1.5" (38 mm) of rainfall, central California received up to 4" (101 mm), and northern California received up to 6" (152 mm) [20,21].

February 2017
The third series of AR events occurred between 02/06 and 02/10. Around 6.5" (165.1 mm) were recorded for areas at higher elevations, as well as reports of river flooding in parts of Sonoma County and the San Francisco Bay area, such as the San Francisquito Creek flooding nearby businesses as well as the US highway 101 bridge. In addition, the Oroville Dam flooding incident [22] that led to the evacuation of over 188,000 people mid-February was partly caused by the heavy rainfall which the lake experienced on 02/07. The fourth series of AR events occurred between 02/15 and 02/21. From the 02/15-17, two storm systems saturated the soil, which easily led to runoff and reservoir overflow during following events, as seen with the Anderson Reservoir spillage on 02/18. On 02/21, the combination of Coyote Creek and Anderson Reservoir overflow caused numerous flooding events in San Jose, inundating several neighborhoods, businesses, and roads [20].
In terms of post flood damages, Santa Clara Valley Water District's 2017 Flooding Report estimated that San Jose identified about $50 million in private property damages, and $23 million in public property damages. Around 14,000 people were evacuated during these events, with no reports of injury or death during rescue operations [20].

Synoptic Conditions and the Observed Characteristics of ARs
The synoptic conditions for the four different AR events are examined to see if a specific synoptic condition can provide favorable conditions for enhancing the long-range aerosol transport, and how the Potential Vorticity associated with the lower-level wind, water vapor, and sea level pressure evolves with time during the AR events. Figure 2 shows the map of ERA-Interim Potential Vorticity (PV) at 250 hPa, specific humidity (q) overlaid by the horizontal wind at 850 hPa, temperature at 850 hPa, and sea level pressure for AR events during January-February 2017. Strong upper-level troughs with large PV were located offshore near the coast of California at this time, which implied that stratospheric air through the tropopause (PVU = 1.5 PVU is regarded as a dynamic tropopause) penetrated downward to the mid-troposphere (~500 hPa, not shown). PV intrusion mostly occurred anticyclonically (i.e., clockwise rotation) over the Pacific and the western U.S., extending toward the West Coast of the U.S (Figure 2a). On 01/19/17, however, PV penetrated the western U.S. cyclonically, mimicking cyclonic wave breakings [23], and leading to stronger sea level low pressure over the northwest of California (Figure 2a,d, second row). The filamentary, northward directing water vapor bands over the Pacific reflected the response of the lower-level water vapor field to the existence of an upper level trough in the developing mid-latitude baroclinic system during AR events (Figure 2b). For all AR events, the high and narrow water vapor band was elongated from southwest to northeast with strong southwesterlies extending toward the California coast for all AR events. As shown in Figure 2a,c, a wavy structure of temperature and its horizontal gradient also appears to be associated with the narrow and deep penetration of elevated upstream PV. Temperature was slightly cooler inland, especially in the 01/07 AR case, but this was also dependent on the location, and extent of the sea level's low and high. There were spatial differences in high and low in sea level pressure for different ARs, but di-pole patterns of the low and high pressures persisted during all AR events (Figure 2d). These strong contrasts between the sea level high and low was shown both in the Pacific and on the West Coast of the U.S., confirming that these are prominent synoptic features of AR events. Low sea level pressure also developed over the eastern North Pacific, associated with the upper level trough for all AR events.
( Figure 2a). On 01/19/17, however, PV penetrated the western U.S. cyclonically, mimicking cyclonic wave breakings [23], and leading to stronger sea level low pressure over the northwest of California (Figure 2a, d, second row). The filamentary, northward directing water vapor bands over the Pacific reflected the response of the lower-level water vapor field to the existence of an upper level trough in the developing mid-latitude baroclinic system during AR events (Figure 2b). For all AR events, the high and narrow water vapor band was elongated from southwest to northeast with strong southwesterlies extending toward the California coast for all AR events. As shown in Figure 2a and 2c, a wavy structure of temperature and its horizontal gradient also appears to be associated with the narrow and deep penetration of elevated upstream PV. Temperature was slightly cooler inland, especially in the 01/07 AR case, but this was also dependent on the location, and extent of the sea level's low and high. There were spatial differences in high and low in sea level pressure for different ARs, but di-pole patterns of the low and high pressures persisted during all AR events ( Figure 2d). These strong contrasts between the sea level high and low was shown both in the Pacific and on the West Coast of the U.S., confirming that these are prominent synoptic features of AR events. Low sea level pressure also developed over the eastern North Pacific, associated with the upper level trough for all AR events. By 2-3 days after each AR event, the higher amount of water vapor has moved inland with cooler temperature as the low pressure moved closer to the western U.S., and the upper level trough was weaker and decayed (not shown). Overall, as the AR approaches the western U.S., a strong upper-level trough with large PV penetrating either anticyclonically or cyclonically is located offshore near the coast of California, leading to strong low-level southerlies and southwesterlies over this same coast of California. This also suggests that the large-scale upper-level trough associated with upper-level jet extending toward northern California may provide a favorable condition for the long-range aerosol transport from Asia during AR events.

Satellite Data
Moderate Resolution Imaging Spectroradiometer (MODIS) Aqua was used to note whether there was Asian long-range transport occurring during the 2017 cases chosen. Aerosol optical depth was measured daily at 500 nm with 1 • resolution. MODIS data was visualized one week before the atmospheric river (AR) event, as this is the usual time it takes for Asian dust to travel across the Northern Pacific Ocean. As shown in Figure 3, there is a visible transport of aerosol across the Northern Pacific. Cloud-cover contributes to most of the missing data, but in some cases MODIS-Aqua was able to pick up Aerosol Optical Depth (AOD) signals along the edges of cloud cover. It is acknowledged that there are limitations to using satellite data because of cloud cover. Alternate methods would be the use of Lidar to visualize from the surface up, however this is not currently possible over the ocean. In the case of CALIPSO data, the swath rarely went over the areas of interest. Another alternative would be to use reanalysis or model data, however this study's scope was to first use observational data.    [24]) provided by NOAA/NESDIS at 0.25 • grid spacing and average twice daily (e.g., 00-12Z and 12-24Z) was employed to visualize the AR event progression off the coast of the Western United States from January and February ( Figure 4). Some selected examples of four case studies (01/07-10, 01/17-19, 02/06-10, 02/15-21) are shown in Figure 4. AR events can be identified when water vapor amounts to above at least 2 cm of integrated water vapor (IWV). Examples of AR events that occurred during these time periods are shown. Water vapor amounts above at least 4 g/cm 2 (>2 cm [7]) represents an AR. Examples of AR events that occurred during these time periods are shown. Water vapor amounts above at least 4 g/cm 2 (>2 cm [7]) represents an AR.

Observation Site Data
The Interagency Monitoring of Protected Visual Environments (IMPROVE) data was used for aerosols measurements in California [10,18,25,26]. Data measurements were taken for twenty-four hours every three days. For the study, a combination of sites were used to represent differences between coastal (Redwood National Park, Point Reyes National Shores, Pinnacles National Park), inland (Trinity, Fresno), and mountainous (Lassen Volcanic National Park, Bliss State Park, Kaiser) areas ( Figure 5). The site Trinity was dropped from the final analysis since measurements from there ended in 2015. Therefore, Fresno was the only site to represent inland conditions. The Interagency Monitoring of Protected Visual Environments (IMPROVE) data was used for aerosols measurements in California [10,18,25,26]. Data measurements were taken for twenty-four hours every three days. For the study, a combination of sites were used to represent differences between coastal (Redwood National Park, Point Reyes National Shores, Pinnacles National Park), inland (Trinity, Fresno), and mountainous (Lassen Volcanic National Park, Bliss State Park, Kaiser) areas ( Figure 5). The site Trinity was dropped from the final analysis since measurements from there ended in 2015. Therefore, Fresno was the only site to represent inland conditions. Iron, calcium, and soil measurements were used to determine Asian dust amounts at the sites. Using a method employed in Creamean et al. (2014), the dust source was separated into Asian or localized, using a ratio of iron (Fe) to calcium (Ca), which is generally accepted as the composition of Asian desert dust [26,27]. Ratios less than a value of one indicated Asian dust presence, hence long-range transport, while ratios greater than one indicated localized dust. Based on Fe/Ca ratios in Equation (1), soil amounts for each site were categorized as either Asian dust or Non-Asian dust for simplicity.

Precipitation Data
Integrated water vapor transport measurements were taken from Modern-Era Retrospective analysis for Research and Applications, version 2.1 (MERRA 2), with a horizontal resolution of 0.5° by 0.66°. 3-hourly reanalysis data was provided by the AR catalog created by Rutz et al. (2014) [7]. Due to the fact that this was a gridded dataset, the integrated water vapor transport (IVT) measurements were extracted to be the closest to the IMPROVE site locations.
To match with IMPROVE dataset dates, daily IVT averages were calculated from the 3-hourly MERRA dataset. IVT measurements were used calculate AR frequency and strength as later discussed in Section 4.2. The Climate Prediction Center Unified Gauge-Based Analysis of Daily Precipitation over CONUS Real-time was used for daily precipitation measurements. The dataset Iron, calcium, and soil measurements were used to determine Asian dust amounts at the sites. Using a method employed in Creamean et al. (2014), the dust source was separated into Asian or localized, using a ratio of iron (Fe) to calcium (Ca), which is generally accepted as the composition of Asian desert dust [26,27]. Ratios less than a value of one indicated Asian dust presence, hence long-range transport, while ratios greater than one indicated localized dust. Based on Fe/Ca ratios in Equation (1), soil amounts for each site were categorized as either Asian dust or Non-Asian dust for simplicity.
Other aerosols analyzed from the sites were black carbon (BC) and sulfate (SO 4 ).

Precipitation Data
Integrated water vapor transport measurements were taken from Modern-Era Retrospective analysis for Research and Applications, version 2.1 (MERRA 2), with a horizontal resolution of 0.5 • by 0.66 • . 3-hourly reanalysis data was provided by the AR catalog created by Rutz et al. (2014) [7]. Due to the fact that this was a gridded dataset, the integrated water vapor transport (IVT) measurements were extracted to be the closest to the IMPROVE site locations.
To match with IMPROVE dataset dates, daily IVT averages were calculated from the 3-hourly MERRA dataset. IVT measurements were used calculate AR frequency and strength as later discussed in Section 4.2. The Climate Prediction Center Unified Gauge-Based Analysis of Daily Precipitation over CONUS Real-time was used for daily precipitation measurements. The dataset has a resolution of 0.25 • × 0.25 • , and was manipulated to match IMPROVE's every 3-day temporal resolution.

Back Trajectory Model Simulation
The NOAA Hybrid Single-Particle Lagrangian Integrated Trajectory model (HYSPLIT) model was used to analyze parcel back trajectories [26]. The Global Data Assimilation System (GDAS) 0.5 degree archive ensembles at Redwood NP (coast), Fresno (inland), and Lassen Volcanic NP (mountain) were simulated at 3, 5, and 7 km for 10 days. These three heights were chosen to observe parcel lifting during the track. During model set-up, heights were chosen in meters, but are displayed in hPa for the trajectory plots. This was to provide an easier way of following parcel heights along the track as well as comparison with synoptic plots.

Dust Measurements
The total dust amounts from 01/01 to 04/28 from the seven sites are investigated. As shown in Figure 6, there is an increase in instances of localized dust presence as the system moves southward among the coastal sites. In terms of the inland site, Fresno dust measurements are mainly attributed to localized dust, which is most likely the result of its surrounding city and/or agricultural land. At the mountainous sites, elevation differences might be a reason for why the mountainous sites did not follow the North-South decrease in Asian dust trends. For example, site measurements could be influenced by whether the site was located toward or away from the windward direction, or influenced by localized terrain differences. Overall, these results are similar to findings in Rutz et al. (2014) [7]. has a resolution of 0.25° × 0.25°, and was manipulated to match IMPROVE's every 3-day temporal resolution.

Back Trajectory Model Simulation
The NOAA Hybrid Single-Particle Lagrangian Integrated Trajectory model (HYSPLIT) model was used to analyze parcel back trajectories [26]. The Global Data Assimilation System (GDAS) 0.5 degree archive ensembles at Redwood NP (coast), Fresno (inland), and Lassen Volcanic NP (mountain) were simulated at 3, 5, and 7 km for 10 days. These three heights were chosen to observe parcel lifting during the track. During model set-up, heights were chosen in meters, but are displayed in hPa for the trajectory plots. This was to provide an easier way of following parcel heights along the track as well as comparison with synoptic plots.

Dust Measurements
The total dust amounts from 01/01 to 04/28 from the seven sites are investigated. As shown in Figure 6, there is an increase in instances of localized dust presence as the system moves southward among the coastal sites. In terms of the inland site, Fresno dust measurements are mainly attributed to localized dust, which is most likely the result of its surrounding city and/or agricultural land. At the mountainous sites, elevation differences might be a reason for why the mountainous sites did not follow the North-South decrease in Asian dust trends. For example, site measurements could be influenced by whether the site was located toward or away from the windward direction, or influenced by localized terrain differences. Overall, these results are similar to findings in Rutz et al., (2014) [7].

2006-2016 Climatological Analysis in Comparison to 2017
The evolution of aerosol for January-April from 2006 to 2016 were constructed. Redwood NP, Fresno, and Lassen Volcanic NP are shown to represent coasts, inlands, and mountains, respectively, Figure 6. 2017 Total dust amounts. The most northern and southern site for the coastal and mountainous regions were selected to show aerosol patterns as the AR system moves along the US West Coast. At each site, dust amounts are totaled from January to April, then separated into Asian-influenced (blue) or non-Asian/localized dust influenced (red).

2006-2016 Climatological Analysis in Comparison to 2017
The evolution of aerosol for January-April from 2006 to 2016 were constructed. Redwood NP, Fresno, and Lassen Volcanic NP are shown to represent coasts, inlands, and mountains, respectively, (Figure 7). It is seen that the majority of dust amounts at coastal sites are attributed to Asia, with some amounts being due to localized dust. Inland at Fresno, localized dust is in the majority with hardly any trace of Asian influenced dust. Again, this is not surprising considering Fresno is located in the valley, where it is expected that localized dust is trapped compared to coastal or mountainous sites. Lastly, the majority of Asian and locally influenced dust varies from year to year at the mountainous regions. This is more so in the case of Lassen Volcanic NP, due the fact that the observation site is located at a lower elevation on the windward side of the mountain.
AR frequency and strength comparisons between the 2006 and 2016 climatology and year 2017 were conducted (Table 1)     In addition to dust, black carbon (BC) and sulfate (SO 4 ) were also measured at each site. A correlation analysis was conducted between these three aerosol types and precipitation amounts.  Table 2 summarizes the correlation coefficient (R-values) for all sites. Overall, R-values vary from being low climatologically, and in the 0.30-40 s for 2017 data. Correlation values with low p-values (α = 0.10) of 0.10 or less are statistically significant and are highlighted. Most of the correlations for both time periods are negative, meaning that aerosols were found to have a suppressing effect on precipitation. This trend was also found when calculating correlations between aerosols and IVT (not shown). Starting with climatological data, statistically significant correlation values appear to be Asian dust and SO 4 for Redwood NP (−0.15 and −0.12). While these values in and of themselves are relatively small, these were the highest at this site that also had p-values < α. BC contributes to the highest statistically significant negative correlations in five out of seven sites at the two most southerly coastal sites, and at each of the mountainous sites (−0.11, −0.18, −0.12, −0.11, and −0.14). SO 4 amounts are also slightly higher at the southerly mountainous sites and Pinnacles NM (−0.09, −0.12, and −0.08). In terms of 2017 data, Asian dust and SO 4 are the most likely contributors to precipitation suppression at Redwood NP with values of −0.44 and −0.43, respectively. At Point Reyes NS, BC and SO 4 have the highest negative correlations (−0.30 and −0.27). Pinnacles NM, Fresno, Lassen Volcanic, and Bliss all had higher negative correlations of locally influenced dust, BC, and SO 4 ( Figure 8). This makes sense for the first three sites as they are more inland, or facing inland (Pinnacles, Fresno, Lassen Volcanic) and located near the southern-most half of northern California.   Figure 9a shows the ten-day back trajectories for Redwood NP, Fresno, and Lassen Volcanic NP starting from 01/07 at 23Z. It appears that from each location, air parcels originate from East Asia, in the middle of China, which is most visibly shown from the Lassen Volcanic NP ensemble. Figure 9b shows parcel heights during the back trajectory. Parcels start near the surface, below 900 hPa, for the first few days before being vertically advected to around 450-300 hPa between 01/02 and 01/03. After two to three days, parcels lower to around 700 hPa when they reach the sites.  Figure 9a shows the ten-day back trajectories for Redwood NP, Fresno, and Lassen Volcanic NP starting from 01/07 at 23Z. It appears that from each location, air parcels originate from East Asia, in the middle of China, which is most visibly shown from the Lassen Volcanic NP ensemble. Figure 9b shows parcel heights during the back trajectory. Parcels start near the surface, below 900 hPa, for the first few days before being vertically advected to around 450-300 hPa between 01/02 and 01/03. After two to three days, parcels lower to around 700 hPa when they reach the sites.

HYSPLIT Back Trajectory Analysis
The 10-day back trajectories and parcel height comparisons at Redwood NP for the ensuing 2017 AR events on 01/07, 01/19, 02/09, and 02/21 are shown in Figure 10. The trajectories of the 3 AR events, as shown in Figure 10a, align with those in Figure 9a suggesting that their origins are to be located between the coast of northern China and Japan. However, it is interesting to note the differences between tracks for each event. For instance, the ensemble trajectories for 01/07 generally show patterns of shorter and steeper ridges and troughs as the parcel is advected over the north Pacific. 01/19 and 02/09, on the other hand, generally display longer and more zonal ridge/trough movement. In Figure 10b, all four events (01/07, 01/19, 02/09, 02/21) show ensemble trajectories where there is steep lifting from the surface to around 450 hPa, with 01/07 and 02/21 being most noticeable. Events with longer/zonal movement generally resulted in lifting from the surface to around 750 hPa, as seen on 01/19 and 02/09, with some air parcels being an exception. The 10-day back trajectories and parcel height comparisons at Redwood NP for the ensuing 2017 AR events on 01/07, 01/19, 02/09, and 02/21 are shown in Figure 10. The trajectories of the 3 AR events, as shown in Figure 10a, align with those in Figure 9a suggesting that their origins are to be located between the coast of northern China and Japan. However, it is interesting to note the differences between tracks for each event. For instance, the ensemble trajectories for 01/07 generally show patterns of shorter and steeper ridges and troughs as the parcel is advected over the north Pacific. 01/19 and 02/09, on the other hand, generally display longer and more zonal ridge/trough movement. In Figure 10b, all four events (01/07, 01/19, 02/09, 02/21) show ensemble trajectories where there is steep lifting from the surface to around 450 hPa, with 01/07 and 02/21 being most noticeable. Events with longer/zonal movement generally resulted in lifting from the surface to around 750 hPa, as seen on 01/19 and 02/09, with some air parcels being an exception. This difference could be a possible explanation for the differences in parcel vertical advection as seen in Figure 10b. In general, the majority of parcels are lifted from the surface to higher levels before settling back to the surface. These track patterns follow and agree with the synoptic conditions described in Figure 2. This difference could be a possible explanation for the differences in parcel vertical advection as seen in Figure 10b. In general, the majority of parcels are lifted from the surface to higher levels before settling back to the surface. These track patterns follow and agree with the synoptic conditions described in Figure 2.

Discussion and Conclusions
Comparing the properties of a 2006-2016 long-term analyses to selected 2017 AR events, the study aimed to understand three main points: (1) The effects of long-range aerosol transport on different terrain regions over Northern California, (2) how aerosols affect precipitation intensity, location, and frequency using IVT as a proxy, and (3) assess how aerosols are correlated with ARassociated precipitation. It is shown that more Asian dust tends to be observed over the coast while more locally produced dust tends to be observed inland. The hypothesis was that aerosols would have enhanced the precipitation from the 2017 AR events (positive correlations). However, based on analysis, dust appears to have a slight suppression of precipitation, particularly farther south in the mountains. Black carbon (BC), on the other hand, shows a slight enhancement of precipitation for sites at higher latitudes than Bliss SP and Kaiser, and closer to the coast, as seen at Pinnacles NP. Sulfate (SO4) also tended to have slight precipitation suppression. In this study our correlation analysis suggests that the long-range transport of aerosols from East Asia had little evidence of effect on either AR events climatologically or in 2017. Reasons for low R-values could be due to wet deposition and scavenging as aerosols travel across the North Pacific. Locally influenced dust, on the other hand, have statistically significant negative correlations with the precipitation during AR events for the majority of the sites. BC and SO4 tend to be associated with precipitation suppression, but whether they are transported or locally generated is not clear and this needs to be investigated further in the future study.
In terms of future work, AR events occurring in March-April, where more dust is present, will be compared to AR events occurring in January-February. Also, in part 2 of the study, numerical

Discussion and Conclusions
Comparing the properties of a 2006-2016 long-term analyses to selected 2017 AR events, the study aimed to understand three main points: (1) The effects of long-range aerosol transport on different terrain regions over Northern California, (2) how aerosols affect precipitation intensity, location, and frequency using IVT as a proxy, and (3) assess how aerosols are correlated with AR-associated precipitation. It is shown that more Asian dust tends to be observed over the coast while more locally produced dust tends to be observed inland. The hypothesis was that aerosols would have enhanced the precipitation from the 2017 AR events (positive correlations). However, based on analysis, dust appears to have a slight suppression of precipitation, particularly farther south in the mountains. Black carbon (BC), on the other hand, shows a slight enhancement of precipitation for sites at higher latitudes than Bliss SP and Kaiser, and closer to the coast, as seen at Pinnacles NP. Sulfate (SO 4 ) also tended to have slight precipitation suppression. In this study our correlation analysis suggests that the long-range transport of aerosols from East Asia had little evidence of effect on either AR events climatologically or in 2017. Reasons for low R-values could be due to wet deposition and scavenging as aerosols travel across the North Pacific. Locally influenced dust, on the other hand, have statistically significant negative correlations with the precipitation during AR events for the majority of the sites. BC and SO 4 tend to be associated with precipitation suppression, but whether they are transported or locally generated is not clear and this needs to be investigated further in the future study.
In terms of future work, AR events occurring in March-April, where more dust is present, will be compared to AR events occurring in January-February. Also, in part 2 of the study, numerical simulations will also be utilized to compare to observation data collected and model aerosol transport effects on local precipitation and cloud microphysics.