Finescale Orographic Precipitation Variability and Gap-Filling Radar Potential in Little Cottonwood Canyon , Utah

Finescale variations in orographic precipitation pose a major challenge for weather prediction, winter road maintenance, and avalanche forecasting andmitigation in mountainous regions. In this investigation, groundbased X-band radar observations collected during intensive observing period 6 (IOP6) of the Storm Chasing Utah Style Study (SCHUSS) are used to provide an example of these variations during a winter storm in the Wasatch Mountains of northern Utah. Emphasis is placed on precipitation features in and around Little Cottonwood Canyon (LCC), which cuts orthogonally eastward into the central Wasatch Mountains. Precipitation during the weakly stratified prefrontal storm stage featured a wavelike barrier-scale reflectivity maximum over theWasatch Crest and upper LCC that extendedweakly westward along the transverse ridges flanking LCC. This precipitation pattern appeared to reflect a veering wind profile, with southwesterly flow over the transverse ridges but cross-barrier westerly flow farther aloft. Sublimation within dry subcloud air further diminished low-level radar reflectivities over lower LCC. In contrast, the cold-frontal stage was associated with stronger reflectivities over lower LCC and the adjoining northto northwest-facing canyon wall, consistent with shallow, northwesterly upslope flow. These results highlight the finescale precipitation variations that can occur during winter storms in complex terrain and demonstrate the potential for improved analysis and forecasting of precipitation in LCC using a gap-filling radar.


Introduction
Finescale (1-10 km) variations in precipitation over areas of highly variable three-dimensional topography pose a significant challenge for weather and climate prediction, hydrologic forecasting, avalanche mitigation, and winter road maintenance in mountainous regions (e.g., Meyers and Steenburgh 2013;Moreno et al. 2013).Interactions between the large-scale flow, regional topography, and local terrain features can create dramatic variations in precipitation rate and structure over small spatial scales during individual storms (Sinclair et al. 1997;Steenburgh 2003;Smith et al. 2003;Minder et al. 2008;Molinié et al. 2012;Moreno et al. 2013), as well as over climatological time periods (e.g., Frei and Schär 1998;Anders et al. 2007).Although much of the research examining orographic precipitation has concentrated on barrier-scale effects, especially windward enhancement processes, the influence of smaller-scale topographic features (e.g., ridge-valley corrugations) and more complex terrain geometries has received less attention (Minder et al. 2008).
The dynamical and microphysical processes that affect the distribution and intensity of orographic precipitation vary with the dynamics and thermodynamics of the incipient airflow, the size and shape of the terrain, and the time scales controlling the growth and fallout of precipitation particles [see Roe (2005), Smith (2006), Houze (2012), Colle et al. (2013), and Stoelinga et al. (2013) for recent reviews].Over the western United States, many orographic storms evolve through stable, transitional (frequently with frontal characteristics), and unstable stages (e.g., Hobbs 1975;Marwitz 1980;Cooper and Saunders 1980;Long et al. 1990;Sassen et al. 1990;Medina et al. 2007).The stable and unstable stages are roughly analogous to the ''flow around'' and ''flow over'' regimes emphasized by Rotunno and Houze (2007).
During the stable stage, the low-level flow is often blocked, resulting in an upstream pool of stagnant air or a barrier jet.Increased low-level stability due to sublimational or evaporative cooling can also contribute to flow blocking, as has been observed upstream of the Wasatch Mountains of northern Utah (Colle et al. 2005;Cox et al. 2005), and create downslope or down-valley flow within canyons and valleys (Hill 1978;Bousquet and Smull 2003;Steiner et al. 2003).The ascent of incident air over the blocked flow can enhance precipitation as much as 100-150 km upstream of the initial mountain slope (Peterson et al. 1991;Houze et al. 2001).In some cases, strong vertical wind shear and turbulence near the top of the blocked air mass may contribute to precipitation generation (Medina et al. 2005;Houze and Medina 2005).
Precipitation processes during the stable stage can also be influenced by barrier-and subbarrier-scale mountain waves, the latter produced by ridge-valley corrugations in the topography.Klimowski et al. (1998), Bruintjes et al. (1994), and Reinking et al. (2000) found that smallscale ridges along Arizona's Mogollon Rim generate gravity waves that modulate the distribution of cloud liquid water and precipitation during winter storms.Colle (2004) showed that the reduced upstream tilt of gravity waves over narrow barriers (;25 km half-width) results in precipitation fallout over the crest and spillover to the lee.Colle (2008) used two-dimensional simulations to show that a series of windward ridges can enhance precipitation over each ridge crest by a factor of 2-3 when the ridge spacing is relatively small (;20 km) and there is strong cross-barrier flow.Over the Cascade Mountains during the second Improvement of Microphysical Parameterization through Observational Verification Experiment (IMPROVE-2), a synergistic effect was found between barrier-scale mountain waves produced by midlevel cross-barrier flow and small-scale mountain waves formed by along-barrier flow over ridge-valley corrugations.This led to greater cloud liquid water production, hydrometeor mixing ratios, and precipitation rates over the subbarrier-scale windward ridges (Garvert et al. 2005(Garvert et al. , 2007)).
The transitional stage can feature a sharp cold-frontal passage (Long et al. 1990;Sassen et al. 1990), a gradual transition in storm structure (Marwitz 1980), or a surge of low equivalent potential temperature u e air aloft that destabilizes the prefrontal environment, leading to convection (Reynolds and Kuciauskas 1988;Steenburgh 2003).The release of potential instability generated by surges of low-u e air aloft or differential u e advection ahead of a surface cold front frequently contributes to orographic precipitation enhancement, especially when coupled with the seeder-feeder process (Browning et al. 1974).Enhanced water vapor fluxes accompanying the prefrontal low-level jet can also contribute to strong orographic precipitation enhancement (Lin et al. 2001;Neiman et al. 2002;White et al. 2003;Neiman et al. 2008).Narrow cold-frontal precipitation bands can exhibit changes in intensity, shape, and orientation as they interact with topographic obstacles (Braun et al. 1997;Colle et al. 1999;Yu and Smull 2000;Colle et al. 2002;Neiman et al. 2004;Bond et al. 2005;Viale et al. 2013).
The unstable storm stage is frequently characterized by orographically initiated or enhanced convection (e.g., Hobbs 1975).Reynolds and Kuciauskas (1988) and Peterson et al. (1991) found that, following cold-frontal passage in the Sierra Nevada, dry air tends to move in aloft with the resulting potential instability released as shallow convection when the flow surmounts the mountain barrier.Over the Oregon Cascades, Medina et al. (2007) showed that the ''late sector'' of extratropical cyclones generally consists of shallow convective echoes, which broaden as they move over the mountain barrier.This was also observed during IMPROVE-2 when postfrontal convective cells strengthened and upstream precipitation became more persistent during flow over the Cascade Mountains (Woods et al. 2005).Although not occurring exclusively in postfrontal or cold-sector environments, a variety of banded, cellular, and mesoscale precipitation features are generated or modified during unstable or potentially unstable flow interaction with complex terrain (e.g., Maddox et al. 1978;Caraceña et al. 1979;Parsons and Hobbs 1983;Sénési et al. 1996;Kirshbaum and Durran 2005;Kirshbaum et al. 2007;Godart et al. 2011).
The rich spectrum of storm environments, incident flow characteristics, topographic shapes and geometries, and dynamical and microphysical processes described above complicates precipitation forecasting and frequently leads to large intra-and interstorm variations in precipitation intensity and distribution (e.g., Steenburgh 2003).Often these variations occur at meso-g (2-20 km) or even micro-a (200 m-2 km) scales (e.g., Rangno 1986;Neiman et al. 2002;Minder et al. 2008;Reuder et al. 2007) and are inadequately observed by operational radars or precipitation-gauge networks (e.g., Westrick et al. 1999;Wood et al. 2003;Beck and Bousquet 2013).Although high-resolution modeling systems are increasingly able to capture climatological precipitation gradients in mountainous regions (e.g., Ikeda et al. 2010), they frequently produce large errors during individual storms (e.g., Minder et al. 2008).
In this paper, we use observations collected by a mobile X-band radar to illustrate the finescale nature of precipitation features produced during a winter storm in an area of highly variable three-dimensional topography: the central Wasatch Mountains of northern Utah (Fig. 1a).The central Wasatch Mountains contain a system of ridges (Alpine, Cottonwood, and Wildcat) and canyons [Little Cottonwood (LCC), Big Cottonwood, and Mill Creek] that are oriented from west to east, orthogonal to the predominantly north-south-oriented Wasatch Crest (Fig. 2a).Although the Wasatch Crest forms a hydrologic divide near the eastern termini of these canyons, the highest terrain is found along the Alpine and Cottonwood Ridges, which flank the deep, glacially carved LCC.Lone Peak (3430 m MSL; all subsequent elevations are MSL unless otherwise specified) and Twin Peaks (3453 m) create the westernmost abutments of the Alpine and Cottonwood Ridges, respectively, and dominate the Salt Lake Valley (SLV) skyline (Fig. 2b).
Dramatic climatological gradients in precipitation are found over the central Wasatch Mountains (Fig. 1b;Dunn 1983).A National Weather Service Cooperative Observer Program (COOP) observing site at Alta (2661 m) in upper LCC averages ;96 cm of liquid equivalent precipitation during the cool season (November-April), whereas Salt Lake City International Airport (KSLC, 1286 m; see Fig. 1a for location) in the northern SLV averages only ;23 cm (Western Regional Climate Center 2013).Estimates produced by the Parameter-Elevation Regression on Independent Slopes Model (PRISM) Climate Group at Oregon State University (Daly et al. 1994) suggest an increase in cool-season liquid equivalent precipitation from ;36 cm at the base of LCC to ;105 cm along the Alpine Ridge (Fig. 1b).
The winter storm examined in this paper was sampled during intensive observing period 6 (IOP6; from 0900 UTC 12 November to 0400 UTC 13 November 2011) of the Storm Chasing Utah Style Study (SCHUSS), a radar education and outreach field campaign involving the University of Utah and the Center for Severe Weather Research (CSWR).During IOP6, University of Utah students operated the CSWR Doppler on Wheels 6 (DOW6) in the southwest corner of the SLV (see Fig. 2a for location) where a mostly unobstructed view of LCC and the western face of the central Wasatch Mountains allowed for detailed observation of precipitation features shaped by multiscale mountain waves and frontmountain interactions.State Route 210 (SR-210), which runs the length of LCC, services two of Utah's most popular ski resorts and has one of the highest snow Avalanche Hazard Indices of any major road in North America (Nepstad et al. 2006).Given that small-scale spatial and temporal gradients in precipitation intensity and hydrometeor type directly affect snow avalanche conditions, this paper also illustrates the potential for using a gapfilling dual-polarized X-band radar to improve the analysis and prediction of precipitation and avalanche hazard within LCC.

a. Doppler on Wheels
High-resolution radar scans from DOW6, an X-band (3.2-cm wavelength) dual-polarized Doppler weather radar developed by the CSWR (Wurman et al. 1997) and operated by students from the University of Utah, provide the cornerstone for this analysis.During IOP6, DOW6 operated near the southwest corner of the SLV at an elevation of 1508 m-about 200 m above the lowest point in the SLV (Fig. 2a).The continuous scanning strategy included plan position indicator (PPI) scans at approximately 18 intervals between 0.58 and 13.78 and range-height indicator (RHI) scans over LCC, the surrounding central Wasatch Mountains, and the eastern SLV.Although DOW6 operated from 1405 UTC 12 November to 0402 UTC 13 November 2011, precipitation after 2130 UTC 13 November 2011 was light in LCC and is not discussed here.
Software programs Solo II and Xltrs II (Oye et al. 1995) were used initially to visualize, process, and translate the DOW6 data.Processing and quality control included the reorientation of PPI and RHI scans relative to true north, the application of a normalized coherent power (a measure of the ''coherence'' of the data in which noisy range gates exhibit values near zero) minimum of 0.25, and the removal of range gates exhibiting both low-velocity and high-reflectivity values as ground clutter.No corrections were made for attenuation, as attenuation rates in snow are very low (Rinehart 1997).A low bias in radar reflectivity may exist down beam of topographic features that generate partial beam blockage (James et al. 2000); however, these range gates were not removed.The maximum unambiguous range was ;48 km, and velocity unfolding was performed at a Nyquist velocity of 19.6 m s 21 .Individual PPI and RHI scans were interpolated onto Cartesian grids with 75-m horizontal and 50-m vertical resolution, factoring in a 4 /3 Earth radius assumption and standard atmospheric refraction conditions (Rinehart 1997).
Time-mean reflectivity and velocity RHIs reflect the average reflectivity or velocity value of each range gate for the specified time period.Composite time-mean reflectivity PPIs are produced by first creating timemean reflectivity PPIs for each scan angle between 3.28 and 13.78, then by calculating the maximum value for each range gate between all scan angles.For both RHIs and PPIs, the averaging process eliminates range gates with undefined values or with fewer than five real values for the time period from the time-mean calculation.
The hydrometeor typing algorithm follows the work of Dolan and Rutledge (2009), which uses a fuzzy-logic approach involving horizontally polarized radar reflectivity, dual-polarization variables (e.g., differential reflectivity, specific differential phase, and correlation coefficient), and vertical temperature profiles (obtained from soundings launched from KSLC and the DOW6 site).Schneebeli et al. (2013) recently used this algorithm to examine the microphysical processes contributing to precipitation in the Swiss Alps.The Dolan and Rutledge (2009) algorithm uses a T-matrix scattering model with seven modeled hydrometeor types (vertically aligned ice, ice crystals, aggregates, low-density graupel, highdensity graupel, drizzle/light rain, and rain).Four of these hydrometeor types were identified by the algorithm during IOP6: 1) ice crystals (assumed to have diameters smaller than 1.5 mm, small axis ratios of 0.1-0.3,densities between 0.4 and 0.9 g cm 23 , and temperatures between 2408 and 2108C), 2) aggregates (assumed to be conglomerates of semispherical, randomly oriented ice crystals with diameters of 1-12 mm, axis ratios of 0.2-0.9,densities between 0.1 and 0.2 g cm 23 , and temperatures between 2158 and 58C), 3) drizzle/light rain (assumed to be small spherical droplets with diameters of 0.3-0.55mm and temperatures above 08C), and 4) low-density graupel (assumed to have axis ratios of 0.5-1.25,densities between 0.25 and 0.55 g cm 23 , and temperatures between 2208 and 58C).The low-density graupel temperature range is a modification from the range (from 2208 to 2108C) used by Dolan and Rutledge (2009) to reflect observations of partially to fully rimed ice crystals and small graupel at higher temperatures during wintertime storms in the central Wasatch Mountains.Specific differential phase, which is dependent on axis ratio and particle concentration and is used most commonly for the detection of heavy rain, was not available for much of IOP6; therefore, the algorithm used in this paper does not include this variable.
The hydrometeor-type RHIs reflect the most frequently identified hydrometeor type for each range gate across the specified time interval, and are presented along with selected images captured by a hydrometeor videosonde (HYVIS; Murakami and Matsuo 1990) snowflake camera system, located near 2990 m at Alta Ski Area in upper LCC, for comparison.The HYVIS system takes twodimensional photographs of snowflakes from the side as they fall through a cylindrical tube, cataloging them by time stamp.While the quality of these images is variable, the shape and size of the silhouettes allow us to subjectively classify the snow crystals.The HYVIS system was located in a sheltered grove of trees ;360 m north of the Up-Canyon South RHI scan, ;1000 m south of the Up-Canyon North RHI scan, and ;340 m north of the Alta-Collins (CLN) weather station (Fig. 2a, locations approximate).The hydrometeor-type algorithm was run using the Up-Canyon South RHI for all time periods except 1430-1700 UTC, when the Up-Canyon South RHI was not available and the Up-Canyon North RHI was used.

b. Supplemental meteorological data
Regional-scale precipitation features were examined using lowest-tilt (0.58) base reflectivity scans from the Salt Lake City, Utah (KMTX), National Weather Service Weather Surveillance Radar-1988 Doppler (WSR-88D), which were obtained in level II format (Crum et al. 1993) from the National Climatic Data Center (NCDC) Next Generation Weather Radar (NEXRAD) archive.KMTX is located on Promontory Point, a major peninsula of the Great Salt Lake, at 2004 m-about 500 m above the DOW6 site and about 700 m above the SLV floor (Fig. 1a).For clarity, KMTX data were smoothed using a moving-average filter with a span of 10 data points.
University of Utah students launched GRAW GPSbased upper-air soundings from the DOW6 site periodically during the event.Profiles of squared dry N 2 d and moist N 2 m Brunt-Väisälä frequencies were calculated from these soundings following Durran and Klemp (1982), and were smoothed using a moving-average filter with a span of 10 data points.Time-height cross sections and plan-view analyses from the Rapid Update Cycle, version 2 (RUC2; Benjamin et al. 1998), available at 13-km grid spacing, were used to assess synoptic conditions.Surface observations, including hourly liquid equivalent precipitation, were obtained from the MesoWest cooperative networks (Horel et al. 2002).

Event overview
During the event sampled by IOP6, a progressive upper-level short-wave trough moved across northern Utah, with the accompanying surface cold front passing KSLC at 1940 UTC 12 November and entering LCC at 2030 UTC 12 November.We divide IOP6 into three stages: prefrontal, frontal (rather than transitional given the existence of a well-defined surface cold front), and postfrontal.Due to the gradual progression of the cold front across the study area, we have defined a later onset of the frontal and postfrontal stages in LCC compared to KSLC.The lengthy prefrontal stage (0900-1940 UTC 12 November at KSLC and 0900-2030 UTC 12 November at LCC) featured a weakly stratified atmosphere with a dry (i.e., relative humidity ,80%) subcloud layer below ;750 hPa (Fig. 3).Over the SLV, winds veered with height from southerly or southeasterly near the surface (see later sounding and surface analysis figures) to westsouthwesterly near crest level (;700 hPa).With the passage of the shallow cold front (frontal stage; 1940-2100 UTC 12 November at KSLC and 2030-2130 UTC 12 November at LCC) the flow below crest level shifted to northwesterly, whereas westerly flow persisted above crest level.During the postfrontal stage (2100-0400 UTC 12-13 November at KSLC and 2130-0400 UTC 12-13 November at LCC), westerly to northwesterly flow prevailed and contained a transient mesoscale snowband that affected KSLC, but weakened before reaching LCC, where scattered convective snow showers were observed.The event concluded at ;0400 UTC 13 November as drier air moved into northern Utah.
The event produced 6 in.(15 cm) of snow and 0.75 in.(1.9 cm) of liquid equivalent at CLN in upper LCC and 0.19 in.(0.47 cm, in the form of rain or wet snow) of liquid equivalent at KSLC in the SLV (see Figs. 1a and 2a for locations).Although CLN recorded measurable precipitation during each of the three storm stages, KSLC did not record any during the prefrontal and frontal stages, with all the measureable precipitation falling during the passage of the mesoscale snowband in the postfrontal stage (Fig. 4).The remainder of this paper describes the structure and evolution of precipitation as observed by DOW6 during the prefrontal and frontal storm stages, with an emphasis on LCC.We elect to not present the postfrontal stage due to the weak precipitation rates and the lack of variation in precipitation structure in and around LCC, opting instead to focus on the richer comparison between the prefrontal and frontal stages.

Orographic precipitation structures during the prefrontal stage
DOW6 operations for IOP6 commenced at 1405 UTC 12 November, near the midpoint of the prefrontal stage.
At 1500 UTC 12 November, the 500-hPa trough axis and vorticity maximum were over Idaho, with the surface cold front and leading edge of the accompanying 700-hPa baroclinic zone positioned over northwest Utah [Figs. 5a,b (surface front not shown)].The prefrontal atmosphere over the SLV was stably stratified with a near moist-adiabatic lapse rate from the surface to ;700 hPa and a dry subcloud layer with dewpoint depressions as large as 9.58C below 630 hPa (Figs.6a,c).Winds veered from southeasterly at the surface to southwesterly near crest level (;700 hPa) and westerly at 600 hPa (Fig. 6c).This veering wind profile likely reflects weak prefrontal warm advection combined with the channeling of flow in the SLV, the latter consistent with the unsaturated conditions and relatively large squared Brunt-Väisälä frequency (i.e., N 2 d ) below crest level (;700 hPa; Fig. 6b).Although it is challenging to effectively utilize the Froude number (U/N d H), with observed data in regions of highly irregular terrain (e.g., Reinecke and Durran 2008), it does help indicate the tendency for air to move either around or over topographic During this period, the KMTX radar showed persistent quasi-stationary echoes in three areas (e.g., Fig. 5d).The first is in the lowlands upstream of the northern Wasatch Mountains where there was an area of windward precipitation enhancement (hereafter the windward precipitation region; labeled A), as frequently occurs during large-scale southwesterly flow (e.g., Cox et al. 2005).The second is over the central Wasatch Mountains east of the SLV, including upper LCC, where there was a barrier-scale precipitation maximum near the Wasatch Crest (hereafter the barrier-scale precipitation maximum; labeled B).The third is over the southern SLV where a west-east-oriented band of higher reflectivity ( 15 The Up-Canyon North RHI scan (see inset for location) from 1543 UTC 12 November illustrates the wavelike radar reflectivity structure that formed over the SLV and LCC by the cross-valley band and barrier-scale precipitation maximum (Fig. 7).High radar reflectivities associated with the cross-valley band descend from west-east across the SLV and then ascend, broaden, and strengthen within the wavelike barrier-scale precipitation maximum over the middle of LCC before sloping downward over upper LCC.During this period CLN recorded light precipitation (1.5-3 mm h 21 liquid equivalent; Fig. 4a) and upper LCC was visually obscured (not shown).
Time-mean reflectivity RHIs for 1430-1700 UTC 12 November (Figs. 8a-c of LCC, respectively.2Note the eastward displacement (downstream relative to the midlevel flow) of the reflectivity maximum relative to Twin Peaks and Lone Peak, despite the fact that these peaks form the windward face of the central Wasatch Mountains, rising ;2000 m above the SLV.The composite time-mean reflectivity PPI for 1430-1700 UTC 12 November also shows this displacement (Fig. 8d), with the barrier-scale reflectivity maximum located along the Wasatch Crest and east (downstream relative to the midlevel flow) of the initial windward peaks.
The Up-Canyon North time-mean Doppler velocity RHI for 1430-1700 UTC 12 November (Fig. 9a) shows a shallow layer of inbound radial velocities over the SLV, consistent with the southeasterly surface flow (Fig. 5c).This is surmounted by a layer of strong Doppler velocity shear that, like the radar reflectivity maximum, slopes upward over lower and middle LCC, reaching a maximum height ;26 km from the radar before sloping downward over upper LCC (cf.Figs.9a and 8b).The Lone Peak time-mean Doppler velocity RHI for 1430-1700 UTC 12 November shows similar characteristics (Fig. 9b), but with a steeper slope over the western face of Lone Peak leading to a maximum height that is about ;500 m higher than is found over LCC (cf.Figs.9a,b).MesoWest observations from Mount Baldy (AMB; 3373 m; see Fig. 2a for location) and surrounding sites show 20-30 m s 21 south-southwesterly flow penetrating across the Alpine Ridge south of LCC during this period (Fig. 5c).This is a stronger flow with a more southerly orientation than is found at a comparable elevation (;650 hPa) in the 1500 UTC 12 November sounding taken over the SLV (Fig. 6).Thus, beneath the sloping shear layer and inferred layer of cross-barrier flow that generated the barrier-scale precipitation maximum, south-southwesterly flow was able to surmount the Alpine Ridge and descend over LCC.A similar structure in Twin Peak RHIs (not shown) suggests a shallow wavelike flow across the ridge-canyon corrugations of the central Wasatch Mountains.Strong directional shear  LCC (Fig. 11a) that is roughly collocated with the wavelike reflectivity maximum noted in radar reflectivity RHIs (Figs. 7 and 8a-c).Ice crystals are shown in the higher reaches of the precipitating cloud, while aggregates are shown west (upstream) and east (downstream) of the region of low-density graupel.The HYVIS camera documented images of hydrometeors in upper LCC during this period that, due to their shape, size, and the temperature within the cloud layer, are most likely rimed crystals and small graupel (Fig. 11b).The presence of graupel and rimed crystals likely reflects higher cloud liquid water concentrations and stronger wavelike ascent into the barrier-scale precipitation maximum.
Between 1700 and 1900 UTC 12 November, the windward precipitation band upstream of the northern Wasatch Mountains dissipated with the approach of the cold front (not shown).Over and around the SLV, the cross-valley band dissipated, the barrier-scale precipitation maximum persisted in a somewhat weakened state (Fig. 12), and precipitation rates at CLN ebbed (Fig. 4).Although the time-mean Doppler velocity RHIs have gaps in coverage over the valley due to the lack of radar returns throughout much of this period, the velocity signatures above mid-to upper LCC and Lone Peak show that the Doppler velocity shear layer is flatter, no longer exhibiting a wavelike shape over LCC (cf.Figs. 9 and 13).In addition, the hydrometeor-type RHI shows aggregates to be the predominant hydrometeor type over the upper canyon, which is consistent with the observed ice crystal types (Figs.11c,d) and reflects a transition from the graupel and heavily rimed crystals observed earlier.
These changes suggest a shift in the dominant precipitation processes compared to earlier in IOP6 when FIG. 12.As in Fig. 8, but for 1700-1900 UTC 12 Nov. the precipitation rates, radar reflectivities, riming, and graupel content were all greater.This is broadly consistent with weaker large-scale ascent and lower surface-725-hPa relative humidities found in the RUC2 analysis at 1800 UTC 12 November (Fig. 3) and the inferred weakening of the wavelike structure in radar reflectivity and Doppler velocity shear from 1700 to 1900 UTC.However, the wavelike reflectivity feature and area of lowdensity graupel redeveloped (not shown) over upper LCC (Figs. 11e,f) from 1900 to 2000 UTC before the cold front entered LCC.These contrasts and variations suggest that the increase in vertical motion during periods of stronger wavelike undulations across the central Wasatch increased the available cloud liquid water and resultant riming over LCC.

Precipitation structures during the frontal stage
At 2000 UTC 12 November, the surface cold front was moving through the SLV (wind shift visible in Fig. 14c) with the 500-hPa trough axis and 700-hPa baroclinic zone just upstream (Figs.14a,b).Just after frontal passage, at 2030 UTC 12 November, a sounding launched from the DOW6 site showed shallow northwesterly flow with two apparent frontal stable layers between 825 and 775 hPa (Fig. 15).
The shallow cold front, which was distorted by the Oquirrh and Wasatch Mountains as it progressed into the SLV, was much weaker at upper elevations (Fig. 16a; front analyzed at lower elevations where it was well defined).For example, the frontal passage produced a dramatic wind shift (from south to north) and large temperature decline (88C) from 1930 to 2030 UTC at KSLC, whereas only a gradual wind shift (from southsouthwest to west-southwest) and temperature fall (1.48C)occurred at AMB from 2000 to 2300 UTC.A band of radar reflectivity ;30 km wide (hereafter the frontal band) accompanied the cold front, and its presence defines the frontal stages at KSLC (1940-2100 UTC 12 November) and LCC (2025-2130 UTC 12 November).As the frontal band moved into the SLV, it was strongly modulated by the terrain with enhancement on the western slopes of the Oquirrh Mountains, shadowing over the western SLV, and enhancement again over the western slopes of the Wasatch Mountains (Fig. 16a).By 2030 UTC 12 November (Fig. 16b), the front had penetrated to the south end of the SLV and was draped across lower LCC.Enhancement of the frontal band continued over the Oquirrh and Wasatch Mountains while the western SLV remained shadowed.By 2100 UTC 12 November (Fig. 16c), the cold front was south of the SLV and LCC and precipitation was confined to very near to and over the Wasatch Mountains.At the same time, a broader area of postfrontal precipitation stalled over the northern end of the SLV.Up-Canyon South radar reflectivity RHIs between 2026 and 2052 UTC 12 November document the precipitation structure in and around LCC during frontal passage (Fig. 17).As the canyon is oriented roughly west-east and the cold front was oriented southwestnortheast, these RHIs cut obliquely through the front.At 2026 UTC 12 November, the RHI slices through the frontal band ;10 km from the radar, where there is an elevated region of 20-30 dBZ reflectivities (Fig. 17a, labeled area 1).Lower reflectivities below this maximum suggest the sublimation of hydrometeors in the subsaturated air below crest level (Fig. 15), and a photograph taken at 1957 UTC 12 November from the DOW6 site (looking toward the west face of the central Wasatch Mountains; Fig. 18) shows the sublimation of precipitation behind the fractus cloud at the leading edge of the cold front.By 2032 UTC 12 November, the frontal band had moved across the Cottonwood Ridge and into LCC, increasing midcanyon radar reflectivities (Fig. 17b, labeled area 2).This area of high reflectivities became shallower by 2039 UTC 12 November (Fig. 17c, labeled area 3) and gradually became concentrated in the lower canyon (Fig. 17d, labeled area 4).After passage of the frontal band, radar reflectivities diminished significantly and became quite shallow in and around LCC.In the lower canyon, however, a small reflectivity maximum averaging 15-20 dBZ persisted until ;2130 UTC 12 November (not shown).
Time-mean reflectivity RHIs for 2025-2100 UTC 12 November show a region of frontal precipitation within LCC, with the strongest reflectivities over the midcanyon (Fig. 19b).Time-mean reflectivities over the Cottonwood and Alpine Ridges to the north and south did not reach a comparable depth or intensity as those within the canyon during this period (Figs.19a,c).These marked contrasts are also visible in the composite time-mean reflectivity PPI scan for the period (Figs. 19d,e), which shows the highest radar reflectivities within lower LCC as well as over the north-to northwest-facing walls of lower Mill Creek Canyon, Big Cottonwood Canyon, and LCC.Reflectivities were lower at high elevations along the Wildcat, Cottonwood, and Alpine Ridges, implying that the storm during this period was relatively shallow with localized, northwesterly upslope flow immediately behind the front creating orographic enhancement along the north-to northwest-facing canyon walls.
As the cold front progressed through LCC, the hydrometeor-type RHI for 2025-2100 UTC 12 November (Fig. 20a) shows a layer of low-density graupel extending over LCC and the SLV with a shallow layer of aggregates at low levels over the SLV.The vigorous vertical motion often associated with surface-based cold fronts (e.g., Sanders 1955;Carbone 1982;Shapiro 1984) likely resulted in riming, and observations of what is most likely graupel and heavily rimed crystals over upper LCC (Fig. 20b) confirm the hydrometeor classification at higher elevations.There is no way to ascertain the actual hydrometeor types over the SLV at this time because of sublimation and a lack of snow crystal identification observations in the SLV; however, it is unlikely that a precipitating cloud would be composed of such a broad layer of aggregates below a layer of graupel.This suggests either a misclassification of the predominant hydrometeor type over the SLV by the Dolan and Rutledge (2009) algorithm or a shift in the actual hydrometeor types over the SLV during passage of the frontal precipitation band due to rapidly changing atmospheric conditions.Houze and Medina (2005) found a similar misclassification between large aggregates and graupel using the hydrometeor typing algorithm of Vivekanandan et al. (1999).

Discussion
The above analysis illustrates the finescale spatiotemporal variability of orographic precipitation observed in the central Wasatch Mountains of northern Utah during IOP6 (12-13 November 2011) of SCHUSS.Although we caution against viewing IOP6 as a ''typical'' event because the characteristics of winter storms over the central Wasatch Mountains vary and feature a wide range of precipitation processes, rates, and distributions (e.g., Dunn 1983;Steenburgh 2003Steenburgh , 2014)), we believe the event provides useful insight into the finescale nature of precipitation features in areas of complex threedimensional topography and the challenges confronting operational meteorologists and forecast consumers during winter storms in mountainous regions.
During IOP6, persistent quasi-stationary radar reflectivity features were observed in and around LCC and the SLV during the prefrontal stage, including 1) a barrier-scale precipitation maximum over the Wasatch Crest and east, or downstream, of the high windward peaks (relative to the midlevel flow); 2) protrusions of higher reflectivity over the transverse ridges flanking LCC; and 3) a cross-valley band of precipitation that extended from the Oquirrh Mountains across the SLV and into the central Wasatch Mountains.RHI scans through these features revealed a wavelike radar reflectivity structure that sloped downward from west to east across the SLV (cross-valley band), while high reflectivities associated with the barrier-scale precipitation maximum ascended, broadened, and strengthened over middle LCC before descending over upper LCC and the Wasatch Crest.Beneath the midlevel, cross-barrier westerly flow, strong southwesterly crest-level winds appeared to generate shallow wavelike flow across the west-east-oriented ridge-canyon corrugations of the central Wasatch Mountains (summarized conceptually in Fig. 21) and protrusions of higher reflectivity over the ridges flanking LCC.Graupel and heavily rimed crystals were identified by the hydrometeor-type algorithm and were observed in upper LCC during periods when the barrier-scale precipitation maximum featured a robust wavelike appearance suggestive of enhanced barrierscale ascent.In contrast, aggregates were identified by the hydrometeor-type algorithm and were observed in upper LCC during periods when the barrier-scale precipitation maximum was weaker and lacked a wavelike structure.
During the frontal stage, the terrain flanking the SLV modulated the cold-frontal precipitation band, producing shadowing over the western SLV to the lee (east) of the Oquirrh Mountains and enhancement over the eastern SLV windward of the central Wasatch  Wasatch Crest eastward (downstream relative to the midlevel flow) of the high windward peaks, with slight protrusions of higher radar reflectivities along the transverse ridges.In contrast, the composite time-mean reflectivity PPI for the frontal stage (2025-2100 UTC 12 November) shows a near reversal of this spatial pattern with the highest radar reflectivities found in lower LCC and at midelevations along the north-to northwest-facing walls of lower LCC and neighboring canyons.
Meteorologists frequently use climatological precipitation-altitude relationships, such as the one presented in Fig. 1b, to make inferences about the distribution of precipitation during winter storms over complex terrain (Meyers and Steenburgh 2013).This includes the interpretation of radar imagery for inferring precipitation rates or nowcasting, and the downscaling or confidence assessment of the numerical forecast guidance available at grid spacings ranging from ;1 to 25 km.For an event like IOP6, however, such an approach is fraught with problems.During the prefrontal stage, it would have resulted in an overprediction of precipitation along the western abutments of the Alpine and Cottonwood Ridges and possibly also in lower LCC where low-level sublimation reduced precipitation rates.During the frontal stage, the strongest radar reflectivities were found within lower to middle LCC and over the north-to northwest-facing walls of lower LCC and neighboring canyons, whereas climatology would suggest an increase in precipitation with elevation, resulting in the heaviest precipitation in upper LCC and along the Alpine and Cottonwood Ridges.Although there is some potential to predict such variations with cloud-permitting or cloud-resolving modeling systems, current real-time modeling systems run at 4-km grid spacing struggle with event-to-event predictions on these finescales (e.g., Minder et al. 2008).LCC provides an excellent example of how these analysis and forecast uncertainties ultimately affect public safety and economic activities in heavily used mountain areas.SR-210 follows the north side of LCC from the SLV to the terminus of LCC and is traveled by an average of more than 7000 vehicles each day during winter.Fifty snow avalanche paths bisect the highway, adjoining roads, and parking lots, with steep road grades and high snowfall rates frequently leading to hazardous driving conditions and bumper-to-bumper traffic during winter storms (Fig. 22).An average of 33 avalanches hit SR-210 each year, and the combination of frequent large avalanches and the volume and speed of traffic gives the highway one of the highest Avalanche Hazard Indices of any major road in North America (Nepstad et al. 2006;Steenburgh 2014).Avalanche forecasters for SR-210 generally close the highway when the avalanche hazard is high and attempt to trigger avalanches artificially using artillery.Based on data from the 1991/92 ski season, the estimated revenue loss for ski resorts in upper LCC during the closure of SR-210 for avalanche hazard is $1.4 million day 21 [$2.3 million day 21 in 2013 dollars (Blattenberger and Fowles 1995)].
Substantial potential exists to use a gap-filling radar, possibly an X-band system positioned at a location similar to that used in this study, to improve weather analysis and forecasting, avalanche control, and winter road maintenance in LCC.The avalanche hazard during a storm is greatly affected by the intensity and amount of new precipitation, as well as by the temporal and spatial patterns of snow crystal.These storm characteristics can vary significantly within LCC, yet avalanche forecasters currently rely on hourly precipitation and wind observations collected primarily in upper LCC, as well as visual evidence of avalanche activity and crystal type, to evaluate changes in avalanche hazard during storms.The KMTX radar, however, is minimally useful during this process due to poor resolution (horizontal and vertical) and beam blockage.For example, the centroid of the lowest-elevation tilt of the KMTX radar clips the Cottonwood Ridge, resulting in ;50% blockage over LCC assuming a standard index of refraction (Wood et al. 2003).The shallow nature of orographic storms within the central Wasatch Mountains often results in overshooting or limited beam filling and coverage at low levels within LCC is nonexistent.A polarimetric gap-filling X-band radar would allow for the improved observation of precipitation features and estimates of precipitation rates along SR-210 and the starting zones of major avalanche paths, including the midcanyon White Pine and Little Pine avalanche paths (see Fig. 22), which are considered the most problematic (Nepstad et al. 2006;Steenburgh 2014).The ability to visualize precipitation intensity along SR-210 in real time, as would be possible with a gap-filling X-band radar, would greatly facilitate the decision-making process for avalanche forecasters.Although work is needed to test and improve classification algorithms, the application of hydrometeor typing in real time during storms might also prove useful for avalanche hazard assessment.The utility of X-band radar to fill gaps in radar coverage over mountainous terrain has been explored in the European Alps, where networks of X-band radars are currently being installed (e.g., Allegretti et al. 2012), and has been considered in mountainous regions within the United States (e.g., Gourley et al. 2009).In northern Utah, such a radar could be sited to also provide weather surveillance in other high-traffic-volume canyons within the central Wasatch (e.g., SR-190 in Big Cottonwood and Interstate Highway 80 in Parleys Canyon), as well as the SLV.

Conclusions
This paper has examined the finescale structure of a winter storm in the central Wasatch Mountains of northern Utah, with emphasis on Little Cottonwood Canyon (LCC).Significant spatiotemporal variations in storm structure occurred during the prefrontal and frontal storm stages.These small-scale variations pose a challenge for operational forecasters and forecast consumers in mountainous regions and are often not well analyzed or anticipated using conventional NEXRAD, existing operational modeling systems, or forecast techniques that employ climatological precipitation-altitude relationships.In addition, the development of numerical forecast models capable of resolving such finescale features is in its infancy, with large errors noted during individual storms in other mountainous regions (e.g., Minder et al. 2008).Significant potential exists to improve the analysis and nowcasting of precipitation within LCC and the surrounding central Wasatch Mountains using a polarimetric gap-filling X-band radar of the type utilized here.Such a system would enable improved estimates of precipitation rates and hydrometeor types along the canyon highway and within the starting zones of major avalanche paths for the benefit of avalanche control work, winter highway maintenance, and the local economy of upper LCC.Further work is needed to evaluate the economic viability of such a radar and the fidelity of X-band estimates of precipitation rates and hydrometeor types.
FIG. 1.(a) Topography and landmarks of the study region.(b) Mean cool-season (November-April) liquid equivalent precipitation (cm; following scale at left) within the inset box in the bottom right-hand corner of (a).[Data source is PRISM Climate Group, Oregon State University (http://prism.oregonstate.edu),created 18 Jun 2013.]

FIG. 2
FIG. 2. (a) Topography of the central Wasatch Mountains, including locations of RHI scans in and around LCC. Dashed line denotes the hydrological divide, or Wasatch Crest.(b) View of LCC and surrounding topography from the DOW6 site.
FIG. 3. RUC2 KSLC time-height section of u e (black contours every 3 K), pressure vertical velocity [contours every 1 Pa s 21 with red (blue) indicating upward (downward) motion], wind (full and half barbs denote 5 and 2.5 m s 21 , respectively), and relative humidity (%; shaded following scale at right).Times of prefrontal, frontal, and postfrontal stages for CLN (red) and KSLC (blue) are annotated at the bottom.
FIG. 6. Upper-air observations at 1500 UTC 12 Nov at DOW6 site: (a) u and u e , (b) N 2 d and N 2 m , and (c) skew T-log p diagram [temperature, dewpoint, and wind barbs (full and half barbs denote 5 and 2.5 m s 21 , respectively)].

FIG. 8 .
FIG. 8. Time-mean radar reflectivity (dBZ; shaded following scale at right) for 1430-1700 UTC 12 Nov. Insets show the location of each RHI scan over the topography in and around LCC. Shown are RHIs for (a) Twin Peaks, (b) Up-Canyon North, and (c) Lone Peak.Time-mean composite radar reflectivity for 1430-1700 UTC 12 Nov (dBZ; shaded following the scale at bottom left) for (d) the SLV region and (e) the central Wasatch Mountains.

FIG. 9 .
FIG. 9. Time-mean Doppler velocity (m s 21 ; shaded following scale at right, positive values denote flow away from radar) for 1430-1700 UTC 12 Nov. Insets show the locations of each RHI scan over the topography in and around LCC. RHIs for (a) Up-Canyon North and (b) Lone Peak are shown.
FIG. 11.Predominant hydrometeor type (shading following the legend at right where HD and LD refers to high density and low density and drizz/lt rain refers to drizzle or light rain) and images of snowflakes observed in upper LCC during the same time period.Insets show the location of each RHI scan over the topography in and around LCC. Shown are (a),(b) Up-Canyon North, 1430-1700 UTC 12 Nov; (c),(d) Up-Canyon South, 1700-1900 UTC 12 Nov; and (e),(f) Up-Canyon South, 1900-2000 UTC 12 Nov.

FIG. 16 .
FIG. 16.KMTX lowest-elevation angle (0.58) base reflectivity (shaded following scale at bottom left) and manual surface-front analysis (heavy black line) at (a) 2002 UTC 12 Nov, (b) 2031 UTC 12 Nov, and (c) 2100 UTC 12 Nov. Black dots mark the locations of MesoWest stations used in the analysis of the surface-front boundary.
928 W E A T H E R A N D F O R E C A S T I N G VOLUME 29

FIG. 18 .
FIG. 18. Photograph of surface-based cold front looking northeast from the DOW6 site.Note the fractus cloud above the leading edge of the cold front, and the sublimation behind it.