Can we quantify controls of snow-albedo feedback using satellite-based snow & albedo estimations ?

1. Can we quantify controls of snow-albedo feedback using satellite-based snow & albedo estimations? Hongxu Zhao and Richard Fernandes Canada Centre for Remote Sensing

2. Outline Part 1: • Motivation & Hall et al’s work • A new snow cover dataset Part 2: • Quantifying Snow-albedo feedback(SAF) • Summary

3. Qu and Hall(2006) calculated the two factors based on outputs of 17 climate models used in IPCC AR4, Atm Factor (atmosphericattenuation) The models agree each other to within 10%. The models also agree with an estimate of the observed value based on ISCCP data. Sur Factor (surface process) exhibits a three-fold spread in these models. Hence this Surface Factor is main source of the divergence in simulations of SAF. Suggestions: “to reduce the divergence, attention should be focused on differing parameterizations of snow process, rather than intermodel variations in the attenuation effect of the atmosphere on surface albedo anomalies.” Snow albedo feedback (SAF) Atmospheric factor Following Qu and Hall (2006, 2007), and Cess and Potter (1988), the strength of SAF can be determined by the product of two terms, one representing the dependence ofplanetary albedo on surface albedo and the other representing the change in surface albedo induced by a unit surface air temperature change. Atmosp surface Surface factor Where Q (constant) and Qnet are the incoming and net shortwave radiation at TOA, αs is the surface albedo, and αpis the planetary albedo.

4. Snow albedo feedback (two loops) 1. Snow cover shrinks:an increase in surface air temperature (Ts) leads to a reduction in snow cover (Sc), which in turn leads to a reduction in surface albedo (αs ). 2. Metamorphosis of remaining snowpack. Reverse effect: The additional absorbed solar radiation results in further warming (K4). The strongest SAF occurs during the spring snow melting season (Qu and Hall, 2007; Hall, 2004). K4 K2 αs K1 Ts Sc K3 (pre melting=0, metamorphosis≠0, post melting≠0) Ideal model: (∆αs/∆Ts )=k3+K2· (∆Sc/ ∆Ts), where (∆Sc/ ∆Ts)=k1 ∆αs/∆Ts (K4-1) controlled by two components, one arising from the metamorphosisof the snowpack in a warmer climate (K3) and the other from the reduction of snow cover (k2*k1). (Qu and Hall 2007)

5. Using Seasonal cycle to simulate climate change Hall and Qu (2006) presented a scatterplot of spring Sur Factor (∆αs/∆Ts) values in climate change (22nd-centery-mean minus 20th-centery-measn) vs. in seasonal cycle (from April to May) averaged over NH continents polarward of 30deg. The observed value is based on ISCCP and ERA40. It shows that intermodel variations in Sur Factor in the seasonal cycle are highly correlated with those in climate change. Hence, they suggest that the Sur Factor (∆αs/∆Ts ) based on the present-day seasonal cycle are excellent predictors of the Sur Factor (∆αs/∆Ts ) in climate change. -1.1

6. A Daily 5km snow cover product extracted from AVHRR Polar Pathfinder (APP) data (APP snow) Current in-situ snow cover datasets have limited spatial coverage while satellite-based snow cover records have either limited historical extent (e.g. MODIS) or limited temporal and spatial resolution (e.g., NOAA weekly snow cover, Robinson) constrained by clouds, specific sensor availability, or processing methodology. Based on a new CCRS snow mapping algorithm, we have produced the new daily 5km APP snow cover (Zhao and Fernandes, 2009 JGR in press), including during cloudy conditions, over Northern Hemisphere land surfaces over 1982-2004. The APP snow cover maps showed an 85% agreement rate or better at 95% of the in-situ sites (at a comparable level of agreement to in-situ snow cover for MODIS 0.05 degree snow cover estimates). The almost continuous spatial and temporal coverage ability of the APP snow product will benefit estimation of spring snowmelt dates and snow-albedo feedback over northern circumpolar regions. Wang and Key (2005) developed all sky APP extended daily 25km albedo product over the period of time. The snow cover (subsampled to 25km) was recalculated into snow fraction based on an approximate 2.5 degree grid cell box while the albedo is simply averaged over the same 2.5 degree grid cell box to match the ERA40 air temperature spatial resolution. Finally, the daily data was averaged to obtain monthly mean fields.

7. Almost continuous spatial coverage of the APP snow maps White=snow; Green=land; Black=not available; Blue=water. Figure 1. Snow cover maps of (a) APP and (b) MODIS on 30 April 2002. Note that white color indicates snow cover, green indicates snow free land, black indicates missing or no data, and blue indicates water bodies or outside the APP domain.

8. Almost continuous temporal coverage of the APP snow maps Numbers of days retrieved from APP and MODIS Figure 2. Comparison of available number of days retrieved per site over Canada between APP (dark) and MODIS (light) snow product in 2002.

9. Similar agreement rates between APP snow and MODIS snow with in-situ measurements Figure 3. Comparison of snow mapping agreement rates between APP snow and MODIS snow for Canadian sites. The top panels represent an estimation for the full period of the year while the bottom ones represent an estimation for a partial period from January 1st to August 31. In-situ sites over Canada & Eurasia

10. Snowmelt dates based on the APP snow Figure 4. Spatial distribution of the mean snowmelt dates (unit: DOY) over Northern Hemisphere, presented (a) by 20 DOY interval (b) continuously but focusing on high latitudes. The value of 250 (white color) is assigned to permanent snow.

11. Climatological mean albedo(%), snow-fraction(%), & temperature(K) during snow melting season (1982-1999) May April Jun APP surface albedo (Wang & Key 2005) APP Snow fraction (Zhao & Fernandes) ERA40 Temperature (Uppala et al 2005)

12. Seasonal changes in albedo(%), snow-fraction(%), & temperature(K) April to May % May to Jun Apr to Jun ∆Albedo (∆α) ∆Snow (∆Sc) ∆Temp (∆Ts)

13. Feedbacks: ∆αs/∆T, ∆Sc/ ∆T, & ∆αs/∆Sc (%/K) April to May May to Jun Apr to Jun ∆αs/∆T time series ∆αs/∆T (K4-1=K3+k2·K1) (ISCCP-1.1) ∆Sc/ ∆T (K1) ∆αs/∆Sc (K2 if k3=0)

14. Approach 1: linear regression Using the method of Ordinary Least Squares, the regression line that minimizes the sum of the squares of errors can be found based on the 1982~1999 satellite observations. (∆αs/∆Ts)=k3+K2·(∆Sc/ ∆Ts) y=a+b · x Robust regression was suggested by Fernandes and Leblanc (2005) when dealing with remote sensing data that are not normally random distributions, say, with outliers.

15. K3 and K2 distribution based on regressions K3 K2 Histogram plot Regular regress Global average K2~= 0.20 K3 =-0.71(%K-1) Robust regress

16. Comparison between regular and Theil’s regressions K3 K2 t-values of t-test Global average K2~= 0.20 K3=-0.71 Regular regress Lower confidence bound K2~= 0.30 K3=-0.52 Thiel’s regress

17. Approach 2: an approximation Qu and Hall (06) developed an approximation method to calculate K2 and K3: K2=1/2·(αfsnow + αpsnow )- αland K3=1/2·(Sfc+Spc)·(αfsnow - αpsnow )/ ∆Ts - αland Where, the superscripts “f” and “p” represent the future and present climate; Snow-free land albedo αland is defined as mean albedo from snowmelt end date to 10 days after; effective snow albedo (completely snow-covered surface) αisnow ={ αis -(1-Sic) αland}/Sic; Sc represents snow cover fraction. K2: the mean contrast between effective snow albedo and snow free land albedo. K3: changes in effective snow albedo with respect to changes in Ts, modulated by mean snow cover.

18. K3 and K2 distribution based on Approach2 K2 K3 Compare Theil’s regress Median=-0.77 K2~=0.30 Regression: -0.71 0.20

19. (unitless) (%/K)

20. Simply averaging (–2.5 ~ 1) rather than considering physical meaning, K3= -0.3 Qu and Hall (2007) Causes of spread in models: SAF∆αs/∆TsK2·(∆Sc/ ∆Ts) K2 K2·(∆Sc/ ∆Ts) -0.76 -0.50 ?

21. Summary • Based on the APP snow and albedo datasets, we can obtain circumpolar quantities of the snow-albedo feedback parameters, and the global average magnitude ∆αs/∆Tsis –1.16 ± 0.06comparable to Qu and Hall ISCCP’s value (-1.1). • k2 (0.2~0.3) is within the range of models’ spread (Qu and Hall, 2006). K2 is of large-scale spatial characteristics, stable for multiple years. • K3 (-0.5~-0.7) is not uniformly distributed and has an apparent regional characteristics and interannual variability.