<?xml-model href='http://www.tei-c.org/release/xml/tei/custom/schema/relaxng/tei_all.rng' schematypens='http://relaxng.org/ns/structure/1.0'?><TEI xmlns="http://www.tei-c.org/ns/1.0">
	<teiHeader>
		<fileDesc>
			<titleStmt><title level='a'>Winter precipitation measurements in New England: results from the Global Precipitation Measurement Ground Validation campaign in Connecticut</title></titleStmt>
			<publicationStmt>
				<publisher>Copernicus Publications</publisher>
				<date>11/03/2025</date>
			</publicationStmt>
			<sourceDesc>
				<bibl> 
					<idno type="par_id">10648528</idno>
					<idno type="doi">10.5194/essd-17-5783-2025</idno>
					<title level='j'>Earth System Science Data</title>
<idno>1866-3516</idno>
<biblScope unit="volume">17</biblScope>
<biblScope unit="issue">11</biblScope>					

					<author>Brian C Filipiak</author><author>David B Wolff</author><author>Aaron Spaulding</author><author>Ali Tokay</author><author>Charles N Helms</author><author>Adrian M Loftus</author><author>Alexey V Chibisov</author><author>Carl Schirtzinger</author><author>Mick J Boulanger</author><author>Charanjit S Pabla</author><author>Larry Bliven</author><author>EunYeol Kim</author><author>Francesc Junyent</author><author>V Chandrasekar</author><author>Hein Thant</author><author>Branislav M Notaros</author><author>Gustavo Britto_Hupsel_de_Azevedo</author><author>Diego Cerrai</author>
				</bibl>
			</sourceDesc>
		</fileDesc>
		<profileDesc>
			<abstract><ab><![CDATA[<p>Abstract. Winter precipitation forecasts of phase and amount are challenging, especially in Northeast United States where mixed precipitation events from various synoptic systems frequently occur. Yet, there are not enough quality observations of winter precipitation, particularly microphysical properties from falling snow or mixed phase precipitation. During the winters of 2021–2022, 2022–2023, and 2023–2024, the NASA Global Precipitation Measurement (GPM) Ground Validation (GV) program conducted a field campaign at the University of Connecticut (UConn). The goal of this campaign was to observe various phases of winter precipitation and winter storm types to validate the GPM satellite precipitation products. Over the three winters at UConn, a total of 40 instruments were deployed across two observing sites that captured 117 precipitation events, including 19 phase transition events as indicated by the PARSIVEL2. These instruments included scanning and vertically pointing radars, along with suites of in-situ sensors. In addition, an unmanned aircraft system has been deployed in 2023–2024. Here, an overview of the different field deployments, instrumentation, and the datasets collected are presented. To showcase the observations, this article features a wide-ranging set of measurements collected from the instrument suite for the 28 February 2023 storm, during which six to eight inches of snow accumulated at the two different observing sites. Also included is a discussion on how these observations can be combined with other datasets to validate ground-based and remote sensing measurements and highlight important atmospheric processes that impact winter precipitation phase and amount. The datasets collected from this GPM GV field campaign are available at https://doi.org/10.5067/GPMGVUCONN/DATA101 (Cerrai et al., 2025).</p>]]></ab></abstract>
		</profileDesc>
	</teiHeader>
	<text><body xmlns="http://www.tei-c.org/ns/1.0" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xmlns:xlink="http://www.w3.org/1999/xlink">
<div xmlns="http://www.tei-c.org/ns/1.0"><head n="1">Introduction</head><p>The Northeast United States is exposed to hazardous winter weather with heavy snowfall, high winds, and freezing rain each year. Owing to its terrain and conducive meteorological conditions, including strong thermal gradients between the land and ocean, the east coast of the United States is susceptible to major winter storms <ref type="bibr">(Maglaras et al., 1995)</ref>. Forecasting when and where precipitation will occur remains one of the key challenges for forecasters. Not only do forecasts have to accurately predict the amount and intensity of precipitation in the correct location but also must identify what phase of precipitation will occur in each location. If there are multiple precipitation phases occurring at once, forecasters must correctly determine the precipitation amount for each precipitation phase <ref type="bibr">(Maglaras et al., 1995;</ref><ref type="bibr">Ralph et al., 2005;</ref><ref type="bibr">Novak et al., 2023)</ref>. Subtle shifts in the environmental conditions over an area can create significant forecast changes and different impacts than forecasted <ref type="bibr">(Novak et al., 2023)</ref>.</p><p>In addition to forecasting challenges, there are issues with the reliability of observations for specific winter weather phenomena. Since the types of observations vary so much, it can be difficult to obtain trustworthy and accurate measurements particularly for snow and ice <ref type="bibr">(Rasmussen et al., 2012;</ref><ref type="bibr">Hurwitz et al., 2020)</ref>. These range from manual observations, in situ sensors, and remote sensing instruments that can be as far away as outer space. Manual observations can vary significantly because some may be reported by trained National Weather Service (NWS) cooperative volunteer observers or other trained observers, such as CoCo-RaHS <ref type="bibr">(Cifelli et al., 2005)</ref>, whereas others can be reported by any person with a mobile app (e.g., mPING, <ref type="bibr">Elmore et al., 2014)</ref> on their phone. This creates a challenge for processing what exactly was measured and how accurate said measurement was. While many in situ sensors can accurately observe winter weather, they still have limitations <ref type="bibr">(Martinaitis et al., 2015)</ref>. This mainly is an issue at Automated Surface Observing System, ASOS <ref type="bibr">(NOAA, 1998)</ref> for precipitation accumulations caused by undercatch, a phenomenon that occurs with precipitation gauges where the orifice size of the measuring device prevents an accurate measurement due to the restrictive collection area <ref type="bibr">(Rasmussen et al., 2012;</ref><ref type="bibr">Martinaitis et al., 2015)</ref>. In addition, non-heated tipping buckets are not effective for frozen precipitation where accumulations will not be measured until the ice melts, potentially days or multiple storms later. Heated tipping buckets can fix some of these issues, but they still are not as accurate for frozen precipitation as weighing bucket precipitation gauges <ref type="bibr">(White et al., 2004;</ref><ref type="bibr">Greeney et al., 2005;</ref><ref type="bibr">Tokay et al., 2010)</ref>. Retrievals from remote sensing instruments are helpful especially where surface observations are not available. However, they have larger uncertainties and require refinement, particularly for frozen winter precipitation <ref type="bibr">(Kulie et al., 2010)</ref>. One significant challenge for frozen precipitation is incorporating the variable particle shape, mass, and density <ref type="bibr">(Brandes et al., 2007;</ref><ref type="bibr">Kulie et al., 2010;</ref><ref type="bibr">You et al., 2017)</ref>, which affects scattering and extinction properties retrieved by these instruments. In general, snow measurements can be impacted by their environment more than other precipitation phases; there is an urgent need to improve spatial and temporal frequency of snowfall observations <ref type="bibr">(Doesken and Robinson, 2009;</ref><ref type="bibr">Hurwitz et al., 2020)</ref>. Many environmental variables, such as wet-bulb temperature, particle size and fall speed, and snow density, contribute to identification of precipitation phase and forecasting snowfall amounts <ref type="bibr">(Rasmussen et al., 2012;</ref><ref type="bibr">Hurwitz et al., 2020)</ref>, but they are not measured by any operational instruments.</p><p>A combination of three observational platforms is ideal to have a robust estimate of precipitation phase and amount with the given uncertainty of each platform. This is feasible with targeted field campaigns, which have contributed to improving forecasters' understanding of thermodynamics and microphysical processes <ref type="bibr">(Novak et al., 2023)</ref>. To this point, a suite of in situ and remote sensing instruments were deployed at the University of Connecticut (UConn) in Northeast United States as a part of a NASA Global Precipitation Measurement (GPM) Ground Validation (GV) field campaign. The timing of this GPM GV campaign allowed for synergistic overlap with the Investigation of Microphysics and Precipitation of Atlantic Coast-Threatening Snowstorms field campaign <ref type="bibr">(IMPACTS, McMurdie et al., 2022)</ref> in the winters of <ref type="bibr">2021-2022 and 2022-2023.</ref> During the overlapping field deployment seasons, three snowfall events were sampled by instruments from both campaigns by having airborne observations of flight legs over or near the UConn site to be combined with remote sensing and ground-based observations at or around UConn. The length of the field campaign, over four months, allowed for the collection of a rich data set from numerous rain, snow, and mixed phase transition events.</p><p>The GPM mission <ref type="bibr">(Hou et al., 2014;</ref><ref type="bibr">Skofronick-Jackson et al., 2017)</ref> was established to unify and advance precipitation measurements from space-borne radar and passive microwave sensors. This is an extremely important mission worldwide because GPM can measure precipitation in areas without observational networks and over the oceans. GPM presented an improvement over its predecessor, the Tropical Rainfall Measurement Mission (TRMM), because of its expansion to include the mid latitude region and ability to detect light precipitation rates and snowfall, albeit not without its challenges <ref type="bibr">(Casella et al., 2017;</ref><ref type="bibr">Milani and Kidd 2023)</ref>. The GPM mission has a core observatory along with a broader constellation of satellites. The core observatory is equipped with the Dual-Frequency Precipitation Radar (DPR) and the GPM Microwave Imager (GMI), a passive microwave radiometer across multiple frequencies. These spaceborne remote sensing instruments can observe a large portion of the globe (&#177;68&#176;of latitude) during both day and night which makes it especially useful for winter weather observations in the midlatitudes.</p><p>While there have been several GPM GV field campaigns across the United States and the globe, only a few have focused on winter weather. These include: LPVEx, the Light Precipitation Validation Experiment <ref type="bibr">(Petersen et al., 2011)</ref> in Finland during the fall of 2010; GCPEx, the GPM Cold-season Precipitation Experiment <ref type="bibr">(Skofronick-Jackson et al., 2015)</ref> in central Ontario, Canada during the winter of 2011-2012; OLYMPEX, the Olympic Mountains Experiment <ref type="bibr">(Houze et al., 2017)</ref> in the Olympic Peninsula of Washington in the winter of 2015-2016; and ICE-POP, International Collaborative Experiment for PyeongChang Olympic and Paralympics <ref type="bibr">(Petersen et al., 2018)</ref> in South Korea during February 2018. There have been a few recent field campaigns (not GPM GV related) that have been undertaken in Northeast United States, (OWLES - <ref type="bibr">Kristovich et al., 2017, IMPACTS, and</ref><ref type="bibr">WINTRE-MIX -Minder et al., 2023)</ref>, but these focused on specific types of winter phenomena: lakeeffect snow, snowstorms, and precipitation phase, respectively.</p><p>Between 1980 and 2023, there have been 22 winter storms that caused at least USD 1 billion of damage in the Continental United States (Consumer Price Index-adjusted); the total estimated loss for the 22 events was USD 97.2 billion (NOAA NCEI 2023). It is important to note that 21 of these 22 events were in the Northeast Climate Region, which has a high concentration and frequency of winter storms compared to the rest of the United States <ref type="bibr">(Changnon 2007;</ref><ref type="bibr">Lazo et al., 2020;</ref><ref type="bibr">NOAA NCEI 2023)</ref>. With the UConn site being located in southern central part of the Northeast United States, it is in an ideal location, due to being far enough inland that colder air is more prolonged when winter storms happen, to sample a wide variety and intensity of winter precipitation that previous winter GPM GV campaigns did not experience.</p><p>This article presents an overview of the GPM GV campaign at UConn for which the main objectives were to collect a rich and consistent set of observations for winter precipitation to validate and improve GPM based retrieval algorithms during winter storms in the Northeast United States. With 117 precipitation events across the three winters, a myriad of winter weather phenomena was captured ranging from heavy snow to freezing rain to ice pellets to events with multiple phase transitions. This unique dataset can be used for examining the GPM's retrieval algorithms for winter weather, validating precipitation phase algorithms or numerical model simulations, and analysing physical processes that can explain key components of winter phenomena, such as snow microphysics. Section 2 presents an overview of the different field deployments and instrumentations. Section 3 details the data processing for each instrument, which includes calibration or validation notes. Section 4 highlights an impactful snowstorm as case study from the 2022-2023 winter to spotlight the variety of observations from the instrumentation. Section 5 details the data availability; Sect. 6 closes with concluding remarks and points to directions of future research.</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head n="2">Deployment Sites and Instrumentation</head><p>A suite of ground-based instruments was assembled with the objective of collecting observations of winter precipitation, improving satellite retrieval algorithms, and understanding the many physical processes and properties of winter weather phenomena. These observations aimed at providing a more complete picture of microphysical and thermodynamic processes that occur near the surface or in the lowest levels of the troposphere. Since it takes significant resources to have such a complete set of instrumentation, including many repetitive measurements across instruments, it was important to obtain as many observations as possible. The location of the UConn sites aligns generally within the swath of many of the winter storms that occur in the Northeast United States (Fig. <ref type="figure">1a</ref>) as shown previous climatological studies <ref type="bibr">(Miller, 1946;</ref><ref type="bibr">Davis et al., 1993)</ref>. UConn's inland location provides the opportunity for prolonged periods of below freezing temperatures for storms that pass to the south and east. Additionally, storms that move from southwest to northeast by UConn allow for the opportunity for phase transitions due to the location of the warm and cold fronts.</p><p>The instruments were deployed over three separate winter seasons, <ref type="bibr">2021-2022, 2022-2023, 2023-2024</ref>. Two separate locations were used during the three years of deployments and will be labeled as GAIL and D3R. The GAIL site (41&#176;48 28 N, 72&#176;17 38 W, elevation of 149 m) and the D3R site (41&#176;49 04 N, 72&#176;15 27 W, elevation of 213 m) are close in distance, approximately 3.2 km apart, but the local terrain could significantly impact the observed precipitation. The GAIL site was located near the base of the Willimantic River valley, whereas the D3R site was located at the top of a hill above the river valley; the 64 m elevation difference could cause differences in precipitation onset, phase, and amount to vary across the two sites (Fig. <ref type="figure">1b</ref>). Because the goal of the observations collected by the instrument suites aims to improve the GPM satellite retrievals, all the instruments provide measurements that can be directly compared to the GPM satellite retrieval algorithms output or be used in the development or refinement of the algorithms. In the subsequent subsections, a description of each instrumentation and location by deployment year.</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head n="2.1">Instrumentation</head></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head n="2.1.1">All-in-One 2 (AIO)</head><p>The All-in-One 2 (AIO, Fig. <ref type="figure">2</ref>) is a sonic weather sensor, produced by Met One, that measures wind speed and direction, ambient air temperature, relative humidity, and barometric pressure, thereby providing the necessary background context for precipitation observations. The temperature and relative humidity sensors are unaspirated and filtered. AIO collects high resolution data at 1 Hz; it has an accuracy of &#177;0.2 &#176;C for ambient air temperature and of &#177;3 % for relative humidity, which is essential for understanding precipitation phase. Observations from the AIO are recorded as 1 min averages. The observations from the AIO can enhance the interpretations and limitations of precipitation observations, such as snow microphysics <ref type="bibr">(King et al., 2024)</ref>. The GPM goal is to observe global precipitation; thus, the ability to constrain the observed meteorological conditions can enhance and validate the end outcome (how much of each precipitation phase fell) and the important quantities remotely sensed by the satellite (reflectivity, precipitation rate, particle size, etc.).</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head n="2.1.2">Micro Rain Radar (MRR)</head><p>The Micro Rain Radar (MRR, Fig. <ref type="figure">3</ref>) is a 24 GHz K-band vertically pointed profiling radar, with a heated radar dish, produced by METEK. It is effective at sampling the lowest levels of the atmosphere and identifying doppler spectra of hydrometeors, which can be used to identify bright brands that indicate falling snow particles melting to rain droplets. Two different METEK MRRs were used during the field deployments.</p><p>In 2021-2022, an MRR-2 was deployed. The MRR-2 has 32 range gates with a minimum detectable reflectivity of 2 dBZ at 1000 m (METEK MRR-2). During its deployment, the vertical resolution was set to 35 m and had a 30 s raw sampling frequency which was averaged to 1 min, which allows for observations of the lowest 2 km of the atmosphere (the two range gates nearest the surface bins are not typically used due to potential ground clutter issues). Raw MRR-2 observations collected were for reflectivity, doppler velocity, and spectrum width with additional processing creating additional quantities such as rain rate and liquid water content. The MRR-2 was designed for liquid precipitation, but a com- monly used post-processing algorithm <ref type="bibr">(Maahn and Kollias, 2012)</ref>, can be applied to improve noise removal and add a de-aliasing component to provide more effective Doppler velocity, reflectivity, and spectral width specifically for frozen and mixed phase precipitation.</p><p>During 2022-2023 and 2023-2024, an MRR-Pro was deployed instead of the MRR-2; this was done due to constraints on the number of available instruments. The MRR-Pro is similar to the MRR-2 and had the same vertical range resolution but has additional range gates, set to 128 for both deployments, and a minimum detectable reflectivity of -8 dBZ at 1000 m (METEK MRR-Pro). Originally in 2022-2023, the range gates were set to 256, but there were issues with the vertical velocity retrievals and was shifted to 128 to avoid this issue. The same sampling frequency and averaging periods were kept between campaigns to create a similar dataset. Unlike for the MRR-2, there is no additional software or algorithms to process MRR-Pro observations for frozen precipitation; METEK has added new internal processing steps to the MRR-Pro that has provided some improvements, but there remains a need for the creation of a processing software similar to <ref type="bibr">Maahn and Kollias (2012)</ref> for the MRR-Pro.</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head n="2.1.3">Pluvio 2 Weighing Gauge (Pluvio)</head><p>The Pluvio 2 Weighing Gauge (Pluvio, Fig. <ref type="figure">4</ref>) is a weighing precipitation gauge manufactured by OTT Hydromet designed to collect all phases of precipitation. Pluvio records precipitation amount and intensity in both real time (output delay of 20 s), and non-real time (output delay of 5 min). With a sampling interval of 1 min, the accuracy of precipitation amount is &#177;0.05 mm over a 60 min measuring window and &#177;0.1 mm min -1 for precipitation intensity. Two types of Pluvio were deployed during the campaign with one being a 200 cm 2 collection area and the other being a 400 cm 2 collection area. The 200 cm 2 model was heated during its deployment. Pluvio's can be subjected to the phenomenon of undercatch due to the orifice limiting accurate totals <ref type="bibr">(Rasmussen et al., 2012;</ref><ref type="bibr">Martinaitis et al., 2015)</ref>, but a double alter wind shield was deployed around it to reduce this phenomenon. In addition, the redundant measurements with other instruments can be used to validate the observations. The measurements from the Pluvio can be used to calibrate the GPM satellite or ground-based radar retrieval algorithms for precipitation amount and rate.</p><p>2.1.4 Platform for In situ Estimation of Rainfall Systems (PIERS+) PARSIVEL 2 PARSIVEL 2 (Particle Size Velocity), a laser-based disdrometer produced by OTT Hydromet, was primarily developed to measure the size and fall velocity of raindrops (Loeffler-Mang and Joss, 2000) but can also measure the size and fall velocity of snowflakes (Loeffler-Mang and Blahak, 2001) with limited accuracy. The PARSIVEL 2 is at the center of a NASA designed Platform for In situ Estimation of Rainfall Systems (PIERS+) platform shown in Fig. <ref type="figure">5</ref>. Particles passing through the laser beam (30 mm wide, 1 mm high, 180 mm long) of the PARSIVEL 2 cause voltage decrease by light extinction. The amplitude and duration of induced voltage signal is related to particle size and fall velocity, respectively. PARSIVEL 2 raw output is 32 by 32 size and fall velocity matrix. Due to sensitivity of laser device to the small particles with minimum detectable size 0.25 mm in equivalent diameter, the first two size bins in the PARSIVEL 2 observations are always empty. In PARSIVEL 2 the bin width increases from 0.125 to 3 mm with size due to the increasing uncertainty of particle size measurement. Therefore, it is not reliable to measure the maximum drop diameter at high accuracy in moderate to heavy rain (&gt; 5 mm h -1 ), and the complex shape of falling snowflakes can be an issue to measure the size accurately. PARSIVEL 2 is also be used to determine present weather and visibility based on fall velocity versus size measurements, which can give synoptic weather codes similar to ASOS stations. Processed PARSIVEL 2 observations can give important microphysical properties like particle size distributions (PSD), liquid water content, and number concentration.</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head>Tipping Bucket Rain Gauges</head><p>Tipping Buckets are a pair of Met One 380 rain gauges deployed on each PIERS+ platform on either side of PARSIVEL 2 . The gauges are made of metal and are used to measure rain only from GPM GV activities as these gauges are not heated. Each gauge utilizes a dual-chambered tipping bucket that is self-emptying via a hole in the bottom of the mechanism. Each tip of the bucket indicates a rain accumulation of 0.01 inch or 0.25 mm. As the bucket tips, it logs the time to the nearest second which can be recorded on a logger or other device. Similar to the Pluvio, the tipping buckets record precipitation amount and intensity and can be subjected to undercatch. However, the pair of buckets allow for verification of the others' observations, which is helpful for identifying instrument errors.</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head>RM Young Mechanical Anemometer</head><p>Attached to the PIERS+ platform is a mechanical RM Young anemometer which is used to measure wind speed and direction. This instrument provides redundancy to the AIO sonic anemometer and allows for verification between the two sensors. The mechanical anemometer has an accuracy of &#177;0.3 m s -1 , with a minimum speed of 1.0 m s -1 , and &#177;3&#176;for direction. The sampling frequency was 1 Hz. Since only one PIERS+ platform was deployed in 2022-2023 and 2023-2024, there was only one mechanical anemometer each deployment year.</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head n="2.1.5">Precipitation Imaging Package (PIP)</head><p>The Precipitation Imaging Package (PIP, Fig. <ref type="figure">6</ref>, <ref type="bibr">Pettersen et al., 2020</ref><ref type="bibr">Pettersen et al., , 2021;;</ref><ref type="bibr">Helms et al., 2022;</ref><ref type="bibr">Tokay et al., 2023;</ref><ref type="bibr">King et al., 2024</ref>) is a NASA-developed video disdrometer consisting of a high-speed camera pointed at a halogen light source. The light source, located 2 m away from the camera lens, acts to backlight the precipitation particles that fall through the sampling volume, which is located 1.33 m from the lens. Critically, PIP uses an open sampling volume that is sufficiently far from the instrument components to prevent contamination of the measurements by particle shattering and other effects of the flow around the instrument under typical wind conditions <ref type="bibr">(Newman et al., 2009)</ref>. The PIP images are 640 by 480 pixels with a calibrated pixel size of 0.1 mm by 0.1 mm and are continuously collected at a rate of 380 frames per second. Due to bandwidth limitations, the images are compressed in the vertical by averaging adjacent pixels resulting in an effective pixel size of 0.1 mm by 0.2 mm. Although the actual sampling volume dimensions are a function of particle size <ref type="bibr">(Newman et al., 2009)</ref>, the sampling volume is no greater than 64 mm by 48 mm in the viewing plane of the camera and has a depth of field of approximately 117 times the maximum particle dimension (0.23 m for a 2 mm).</p><p>The PIP analysis software generates a wide variety of direct particle measurements, although the two key measurements are the equivalent diameter and vertical motion as these are used to produce the PSD and the particle density estimates. PIP produces multiple levels of data output, as described in King and Coauthors 2024, and its highest order processed observations include particle count, rain and non-rain rate, particle density, and particle fall speed.</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head n="2.1.6">Aerosol, Cloud, Humidity Interactions Exploring and Validation Enterprise (ACHIEVE)</head><p>The NASA-Goddard Space Flight Center Aerosol, Cloud, Humidity Interactions Exploring and Validation Enterprise (ACHIEVE, Fig. <ref type="figure">7</ref>) is a transportable ground-based laboratory for obtaining remotely sensed measurements of cloud and precipitation properties <ref type="bibr">(Loftus et al., 2016)</ref>. ACHIEVE's primary instrument is the scanning W-band (94 GHz) Doppler cloud radar, manufactured by ProSensing Inc., and provides high-resolution measurements of reflectivity, Doppler velocities, Linear Depolarization Ratio (LDR), Doppler spectrum width, as well as Doppler power spectra. ACHIEVE was deployed in the 2022-2023 campaign at the D3R site. During the deployment, the W-band radar oper- ated mostly in vertically pointing mode in order to capture the evolution of cloud and precipitation system vertical structures as they passed overhead. During active precipitation periods, additional Range Height Indicator (RHI) scans were performed every 10 min over the GAIL site to the west as well as N-S oriented in an effort to sample spatial variability of the systems and for obtaining near-coincident data with the D3R scanning radars. Coincident RHI scans along the ER-2 flight paths were also conducted during overflights on 28 February 2023. The W-band radar collected data from 11 precipitation events, six of which were classified as phase transition events (rain to snow or snow to rain), and one mixed phase event with graupel and freezing rain.  (shown in Fig. <ref type="figure">1a</ref>) with 150 m range resolution. With this configuration, the minimum detectable signal at 15 km is -8 dBZ for the Ku band and -3 dBZ for the Ka band. The D3R operates like a more traditional radar with both plan position indicator (PPI) and RHI scans compared to the MRRs and ACHIEVE; the PPI scans covered the areas WSW to ENE of the radar location, which includes the area over the GAIL site. Individual RHIs in the direction of the GAIL site, as well as true north, were also conducted to match the ACHIEVE. Besides its observations of reflectivity and doppler velocity, the dual-polarization nature of the D3R allows for the hydrometeors to be observed with greater detail with a suite of dual-polarization products. The D3R's ability to provide dual-polarized derived products is necessary for understanding the microphysical characteristics of the falling particles.</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head n="2.1.8">Snowflake Measurement and Analysis System (SMAS)</head><p>The Snowflake Measurement and Analysis System (SMAS, Fig. <ref type="figure">9</ref>) is a novel system for measurement and analysis of snow particles in freefall, designed, developed, and built in the Electromagnetics Laboratory at Colorado State University (CSU). The main features of the SMAS are: seven highresolution cameras strategically placed in a 3D fashion for reconstruction of shapes of snowflakes based on photographs of particles in freefall from multiple views <ref type="bibr">(Kleinkort et al., 2017)</ref>, capability for measurement and analysis of multiple snowflakes at once, namely, processing of images with multiple snowflakes per frame, and comprehensive characterization of microphysical properties of snowflakes. Some parameters of the SMAS are 1440 cm 3 imaging volume, 5-Megapixel cameras, and 31 &#181;m pixel resolution (Thant et al., 2022). The instrument is also equipped with six flashes, five sensor boards, two lasers, a printed circuit board, which triggers the cameras and flashes when a snowflake falls through the system and the laser plane gets blocked, and a built-in computer, which organizes and stores the images captured by the cameras. Measured and estimated particle properties are fall speed, density, effective dielectric constant <ref type="bibr">(Notaro&#353;, 2021)</ref>, and scattering observables <ref type="bibr">(Chobanyan et al., 2015)</ref>. The SMAS is also capable of calculating the PSD <ref type="bibr">(Huang et al., 2017)</ref>, e.g., the particle cross section area and diameter, respectively, in real time. Multiple studies have shown that neural networks, can be applied to classify snowflakes into six geometrical categories (small particle, planar crystal, columnar crystal, combination of columnar and planar crystal, aggregate, and graupel), as well as five separate degrees of riming can be observed, which is essential for understanding the formation and meteorological environment of snow <ref type="bibr">(Hicks and Notaro&#353;, 2019;</ref><ref type="bibr">Key et al., 2021;</ref><ref type="bibr">Thant et al., 2023a)</ref>.</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head n="2.1.9">CSU-Modified Multi-Angle Snowflake Camera (MASC)</head><p>The Multi-Angle Snowflake Camera (MASC, Fig. <ref type="figure">10</ref>; <ref type="bibr">Garrett et al., 2012)</ref>, a commercially available instrument, uses three cameras in the horizontal plane separated by 36&#176;to capture high-resolution photographs of snowflakes or other frozen hydrometeors in freefall from three views, while simultaneously measuring their fall speed. The pixel resolution is 50 &#181;m, the virtual measurement area is 30 cm 2 , and the measurement volume rate is &#8764; 200 cm 3 s -1 . The instrument includes two near-IR emitter-receiver pair arrays positioned one above the other, and as a particle falls through the lower array, it triggers all cameras. Additionally, the particle fall speed is obtained from the fall time between the two triggers.</p><p>To enable 3D shape reconstruction, we added two cameras to the CSU-MASC, "externally", above the original cameras and at a 55&#176;angle with respect to the horizon, to provide additional views <ref type="bibr">(Kleinkort et al., 2017;</ref><ref type="bibr">Notaro&#353;, 2021)</ref>. Additionally, we can estimate PSD <ref type="bibr">(Huang et al., 2017)</ref> and classify snowflake geometry, riming degree, and melt/dry state <ref type="bibr">(Hicks and Notaro&#353;, 2019;</ref><ref type="bibr">Key et al., 2021;</ref><ref type="bibr">Thant et al., 2023b)</ref>.</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head n="2.1.10">WxUAS</head><p>The Oklahoma StormTrooper is an Unmanned Aerial Weather Measurement System (WxUAS, Fig. <ref type="figure">11</ref>) capable of sampling the low-altitude wintry precipitation environment, including mix-phased and freezing precipitation (Britto Hupsel de Azevedo 2024; <ref type="bibr">Azevedo et al., 2024</ref><ref type="bibr">Azevedo et al., , 2025))</ref>. This WxUAS provides in-situ samples of pressure, temperature, relative humidity, and particle size counts. Additionally, it provides remote samples of reflectivity and velocity from an onboard, vertically pointing 74 GHz millimeter wave (mmWave) radar with 0.4887 m resolution (see Table <ref type="table">C</ref>.1 of Britto Hupsel de Azevedo, 2024 for radar specifications). Leveraging the mobility of UAVs, the proposed WxUAS creates a spatial distribution of the measured atmospheric parameters. Via repetition of the flight pattern at regular intervals, the StormTrooper captures the temporal evolution of the measured parameters' spatial distribution. Together, the spatiotemporal aspects of this WxUAS' measurements produce insight into atmospheric mechanisms that govern the low-altitude wintry precipitation environment. Temperature, relative humidity measurements and particle counts are produced in an actively aspirated sensor shield at 1 Hz. Due to its shape, this shield design can leverage gravity to mechanically separate large precipitation drops from the air. This separation occurs because larger droplets are dominated by their inertia while small droplets will follow the streamlines on the intake system according to Stokes flow (Britto Hupsel <ref type="bibr">de Azevedo, 2024;</ref><ref type="bibr">Azevedo et al., 2024</ref><ref type="bibr">Azevedo et al., , 2025))</ref>. This flow characteristic allows the system to sample the atmospheric conditions without exposure to harmful freezing precipitation. The 74 GHz, vertically pointing, phased array radar can use beam forming techniques to compensate aircraft motion, making Doppler velocity an approximation of hydrometeor fall velocity. WxUAS provides a unique perspective that connects the ground instrumentation with remote sensing instruments collecting observations above ground level. Its utility in capturing thermodynamic changes in the lowest levels of the atmosphere is essential for understanding precipitation events in near freezing conditions.</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head n="2.1.11">Ceilometer</head><p>Ceilometers are useful tools for observing overhead cloud and visibility conditions and can provide additional information about the atmospheric boundary layer. This ceilometer, produced by Vaisala (Fig. <ref type="figure">12</ref>), uses pulse diode laser Light Detection And Ranging, LIDAR, technology at 910 nm wavelengths; the vertical observed profiles are able to measure aerosol backscatter up to 15 km, and retrieve cloud base heights and boundary layer structure up to 4.5 km. These measurements can be compared to estimates of cloud base levels from the WxUAS, which provides validation and context of the vertical profiles it has collected. Observations of cloud layers can also be compared to ground-based radar observations for additional insight.</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head n="2.2">Deployment Sites</head></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head n="2.2.1">2021-2022</head><p>During the first winter of the campaign, the GAIL site was used to host all the instrumentation. A core set of instrumentation was identified to be located at UConn during this and all future winters. This set consisted of: (i) an MRR (MRR-2 was used in 2021-2022 only), (ii) a Pluvio, (iii) an AIO, (iv) a PARSIVEL 2 , and (v) a PIP. This core instrumentation provided comparable and consistent measurements year to year and collected a wide array of observations from microphysical observations, precipitation amounts, meteorological conditions, and ground-based radar. The instruments were deployed from 1 December 2021 to 26 April 2022.</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head n="2.2.2">2022-2023</head><p>For the second deployment, the goal was to increase the diversity of instrumentation as well as create an additional deployment site at a higher elevation. This led to instruments being located at both the GAIL and D3R sites with the same core instruments at both sites, except for the MRR-2 being replaced by an MRR-Pro. PIERS+ were also added to the core instrumentation set since there is a lack of qualitycontrolled rain gauge observations near UConn. The additional precipitation measurements were also valuable during the warm months of the deployment because they could be compared with the Pluvio for further validation. At the D3R site, the D3R and the 94 GHz W-band radar mounted on the ACHIEVE trailer were co-located to allow for comparable radar observations. At the GAIL site, additional instruments to measure microphysical properties of snow were included. The SMAS and MASC collected observations to help verify and enhance the observations from the PIP. The instruments were deployed at UConn from 23 December 2022 to 11 April 2023. </p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head n="2.2.3">2023-2024</head><p>During the 2023-2024 deployment, the goal was to further validate the instrumentation used in previous deployments, as well as to expand the vertical atmospheric profile measurements. To do this, only the GAIL site was used, and two sets of each core instrumentation was deployed in close proximity to one another; this technique has been used before in other field campaigns and classified as a "supersite" <ref type="bibr">(Gultepe et al., 2019;</ref><ref type="bibr">Kim et al., 2021)</ref>. In addition to the extra instrumentation, a Vaisala CL-51 ceilometer was deployed to provide measurements essential for comparison with the low-altitude WxUAS. The WxUAS is a novel tool to sample the lowest levels of the atmosphere at quicker intervals than would be possible with radiosondes. The WxUAS was deployed in the style of Intense Observing Periods (IOPs) as opposed to the instruments, which were collecting observations continuously. The instrumentation was deployed from 15 December 2023 to 21 May 2024. For an overall summary, Table <ref type="table">1</ref> summarizes the instrumentation deployed during each campaign; Fig. <ref type="figure">13</ref> showed photographs of how the deployment sites were setup in 2022-2023 and 2023-2024.</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head n="3">Data Processing and Quality Control</head><p>In this section, a list of precipitation events across the three years of deployments are derived from the PARSIVEL 2 to provide a starting point for specific research applications (Appendix A). To accompany the event list, additional indepth descriptions are given for the quality control or data collection on certain instrumentation. All datasets come with a base level documentation about the quality control performed, but here, more in depth analysis and explanations are provided for non-standard quality control or data collection practices.</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head n="3.1">Precipitation Events</head><p>Over the three-year deployment, 117 unique precipitation events were observed at the UCONN site which includes 19 snow only events, 79 rain only events, and 19 mixed precipitation phase events. Event start and end times were identified based on the PARSIVEL 2 observations; the event phase is determined based on the present weather algorithm. Events are defined as periods of precipitation separated by greater than 60 min of no precipitation. Table <ref type="table">2</ref> lists out the number of events by phase and year; Appendix A contains the full list of events and times. Mixed phase events can have one or more phase transitions between rain and snow or other frozen precipitation phases.</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head n="3.2">Data Quality Control and Processing</head></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head n="3.2.1">All-in-One 2 (AIO)</head><p>The AIO provides a standard collection of surface meteorological observations to help provide context for the precipitation and microphysical measurements from other instruments. Of the five key measurements it records (temperature, relative humidity, surface pressure, wind speed and direction), temperature, surface pressure and wind speed required additional quality control based on initial comparisons to other nearby surface stations. Relative humidity was too difficult to quality control because of the distance to the nearest reference surface observations; wind direction quality control could be connected to the wind speed since invalid wind speed observations likely would have invalid direction measurements. In addition to the quality control of the observations, each year of data collected was also checked from missing observation minutes; these were filled in during the quality control with -9999 and a summary table of the AIO deployments is available in Table <ref type="table">3</ref>.</p><p>To verify the AIO observations, the reference data used for comparisons were from the 1 min ASOS observations at Willimantic, CT (KIJD), daily climate records from the  Temperature and pressure were validated using the ASOS observations as that is the nearest surface observations to the deployment location; due to the chaotic nature of wind measurements, the RM Young anemometer was used to quality control the wind speed observations due to it being colocated with the AIO. The KIJD ASOS site is located to the SE of the deployment sites at 76 m elevation and 10.5 km away from D3R and 11.7 km away from GAIL. The Storrs GHCN site is located to the SE of the deployment sites at 202 m elevation and 3.5 km from D3R and 5.6 km away from GAIL.</p><p>For the quality control of the surface pressure, all observations were converted to sea level pressure to remove the bias due to the elevation differences. Outlier sea level pressure values, outside the range of 920 to 1050 hPa, were removed and marked as erroneous; the remaining values were used to calculate the bias and standard deviation for each respective year. Additional outliers were identified by identifying periods where the differences between the AIO and KIJD outside the bias &#177;2 standard deviations. All outliers were manually checked against the KIJD observations. If the outlier occurred for only a 1 min period, then the value was corrected using linear interpolation of the observations 5 minutes on either side of the outlier. All outliers that were erroneous values were changed to -9999.</p><p>For the wind speed, initial plots were made of the distribution of wind speeds from the AIO and RM Young; there was a clear high bias in the AIO and produced numerous outlier observations. All wind speeds more than 40 m s -1 were marked as erroneous. Based on these distributions, two thresholds were set based on whether the AIO was located on the ground, or on the roof of a trailer with the threshold being 20 m s -1 for the ground and 25 m s -1 for the roof locations. These were determined because both the RM Young and ASOS did not report 1 min wind speeds in excess of 20 m s -1 over the three-year period, and an extra 5 m s -1 was added for the roof locations due to the increase in wind speed with height above the ground level. All suspect and erroneous values were checked manually for consistency as it is possible to get higher wind gusts. Upon completion of this process, suspect values that were determined to be erroneous were changed and marked as so. All erroneous values were set to -9999, and the suspect values were flagged with the original values retained.</p><p>For the quality control of the temperature observations, most of it was done manually due to the changes of the different elevations of the sites, the local terrain, and the diurnal nature of the temperature cycle. With the local topography of the Willimantic Valley, it was possible to have a wide range of temperatures particularly at night with the ASOS site at the bottom of the valley, GAIL site on the valley slope, and the D3R site at the top of the valley. To provide some baseline for outlier observations, the mean bias and standard deviation were calculated between the AIO and ASOS, as was done for pressure. The bias plus two standard deviations was then used to determine two thresholds based on the site location; the GAIL site had a threshold set of a difference of 3 &#176;C, and the D3R had a threshold set for 4 &#176;C. To aid the manual quality control, weekly timeseries were produced using the AIO and ASOS observations. Figure <ref type="figure">14a</ref> shows an example timeseries with the actual observations over a two-week period. A second plot was made (not shown) of the difference between the two observations. Any periods with differences greater than the threshold were marked down to be investigated further. The GHCN site in Storrs, CT, was used to help determine did the diurnal cycle seem to match in terms of maximum and minimum temperatures. This was especially helpful if the sites had significant differences between them as the GHCN site is more closely located and at a closer elevation to the deployment sites.</p><p>One issue to note with the temperature observations was in periods after the computer connected to the AIO was restarted, the temperature would not produce negative temperature observations right away. Figure <ref type="figure">14a</ref> shows an example of this during 2023-2024, where the KIJD records temperatures well below freezing but the AIO is reporting temperatures of 5-10 &#176;C above freezing. The GHCN data was useful in verifying if the minimum temperature went below freezing or not on a given day. Periods when this occurred were manually identified and corrected by multiplying the AIO observations by a value of -1. Figure <ref type="figure">14b</ref> shows the absolute values of both temperature time series over the same time window, and the magnitude from both sets of observations match further confirming the issue with the AIO.</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head n="3.2.2">Micro Rain Radar (MRR)</head><p>Two varieties of MRRs, as described in Sect. 2.1.2, were used during the three deployment winters. The MRR-2 has a commonly used post-processing algorithm <ref type="bibr">(Maahn and Kollias, 2012)</ref>, which enhances the radar returns in frozen and mixed precipitation. This can be done to improve noise removal and add a de-aliasing component to provide more effective doppler velocity, reflectivity, and spectral width. This algorithm could be applied to the MRR-2 data in 2021-2022. However, the MRR-Pro, which was used in the following two deployments, does not have a generally accepted method to improve the radar returns in frozen or mixed precipitation akin to <ref type="bibr">Maahn and Kollias (2012)</ref>. New approaches to provide improved radar data are being developed, and here, the MRR-Pro data was processed through a new method detailed in <ref type="bibr">Williams (2023)</ref>. This approach has two versions: a calibrated version with velocity de-aliasing and recalculated moments and a disdrometer-calibrated version using a surface disdrometer, which follows a similar approach. The general approach removes range dependent noise and interference to find the lowest clutter free range gate and then applies a clutter filter adjustment to the clutter free range gates; from there, the spectra is reprocessed so velocity de-aliasing can be applied. Here is where the two methods differ as the calibration can be done compared to the originally observed spectra or against a surface disdrometer. In Sect. 4, we will present MRR-Pro data using this method with the calibration against the original spectra.</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head n="3.2.3">WxUAS</head><p>As detailed in Sect. 2.2.3, the StormTrooper WxUAS was deployed in two IOPs in March 2024. In these deployments, the WxUAS sampled two storms targeting mixed-phased precipitation often seen in the temperature and pressure gradients caused by the leading and trailing edges of these storms. To increase the chances of capturing the low-altitude phase transition periods while accommodating the limited flight crew size and the FAA crew resource management practices, the WxUAS deployment plan was based on an 8 h window. The intention was to start sampling one to two hours before the target storm effects were forecasted to occur and end it another one to two hours after. For the approaching storm cases, the targeted deployment would characterize the local conditions before the arrival of the storm's frontal boundaries, its arrival, and its transition across boundaries. For departing storms, this same strategy would characterize the local conditions during the storm, its departure, and the conditions behind its trailing edge. Considering the timescale in which precipitation phase changes had been reported by the surface observations in the past three winters and the time necessary to perform the pre-and post-flight safety procedures, the 8 h WxUAS data collection window was broken down into flights every 15 min. For these profiles, the WxUAS maintained a 3 m s -1 ascent rate while countering the winds to hold latitude and longitude over its take-off spot near the surface instruments. The targeted altitude for the WxUAS profiles was 900 m (approximately 3000 ft.). However, given the FAA's 120 m flight ceiling for commercial UAV flights, an altitude restriction waiver request was submitted. Even though the study site is in uncontrolled airspace with low traffic, particularly when hazardous winter storms are present, and the surrounding surface near the campus is sparsely populated, the FAA rejected the request without much justification. Nonetheless, leveraging nearby terrain elevation change and structures, the FAA Part 107 rules allowed the aircraft to reach up to 200 m (650 ft.) above ground level.</p><p>As discussed in Sect. 2.1.10, the WxUAS collected thousands of samples per variable, per profile. These samples produced high-resolution snapshots of the atmosphere, allowing researchers to tease out minor details in each variable for each flight. Figure <ref type="figure">15</ref> showed a few visualization examples of the data from the 9-10 March 2024, IOP. Figure <ref type="figure">15a</ref> exemplified the cloud droplet size distribution profile where the yaxis represents altitude, and the x-axis the total particle count for that altitude. In the cloud droplet size distribution profile, the outer color represents the distribution's mode at each altitude, that is, the cloud droplet size bin with the highest count while the inner color represents the largest particle size detected at that altitude. Here, the cloud droplet distribution explained what type of cloud region the WxUAS sampled. Near the surface, there were low particle counts and both the inside and outside coloring were blue, so the WxUAS was in clear air (Fig. <ref type="figure">15a</ref>); as the WxUAS approached 100 m, the inner color was red and the outer was blue, which indicates a mix of small and large particles see in cloud transition layer, so the WxUAS likely was starting to enter clouds. Figure <ref type="figure">15b</ref> presented a micro-range water content which indicates the vertical variations of water content in the cloud droplet size range. Figure <ref type="figure">15c</ref>, d represented the vertical profiles for temperature and relative humidity. Finally, Fig. <ref type="figure">15e</ref>, f highlighted time series for the radar derived approximated hydrometeor fall velocity and average reflectivity.</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head n="4">Highlights and Comparisons of the</head><p>Instrumentation during the 28 February 2023, Nor'easter 28 February was one of the largest snow events to occur in Southern New England during the winter of 2022-2023. Large portions of Connecticut, Rhode Island, and Western Massachusetts received 152.4-203.2 mm of snow with some areas reaching up to 254 mm. This winter storm was forced by a strong mid-level low pressure moving through the Ohio River Valley. There were double surface lows associated with the storm with the main low crossing the Great Lakes and the secondary low developing over the Southern Chesapeake Bay; the southern low-pressure center was closest to UConn as it moved along the coast of the Northeast United States. Significant moisture and forcing was in the region, along with a cold thermodynamic profile, which led to a long duration snow event. The UConn GPM GV sites were in the path of some of the heaviest snow rates, which occurred in the early overnight hours of the 28th. Official NWS Boston snow reports showed that observing locations around UConn received 177.8 mm of snow as of 07:00 a.m. EST (12:00 UTC). Despite the GAIL site's location being further west, the first precipitation occurred, as detected by the PARSIVEL 2 , at 01:41 UTC at the D3R site. Precipitation at the GAIL site started slightly later at 02:02 UTC. While one would expect the western site, GAIL, to experience precipitation first due to the storm's general west to east motion, this was not the case due to the higher elevation at the D3R site. There are many possible contributors to this occurring; with the relative humidity was slightly lower at the GAIL site and combining that with a small dry slot aloft, it is likely that the terrain and elevation change of 64 m caused snowflakes to sublimate before reaching the ground. The precipitation lasted until 22:18 UTC at the D3R site and ended shortly after at the GAIL site at 22:37 UTC. This event registered completely as snow to the PARSIVEL 2 's present weather algorithm at both observing locations.</p><p>At the GAIL site, the AIO (Fig. <ref type="figure">16a-e</ref>) registered significant changes before and during the storm. Fig. <ref type="figure">16a</ref> highlighted that during the second half of the snow event the air temperature and dewpoint had risen above 0 &#176;C as opposed to being below freezing during the heaviest snowfall; this was important to note because some precipitation phase classifications may misclassify this (e.g., MRMS precipitation phase would classify periods of snow during this event as periods of rain based on NEXRAD radar data, see <ref type="bibr">Zhang et al., 2016, Fig. 10</ref>). Figure <ref type="figure">16c</ref> showed surface pressure values dropping until between 07:00 and 08:00 UTC, where they remained relatively constant until after 18 UTC. The wind speed also dropped off significantly during the event after precipitation onset (Fig. <ref type="figure">16d</ref>, <ref type="figure">e</ref>). The winds were strongest around 8 m s -1 out of the south before onset of precipitation and dropped to about 3 m s -1 for most of the day; during the event, the wind shifted from the south to the east to the north as the low-pressure center moved northeast along the coast.</p><p>The PARSIVEL 2 and PIP both highlight the heaviest observed precipitation rates occurring between 03:00 and 07:00 UTC (Fig. <ref type="figure">16g</ref>, <ref type="figure">h</ref>). The PARSIVEL 2 snow rate (Fig. <ref type="figure">16g</ref>) estimates the amount of snowfall accumulation on the ground, hence why it is around a factor of 10 higher than the PIP precipitation rate (Fig. <ref type="figure">16h</ref>). Most of the precipitation occurred before 13:00 UTC and became more intermittent as temperatures increased (Fig. <ref type="figure">16a</ref>). Towards the end of the precipitation, larger snowflakes were forming, according to the maximum and mass-weighted mean diameter (Fig. <ref type="figure">16f</ref>). This could be an indication of particles aggregating together or riming more as surface temperatures rose. The PIP's equivalent bulk density (eDen) measurements (Fig. <ref type="figure">16i</ref>) also pointed to this fact (a typical 10 : 1 snow to liquid ratio appears as 0.1 g cm -3 ); after 16:00 UTC, eDen values rose to being significantly above 0.1 g cm -3 indicating heavier, rimed snowflakes. The PSD for both PARSIVEL 2 and PIP (Fig. <ref type="figure">16l</ref>, <ref type="figure">m</ref>) matched well and highlighted the heavier precipitation at the beginning of the snow event with more numerous larger particles and the largest amounts of the small particles as well. The PARSIVEL 2 tended to observe more small particles than the PIP later in the event, which was important to note when checking validity of the observations. If only one of these instruments was present at GAIL, this difference would not have been noticeable, reinforcing the need for redundant observations.</p><p>From a radar perspective, the MRR observations were consistent with the microphysical instrument observations. The highest reflectivity values occurred just after precipitation onset throughout the periods of heaviest precipitation (Fig. <ref type="figure">16j</ref>), and after 13:00 UTC, the precipitation was lighter, as previously seen. The Doppler velocities (Fig. <ref type="figure">16k</ref>) seen by the MRR also matched what would be expected for a snow event; this was because the velocities are slower, typical velocities for snow is 1-3 m s -1 , and there was no bright band, a sharp gradient in doppler velocities from slow to fast, to indicate a change from ice or snow to rain (Fig. <ref type="figure">16j</ref>). Greater downward velocities highlighted the heaviest snow rates at the beginning of the event, whereas smaller velocities dominated at low levels as the day progressed.</p><p>The MASC and SMAS were able to observe a variety of snowflake types due to the variety of conditions from the heavy snow bands to the lighter, denser snow later in the day. Graupel and aggregates were the most common types of snowflakes seen out of the six available geometric categories for classification (Fig. <ref type="figure">16n</ref>). In addition, 79.7 % of snowflakes classified had some form of riming during this event (Fig. <ref type="figure">16o</ref>). The rising temperatures towards the end of the event, visible in Fig. <ref type="figure">16a</ref>, raised wet-bulb temperatures above freezing, which would induce more riming or parti-cle aggregates as the edges of the ice or snow particles begin to melt increase the water content and density, confirmed by Fig. <ref type="figure">16i</ref>.</p><p>At the D3R site, the Pluvio told a similar precipitation story as at the GAIL site. The total accumulation from the Pluvio (Fig. <ref type="figure">17a</ref>) showed a steady increase in precipitation during the early morning hours matching the high precipitation rates. Later in the day, precipitation rates decreased with less overall accumulating precipitation. The PARSIVEL 2 showed similar trends as its counterpart at the GAIL site. The largest snowflakes at the D3R site occurred during the period of heaviest precipitation at the beginning of the event, and the D3R site tended to have slightly larger flakes than the GAIL Site throughout the second half of the event (Fig. <ref type="figure">17b</ref>). Despite the larger snowflakes at the D3R site later in the event, the snow rates did not differ significantly from the GAIL site (Fig. <ref type="figure">17c</ref>). It was not surprising that the two sites return similar observations due to the storm motion.</p><p>For the PIP, there were some notable differences between the two sites. The D3R PIP recorded some heavier precipitation rates particularly at the beginning of the event with rates reaching over 5 mm h -1 (Fig. <ref type="figure">17d</ref>), whereas the GAIL site's max precipitation rate was around 2.5 mm h -1 (Fig. <ref type="figure">16h</ref>). In addition, precipitation rates tended to be greater during the second half of the event due to the D3R site being slightly colder than the GAIL site because of the elevation difference. The colder temperatures may have prevented the snowflakes from beginning to rime or partially melt. This was reflected by the lower eDen values captured by the PIP (Fig. <ref type="figure">17e</ref>). The PSDs for both the PARSIVEL 2 and PIP (Fig. <ref type="figure">17f</ref>, <ref type="figure">g</ref>) are very similar during the entire event with the main exception being that the PIP saw more numerous very small and very large particles throughout most of the event.</p><p>As for the radar perspective, there were some interesting features visible from the three different radars at the D3R site. For the MRR, there was a similar pattern to what was exhibited at the GAIL site. The reflectivity from the MRR showed the highest values very early in the event (Fig. <ref type="figure">17h</ref>); it also had some higher values later in the event which were not seen at the GAIL site. The velocity values (Fig. <ref type="figure">17i</ref>) were also similar to the GAIL site, but they tended to be slightly slower across the entire event. The GAIL had velocities between 2-6 m s -1 between 1 and 3 km, whereas the D3R site velocity barely reached 3 m s -1 consistently. The ACHIEVE W-band radar is a nice complement to the MRR as it was generally pointed vertically for the duration of the event and scans to higher altitudes than the MRR. Between 00:00-09:00 UTC, the reflectivity peaks from the W-band (Fig. <ref type="figure">17j</ref>) generally matched the MRR, even if they had differing intensities. Because the W-band transmitter was more powerful, it better highlighted higher cloud structures during the beginning of event, although the W-band signal attenuates more as precipitation rates increase. Figure <ref type="figure">17j</ref> clearly shows the descending clouds approaching at the beginning, and it shows more intense periods of reflectivity associated with the tim-  <ref type="bibr">(Thant et al., 2023b)</ref>, and a machine learning classification <ref type="bibr">(Thant et al., 2023b)</ref> of riming degrees of snowflakes (o) based on all MASC and SMAS images during 28 February event. ing of heavier precipitation rates. As a complement, the Wband's LDR (Fig. <ref type="figure">17k</ref>) indicates departures from spherical particles. During the period of heaviest precipitation, pockets of larger LDR values (light blue shading) are most likely associated with larger aggregates and rimed snow particles.</p><p>The D3R radar performed PPI scans, as well as RHI scans over the GAIL site (315&#176;azimuth) and along the ER-2 aircraft track for this snow event. Fig. <ref type="figure">17l</ref> showed the reflectivity as measured by the Ka band radar on the D3R just after precipitation started. The highest reflectivity values at the lowest height levels indicated the heaviest precipitation was still further away (over 15 km) from the D3R site. The RHI scan NW of the radar showed light reflectivity values near the surface. This correlated with the precipitation not starting at the GAIL site until 02:02 UTC with most of this reflectivity tied to precipitation aloft. Figure <ref type="figure">17m</ref> highlighted the doppler velocity from the D3R; this indicated the precipitation was moving away from the D3R (purple and red colors) at the lowest level, which matches the AIO observations recording wind direction from the south southeast. For the most part, the doppler velocities were fairly uniform except at the lowest scanning levels, where there was variation of higher and lower speeds.</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head n="5">Data availability</head><p>The dataset presented in this article is available through NASA's Global Hydrometeorology Resource Center (GHRC) at the following campaign wide DOI: <ref type="url">https://doi.org/10.5067/GPMGVUCONN/DATA101</ref>  <ref type="bibr">(Cerrai et al., 2025)</ref>. Table <ref type="table">4</ref> below includes the individual dataset DOIs. All datasets processed by GHRC goes through standardized quality control checks and documentation as certified by the World Data System CoreTrustSeal <ref type="bibr">(Global Hydrometeorology Resource Center, 2024)</ref>. This data processing done at GHRC includes all steps of a dataset's life cycle from ingestion to processing, archiving, cataloguing, and documenting; their work also supports the distribution on final, published datasets. Accompanying each individual dataset is a user guide that provides additional information and metadata about the dataset including a description of the instrument, file naming conventions, the variables stored within each file, and the data format of the files.</p><p>For the external data, the KIJD ASOS observations, used to validate the AIO, can be downloaded from Iowa State University's Environmental Mesonet here: <ref type="url">https://mesonet</ref>. agron.iastate.edu/request/asos/1min.phtml (Iowa Environmental Mesonet, 2025). The Storrs, CT GHCN site data, used to validate the AIO, can be viewed here: <ref type="url">https:  //www.nrcc.cornell.edu/wxstation/nowdata.html</ref> (Northeast Regional Climate Center, 2025).</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head n="6">Conclusions</head><p>These GPM GV campaigns collected high-resolution datasets, with multiple complimentary instruments, across three winters in Connecticut, which produced 117 precipitation events of varying characteristics. These observations captured a wide variety of winter phenomena during their deployment and provided a rich set of data that can be used to conduct research in combination with multiple datasets including other observational field campaigns, numerical weather models, and satellites and other remote sensing instruments. The three years of deployments allowed for variations in instrumentation and site locations, which enhances the quality of observations collected; this diverse collection provides a variety of different ways to combine this observational dataset with other relevant datasets.</p><p>A natural pairing of this data is with the GPM satellite datasets as this was the main reason behind the field deployments. Combining satellite products with observations collected at UCONN allows for verification of satellite derived products and evaluation of algorithms in a variety of winter conditions. The wide variety of precipitation phases and weather conditions allow for a substantial evaluation of satellite products in the Northeast United States, where strong and impactful winter storms occur regularly.</p><p>Another logical pairing of data is with the IMPACTS field campaign data as there were three intensive observing periods where aircraft passed over or near UCONN sites or mobile sounding teams launched radiosondes near the sites. This dataset presents the opportunity to have a more expanded view of winter storms in the Northeast United States with a broader set of observations. This could include verification of numerical models, comparisons between aircraft, groundbased and remote sensing instruments, which can further refine understanding of winter precipitation processes.</p><p>The high-resolution nature of the data collected allows for it to be used in experiments seeking to better understand the physical processes of precipitation. Experiments involving intercomparisons between the ground-based instruments would allow for understanding of how the different types of instrumentation observe the same precipitation. Comparisons between the remote sensing instruments and in-situ instruments are also possible <ref type="bibr">(Billault-Roux et al., 2023;</ref><ref type="bibr">Shates et al., 2025)</ref>. Other research using this dataset has already been started including improving the PIP's algorithm to accurately identify rain versus snow, refining SWER(ZE) relationships <ref type="bibr">(Tokay et al., 2023;</ref><ref type="bibr">Inglis et al., 2024)</ref>, used to estimate snow water equivalent rates from radar reflectivity, between the PIP and MRR, understanding microphysical properties of snowfall <ref type="bibr">(Chang et al., 2024;</ref><ref type="bibr">King et al., 2024)</ref>, and intercomparison of precipitation accumulations from different instrumentation.</p><p>Additional high-resolution observations of winter phenomena are an essential part of furthering research to improve our prior knowledge and ability to forecast them. Future work to improve the measurements and understanding of precipitation across the globe relies on high-quality observations, like the ones collected in this NASA GPM GV campaign.</p><p>Table A1. Continued. Site Start Date Start End Date End Precipitation Event Total precipitation HH:MM HH:MM Minutes Phase (mm) GAIL 2/13/2022 08:49 2/14/2022 02:55 686 S 137.95 GAIL 2/18/2022 03:26 2/18/2022 13:10 496 R 21.54 GAIL 2/19/2022 17:42 2/19/2022 23:23 129 S 41.43 GAIL 2/22/2022 17:48 2/23/2022 03:33 450 R 26.39 GAIL 2/25/2022 07:22 2/25/2022 22:17 870 M 116.2 GAIL 3/2/2022 00:54 3/2/2022 04:33 174 R 3.96 GAIL 3/3/2022 06:48 3/3/2022 09:31 158 R 6.68 GAIL 3/6/2022 15:26 3/6/2022 16:28 56 R 3 GAIL 3/8/2022 03:09 3/8/2022 03:52 23 R 2.46 GAIL 3/9/2022 15:47 3/10/2022 01:17 563 S 178.06 GAIL 3/12/2022 06:25 3/13/2022 02:56 799 M 77.46 GAIL 3/16/2022 02:04 3/16/2022 03:52 109 R 1.35 GAIL 3/17/2022 13:05 3/17/2022 20:47 167 R 1.16 GAIL 3/20/2022 02:18 3/20/2022 05:03 36 R 2.12 GAIL 3/24/2022 01:50 3/24/2022 19:30 932 R 22.46 GAIL 3/25/2022 01:27 3/25/2022 07:58 333 R 26.76 GAIL 3/27/2022 01:46 3/27/2022 03:30 37 R 1.22 GAIL 4/1/2022 00:00 4/1/2022 12:05 481 R 31.7 GAIL 4/3/2022 16:08 4/4/2022 03:40 474 R 5.09 GAIL 4/6/2022 06:31 4/6/2022 18:03 665 R 10.77 GAIL 4/7/2022 20:22 4/8/2022 13:16 901 R 33.55 GAIL 4/9/2022 15:50 4/9/2022 22:19 94 R 5.03 GAIL 4/12/2022 13:07 4/12/2022 15:32 146 R 5.03 GAIL 4/14/2022 23:02 4/15/2022 03:36 107 R 4.98 GAIL 4/16/2022 22:35 4/17/2022 04:04 165 R 10.03 GAIL 4/19/2022 04:40 4/19/2022 18:22 650 R 48.75 GAIL 4/26/2022 13:40 4/27/2022 05:09 517 R 4.95 GAIL 12/22/2022 20:41 12/23/2022 21:20 1104 R 38.97 D3R 12/22/2022 21:56 12/23/2022 21:19 1053 R 40 GAIL 12/31/2022 18:58 1/1/2023 07:43 620 R 11.71 D3R 12/31/2022 19:00 1/1/2023 07:45 621 R 11.48 GAIL 1/3/2023 13:43 1/4/2023 09:47 758 R 16.87 D3R 1/3/2023 13:15 1/4/2023 09:48 778 R 17.21 GAIL 1/4/2023 23:32 1/5/2023 08:55 274 R 3.67 D3R 1/4/2023 23:34 1/5/2023 08:58 272 R 3.26 GAIL 1/6/2023 09:55 1/6/2023 20:51 555 R 8.25 D3R 1/6/2023 10:18 1/6/2023 21:03 534 R 7.33 GAIL 1/12/2023 12:21 1/13/2023 15:01 872 R 22.52 D3R 1/12/2023 12:26 1/13/2023 14:59 790 R 20.59 GAIL 1/14/2023 18:21 1/14/2023 20:22 115 S 5.48 GAIL 1/19/2023 18:05 1/20/2023 05:40 659 R 29.75 D3R 1/19/2023 18:09 1/20/2023 05:40 654 R 30.08 GAIL 1/22/2023 22:15 1/23/2023 23:53 1451 M 111.27 D3R 1/22/2023 22:22 1/23/2023 23:48 1424 M 126.49 GAIL 1/25/2023 20:36 1/25/2023 21:33 57 S 5.71 GAIL 1/25/2023 23:55 1/26/2023 10:38 642 R 38.59 D3R 1/25/2023 20:38 1/25/2023 21:33 55 S 4.61 D3R 1/25/2023 23:56 1/26/2023 10:40 644 R 34.46 GAIL 1/31/2023 09:22 1/31/2023 13:10 122 S 3.04 D3R 1/31/2023 09:22 1/31/2023 13:12 155 S 4.72 GAIL 2/8/2023 01:56 2/8/2023 03:57 97 R 3.94 D3R 2/8/2023 01:37 2/8/2023 03:57 96 R 3.84 GAIL 2/17/2023 14:15 2/17/2023 23:11 234 R 4.2 Table A1. Continued. Site Start Date Start End Date End Precipitation Event Total precipitation HH:MM HH:MM Minutes Phase (mm) Table A1. Continued. Site Start Date Start End Date End Precipitation Event Total precipitation HH:MM HH:MM Minutes Phase (mm) GAIL 4/3/2024 11:41 4/4/2024 18:17 1448 S 30.92 GAIL 4/11/2024 12:04 4/12/2024 19:56 890 R 28.37 GAIL 4/17/2024 23:49 4/18/2024 21:03 794 R 18.37 GAIL 4/20/2024 07:48 4/20/2024 14:18 234 R 5.83 GAIL 4/28/2024 01:43 4/28/2024 08:09 181 R 4.07 GAIL 4/28/2024 21:40 4/29/2024 03:15 155 R 2.83 GAIL 5/1/2024 02:10 5/1/2024 09:05 236 R 3.14 GAIL 5/5/2024 18:05 5/6/2024 07:28 650 R 24.14 GAIL 5/8/2024 11:48 5/8/2024 13:43 115 R 15.02 GAIL 5/9/2024 00:45 5/9/2024 00:56 12 R 4.05 GAIL 5/10/2024 00:45 5/10/2024 11:49 320 R 2.44 GAIL 5/15/2024 14:01 5/16/2024 16:50 1473 R 45.3 GAIL 5/18/2024 10:12 5/19/2024 02:18 597 R 9.9</p><p>Author contributions. BF wrote the original first draft, conducted some of the analysis, and produced most of the visualizations. DW and DC were responsible for the conceptualization, supervision, and funding of the field campaigns. AT, CH, AL, LB, VC, and GA provided input on the instrumentation used and deployed during the field campaigns. BF, AS, AC, CS, MB, EK, FJ, HT, and GA were responsible for installing and maintaining the instrumentation. AT, CP, AL, EK, HT, GA assisted with the development and production of some of the visualizations and analysis. BF, DW, AT, CH, AL, CP, VC, BN, GH, and DC all contributed to the writing and editing of the manuscript.</p><p>Competing interests. The contact author has declared that none of the authors has any competing interests.</p><p>Disclaimer. Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims made in the text, published maps, institutional affiliations, or any other geographical representation in this paper. While Copernicus Publications makes every effort to include appropriate place names, the final responsibility lies with the authors. Views expressed in the text are those of the authors and do not necessarily reflect the views of the publisher.</p><p>Acknowledgements. This field campaign was funded by NASA grant 80NSSC20K0577. We would like to thank the UConn Facilities staff for all their assistance in preparing the sites and helping maintain operations throughout the duration of the campaign, and NASA's Global Hydrometeorology Resource Center staff for their assistance with archiving the data. We would also like to thank George Huffman, GPM Project Scientist, and Will McCarty, GPM Program Scientist, for their support. The participation of V. Chandrasekar, EunYeol Kim, and Francesc Junyent was supported by NASA grant (80NSSC19K0877).</p><p>Financial support. This research has been supported by the Goddard Earth Sciences (grant no. 80NSSC20K0577, 80NSSC19K0877).</p><p>Review statement. This paper was edited by Graciela Raga and reviewed by two anonymous referees.</p></div><note xmlns="http://www.tei-c.org/ns/1.0" place="foot" xml:id="foot_0"><p>Earth Syst. Sci.Data, 17, 5783-5810, 2025   https://doi.org/10.5194/essd-17-5783-2025</p></note>
			<note xmlns="http://www.tei-c.org/ns/1.0" place="foot" xml:id="foot_1"><p>https://doi.org/10.5194/essd-17-5783-2025Earth Syst. Sci. Data, 17, 5783-5810, 2025</p></note>
		</body>
		</text>
</TEI>
