SEARCH · Search NASA

Results for “Snowfall”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 127 records · Page 7

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Permafrost thermal conditions are sensitive to shifts in snow timing

Abstract Changes in snow precipitation at high latitudes can significantly affect permafrost thermal conditions and thaw depth, potentially exposing more carbon-laden soil to microbial decomposition. A fully coupled process-based surface/subsurface thermal hydrology model with surface energy balance is used to analyze the impact of intra-annual variability in snow on permafrost thermal regime and the active layer thickness. In the four numerical scenarios considered, simulations were forced by the same meteorological data, except the snow precipitation, which was systematically altered to change timing of snowfall. The scenarios represent subtle shifts in snow timing, but the snow onset/melt days, the end of winter snowpack depth, and total annual snow precipitation are unchanged among scenarios. The simulations show that small shifts in the timing of snow accumulation can have significant effects on subsurface thermal conditions leading to active layer deepening and even talik formation when snowfall arrives earlier in the winter. The shifts in snow timing have a stronger impact on wetter regions, especially soil underneath small ponds, as compared to drained regions. This study highlights the importance of understanding potential changes in winter precipitation patterns for reliable projections of active-layer thickness in a changing Arctic climate.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

The Extraordinary March 2022 East Antarctica “Heat” Wave. Part II: Impacts on the Antarctic Ice Sheet

Abstract Between 15 and 19 March 2022, East Antarctica experienced an exceptional heat wave with widespread 30°–40°C temperature anomalies across the ice sheet. In Part I, we assessed the meteorological drivers that generated an intense atmospheric river (AR) that caused these record-shattering temperature anomalies. Here, we continue our large collaborative study by analyzing the widespread and diverse impacts driven by the AR landfall. These impacts included widespread rain and surface melt that was recorded along coastal areas, but this was outweighed by widespread high snowfall accumulations resulting in a largely positive surface mass balance contribution to the East Antarctic region. An analysis of the surface energy budget indicated that widespread downward longwave radiation anomalies caused by large cloud-liquid water contents along with some scattered solar radiation produced intense surface warming. Isotope measurements of the moisture were highly elevated, likely imprinting a strong signal for past climate reconstructions. The AR event attenuated cosmic ray measurements at Concordia, something previously never observed. Last, an extratropical cyclone west of the AR landfall likely triggered the final collapse of the critically unstable Conger Ice Shelf while further reducing an already record low sea ice extent. Significance Statement Using our diverse collective expertise, we explored the impacts from the March 2022 heat wave and atmospheric river across East Antarctica. One key takeaway is that the Antarctic cryosphere is highly sensitive to meteorological extremes originating from the midlatitudes and subtropics. Despite the large positive temperature anomalies driven from strong downward longwave radiation, this event led to huge amounts of snowfall across the Antarctic interior desert. The isotopes in this snow of warm airmass origin will likely be detectable in future ice cores and potentially distort past climate reconstructions. Even measurements of space activity were affected. Also, the swells generated from this storm helped to trigger the final collapse of an already critically unstable Conger Ice Shelf while further degrading sea ice coverage.

Wille, Jonathan D.↗

DOE OSTI2026

Seasonal Precipitation Classification during Surface Atmosphere Integrated Field Laboratory Campaign

The Surface Atmosphere Integrated Field Laboratory (SAIL) campaign, conducted from September 2021 to June 2023 in Crested Butte, Colorado, aimed to characterize precipitation processes in the Upper Colorado River Basin (UCRB). This increased observations of snowfall accumulation in this hydrologically significant watershed would be useful for quantitative precipitation estimates (QPE). Therefore, the Surface Quantitative Precipitation Estimate (SQUIRE) product was developed using the ARM-supported Colorado State University (CSU) X-band Precipitation Radar. Although SQUIRE will be only released for snowfall, by categorizing precipitation types, users can effectively utilize relevant datasets under diverse meteorological conditions. Moreover, the dataset facilitates validation of the QPE product and the analysis of seasonal variations in precipitation types at the surface. Hydrometeors classes are organized based on their phase and physical characteristics mapping the CSU (both winter Summer) and Py-ART classifications into four groups. 1. Liquid Precipitation: includs drizzle, rain, and large raindrops. 2. Frozen Snow and Ice : Pure Snow, combining ice crystals, aggregates, and vertically oriented ice structures. 3. Dense and Large frozen hydrometeors: including low- and high-density graupel and dry hail. 4.Melting: Wet Snow and Melting Hail, hydrometeors exhibiting both liquid and frozen characteristics.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2026

Xanthos-Lake Dataset

The Xanthos-Lake v1.0 dataset provides the input data, trained machine-learning models, and simulation outputs needed to characterize lake water balance, snow and ice conditions, and mixing-layer temperature within the Xanthos global hydrological modeling framework. The dataset supports lake representation across a wide range of lake sizes and hydroclimatic conditions by combining xLSIM, a basin-specific machine-learning emulator of lake snow, ice, ice-cover fraction, and mixing-layer temperature, with the Xanthos-Lake water-balance model. The archive contains NetCDF datasets used to train and evaluate xLSIM, trained model weights, processed meteorological and lake-property inputs, and basin- and lake-category-specific simulation outputs. These materials are organized into four primary data groups, described below. Snowice_model_inputs: Contains the NetCDF input data used to train xLSIM. The xLSIM machine-learning framework uses three lake-based datasets. The meteorological forcing dataset provides monthly relative humidity, specific humidity, surface wind speed, maximum and minimum air temperature, downward longwave and shortwave radiation, snowfall, surface air pressure, and total precipitation. Lake surface area is included as an additional static predictor. The target-state dataset provides lake ice thickness, snow depth, snow cover, and lake mixing-layer temperature, while a companion lake-surface dataset provides the lake ice-cover fraction. Before training, ice thickness and snow depth are converted from meters to centimeters, mixing-layer temperature is converted from kelvin to degrees Celsius and constrained to nonnegative values, and ice-cover fraction is converted from a fraction to a percentage. The predictor variables are normalized using statistics calculated across the selected lakes and time steps. Snowice_model_outputs: Contains the NetCDF outputs generated by xLSIM. For each basin, xLSIM produces a file containing observed and predicted lake-state variables for the training, validation, and testing periods. The modeled variables include lake ice thickness, snow depth, snow cover, mixing-layer temperature, and lake ice-cover fraction. For basins without a sufficiently persistent snow-and-ice signal, the emulator predicts only mixing-layer temperature. The outputs also include training and validation loss histories, the selected model configuration, identifiers of the lakes used in training, and SHAP-based feature-importance information at the global, lake, and seasonal-regime levels. The trained machine-learning model weights are provided separately within the dataset archive. Together, these files support model evaluation and subsequent coupling with the Xanthos-Lake water-balance framework. XanthosLAKES: Contains the NetCDF input data used by the Xanthos-Lake framework. Monthly meteorological inputs include relative and specific humidity, downward shortwave and longwave radiation, mean, maximum, and minimum air temperature, wind speed, precipitation, snowfall, and surface air pressure. Static lake-property datasets provide lake identifiers, geographic locations, surface area, volume, mean depth, elevation, drainage area, fetch, outlet-routing information, and associated Xanthos grid-cell attributes. Separate bathymetric datasets provide the coefficients of the area–depth and volume–depth relationships for each aggregated lake unit. GLEV-based records provide observed lake surface area and evaporation data used to initialize lake states, define reference conditions, and calibrate and evaluate the model. Xanthos-Lake Outputs: Contains the basin- and lake-category-specific NetCDF outputs generated by Xanthos-Lake. Monthly variables include lake surface area, storage volume, outlet discharge, evaporation rate, evaporation volume, lake–groundwater exchange, lake inflow, ice thickness, snow depth, snow-cover fraction, ice-cover fraction, and mixing-layer temperature. The files also contain lake-specific calibration and validation statistics, including normalized root-mean-square error, mean absolute error, Nash–Sutcliffe efficiency, Kling–Gupta efficiency, and percent bias. Stored calibrated and derived parameters include the weir discharge coefficient, fractional freeboard, groundwater exchange coefficient, reference water level, corresponding reference surface area and storage volume, weir-width adjustment factor, and the fraction of routed inflow entering the lake. Basin identifiers, lake category, simulation period, calibration and validation periods, and parameter-schema information are retained as NetCDF metadata.

Abeshu, Guta [Pacific Northwest National Laborator↗

NASA NTRS1978

Remote sensing applications to hydrologic modeling in the southern Sierra Nevada and portions of the San Joaquin Valley, volume 1

The author has identified the following significant results. Characteristics of LANDSAT MSS imagery present problems in using satellite radiation measurements to estimate the shortwave albedo of an alpine snow cover. Every 15 minute USGS quadrangle contains over 100,000 pixels which poses a computation problem if each pixel is to be evaluated individually. The sampling interval may be sufficiently great to mask some effects of terrain and vegetation on reflectance. Three frames of LANDSAT imagery are needed for complete coverage of the study area, yet less than one third of the area coverage from each frame covers an area of interest. Because of distortions inherent in the imagery, information regarding spacecraft altitude, attitude, and position must be statistically derived with respect to ground control points in the image whose geodetic locations are known. An inspection of shade points indicates that up to one third of the most heavily snow covered areas may saturate in bands 4 through 6. LANDSAT's 9 day repeat cycle is not optimum for snow cover reflectance modeling because the most pronounced changes in albedo occur most nearly following a new snowfall. Such a snowfall, occurring between overpasses, is inadequately represented by extrapolation from the previous overpasses.

Dozier, J.↗

NASA NTRS1986

An examination of vertical circulations associated with heavy snow events along the East Coast

Observational data are presented to show how a particular orientation of upper-level jet streaks (ULJS) associated with a diffluent trough near the East Coast of the U.S. can lead to heavy snowfalls. Surface fronts, precipitation and isobar maps, vertical motion vectors, and geostrophic wind streamlines preceding and during snowstorms on February 11, 1983, February 19, 1979, January 20, 1978, and January 13, 1964 are examined. The interaction of a direct vertical ageostrophic circulation within the confluent entrance of a ULJS over the northeastern U.S. or southeastern Canada and an indirect circulation within the exit region of a ULJS accompanying a trough nearing the East Coast are demonstrated to be conducive to heavy snowfall on the East Coast.

Kocin, Paul J.↗

NASA NTRS1990

Scatterers in Triton's atmosphere - Implications for the seasonal volatile cycle

Nitrogen and methane ices on the surface of Triton, Neptune's largest satellite, are exchanged between the summer and winter hemispheres on a seasonal time scale. Images of the satellite's sky obtained by the Voyager 2 spacecraft show the presence of several types of scattering materials that provide insights into this seasonal cycle of volatiles. Discrete clouds, probably composed of N2 ice particles, arise in regions of active sublimation. They are found chiefly poloward of 30 deg S in the southern, summer hemisphere. Haze particles, probably made of hydrocarbon ices, are present above most, but not all places. Recent snowfall may have occurred at low southern latitudes in places where they are absent. The latent heat released in the formation of the discrete clouds may have a major impact on the thermal balance of the lower atmosphere. Triton may have been less red at the time of the Voyager flyby than 12 years earlier due to recent N2 snowfall at a wide range of latitudes.

Pollack, James B.↗

NASA NTRS1989

A comparison of two exceptional snowstorms during the early 20th century

Two East-Coast snowstorms of the early twentieth century, which generated massive accumulations over portions of the Atlantic coastal plain, are described. The 'Christmas Storm' of December 25-26, 1909 produced snowfalls of greater than 25 cm over most of the Northeast Corridor from Maryland to southern New England, while the 'Knickerbocker Storm' of January 27-29, 1922 produced a total of 71 cm at Washington, DC. The storms are examined by constructing surface-weather analyses based on original data; surface observations for the continental United States at 0100 and 1300 UTC are obtained from the National Climatic Data Center (NCDC), ship reports from American vessels are taken from the NCDC archives, and the 1300 UTC data is supplemented by Canadian observations and European ship reports. In spite of the vast differences in the evolutions of the storms, it is found that both produced similarly massive snowfalls, therefore, it is concluded that a wide range of synoptic conditions can give rise to major snowstorms along the East Coast.

Weiss, Alan D.↗

NASA NTRS1992

Polarization Lidar Liquid Cloud Detection Algorithm for Winter Mountain Storms

We have collected an extensive polarization lidar dataset from elevated sites in the Tushar Mountains of Utah in support of winter storm cloud seeding research and experiments. Our truck-mounted ruby lidar collected zenith, dual-polarization lidar data through a roof window equipped with a wiper system to prevent snowfall accumulation. Lidar returns were collected at a rate of one shot every 1 to 5 min during declared storm periods over the 1985 and 1987 mid-Jan. to mid-Mar. Field seasons. The mid-barrier remote sensor field site was located at 2.57 km MSL. Of chief interest to weather modification efforts are the heights of supercooled liquid water (SLW) clouds, which must be known to assess their 'seedability' (i.e., temperature and height suitability for artificially increasing snowfall). We are currently re-examining out entire dataset to determine the climatological properties of SLW clouds in winter storms using an autonomous computer algorithm.

Sassen, Kenneth↗

NASA NTRS2004

Investigating Downscaling Methods and Evaluating Climate Models for Use in Estimating Regional Water Resources in Mountainous Regions under Changing Climatic Conditions

The purpose of this three-year study is to develop and evaluate techniques to estimate the range of potential hydrological impacts of climate change in mountainous areas. Three main objectives are set out in the proposal. (1) To develop and evaluate transfer functions to link tropospheric circulation to regional snowfall. (2) To evaluate a suite of General Circulation Models (GCMs) for use in estimating synoptic scale circulation and the resultant regional snowfall. And (3) to estimate the range of potential hydrological impacts of changing climate in the two case study areas: the Upper Colorado River basin, and the Catskill Mountains of southeastern New York State. Both regions provide water to large populations.

Frei, Allan↗

NASA NTRS2005

Melting in Martian Snowbanks

Precipitation as snow is an emerging paradigm for understanding water flow on Mars, which gracefully resolves many outstanding uncertainties in climatic and geomorphic interpretation. Snowfall does not require a powerful global greenhouse to effect global precipitation. It has long been assumed that global average temperatures greater than 273K are required to sustain liquid water at the surface via rainfall and runoff. Unfortunately, the best greenhouse models to date predict global mean surface temperatures early in Mars' history that differ little from today's, unless exceptional conditions are invoked. Snowfall however, can occur at temperatures less than 273K; all that is required is saturation of the atmosphere. At global temperatures lower than 273K, H2O would have been injected into the atmosphere by impacts and volcanic eruptions during the Noachian, and by obliquity-driven climate oscillations more recently. Snow cover can accumulate for a considerable period, and be available for melting during local spring and summer, unless sublimation rates are sufficient to remove the entire snowpack. We decided to explore the physics that controls the melting of snow in the high-latitude regions of Mars to understand the frequency and drainage of snowmelt in the high martian latitudes.

Zent, A. P.↗

NASA NTRS2006

GPM Constellation Reconfiguration and Mission Status

The Global Precipitation Measurement (GPM) Mission is an international satellite mission that uses combined active and passive remote sensing techniques to improve global precipitation measurements derived from dedicated and operational passive microwave sensors. GPM is a science mission with integrated applications goals to (1) advance the knowledge of precipitation physics and the global water cycle variability, and (2) improve weather, climate, and hydrological prediction capabilities through more accurate and frequent measurements of global precipitation and innovative application methods. The GPM Mission is currently a partnership between NASA and the Japan Aerospace Exploration Agency (JAXA), with opportunities for additional domestic and international partners in satellite constellation buildup and ground validation activities. The GPM concept is centered on employing a dualfrequency precipitation radar and a microwave radiometer with high-frequency capabilities on a core satellite to serve as a physics observatory and calibration standard to provide a consistent framework for unifying precipitation measurements from a heterogeneous constellation of passive microwave sensors. Building on the siccess of TRMM, GPM extends combined radadradiometer precipitation measurements into the mid and high latitudes, with new science foci on improved capabilities for light-rain and snowfall measurements, as well as more accurate precipitation retrievals over land. With recent studies indicating that AMSU-B rainfall estimates are comparable in quality to those derived from conically-scanning radiometers over land, it is envisioned that cross-track microwave sounders with high-frequency channels on operational satellites such as the National Polar-orbiting Operational Environmental Satellite System (NPOESS), NPOESS Preparatory Project (NPP), NOAA-N', and MetOp satellites can play a significant role in augmenting conically-scanning microwave radiometers to achieve better sampling and coverage over land. Plans to reconfigure the baseline GPM constellation to include cross-track microwave sounders over land and the status of other GPM mission elements such as snowfall algorithm development and ground validation plans will be presented

Hou, Arthur Y.↗

NASA NTRS2007

Antarctic Sea Ice Thickness and Snow-to-Ice Conversion from Atmospheric Reanalysis and Passive Microwave Snow Depth

Passive microwave snow depth, ice concentration, and ice motion estimates are combined with snowfall from the European Centre for Medium Range Weather Forecasting (ECMWF) reanalysis (ERA-40) from 1979-200 1 to estimate the prevalence of snow-to-ice conversion (snow-ice formation) on level sea ice in the Antarctic for April-October. Snow ice is ubiquitous in all regions throughout the growth season. Calculated snow- ice thicknesses fall within the range of estimates from ice core analysis for most regions. However, uncertainties in both this analysis and in situ data limit the usefulness of snow depth and snow-ice production to evaluate the accuracy of ERA-40 snowfall. The East Antarctic is an exception, where calculated snow-ice production exceeds observed ice thickness over wide areas, suggesting that ERA-40 precipitation is too high there. Snow-ice thickness variability is strongly controlled not just by snow accumulation rates, but also by ice divergence. Surprisingly, snow-ice production is largely independent of snow depth, indicating that the latter may be a poor indicator of total snow accumulation. Using the presence of snow-ice formation as a proxy indicator for near-zero freeboard, we examine the possibility of estimating level ice thickness from satellite snow depths. A best estimate for the mean level ice thickness in September is 53 cm, comparing well with 51 cm from ship-based observations. The error is estimated to be 10-20 cm, which is similar to the observed interannual and regional variability. Nevertheless, this is comparable to expected errors for ice thickness determined by satellite altimeters. Improvement in satellite snow depth retrievals would benefit both of these methods.

Markus, Thorsten↗

NASA NTRS2009

NASA Global Precipitation Mission Ground Validation Implementation

The Global Precipitation Mission (GPM; core-satellite launch 2013) will provide Ka/Ku-band dual-frequency precipitation radar (DPR) and accompanying passive microwave radiometer-diagnosed precipitation estimates over a latitude range of 65 N to 65 S. The extended latitudinal domain of GPM coverage combined with requirements to detect (and in the case of liquid, estimate) liquid and frozen precipitation rates for values ranging from several hundred to just a few tenths of a millimeter per hour present new challenges to the development of physically-based satellite precipitation retrieval algorithms. On regional scales select national and international resources such as existing calibrated radar and rain gauge networks can provide basic datasets that enable direct statistical validation of GPM core-satellite reflectivitys and core/constellation rain rate measurements. Near-term planned field campaign involvements include Finland/Baltic Sea (fall 2010; joint CloudSat,GPM, and European study of precipitation in low-altitude melting layers and snowfall in the vicinity of the Helsinki testbed), central Oklahoma (spring 2011; joint with DOE ARM- precipitation retrievals over a mid-latitude continental land surface), and the Great Lakes region (winter 2011-12, snowfall retrieval).

Petersen, Walter A.↗

NASA NTRS2011

Performance Tests of Snow-Related Variables Over the Tibetan Plateau and Himalayas Using a New Version of NASA GEOS-5 Land Surface Model that Includes the Snow Darkening Effect

The snow darkening effect, i.e. the reduction of snow albedo, is caused by absorption of solar radiation by absorbing aerosols (dust, black carbon, and organic carbon) deposited on the snow surface. This process is probably important over Himalayan and Tibetan glaciers due to the transport of highly polluted Atmospheric Brown Cloud (ABC) from the Indo-Gangetic Plain (IGP). This effect has been incorporated into the NASA Goddard Earth Observing System model, version 5 (GEOS-5) atmospheric transport model. The Catchment land surface model (LSM) used in GEOS-5 considers 3 snow layers. Code was developed to track the mass concentration of aerosols in the three layers, taking into account such processes as the flushing of the compounds as liquid water percolates through the snowpack. In GEOS-5, aerosol emissions, transports, and depositions are well simulated in the Goddard Chemistry Aerosol Radiation and Transport (GO CART) module; we recently made the connection between GOCART and the GEOS-5 system fitted with the revised LSM. Preliminary simulations were performed with this new system in "replay" mode (i.e., with atmospheric dynamics guided by reanalysis) at 2x2.5 degree horizontal resolution, covering the period 1 November 2005 - 31 December 2009; we consider the final three years of simulation here. The three simulations used the following variants of the LSM: (1) the original Catchment LSM with a fixed fresh snowfall density of 150 kg m-3 ; (2) the LSM fitted with the new snow albedo code, used here without aerosol deposition but with changes in density formulation and melting water effect on snow specific surface area, (3) the LSM fitted with the new snow albedo code as same as (2) but with fixed aerosol deposition rates (computed from GOCART values averaged over the Tibetan Plateau domain [Ion.: 60-120E; lat.: 20-50N] during March-May 2008) applied to all grid points at every time step. For (2) and (3), the same setting on the fresh snowfall density as in (1) was used.

Yasunari, Tppei J.↗

NASA NTRS2012

Satellite and Surface Perspectives of Snow Extent in the Southern Appalachian Mountains

Assessing snow cover patterns in mountain regions remains a challenge for a variety of reasons. Topography (e.g., elevation, exposure, aspect, and slope) strongly influences snowfall accumulation and subsequent ablation processes, leading to pronounced spatial variability of snow cover. In-situ observations are typically limited to open areas at lower elevations (<1000 m). In this paper, we use several products from the Moderate Resolution Imaging Spectroradiometer (MODIS) to assess snow cover extent in the Southern Appalachian Mountains (SAM). MODIS daily snow cover maps and true color imagery are analyzed after selected snow events (e.g., Gulf/Atlantic Lows, Alberta Clippers, and Northwest Upslope Flow) from 2006 to 2012 to assess the spatial patterns of snowfall across the SAM. For each event, we calculate snow cover area across the SAM using MODIS data and compare with the Interactive Multi-sensor Snow and ice mapping system (IMS) and available in-situ observations. Results indicate that Gulf/Atlantic Lows are typically responsible for greater snow extent across the entire SAM region due to intensified cyclogenesis associated with these events. Northwest Upslope Flow events result in snow cover extent that is limited to higher elevations (>1000 m) across the SAM, but also more pronounced along NW aspects. Despite some limitations related to the presence of ephemeral snow or cloud cover immediately after each event, we conclude that MODIS products are useful for assessing the spatial variability of snow cover in heavily forested mountain regions such as the SAM.

Sugg, Johnathan W.↗

NASA NTRS2017

So, How Much of the Earth's Surface Is Covered by Rain Gauges?

The measurement of global precipitation, both rainfall and snowfall, is critical to a wide range of users and applications. Rain gauges are indispensable in the measurement of precipitation, remaining the de facto standard for precipitation information across Earths surface for hydrometeorological purposes. However, their distribution across the globe is limited: over land their distribution and density is variable, while over oceans very few gauges exist and where measurements are made, they may not adequately reflect the rainfall amounts of the broader area. Critically, the number of gauges available, or appropriate for a particular study, varies greatly across the Earth owing to temporal sampling resolutions, periods of operation, data latency, and data access. Numbers of gauges range from a few thousand available in near real time to about 100,000 for all official gauges, and to possibly hundreds of thousands if all possible gauges are included. Gauges routinely used in the generation of global precipitation products cover an equivalent area of between about 250 and 3,000 sq m. For comparison, the center circle of a soccer pitch or tennis court is about 260 sq m. Although each gauge should represent more than just the gauge orifice, autocorrelation distances of precipitation vary greatly with regime and the integration period. Assuming each Global Precipitation Climatology Centre (GPCC) available gauge is independent and represents a surrounding area of 5-km radius, this represents only about 1% of Earths surface. The situation is further confounded for snowfall, which has a greater measurement uncertainty.

Rain gauge↗

NASA NTRS2017

Contribution of Lake-Effect Snow to the Catskill Mountains Snowpack

Meltwater from snow that falls in the Catskill Mountains in southern New York contributes to reservoirs that supply drinking water to approximately nine million people in New York City. Using the NOAA National Ice Centers Interactive Multisensor Snow and Ice Mapping System (IMS) 4km snow maps, we have identified at least 32 lake-effect (LE) storms emanating from Lake Erie andor Lake Ontario that deposited snow in the CatskillDelaware Watershed in the Catskill Mountains of southern New York State between 2004 and 2017. This represents a large underestimate of the contribution of LE snow to the Catskills snowpack because many of the LE snowstorms are not visible in the IMS snow maps when they travel over snow-covered terrain. Most of the LE snowstorms that we identified originate from Lake Ontario but quite a few originate from both Erie and Ontario, and a few from Lake Erie alone. Using satellite, meteorological and reanalysis data we identify conditions that contributed to LE snowfall in the Catskills. Clear skies following some of the storms permitted measurement of the extent of snow cover in the watershed using multiple satellite sensors. IMS maps tend to overestimate the extent of snow compared to MODerate resolution Imaging Spectroradiometer (MODIS) and Landsat-derived snow-cover extent maps. Using this combination of satellite and meteorological data, we can begin to quantify the important contribution of LE snow to the Catskills Mountain snowpack. Changes that are predicted in LE snowfall from the Great Lakes could impact the distribution of rain vs snow in the Catskills which may affect future reservoir operations in the NYC Water Supply System.

Lake Ontario↗