Evaluation of SMAP Freeze/Thaw Retrieval Accuracy at Core Validation Sites in the Contiguous United States

Seasonal freeze-thaw (FT) impacts much of the northern hemisphere and is an important control on its water, energy, and carbon cycle. Although FT in natural environments extends south of 45°N, FT studies using the L-band have so far been restricted to boreal or greater latitudes. This study addresses this gap by applying a seasonal threshold algorithm to Soil Moisture Active Passive (SMAP) data (L3_SM_P) to obtain a FT product south of 45°N (‘SMAP FT’), which is then evaluated at SMAP core validation sites (CVS) located in the contiguous United States (CONUS). SMAP landscape FT retrievals are usually in good agreement with 0–5 cm soil temperature at SMAP grids containing CVS stations (\u3e70%). The accuracy could be further improved by taking into account specific overpass time (PM), the grid-specific seasonal scaling factor, the data aggregation method, and the sampling error. Annual SMAP FT extent maps compared to modeled soil temperatures derived from the Goddard Earth Observing System Model Version 5 (GEOS-5) show that seasonal FT in CONUS extends to latitudes of about 35–40°N, and that FT varies substantially in space and by year. In general, spatial and temporal trends between SMAP and modeled FT were similar

Retrieving landscape freeze/thaw state fromSoil Moisture Active Passive (SMAP) radar and radiometer measurements

Over one-third of the global land area undergoes a seasonal transition between predominantly frozen and non-frozen conditions each year, with the land surface freeze/thaw (FT) state a significant control on hydrological and biospheric processes over northern land areas and at high elevations. The NASA Soil Moisture Active Passive (SMAP) mission produced a daily landscape FT product at 3-km spatial resolution derived from ascending and descending orbits of SMAP high-resolution L-band (1.4 GHz) radar measurements. Following the failure of the SMAP radar in July 2015, coarser (36-km) footprint SMAP radiometer inputs were used to develop an alternative daily passive microwave freeze/thaw product. In this study, in situ observations are used to examine differences in the sensitivity of the 3-km radar versus the 36-km radiometer measurements to the landscape freeze/thaw state during the period of overlapping instrument operation. Assessment of the retrievals at high-latitude SMAP core validation sites showed excellent agreement with in situ flags, exceeding the 80% SMAP mission accuracy requirement. Similar performance was found for the radar and radiometer products using both air temperature and soil temperature derived FT reference flags. There was a tendency for SMAP thaw retrievals to lead the surface flags due to the influence of wet snow cover conditions on both the radar and radiometer signal. Comparison with other satellite derived FT products showed those derived from passive measurements (SMAP radiometer; Aquarius radiometer; Advanced Microwave Scanning Radiometer - 2) retrieved less frozen area than the active products (SMAP radar; Aquarius radar)

Amélioration de la caractérisation de la neige et du sol arctique afin d’améliorer la prédiction de l’équivalent en eau de la neige en télédétection micro-ondes

Le phénomène de l’amplification arctique consiste en une augmentation plus prononcée des températures de surface dans cette région que sur le reste du globe. Ce phénomène est notamment dû à la diminution marquée du couvert nival provoquant un déséquilibre dans le bilan d’énergie de surface via une réduction généralisée de l’albédo (rétroaction positive). L’accélération du réchauffement est jusqu’à trois fois plus élevée dans ces régions. Il est donc primordial, dans un contexte de changement climatique arctique, de poursuivre et d’améliorer le suivi à grande échelle du couvert nival afin de mieux comprendre les processus gouvernant la variabilité spatio-temporelle du manteau neigeux. Plus spécifiquement, l’Équivalent en Eau de la Neige (EEN) est généralement utilisé pour quantifier deux propriétés (hauteur et densité) de la neige. Son estimation à grande échelle dans les régions éloignées tel que l’Arctique provient actuellement essentiellement de produits en micro-ondes passives satellitaires. Cependant, il existe encore beaucoup d’incertitudes sur les techniques d’assimilation de l’ÉEN par satellite et ce projet vise une réduction de l’erreur liée à l’estimation de l’ÉEN en explorant deux des principales sources de biais tels que : 1) la variabilité spatiale de l’épaisseur et des différentes couches du manteau neigeux arctique liées à la topographie et la végétation au sol influençant l’estimation de l’ÉEN; et 2) les modèles de transfert radiatif micro-ondes de la neige et du sol ne bénéficient pas actuellement d’une bonne paramétrisation en conditions arctiques, là où les erreurs liées à l’assimilation de l’ÉEN sont les plus importantes. L’objectif global est donc d’analyser les propriétés géophysiques du couvert nival en utilisant des outils de télédétection et de modélisation pour diminuer l’erreur liée à la variabilité spatiale locale dans l’estimation du ÉEN à grande échelle, tout en améliorant la compréhension des processus locaux qui affectent cette variabilité. Premièrement, une analyse haute résolution à l’aide de l’algorithme Random Forest a permis de prédire la hauteur de neige à une résolution spatiale de 10 m avec une RMSE de 8 cm (23%) et d’en apprendre davantage sur les processus de distribution de la neige en Arctique. Deuxièmement, la variabilité du manteaux neigeux arctique (hauteur et microstructure) a été incorporée dans des simulations en transfert radiatif micro-ondes de la neige et comparée au capteur satellitaire SSMIS. L’ajout de variabilité améliore la RMSE des simulations de 8K par rapport à un manteau neigeux uniforme. Finalement, une paramétrisation du sol gelé est présentée à l’aide de mesures de rugosité provenant de photogrammétrie (Structure-from-Motion). Cela a permis d’investiguer trois modèles de réflectivité micro-ondes du sol ainsi que la permittivité effective du sol gelé avec une rugosité SfM d’une précision de 0.1 mm. Ces données de rugosité SfM avec une permittivité optimisée (ε'_19 = 3.3, ε'_37 = 3.6) réduisent significativement l’erreur des températures de brillance simulées par rapport à des mesures au sols (RMSE = 3.1K, R^2 = 0.71) pour toutes les fréquences et polarisations. Cette thèse offre une caractérisation des variables de surface (neige et sol) en Arctique en transfert radiatif micro-ondes qui bénéficie aux multiples modélisations (climatiques et hydrologiques) de la cryosphère

Analyse des cycles gel/dégel des régions nordiques par télédétection micro-ondes passives en bande L

Le réchauffement climatique dans les régions nordiques, fort important depuis le milieu du siècle dernier, a de multiples impacts sur la dynamique des écosystèmes, notamment sur les cycles gel/dégel de surface qui influencent les flux de carbone, l'activité biogéochimique des sols, l'hydrologie et le pergélisol aux hautes latitudes. La télédétection satellitaire du gel/dégel par micro-ondes passives est un outil très prometteur permettant un suivi continu et global, mais comporte des difficultés souvent reliées à l’effet d’hétérogénéité spatiale intra-pixel relié aux résolutions grossières des capteurs micro-ondes passives à basse fréquence. L’objectif principal du projet est d’évaluer l’utilisation de la télédétection micro-onde passive en bande L (1.4 GHz) pour le suivi de l’état de gel/dégel de la surface en forêt boréale. Un premier objectif spécifique est d’évaluer un nouveau produit des cycles de gel/dégel de surface estimée à partir des radiomètres bande L satellitaires Aquarius. Cette base de données de 3.5 années a été mise en ligne au National Snow and Ice Data Center (NSIDC). Le deuxième objectif spécifique est d’analyser l’effet de la variabilité spatiale intrapixel de l’état de gel du sol et de son impact sur les températures de brillance (TB) mesurées par le radiomètre de la mission Soil Moisture Active Passive (SMAP) en période de transition afin de quantifier la fraction de sol gelé. Les résultats pour le premier objectif montrent que la nouvelle base de données possède une bonne capacité à estimer l’état de gel/dégel de la surface sur l’ensemble de l’Hémisphère Nord (> 50°N). Cette recherche offre également une rare intercomparaison entre produits de gel/dégel satellitaires en comparant le produit Aquarius au Freeze/Thaw-Earth System Data Record (FT-ESDR) développé avec les données à plus hautes fréquences du capteur SSM/I. Pour le deuxième objectif, des capteurs de température distribués le long de transects de plusieurs kilomètres sur deux différents sites de taïga montrent que la variabilité spatiale du gel à l’automne peut être de 7.5 à 9.5 semaines. Il est également démontré que les mesures de SMAP sont sensibles à cette variabilité et un algorithme développé permet d’estimer le pourcentage intrapixel de sol gelé avec des coefficients de détermination (R2) entre 0.63 et 0.88 lorsque comparé aux mesures in situ. Ces résultats offrent de nouveaux outils pour mieux comprendre et quantifier les cycles de gel/dégel de l’environnement boréal et leurs impacts sur les processus biogéophysiques, hydrologiques et sur le pergélisol.Abstract: Climate change in nordic regions, which has been of growing significance over the past century has multiple impacts on the dynamic of ecosystems, notably on the surface freeze/thaw cycles, which influences carbon flux, soil biogeochemical activity, hydrology and permafrost at high latitudes. Satellite remote sensing of freeze/thaw with passive microwaves is a promising tool to offer continuous and global monitoring, but can also entail some difficulties due to intra-pixel spatial variability effects coming from the low resolution of low-frequency passive microwave sensors. The primary objective of the project is to evaluate the use of passive microwave remote sensing in L-band (1.4 GHz) for monitoring of the surface freeze/thaw in the boreal forest. A first specific objective is to evaluate a new surface freeze/thaw product estimated by the Aquarius satellite L-band radiometers. This 3.5 year-old database has been put online at the National Snow and Ice Data Center (NSIDC) website. The second specific objective is to analyse the effect of intra-pixel spatial variability of freeze/thaw and its impact on brightness temperatures (TB) measured by the Soil Moisture Active Passive (SMAP) radiometer during transition periods in order to quantify the frozen soil fraction. Results for the first objective show that the new database possesses a good capacity to estimate the surface freeze/thaw state for the entirety of the Northern Hemisphere (>50°N). This research also offers a rare intercomparison between freeze/thaw satellite products by comparing the Aquarius product to the Freeze/Thaw-Earth System Data Record (FT-ESDR) product developed with higher frequencies data of the SSM/I sensor. For the second objective, temperature sensors distributed along transects of several kilometers on two different taiga sites show that the spatial variability of autumn soil freeze onset can be between 7.5 and 9.5 weeks. It demonstrates that SMAP measurements are sensitive to this variability and a developed algorithm offers estimations of the intrapixel soil frozen fraction with coefficients of determination (R2) between 0.63 and 0.88 when compared to in situ measurements. These results offer new tools for a better understanding and quantification of freeze/thaw cycles in boreal environments and their impacts on biogeochemical and hydrologic processes and on permafrost

An Investigation into the Effects of Variable Lake Ice Properties on Passive and Active Microwave Measurements Over Tundra Lakes Near Inuvik, N.W.T.

The accurate estimation of snow water equivalent (SWE) in the Canadian sub-arctic is integral to climate variability studies and water availability forecasts for economic considerations (drinking water, hydroelectric power generation). Common passive microwave (PM) snow water equivalent (SWE) algorithms that utilize the differences in brightness temperature (Tb) at 37 GHz – 19 GHz falter in lake-rich tundra environments because of the inclusion of lakes within PM pixels. The overarching goal of this research was to investigate the use of multiple platforms and methodologies to observe and quantify the effects of lake ice and sub-ice water on passive microwave emission for the purpose of improving snow water equivalent (SWE) retrieval algorithms. Using in situ snow and ice measurements as input, the Helsinki University of Technology (HUT) multi-layer snow emission model was modified to include an ice layer below the snow layer. Emission for 6.9, 19, 37 and 89 GHz were simulated at horizontal and vertical polarizations, and were validated by high resolution airborne passive microwave measurements coincident with in situ sampling sites over two lakes near Inuvik, Northwest Territories (NWT). Overall, the general magnitude of brightness temperatures were estimated by the HUT model for 6.9 and 19 GHz H/V, however the variability was not. Simulations produced at 37 GHz exhibited the best agreement relative to observed temperatures. However, emission at 37 GHz does not interact with the radiometrically cold water, indicating that ice properties controlling microwave emission are not fully captured by the HUT model. Alternatively, active microwave synthetic aperture radar (SAR) measurements can be used to identify ice properties that affect passive microwave emission. Dual polarized X-band SAR backscatter was utilized to identify ice types by the segmentation program MAGIC (MAp Guided Ice Classification). Airborne passive microwave transects were grouped by ice type classes and compared to backscatter measurements. In freshwater, where there were few areas of high bubble concentration at the ice/water interface Tbs exhibited positive correlations with cross-polarized backscatter, corresponding to ice types (from low to high emission/backscatter: clear ice, transition zone between clear and grey ice, grey ice and rafted ice). SWE algorithms were applied to emission within each ice type producing negative or near zero values in areas of low 19 GHz Tbs (clear ice, transition zone), but also produced positive values that were closer to the range of in situ measurements in areas of high 19 GHz Tbs (grey and rafted ice). Therefore, cross-polarized X-band SAR measurements can be used as a priori ice type information for spaceborne PM algorithms, providing information on ice types and ice characteristics (floating, frozen to bed), integral to future tundra-specific SWE retrieval algorithms

POTENTIAL CONTRASTS IN CO2 AND CH4 FLUX RESPONSE UNDER CHANGING CLIMATE CONDITIONS: A SATELLITE REMOTE SENSING DRIVEN ANALYSIS OF THE NET ECOSYSTEM CARBON BUDGET FOR ARCTIC AND BOREAL REGIONS

The impact of warming on the net ecosystem carbon budget (NECB) in Arctic-boreal regions remains highly uncertain. Heightened CH4 emissions from Arctic-boreal ecosystems could shift the northern NECB from an annual carbon sink further towards net carbon source. Northern wetland CH4 fluxes may be particularly sensitive to climate warming, increased soil temperatures and duration of the soil non-frozen period. Changes in northern high latitude surface hydrology will also impact the NECB, with surface and soil wetting resulting from thawing permafrost landscapes and shifts in precipitation patterns; summer drought conditions can potentially reduce vegetation productivity and land sink of atmospheric CO2 but also moderate the magnitude of CH4 increase. The first component of this work develops methods to assess seasonal variability and longer term trends in Arctic-boreal surface water inundation from satellite microwave observations, and quantifies estimate uncertainty. The second component of this work uses this information to improve understanding of impacts associated with changing environmental conditions on high latitude wetland CH4 emissions. The third component focuses on the development of a satellite remote sensing data informed Terrestrial Carbon Flux (TCF) model for northern wetland regions to quantify daily CH4 emissions and the NECB, in addition to vegetation productivity and landscape CO2 respiration loss. Finally, the fourth component of this work features further enhancement of the TCF model by improving representation of diverse tundra and boreal wetland ecosystem land cover types. A comprehensive database for tower eddy covariance CO2 and CH4 flux observations for the Arctic-boreal region was developed to support these efforts, providing an assessment of the TCF model ability to accurately quantify contemporary changes in regional terrestrial carbon sink/source strength

Tundra Snow Cover Properties from \u3cem\u3eIn-Situ\u3c/em\u3e Observation and Multi-Scale Passive Microwave Remote Sensing

Tundra snow cover is important to monitor as it influences local, regional, and global scale surface water balance, energy fluxes, and ecosystem and permafrost dynamics. Moreover, recent global circulation models (GCM) predict a pronounced shift in high latitude winter precipitation and mean annual air temperature due to the feedback between air temperature and snow extent. At regional and hemispheric scales, the estimation of snow extent, snow depth and, snow water equivalent (SWE) is important because high latitude snow cover both forces and reacts to atmospheric circulation patterns. Moreover, snow cover has implications on soil moisture dynamics, the depth, formation and growth of the permafrost active layer, the vegetation seasonality, and the respiration of C02. In Canada, daily snow depth observations are available from 1955 to present for most meteorological stations. Moreover, despite the abundance and dominance of a northern snow cover, most, if not all, long term snow monitoring sites are located south of 550N. Stations in high latitudes are extremely sparse and coastally biased. In Arctic regions, it can be logistically difficult and very expensive to acquire both spatially and temporally extensive in-situ snow data. Thus, the possibility of using satellite remote sensing to estimate snow cover properties is appealing for research in remote northern regions. Remote sensing techniques have been employed to monitor the snow since the 1960s when the visible light channels were used to map snow extent. Since then, satellite remote sensing has expanded to provide information on snow extent, depth, wetness, and SWE. However, the utility of satellite sensors to provide useful, operational tundra snow cover data depends on sensor parameters and data resolution. Passive microwave data are the only currently operational sources for providing estimates of dry snow extent, SWE and snow depth. Currently, no operational passive microwave algorithms exist for the spatially expansive tundra and high Arctic regions. The heterogeneity of sub-satellite grid tundra snow and terrain are the main limiting factors in using conventional SWE retrieval algorithm techniques. Moreover, there is a lack of in-situ data for algorithm development and testing. The overall objective of this research is to improve operational capabilities for estimating end of winter, pre-melt tundra SWE in a representative tundra study area using satellite passive microwave data. The study area for the project is located in the Daring-Exeter-Yamba portion of the Upper-Coppermine River Basin in the Northwest Territories. The size, orientation and boundaries of the study area were defined based on the satellite EASE grid (25 x 25 km) centroid located closest to the Tundra Ecosystem Research Station operated by the Government of the Northwest Territories. Data were collected during intensive late winter field campaigns in 2004, 2005, 2006, 2007, 2008, and 2009. During each field campaign, snow depth, density and stratigraphy were recorded at sites throughout the study area. During the 2005 and 2008 seasons, multi-scale airborne passive microwave radiometer data were also acquired. During the 2007 season, ground based passive microwave radiometer data were acquired. For each year, temporally coincident AMSR-E satellite Tb were obtained. The spatial distribution of snow depth, density and SWE in the study area is controlled by the interaction of blowing snow with terrain and land cover. Despite the spatial heterogeneity of snow cover, several inter-annual consistencies were identified. Tundra snow density is consistent when considered on a site-by-site basis and among different terrain types. A regional average density of 0.294 g/cm3 was derived from the six years of measurements. When applied to site snow depths, there is little difference in SWE derived from either the site or the regional average density. SWE is more variable from site to site and year to year than density which requires the use of a terrain based Classification to better quantify regional SWE. The variability in SWE was least on lakes and flat tundra, while greater on slopes and plateaus. Despite the variability, the interannual ratios of SWE among different terrain types does not change that much. The variability (CV) in among terrain categories was quite similar. The overall weighted mean CV for the study area was 0.40, which is a useful regional generalization. The terrain and landscape based classification scheme was used to generalize and extrapolate tundra SWE. Deriving a weighted mean SWE based on the spatial proportion of landscape and terrain features was shown as a method for generalizing the regional distribution of tundra SWE. The SWE data from each year were compared to AMSR-E satellite Tb. Within each season and among each of the seasons, there was little difference in 19 GHz Tb. However, there was always a large decrease in 37 GHz Tb from early November through April. The change in ΔTb37-19 throughout each season showed that the Tb at 37 GHz is sensitive to parameters which evolve over a winter season. A principal component analysis (PCA) showed that there are differences in ΔTb37-19 among different EASE grids and that land cover may have an influence on regional Tb. However, the PCA showed little relationship between end of season ΔTb37-19 and lake fraction. A good relationship was found between ΔTb37-19 and in-situ SWE. A quadratic function was fitted to explain 89 percent of the variance in SWE from the ΔTb37-19. The quadratic relationship provides a good fit between the data; however, the nature of the relationship is opposite to the expected linear relationship between ΔTb37-19 and SWE. Airborne Tb data were used to examine how different snow, land cover and terrain properties influence microwave emission. In flat tundra, there was a significant relationship between SWE and high resolution ΔTb37-19. On lakes and slopes, no strong relationships were found between SWE and high resolution ΔTb37-19. Due to the complexity of snow and terrain in high resolution footprints, it was a challenge to isolate a relationship between SWE and Tb. However, as the airborne footprint size increased the amplitude of variability in Tb decrease considerably to the point that Tb in large footprints is not sensitive to local scale variability in SWE. As such, most of the variability evident in the high and mid resolution airborne data will not persist at the EASE grid scale. Despite the many challenges, algorithm development should be possible at the satellite scale. The AMSR-E ΔTb37-19 changes from year to year in response to differences in snow cover properties. However, the multiple years of in-situ snow data remain the most important contribution in linking Tb with SWE

Theoretical Modeling and Analysis of L- and P-band Radar Backscatter Sensitivity to Soil Active Layer Dielectric Variations

Freeze-thaw (FT) and moisture dynamics within the soil active layer are critical elements of boreal, arctic and alpine ecosystems, and environmental change assessments. We evaluated the potential for detecting dielectric changes within different soil layers using combined L- and P-band radar remote sensing as a prerequisite for detecting FT and moisture profile changes within the soil active layer. A two-layer scattering model was developed and validated for simulating radar responses from vertically inhomogeneous soil. The model simulations indicated that inhomogeneity in the soil dielectric profile contributes to both L- and P-band backscatter, but with greater P-band sensitivity at depth. The difference in L- and P-band responses to soil dielectric profile inhomogeneity appears suitable for detecting associated changes in soil active layer conditions. Additional evaluation using collocated airborne radar (AIRSAR) observations and in situ soil moisture measurements over alpine tundra indicates that combined L- and P-band SAR observations are sensitive to soil dielectric profile heterogeneity associated with variations in soil moisture and FT conditions

Modélisation de l’émission micro-onde hivernale en forêt boréale canadienne

La caractérisation du couvert nival en forêt boréale est un élément important pour la compréhension des régimes climatiques et hydrologiques. Depuis plusieurs années, l’utilisation des micro-ondes passives est étudiée pour l’estimation de l’équivalent en eau de la neige (SWE : Snow Water Equivalent) à partir de capteurs satellitaires. Les algorithmes empiriques traditionnels étant limités en forêt boréale, le couplage d’un modèle de transfert radiatif (MTR) micro-onde passive (qui prend en compte les contributions du sol, de la neige, de la végétation et de l’atmosphère) avec un modèle de neige pour l’inversion du SWE semble une avenue prometteuse. La thèse vise donc à coupler un MTR avec le schéma de surface du modèle climatique canadien (CLASS) dans une perspective d’application opérationnelle pour les estimations de SWE à partir de données satellitaires micro-onde à 10.7, 19 et 37 GHz. Dans ce contexte, certains aspects centraux du MTR, dont l’effet de la taille des grains ainsi que la contribution de la végétation sont développés et quantifiés. Le premier aspect étudié dans la thèse concerne l’adaptation du modèle d’émission micro-onde passive DMRT-ML (Dense media radiative transfer theory – multi layer) pour l’intégration d’une nouvelle métrique représentant la taille des grains (surface spécifique des grains de neige: SSA). L’étude basée sur des mesures radiométriques et de neige in situ, montre la pertinence de l’utilisation de la SSA dans DMRT-ML et permet d’analyser le sens physique de l’adaptation nécessaire pour amener le modèle à simuler les températures de brillance (T[indice inférieur B) de la neige avec une erreur quadratique moyenne minimale de l’ordre de 13 K. Dans un contexte du couplage entre le modèle de neige de CLASS et DMRT-ML, un modèle d’évolution de la SSA est ensuite implémenté dans CLASS. Les SSA simulées par le module développé sont validées avec des données in situ basées sur la réflectance de la neige dans l’infrarouge à courte longueur d’onde pour différents types d’environnement. Au niveau de la contribution de la végétation, le modèle γ-ω a été étudié à partir de différentes bases de données (satellite, avion et au sol) en forêt boréale dense. L’étude montre l’importance de la considération de la diffusion (ω) pour l’estimation de l’émission de la végétation, paramètre auparavant généralement négligé aux hautes fréquences. Ensuite, des relations entre les transmissivités et certains paramètres structuraux de la forêt, dont l’indice de surface foliaire (LAI), ont été établies pour des forêts boréales en été. Des valeurs d’albédo de diffusion (ω) ainsi que les paramètres définissant la réflectivité du sol (QH) en forêt boréale ont aussi été inversées. Finalement, les simulations de T [indice inférieur] B issues du couplage du MTR (DMRT-ML, modèle γ-ω, et modèle atmosphérique) avec CLASS (dont les SSA simulées) ont été comparées avec les données AMSR-E sur une série temporelle continue de sept ans. Les premières comparaisons montrent une différence entre les paramètres de végétation (γ-ω) d’été et d’hiver, ainsi qu’une importante contribution des croûtes de glace dans la neige au signal. Les simulations du modèle ajusté montrent une bonne correspondance avec les observations d’AMSR-E (de l’ordre de 3 à 7 K selon la fréquence et la polarisation). Des tests de sensibilité montrent par contre une faible sensibilité du MTR/CLASS au SWE pour des forêts denses et des couverts nivaux épais. Le MTR-CLASS développé pourrait permettre l’assimilation de températures de brillance satellitaires en forêt boréale dans des systèmes opérationnels pour l’amélioration de paramètres de surface, dont la neige, dans les modèles météorologiques et climatiques