Mapping Snow Grain Size over Greenland from MODIS

This paper presents a new automatic algorithm to derive optical snow grain size (SGS) at 1 km resolution using Moderate Resolution Imaging Spectroradiometer (MODIS) measurements. Differently from previous approaches, snow grains are not assumed to be spherical but a fractal approach is used to account for their irregular shape. The retrieval is conceptually based on an analytical asymptotic radiative transfer model which predicts spectral bidirectional snow reflectance as a function of the grain size and ice absorption. The analytical form of solution leads to an explicit and fast retrieval algorithm. The time series analysis of derived SGS shows a good sensitivity to snow metamorphism, including melting and snow precipitation events. Preprocessing is performed by a Multi-Angle Implementation of Atmospheric Correction (MAIAC) algorithm, which includes gridding MODIS data to 1 km resolution, water vapor retrieval, cloud masking and an atmospheric correction. MAIAC cloud mask (CM) is a new algorithm based on a time series of gridded MODIS measurements and an image-based rather than pixel-based processing. Extensive processing of MODIS TERRA data over Greenland shows a robust performance of CM algorithm in discrimination of clouds over bright snow and ice. As part of the validation analysis, SGS derived from MODIS over selected sites in 2004 was compared to the microwave brightness temperature measurements of SSM\I radiometer, which is sensitive to the amount of liquid water in the snowpack. The comparison showed a good qualitative agreement, with both datasets detecting two main periods of snowmelt. Additionally, MODIS SGS was compared with predictions of the snow model CROCUS driven by measurements of the automatic whether stations of the Greenland Climate Network. We found that CROCUS grain size is on average a factor of two larger than MODIS-derived SGS. Overall, the agreement between CROCUS and MODIS results was satisfactory, in particular before and during the first melting period in mid-June. Following detailed time series analysis of SGS for four permanent sites, the paper presents SGS maps over the Greenland ice sheet for the March-September period of 2004

Aerosol Remote Sensing in Polar Regions

Multi-year sets of ground-based sun-photometer measurements conducted at 12 Arctic sites and 9 Antarctic sites were examined to determine daily mean values of aerosol optical thickness tau(lambda) at visible and near-infrared wavelengths, from which best-fit values of ngstrm's exponent alpha were calculated. Analyzing these data, the monthly mean values of tau(0.50 micrometers) and alpha and the relative frequency histograms of the daily mean values of both parameters were determined for winter-spring and summer-autumn in the Arctic and for austral summer in Antarctica. The Arctic and Antarctic covariance plots of the seasonal median values of alpha versus tau(0.50 micrometers) showed: (i) a considerable increase in tau(0.50 micrometers) for the Arctic aerosol from summer to winter-spring, without marked changes in alpha; and (ii) a marked increase in tau(0.50 micrometer) passing from the Antarctic Plateau to coastal sites, whereas alpha decreased considerably due to the larger fraction of sea-salt aerosol. Good agreement was found when comparing ground-based sun-photometer measurements of tau(lambda) and alpha at Arctic and Antarctic coastal sites with Microtops measurements conducted during numerous AERONET/MAN cruises from 2006 to 2013 in three Arctic Ocean sectors and in coastal and off-shore regions of the Southern Atlantic, Pacific, and Indian Oceans, and the Antarctic Peninsula. Lidar measurements were also examined to characterize vertical profiles of the aerosol backscattering coefficient measured throughout the year at Ny-lesund. Satellite-based MODIS, MISR, and AATSR retrievals of tau(lambda) over large parts of the oceanic polar regions during spring and summer were in close agreement with ship-borne and coastal ground-based sun-photometer measurements. An overview of the chemical composition of mode particles is also presented, based on in-situ measurements at Arctic and Antarctic sites. Fourteen log-normal aerosol number size-distributions were defined to represent the average features of nuclei, accumulation and coarse mode particles for Arctic haze, summer background aerosol, Asian dust and boreal forest fire smoke, and for various background austral summer aerosol types at coastal and high-altitude Antarctic sites. The main columnar aerosol optical characteristics were determined for all 14 particle modes, based on in-situ measurements of the scattering and absorption coefficients. Diurnally averaged direct aerosol-induced radiative forcing and efficiency were calculated for a set of multimodal aerosol extinction models, using various Bidirectional Reflectance Distribution Function models over vegetation-covered, oceanic and snow-covered surfaces. These gave a reliable measure of the pronounced effects of aerosols on the radiation balance of the surface-atmosphere system over polar regions

Analysis of the Discrete Theory of Radiative Transfer in the Coupled "Ocean-Atmosphere" System: Current Status, Problems and Development Prospects

In this paper, we analyze the current state of the discrete theory of radiative transfer. One-dimensional, three-dimensional and stochastic radiative transfer models are considered. It is shown that the discrete theory provides a unique solution to the one-dimensional radiative transfer equation. All approximate solution techniques based on the discrete ordinate formalism can be derived based on the synthetic iterations, the small-angle approximation, and the matrix operator method. The possible directions for the perspective development of radiative transfer are outlined

The Determination of the Snow Optical Grain Diameter and Snowmelt Area on the Greenland Ice Sheet Using Spaceborne Optical Observations

The optical diameter of the surface snow grains impacts the amount of energy absorbed by the surface and therefore the onset and magnitude of surface melt. Snow grains respond to surface heating through grain metamorphism and growth. During melt, liquid water between the grains markedly increases the optical grain size, as wet snow grain clusters are optically equivalent to large grains. We present daily surface snow grain optical diameters (dopt) retrieved from the Greenland ice sheet at 1 km resolution for 2017–2019 using observations from Ocean and Land Colour Instrument (OLCI) onboard Sentinel-3A. The retrieved dopt are evaluated against 3 years of in situ measurements in Northeast Greenland. We show that higher dopt are indicative of surface melt as calculated from meteorological measurements at four PROMICE automatic weather stations. We deduce a threshold value of 0.64 mm in dopt allowing categorization of the days either as melting or nonmelting. We apply this simple melt detection technique in Northeast Greenland and compare the derived melting areas with the conventional passive microwave MEaSUREs melt flag for June 2019. The two flags show generally consistent evolution of the melt extent although we highlight areas where large grain diameters are strong indicators of melt but are missed by the MEaSUREs melt flag. While spatial resolution of the optical grain diameter-based melt flag is higher than passive microwave, it is hampered by clouds. Our retrieval remains suitable to study melt at a local to regional scales and could be in the future combined with passive microwave melt flags for increased coveragepublishedVersio

Retrieval of snow properties from the Sentinel-3 Ocean and Land Colour Instrument

The Sentinel Application Platform (SNAP) architecture facilitates Earth Observation data processing. In this work, we present results from a new Snow Processor for SNAP. We also describe physical principles behind the developed snow property retrieval technique based on the analysis of Ocean and Land Colour Instrument (OLCI) onboard Sentinel-3A/B measurements over clean and polluted snow fields. Using OLCI spectral reflectance measurements in the range 400–1020 nm, we derived important snow properties such as spectral and broadband albedo, snow specific surface area, snow extent and grain size on a spatial grid of 300 m. The algorithm also incorporated cloud screening and atmospheric correction procedures over snow surfaces. We present validation results using ground measurements from Antarctica, the Greenland ice sheet and the French Alps. We find the spectral albedo retrieved with accuracy of better than 3% on average, making our retrievals sufficient for a variety of applications. Broadband albedo is retrieved with the average accuracy of about 5% over snow. Therefore, the uncertainties of satellite retrievals are close to experimental errors of ground measurements. The retrieved surface grain size shows good agreement with ground observations. Snow specific surface area observations are also consistent with our OLCI retrievals. We present snow albedo and grain size mapping over the inland ice sheet of Greenland for areas including dry snow, melted/melting snow and impurity rich bare ice. The algorithm can be applied to OLCI Sentinel-3 measurements providing an opportunity for creation of long-term snow property records essential for climate monitoring and data assimilation studies—especially in the Arctic region, where we face rapid environmental changes including reduction of snow/ice extent and, therefore, planetary albedo.publishedVersio