Remote sensing of snow and ice / La télédétection de neige et de glace

Abstract

Monitoring of snow and ice on the Earth's surface will require increasing use of satellite remote sensing techniques. These techniques are evolving rapidly. Active and passive sensors operating in the visible, near infrared, thermal infrared, and microwave wavelengths are described in regard to general applications and in regard to specific USA or USSR satellites. Meteorological satellites (frequent images of relatively crude resolution) and Earth resources satellites such as Landsat (less frequent images of higher resolution) have been used to monitor the areal extent of seasonal snow, but problems exist with cloud cover or dense forest canopies. Snow mass (water equivalent) can be measured from a low-flying aircraft using natural radioactivity, but cannot yet be measured from satellite altitudes. A combination of active and passive microwave sensors may permit this kind of measurement, but not until more is known about radiation scattering in snow. Satellite observations are very useful in glacier inventories, correcting maps of glacier extent, estimating certain mass balance parameters, and monitoring calving or surging glaciers. Ground ice is virtually impossible to monitor from satellites; ice on rivers and lakes can be monitored only with very high-resolution sensors. Microwave sensors, due to their all-weather capability (the ability to see through clouds) provide exciting data on sea ice distribution. Analysis of digital tapes of satellite data requires the archiving and scanning of huge amounts of data. Simple methods for extracting quantitative data from satellite images are described.     

Snow distribution from SSM/I and its relationships to the hydroclimatology of the Mackenzie River Basin,Canada

The spatial and temporal distribution of snow cover extent (SCE) and snow water equivalent (SWE) play vital roles in the hydrology of northern watersheds. We apply remotely sensed Special Sensor Microwave Imager (SSM/I) data from 1988 to 2007 to explore the relationships between snow distribution and the hydroclimatology of the Mackenzie River Basin (MRB) of Canada and its major sub-basins. The Environment Canada (EC) algorithm is adopted to retrieve the SWE from SSM/I data. Moderate Resolution Imaging Spectroradiometer (MODIS) 8-day maximum snow cover extent products (MOD10A2) are used to estimate the different thresholds of retrieved SWE from SSM/I to classify the land cover as snow or no snow for various sub-basins in the MRB. The sub-basins have varying topography and hence different thresholds that range from 10 mm to 30 mm SWE. The accuracy of snow cover mapping based on the combination of several thresholds for the different sub-basins reaches ≈ 90%. The northern basins are found to have stronger linear relationships between the date on which snow cover fraction (SCF) reaches 50% or when SWE reaches 50% and mean air temperatures, than the southern basins. Correlation coefficients between SCF, SWE, and hydroclimatological variables show the new SCF products from SSM/I perform better than SWE from SSM/I to analyze the relationships with the regional hydroclimatology. Statistical models relating SCF and SWE to runoff indicate that the SCF and SWE from EC algorithms are able to predict the discharge in the early snow ablation seasons in northern watersheds.     

Pixel‐scale evaluation of SSM/I sea‐ice algorithms in the marginal ice zone during early fall freeze‐up

Observed reduction in recent sea ice areal extent and thickness has focused attention on the fact that the Arctic marine system appears to be responding to global‐scale climate variability and change. Passive microwave remote‐sensing data are the primary source underpinning these reports, yet problems remain in geophysical inversion of information on ice type and concentration. Uncertainty in sea‐ice concentration (SIC) retrievals is highest in the summer and fall, when water occurs in liquid phase within the snow–sea‐ice system. Of particular scientific interest is the timing and rate of new ice formation due to the control that this form of sea ice has on mass, energy and gas fluxes across the ocean–sea‐ice–atmosphere interface. In this paper we examine the critical fall freeze‐up period using in situ data from a ship‐based and aerial survey programme known as the Canadian Arctic Shelf Exchange study combined with microwave and optical Earth observations data. Results show that: (1) the overall physical conditions observed from aerial survey photography were well matched with coincident moderate‐resolution imaging spectroradiometer data and Radarsat ScanSAR imagery; (2) the shortwave albedo was linearly related to old ice concentration derived from survey photography; (3) the three SSM/I SIC algorithms (NASA Team (NT), NASA Team 2 (NT2), and Bootstrap (BT)) showed considerable discrepancies in pixel‐scale comparison with the Radarsat ScanSAR SICs well calibrated by the aerial survey data. The major causes of the discrepancies are attributed to (1) the inherent inability to detect the new thin ice in the NT and BT algorithms, (2) mismatches of the thin‐ice tie point of the NT2 algorithm, and (3) sub‐pixel ambiguity between the thin ice and the mixture of open water and sea ice. These results suggest the need for finer resolution of passive microwave sensors, such as AMSR‐E, to improve the precision of the SSM/I SIC algorithms in the marginal ice zone during early fall freeze‐up. Copyright © 2006 John Wiley & Sons, Ltd.     

Passive microwave remote sensing of the historic February 2010 snowstorms in the Middle Atlantic region of the USA

The snowfall in the Baltimore/Washington metropolitan area during the winter of 2009/2010 was unprecedented and caused serious snow‐related disruptions. In February 2010, snowfall totals approached 2 m, and because maximum temperatures were consistently below normal, snow remained on the ground the entire month. One of the biggest contributing factors to the unusually severe winter weather in 2009/2010, throughout much of the middle latitudes, was the Arctic Oscillation. Unusually high pressure at high latitudes and low pressure at middle latitudes forced a persistent exchange of mass from north to south. In this investigation, a concerted effort was made to link remotely sensed falling snow observations to remotely sensed snow cover and snowpack observations in the Baltimore/Washington area. Specifically, the Advanced Microwave Scanning Radiometer onboard the Aqua satellite was used to assess snow water equivalent, and the Advanced Microwave Sounding Unit‐B and Microwave Humidity Sounder were employed to detect falling snow. Advanced Microwave Scanning Radiometer passive microwave signatures in this study are related to both snow on the ground and surface ice layers. In regard to falling snow, signatures indicative of snowfall can be observed in high frequency brightness temperatures of Advanced Microwave Sounding Unit‐B and Microwave Humidity Sounder. Indeed, retrievals show an increase in snow water equivalent after the detection of falling snow. Yet, this work also shows that falling snow intensity and/or the presence of liquid water clouds impacts the ability to reliably detect snow water equivalent. Moreover, changes in the condition of the snowpack, especially in the surface features, negatively affect retrieval performance. Copyright © 2011. This article is a U.S. Government work and is in the public domain in the USA.     

Overview of Overland Satellite Rainfall Estimation for Hydro-Meteorological Applications

This paper reviews current techniques on rainfall estimation from satellite sensor observations. The sensors considered in this study are the Precipitation Radar (PR) and radiometer (TMI) onboard TRMM (Tropical Rainfall Measuring Missio) satellite, the Special Sensor Microwave/Imager (SSM/I) onboard Defense Meteorological Satellite Program (DMSP) platforms, and infrared (IR) sensors onboard geostationary satellites. We present the physical basis and mathematical formulation of a newly developed combined radar-radiometer (PR/TMI) retrieval for TRMM and its application for overland rain estimation. Subsequently we discuss the current state-of-the-art in overland passive microwave (TMI and SSM/I) rain estimation techniques, and outstanding issues associated with the inverse problem. The significance of lightning information in advancing high-frequency rainfall estimation from passive microwave-calibrated IR retrieval techniques is discussed on the basis of newly developed techniques. Finally, current approaches are presented on merging the infrequent passive microwave-based rainfall estimates with the high-frequency, but lower accuracy, rainfall fields derived from proxy parameters (e.g., lightning and IR). The paper provides useful insights on satellite rainfall estimation and discusses issues and applications.     

II. Snow hydrology processes and remote sensing

In the last 20 years remote sensing research has led to significant progress in monitoring and measuring certain snow hydrology processes. Snow distribution in a drainage basin can be adequately assessed by visible sensors. Although there are still some interpretation problems, the NOAA-AVHRR sensor can provide frequent views of the areal snow cover in a basin, and snow cover maps are produced operationally by the National Weather Service on about 3000 drainage basins in North America. Measurement of snow accumulation or snow water equivalent with microwave remote sensing has great potential because of the capabilities for depth penetration, all-weather observation and night-time viewing. Several critical areas of research remain, namely, the acquisition of snow grain size information for input to microwave models and improvement in passive microwave resolution from space. Methods that combine both airborne gamma ray and visible satellite remote sensing of the snowpack with field measurements also hold promise for determining areal snow water equivalent. Some remote sensing techniques can also be used to detect different stages of snow metamorphism. Various aspects of snowpack ripening can be detected using microwave and thermal infra-red capabilities. The capabilities for measurement of snow albedo and surface temperature have direct application in both snow metamorphism and snowpack energy balance studies. The potentially most profitable research area here is the study of the bidirectional reflectance distribution function to improve snow albedo measurements. Most of the remote sensing capabilities in snow hydrology have been developed for improving snowmelt-run-off forecasting. Most applications have used the input of snow cover extent to deterministic models, both of the degree day and energy balance types. Snowmelt-run-off forecasts using satellite derived snow cover depletion curves and the models have been successfully made. As the extraction of additional snow cover characteristics becomes possible, remote sensing will have an even greater impact on snow hydrology. Important remote sensing capabilities will become available in the next 20 years through space platform observing systems that will improve our capability to observe the snowpack on an operational basis.     

Spatial variability of snowmelt timing from AMSR‐E and SSM/I passive microwave sensors,Pelly River,Yukon Territory,Canada

Spring snow melt run‐off in high latitude and snow‐dominated drainage basins is generally the most significant annual hydrological event. Melt timing, duration, and flow magnitude are highly variable and influence regional climate, geomorphology, and hydrology. Arctic and sub‐arctic regions have sparse long‐term ground observations and these snow‐dominated hydrologic regimes are sensitive to the rapidly warming climate trends that characterize much of the northern latitudes. Passive microwave brightness temperatures are sensitive to changes in the liquid water content of the snow pack and make it possible to detect incipient melt, diurnal melt‐refreeze cycles, and the approximate end of snow cover on the ground over large regions. Special Sensor Microwave Imager (SSM/I) and Advanced Microwave Scanning Radiometer for EOS (AMSR‐E) passive microwave brightness temperatures (T_b) and diurnal amplitude variations (DAV) are used to investigate the spatial variability of snowmelt onset timing (in two stages, ‘DAV onset’ and ‘melt onset’) and duration for a complex sub‐arctic landscape during 2005. The satellites are sensitive to small percentages of liquid water, and therefore represent ‘incipient melt’, a condition somewhat earlier than a traditional definition of a melting snowpack. Incipient melt dates and duration are compared to topography, land cover, and hydrology to investigate the strength and significance of melt timing in heterogeneous landscapes in the Pelly River, a major tributary to the Yukon River. Microwave‐derived melt onset in this region in 2005 occurred from late February to late April. Upland areas melt 1–2 weeks later than lowland areas and have shorter transition periods. Melt timing and duration appear to be influenced by pixel elevation, aspect, and uniformity as well as other factors such as weather and snow mass distribution. The end of the transition season is uniform across sensors and across the basin in spite of a wide variety of pixel characteristics. Copyright © 2007 John Wiley & Sons, Ltd.     

Improvement of snow depth retrieval for FY3B-MWRI in China

The primary objective of this work is to develop an operational snow depth retrieval algorithm for the FengYun3B Microwave Radiation Imager(FY3B-MWRI)in China.Based on 7-year(2002–2009)observations of brightness temperature by the Advanced Microwave Scanning Radiometer-EOS(AMSR-E)and snow depth from Chinese meteorological stations,we develop a semi-empirical snow depth retrieval algorithm.When its land cover fraction is larger than 85%,we regard a pixel as pure at the satellite passive microwave remote-sensing scale.A 1-km resolution land use/land cover(LULC)map from the Data Center for Resources and Environmental Sciences,Chinese Academy of Sciences,is used to determine fractions of four main land cover types(grass,farmland,bare soil,and forest).Land cover sensitivity snow depth retrieval algorithms are initially developed using AMSR-E brightness temperature data.Each grid-cell snow depth was estimated as the sum of snow depths from each land cover algorithm weighted by percentages of land cover types within each grid cell.Through evaluation of this algorithm using station measurements from 2006,the root mean square error(RMSE)of snow depth retrieval is about 5.6 cm.In forest regions,snow depth is underestimated relative to ground observation,because stem volume and canopy closure are ignored in current algorithms.In addition,comparison between snow cover derived from AMSR-E and FY3B-MWRI with Moderate-resolution Imaging Spectroradiometer(MODIS)snow cover products(MYD10C1)in January 2010 showed that algorithm accuracy in snow cover monitoring can reach 84%.Finally,we compared snow water equivalence(SWE)derived using FY3B-MWRI with AMSR-E SWE products in the Northern Hemisphere.The results show that AMSR-E overestimated SWE in China,which agrees with other validations.     

Evaluating the potential of a blended passive microwave‐interactive multi‐sensor product for improved mapping of snow cover and estimations of snow water equivalent

The retrieval of Snow Water Equivalent (SWE) from remote sensing satellites continues to be a very challenging problem. In this paper, we evaluate the accuracy of a new SWE product derived from the blending of a passive microwave SWE product based on the Advanced Microwave Sounding Unit (AMSU) with a multi‐sensor snow cover extent product based on the Interactive Multi‐sensor Snow and Ice Mapping System (IMS). The microwave measurements have the ability to penetrate the snow pack, and thus, the retrieval of SWE is best accomplished using the AMSU. On the other hand, the IMS maps snow cover more reliably due to the use of multiple satellite and ground observations. The evolution of global snow cover from the blended, the AMSU and the IMS products was examined during the 2006 snow season. Despite the overall good inter‐product agreement, it was shown that the retrievals of snow cover extent in the blended product are improved when using IMS, with implications for improved microwave retrievals of SWE. In a separate investigation, the skill of the microwave SWE product was also examined for its ability to correctly estimate SWE globally and regionally. Qualitative evaluation of global SWE retrievals suggested dependence on land surface temperature: the lower the temperature, the higher the SWE retrieved. This temperature bias was attributed in part to temperature effects on those snow properties that impact microwave response. Therefore, algorithm modifications are needed with more dynamical adjustments to account for changing snow cover. Quantitative evaluation over Slovakia in central Europe, for a limited period in 2006, showed reasonably good performance for SWE less than 100 mm. Sensitivity to deeper snow decreased significantly. Copyright © 2007 John Wiley & Sons, Ltd.     

Satellite remote sensing of river inundation area,stage, and discharge: a review

The growing availability of multi-temporal satellite data has increased opportunities for monitoring large rivers from space. A variety of passive and active sensors operating in the visible and microwave range are currently operating, or planned, which can estimate inundation area and delineate flood boundaries. Radar altimeters show great promise for directly measuring stage variation in large rivers. It also appears to be possible to obtain estimates of river discharge from space, using ground measurements and satellite data to construct empirical curves that relate water surface area to discharge. Extrapolation of these curves to ungauged sites may be possible for the special case of braided rivers. Where clouds, trees and floating vegetation do not obscure the water surface, high-resolution visible/infrared sensors provide good delineation of inundated areas. Synthetic aperture radar (SAR) sensors can penetrate clouds and can also detect standing water through emergent aquatic plants and forest canopies. However, multiple frequencies and polarizations are required for optimal discrimination of various inundated vegetation cover types. Existing single-polarization, fixed-frequency SARs are not sufficient for mapping inundation area in all riverine environments. In the absence of a space-borne multi-parameter SAR, a synergistic approach using single-frequency, fixed-polarization SAR and visible/infrared data will provide the best results over densely vegetated river floodplains. © 1997 John Wiley & Sons, Ltd.     

Ice‐cover variability on shallow lakes at high latitudes: model simulations and observations

A one‐dimensional thermodynamic model for simulating lake‐ice phenology is presented and evaluated. The model can be driven with observed daily or hourly atmospheric forcing of air temperature, relative humidity, wind speed, cloud amount and snowfall. In addition to computing the energy balance components, key model output includes the temperature profile at an arbitrary number of levels within the ice/snow (or the water temperature if there is no ice) and ice thickness (clear ice and snow‐ice) on a daily basis, as well as freeze‐up and break‐up dates. The lake‐ice model is used to simulate ice‐growth processes on shallow lakes in arctic, sub‐arctic, and high‐boreal forest environments. Model output is compared with field and remote sensing observations gathered over several ice seasons. Simulated ice thickness, including snow‐ice formation, compares favourably with field measurements. Ice‐on and ice‐off dates are also well simulated when compared with field and satellite observations, with a mean absolute difference of 2 days. Model simulations and observations illustrate the key role that snow cover plays on the seasonal evolution of ice thickness and the timing of spring break‐up. It is also shown that lake morphometry, depth in particular, is a determinant of ice‐off dates for shallow lakes at high latitudes. Copyright © 2003 John Wiley & Sons, Ltd.     

Characteristics of thermal radiation of ice covers on water bodies with different mineralization

The microwave brightness temperature of ice cover on fresh and slightly saline noneutrophic water bodies with mineralization varying from 0.1 to 3 g/l is examined. The optimal conditions for the determination of salt concentration in ice cover by radio thermal radiation are identified.     

星载微波ERS-1和多通道SSM/I观测海洋的相关性和随机粗糙面复合模型的数值模拟

欧洲遥感卫星（ERS）和美国防卫气象卫星计划（DMSP）空对地微波遥感是当今研究全球大气地表微波散射辐射和反演地球物理与水文参数的主要数据来源．本文研究了ERS－1散射计和SSM／I多通道辐射计在中国海域观测到的后向散射和热辐射数据,论证了同一地区同一时间段内ERS主动散射计和SSM／I被动辐射计观测数据的相关性．用海域主、被动遥感数据的比较,阐述了主、被动联合多通道分析方法有利于监视和分析复杂地表和海面在时间和空间尺度上的变化．用带泡沫散射层的双尺度随机粗糙面的复合模型计算后向散射和热发射,用以数值模拟ERS和SSM／I数据．并讨论了后向散射与热辐射数值模拟结果的相关性,以及与星载微波遥感器实际观测结果的比较．     

Observations of late winter Canadian tundra snow cover properties

Tundra snow cover is important to monitor as it influences local, regional, and global‐scale surface water balance, energy fluxes, as well as ecosystem and permafrost dynamics. Observations are already showing a decrease in spring snow cover duration at high latitudes, but the impact of changing winter season temperature and precipitation on variables such as snow water equivalent (SWE) is less clear. A multi‐year project was initiated in 2004 with the objective to quantify tundra snow cover properties over multiple years at a scale appropriate for comparison with satellite passive microwave remote sensing data and regional climate and hydrological models. Data collected over seven late winter field campaigns (2004 to 2010) show the patterns of snow depth and SWE are strongly influenced by terrain characteristics. Despite the spatial heterogeneity of snow cover, several inter‐annual consistencies were identified. A regional average density of 0.293 g/cm³ was derived and shown to have little difference with individual site densities when deriving SWE from snow depth measurements. The inter‐annual patterns of SWE show that despite variability in meteorological forcing, there were many consistent ratios between the SWE on flat tundra and the SWE on lakes, plateaus, and slopes. A summary of representative inter‐annual snow stratigraphy from different terrain categories is also presented. © 2013 Her Majesty the Queen in Right of Canada. Hydrological Processes. © 2013 John Wiley & Sons, Ltd.     

Detection and classification of surface roughness in an Arctic marginal sea ice zone

Sea ice dynamic and thermodynamic processes are important and highly variable elements of the marginal ice zone (MIZ). This study examines the detection and classification of statistically separable sea ice classes in the MIZ through a range of temporal and spatial scales. A helicopter‐based laser system was used to obtain large‐scale and a ship‐based laser profiler to identify small‐scale roughness types, respectively. The analysis of variance of surface height data from helicopter‐ and ship‐based laser systems, active microwave (AMW) C‐band backscattering data and passive microwave (PMW) (37 and 89 GHz) brightness temperature data reveal different classes that statistically differ from one another. We found significant statistical difference in variances in AMW data with six classes that differ in VV polarization, three classes in VH polarization, and five classes in HH polarization in the MIZ (e.g. snow‐covered first‐year ice, ice rubble, pancake ice, frost flowers, melt pond, flooded ice, and ice edge) of southeastern Beaufort Sea. The PMW emission was not as effective at discrimination, yielding only one statistically separable class. The results can potentially be extended to satellite‐based investigations of the MIZ at regional scales. Copyright © 2012 John Wiley & Sons, Ltd.     

New satellite derived sea ice motion tracks Arctic contamination

Sea ice has been reported to contain contaminants from atmospheric and nearshore sediment resuspension processes. In this study successive passive microwave images from the 85.5 GHz channels on the Special Sensor Microwave Imager (SSM/I) were merged with drifting buoy trajectories from the International Arctic Buoy Program to compute Arctic sea ice motion in the Russian Arctic between 1988 and 1994. Smooth daily motion fields were averaged to prepare monthly maps making it possible to compute the 7-year mean and mean seasonal ice motions as well as principal components of directional variability of sea ice motion for the entire Arctic and surrounding basins. These mean motion vectors are used to simulate the advection of contaminants deposited on or contained within the sea ice and subsequently transported into the Arctic Ocean in order to predict both their mean trajectories and dispersal over time. The 3-year displacement of contaminants from a number of Russian sites and one American site display various behaviours from substantial displacement and dispersal to almost no movement. This computational procedure could be applied to realtime SSM/I and ice buoy data to provide detailed, all-weather, vector motion maps of ice circulation to predict the path and dispersal of any new substance introduced to the sea ice and transported into the Arctic or Antarctic ocean surface.     

卫星遥感藏北积雪分布及影响因子分析

利用1993~2004年SSM/I被动微波辐射仪反演的雪深资料,1996~2004年NOAA/AVHRR可见光和红外反演的积雪覆盖面积资料,1966~2003年藏北地区6个地面台站的积雪观测资料来检验卫星资料的可用性,并研究近年来藏北积雪的时空分布和影响因素.结果表明,SSM/I, NOAA/AVHRR和实际观测的积雪资料具一致性.从积雪时间变化看：季节尺度上,藏北地区秋冬季积雪迅速增加,但春季(3~5月)融雪速度不快,呈现正反馈特征;年际尺度上,藏北地区20世纪60年代末期起积雪开始减少,80年代积雪增加,90年代起到2003年积雪总体上减少,呈现出减少—增加—减少趋势.采用小波分析发现积雪振荡周期存在着一个准2~3年,准9年和13年的周期,从20世纪70年代初到90年代中期还有一个5年的周期.积雪空间上看,藏北地区积雪主要集中在东部地区,该区每个冬春年积雪覆盖旬数超过15旬,显著高于西部少雪区,大部分积雪集中在4900~5600 m的高度左右;藏北高原积雪变动的显著区位于藏北中东部的安多和聂荣地区.利用藏北地区1966~2003年的地面温度和降水资料建立回归方程模拟年累积雪日,结果表明模拟值与实测值之间的相关系数达0.74.积雪时空分布受温度、降水等因子影响明显.1998~2003年藏北积雪的减少与全球变暖有关,但降水的减少可能是导致近年来藏北积雪减少的更主要因素.     

Improved sea ice parcel trajectories in the Arctic via data assimilation

An assimilated sea ice motion product is used to track ice parcels in several regions of the Arctic over time periods of one day to several weeks during 1992-1993. Motions simulated using a two-dimensional, dynamic-thermodynamic sea ice model are combined with motions derived from daily 85 GHz special sensor microwave/imager (SSM/I) imagery using an optimal interpolation method that minimizes error covariance. Assimilation attenuates the tracking error over the stand-alone model in comparison to buoy trajectories with the same starting location and time. The average 14-day assimilated trajectory's displacement error is as much as 34% lower than the model trajectory, while the RMS direction error is decreased by up to 10 degrees (24%). Assimilation can also yield an estimate of dispersion, which is not retrievable by point buoy observations. An assimilation approach improves estimates of ice drift and has the potential to further the understanding of ice mass flux, freshwater flux, and pollutant transport in the polar regions.     

查干湖冰封期光谱特征及影响因素

湖冰光谱特征是湖冰遥感反演的物理基础,是研究湖冰光学特性和空间分布的理论依据。本文以查干湖为例,使用ASD Field Spec 4便携式地物光谱仪采集冰封期不同类型湖冰、积雪和水体光谱,利用Savitzky-Golay滤波法和包络线去除法分析白冰、灰冰、黑冰、雪冰、积雪和水体的反射光谱特征,探索气泡对湖冰反射光谱特征的影响。积雪和雪冰、白冰和灰冰、黑冰和水体的反射特征随着波长的变化特征基本一致,冰的反射率介于积雪和水体之间,其中白冰的反射率高于灰冰和黑冰,在包络线去除结果中,黑冰和水体在440 nm吸收谷处的吸收面积为5.184和10.878、吸收深度为0.052和0.106,雪、雪冰、白冰、灰冰在800和1030 nm吸收谷处的吸收面积和吸收深度的变化表现为雪<雪冰<灰冰<白冰。气泡是影响湖冰光谱特征的重要因素,气泡使白冰反射率减小和黑冰反射率增大,并且气泡使得白冰在800/1030nm和黑冰在440 nm处的吸收面积和吸收深度减小,其中气泡大小和疏密程度的不同会导致湖冰反射率的影响程度存在差异。同时,本文选取时间同步的Landsat 8 OLI遥感影像,在完成辐...     

Assimilation of ice concentration in a coupled ice–ocean model,using the Ensemble Kalman filter

An implementation of the Ensemble Kalman filter (EnKF) with a coupled ice–ocean model is presented. The model system consists of a dynamic–thermodynamic ice model using the elastic-viscous-plastic (EVP) rheology coupled with the HYbrid Coordinate Ocean Model (HYCOM). The observed variable is ice concentration from passive microwave sensor data (SSM/I). The assimilation of ice concentration has the desired effect of reducing the difference between observations and model. Comparison of the assimilation experiment with a free-run experiment shows that there are large differences, especially in summer. In winter the differences are relatively small, partly because the atmospheric forcing used to run the model depends upon SSM/I data. The assimilation has the strongest impact close to the ice edge, where it ensures a correct location of the ice edge throughout the simulation. An inspection of the model ensemble statistics reveals that the error estimates of the model are too small in winter, partly a result of too low model ice-concentration variance in the central ice pack. It is found that the ensemble covariance between ice concentration and sea-surface temperature in the same grid cell is of the same sign (negative) throughout the year. The ensemble covariance between ice concentration and salinity is more dependent upon the physical mechanisms involved, with ice transport and freeze/melt giving different signs of the covariances. The ice-transport and ice-melt mechanisms also impact the ice-concentration variance and the covariance between ice concentration and ice thickness. The ensemble statistics show a high degree of complexity, which to some extent merits the use of computationally expensive assimilation methods, such as the Ensemble Kalman filter. The present study focuses on the assimilation of ice concentration, but it is understood that assimilation of other datasets, such as sea-surface temperature, would be beneficial.Responsible Editor: Jin-Song von Storch