Atmospheric remote sensing and radiopropagation: from numerical modeling to spaceborne and terrestrial applications

The remote sensing of electromagnetic wave properties is probably the most viable and fascinating way to observe and study physical media, comprising our planet and its atmosphere, at the same time ensuring a proper continuity in the observations. Applications are manifold and the scientific community has been importantly studying and investing on new technologies, which would let us widen our knowledge of what surrounds us. This thesis aims at showing some novel techniques and corresponding applications in the field of the atmospheric remote sensing and radio-propagation, at both microwave and optical wavelengths. The novel Sun-tracking microwave radiometry technique is shown. The antenna noise temperature of a ground-based microwave radiometer is measured by alternately pointing toward-the-Sun and off-the-Sun while tracking it along its diurnal ecliptic. During clear sky the brightness temperature of the Sun disk emission at K and Ka frequency bands and in the under-explored millimeter-wave V and W bands can be estimated by adopting different techniques. Parametric prediction models for retrieving all-weather atmospheric extinction from ground-based microwave radiometers are tested and their accuracy evaluated. Moreover, a characterization of suspended clouds in terms of atmospheric path attenuation is presented, by exploiting a stochastic approach used to model the time evolution of the cloud contribution. A model chain for the prediction of the tropospheric channel for the downlink of interplanetary missions operating above Ku band is proposed. On top of a detailed description of the approach, the chapter presents the validation results and examples of the model-chain online operation. Online operation has already been tested within a feasibility study applied to the BepiColombo mission to Mercury operated by the European Space Agency (ESA) and by exploiting the Hayabusa-2 mission Ka-band data by the Japan Aerospace Exploration Agency (JAXA), thanks to the ESA cross-support service. A preliminary (and successful) validation of the model-chain has been carried out by comparing the simulated signal-to-noise ratio with the one received from Hayabusa-2. At the next ITU World Radiocommunication Conference 2019, Agenda Item 1.13 will address the identification and the possible additional allocation of radio-frequency spectrum to serve the future development of systems supporting the fifth generation of cellular mobile communications (5G). The potential impact of International Mobile Telecommunications (IMT) deployments is shown in terms of received radio frequency interference by ESA’s telecommunication links. Received interference can derive from several radio-propagation mechanisms, which strongly depend on atmospheric conditions, radio frequency, link availability, distance and path topography; at any time a single mechanism, or more than one may be present. Results are shown in terms of required separation distances, i.e. the minimum distance between the earth station and the IMT station ensuring that the protection criteria for the earth station are met

Deriving vertical total electron content maps from SMOS full polarimetric data to compensate the Faraday rotation effect

The Faraday rotation is a geophysical effect that causes a rotation of the electromagnetic field components emitted by the Earth when it propagates through the ionosphere. It depends on the vertical total electron content (VTEC) of the ionosphere, the geomagnetic field, and the frequency. For satellite measurements at the L band, this effect is not negligible and must be compensated for. This is the case of the Soil Moisture and Ocean Salinity (SMOS) mission, where the measured polarimetric brightness temperature must be corrected from the Faraday rotation effect before the retrieval of the geophysical parameters. The Faraday rotation angle (FRA) can be estimated using a theoretical formulation that makes use of external sources for the VTEC and the geomagnetic field. Alternatively, it can be continuously retrieved from the SMOS full-polarimetric data. However, this is not straightforward due to the relatively poor radiometric sensitivity (thermal noise) and accuracy (spatial bias) of its payload MIRAS (Microwave Interferometer Radiometer by Aperture Synthesis). In this thesis, a methodology for estimating the total electron content of the ionosphere by using an inversion procedure from the measured rotation angle has been developed. These SMOS VTEC maps are derived from SMOS measurements in the Extended Alias-Free Field of View (EAF-FoV) by applying spatio-temporal filtering techniques to mitigate the radiometric errors present in the full-polarimetric measured brightness temperatures. Systematic error patterns found in the Faraday rotation angle retrieval have been characterized along the mission and corrected. The methodology is independent, not only of external databases and forward models, but also of the target that is being measured. Eventually, these SMOS-derived VTEC maps can then be used in the SMOS level 2 processors to improve the geophysical retrievals. The impact of using these SMOS VTEC maps to correct the FRA in the SMOS mission instead of the commonly used VTEC data from GPS has also been assessed, particularly over ocean, where the ionospheric effect is stronger. This assessment has demonstrated improvements in the spatial biases, in the stability of the brightness temperatures (especially in the third Stokes parameter), and in the reduction of the latitudinal gradient present in the third Stokes parameters. All these quality indicators point to a better quality of the geophysical retrievals.La rotación de Faraday es un efecto geofísico que causa un giro en las componentes del campo electromagnético emitido por la Tierra cuando éste se propaga a través de la ionosfera. Ésta depende del contenido vertical total de electrones (VTEC) en la ionosfera, el campo geomagnético y la frecuencia. En las medidas de los satélites que operan en banda L, este efecto no es despreciable y se debe compensar. Este es el caso de la misión SMOS (Soil Moisture and Ocean Salinity), por lo que el efecto de Faraday se tiene que corregir en las medidas polarimétricas captadas por el instrumento antes de obtener parámetros geofísicos. El ángulo de rotación de Faraday (FRA) se puede estimar con una fórmula teórica que usa bases de datos externas para el VTEC y el campo geomagnético. Alternativamente, se puede obtener de una manera continua a partir de los datos polarimétricos de SMOS. Sin embargo, esto no se logra con un cálculo directo debido a la pobre sensibilidad radiométrica (ruido térmico) y a la baja precisión (sesgos espaciales) que presenta el instrumento MIRAS (Microwave Interferometer Radiometer by apertura Synthesis), que se encuentra a bordo del satélite. En esta tesis, se desarrolla una metodología para estimar el VTEC de la ionosfera usando un proceso inverso a partir del ángulo de rotación medido. Estos mapas de VTEC se derivan de medidas en todo el campo de visión extendido en donde no hay aliasing. Para mitigar los errores radiométricos en las temperaturas de brillo polarimétricas, se aplican técnicas de filtrados temporales y espaciales. En el ángulo de rotación de Faraday recuperado se detectaron errores sistemáticos. Estos se caracterizaron a lo largo de la misión y se corrigieron. La metodología es independiente, no solo de bases de datos externas y modelos de océano, sino también de la superficie medida. Estos mapas de VTEC derivados de los datos SMOS se pueden usar en el procesador de nivel 2 para mejorar las recuperaciones geofísicas. Se ha evaluado el impacto de usar estos mapas para corregir el FRA en la misión, en vez de los datos de VTEC que comúnmente se emplean (mapas provenientes de datos de GPS), particularmente sobre océano, en donde los efectos de la ionosfera son más críticos. Esta verificación ha demostrado mejoras en el sesgo espacial, en la estabilidad de las temperaturas de brillo (especialmente en el tercer parámetro de Stokes) y en la reducción del gradiente latitudinal presente en el tercer parámetro de Stokes. Todos estos indicadores de calidad apuntan a la obtención de parámetros geofísicos de mejor calidad.Postprint (published version

Contributions to radio frequency interference detection and mitigation in Earth observation

Radio Frequency Interference (RFI) is the most common problem for electronic measuring systems. The presence of those electromagnetic waves can harm the measurements taken from very sensitive instruments, like microwave radiometry or navigation systems. The accuracy and precision are compromised. A first step to mitigate those unwanted effects is to study the RFI properties. Different algorithms have been proposed to detect the interferences, but there is no method that works in all cases. The scope of this dissertation is the design, implementation and testing of different detection and mitigation methods in real-time. Performed surveys and characterization of RFI sources provide a great contribution to optimize the current mitigation techniques. In the mitigation area, two real-time hardware systems have been implemented: a wavelet denoise system to model the RFI and mitigate it, and a circuit to allow a navigation system to continue operational under the effects of a jammer.El problema més comú en els sistemes electrònics de mesura són les interferències electromagnètiques. La presència d'aquests senyals pot danyar les mesures preses per instruments molt sensibles, com radiòmetres de microones o sistemes de navegació. L'exactitud i precisió es veuen compromeses. El primer pas per mitigar aquests efectes no desitjats és estudiar les propietats de les interferències electromagnètiques. Diversos algoritmes han estat proposats per detectar interferències, però no hi ha mètode que funcioni bé en tots els casos . Aquest treball comprèn el disseny, implementació i comprovació de diferents mètodes de detecció i mitigació en temps real. Els estudis i caracterització de les fonts d'interferències són una gran contribució per a optimitzar les tècniques de mitigació actuals. En el tema de mitigació, dos sistemes en temps real han estat implementats en hardware: un sistema que utilitza wavelets per modelar la interferència i mitigar-la, i un circuit que permet a un sistema de navegació continuar funcionant sota els efectes d'un interferidor comercial ( jammer )

Radio frequency interference detection and mitigation techniques for navigation and Earth observation

Radio-Frequency Interference (RFI) signals are undesired signals that degrade or disrupt the performance of a wireless receiver. RFI signals can be troublesome for any receiver, but they are especially threatening for applications that use very low power signals. This is the case of applications that rely on the Global Navigation Satellite Systems (GNSS), or passive microwave remote sensing applications such as Microwave Radiometry (MWR) and GNSS-Reflectometry (GNSS-R). In order to solve the problem of RFI, RFI-countermeasures are under development. This PhD thesis is devoted to the design, implementation and test of innovative RFI-countermeasures in the fields of MWR and GNSS. In the part devoted to RFI-countermeasures for MWR applications, first, this PhD thesis completes the development of the MERITXELL instrument. The MERITXELL is a multi-frequency total-power radiometer conceived to be an outstanding platform to perform detection, characterization, and localization of RFI signals at the most common MWR imaging bands up to 92 GHz. Moreover, a novel RFI mitigation technique is proposed for MWR: the Multiresolution Fourier Transform (MFT). An assessment of the performance of the MFT has been carried out by comparison with other time-frequency mitigation techniques. According to the results, the MFT technique is a good trade-off solution among all other techniques since it can mitigate efficiently all kinds of RFI signals under evaluation. In the part devoted to RFI-countermeasures for GNSS and GNSS-R applications, first, a system for RFI detection and localization at GNSS bands is proposed. This system is able to detect RFI signals at the L1 band with a sensitivity of -108 dBm at full-band, and of -135 dBm for continuous wave and chirp-like signals when using the averaged spectrum technique. Besides, the Generalized Spectral Separation Coefficient (GSSC) is proposed as a figure of merit to evaluate the Signal-to-Noise Ratio (SNR) degradation in the Delay-Doppler Maps (DDMs) due to the external RFI effect. Furthermore, the FENIX system has been conceived as an innovative system for RFI detection and mitigation and anti-jamming for GNSS and GNSS-R applications. FENIX uses the MFT blanking as a pre-correlation excision tool to perform the mitigation. In addition, FENIX has been designed to be cross-GNSS compatible and RFI-independent. The principles of operation of the MFT blanking algorithm are assessed and compared with other techniques for GNSS signals. Its performance as a mitigation tool is proven using GNSS-R data samples from a real airborne campaign. After that, the main building blocks of the patented architecture of FENIX have been described. The FENIX architecture has been implemented in three real-time prototypes. Moreover, a simulator named FENIX-Sim allows for testing its performance under different jamming scenarios. The real-time performance of FENIX prototype has been tested using different setups. First, a customized VNA has been built in order to measure the transfer function of FENIX in the presence of several representative RFI/jamming signals. The results show how the power transfer function adapts itself to mitigate the RFI/jamming signal. Moreover, several real-time tests with GNSS receivers have been performed using GPS L1 C/A, GPS L2C, and Galileo E1OS. The results show that FENIX provides an extra resilience against RFI and jamming signals up to 30 dB. Furthermore, FENIX is tested using a real GNSS timing setup. Under nominal conditions, when no RFI/jamming signal is present, a small additional jitter on the order of 2-4 ns is introduced in the system. Besides, a maximum bias of 45 ns has been measured under strong jamming conditions (-30 dBm), which is acceptable for current timing systems requiring accuracy levels of 100 ns. Finally, the design of a backup system for GNSS in tracking applications that require high reliability against RFI and jamming attacks is proposed.Les interferències de radiofreqüència (RFI) són senyals no desitjades que degraden o interrompen el funcionament dels receptors sense fils. Les RFI poden suposar un problema per qualsevol receptor, però són especialment amenaçadores per les a aplicacions que fan servir senyals de molt baixa potència. Aquest és el cas de les aplicacions que depenen dels sistemes mundials de navegació per satèl·lit (GNSS) o de les aplicacions de teledetecció passiva de microones, com la radiometria de microones (MWR) i la reflectometria GNSS (GNSS-R). Per combatre aquest problema, sistemes anti-RFI s'estan desenvolupament actualment. Aquesta tesi doctoral està dedicada al disseny, la implementació i el test de sistemes anti-RFI innovadors en els camps de MWR i GNSS. A la part dedicada als sistemes anti-RFI en MWR, aquesta tesi doctoral completa el desenvolupament de l'instrument MERITXELL. El MERITXELL és un radiòmetre multifreqüència concebut com una plataforma excepcional per la detecció, caracterització i localització de RFI a les bandes de MWR més utilitzades per sota dels 92 GHz. A més a més, es proposa una nova tècnica de mitigació de RFI per MWR: la Transformada de Fourier amb Multiresolució (MFT). El funcionament de la MFT s'ha comparat amb el d'altres tècniques de mitigació en els dominis del temps i la freqüència. D'acord amb els resultats obtinguts, la MFT és una bona solució de compromís entre les altres tècniques, ja que pot mitigar de manera eficient tots els tipus de senyals RFI considerats. A la part dedicada als sistemes anti-RFI en GNSS i GNSS-R, primer es proposa un sistema per a la detecció i localització de RFI a les bandes GNSS. Aquest sistema és capaç de detectar senyals RFI a la banda L1 amb una sensibilitat de -108 dBm a tota la banda, i de -135 dBm per a senyals d'ona contínua i chirp fen un mitjana de l'espectre. A més a més, el Coeficient de Separació Espectral Generalitzada (GSSC) es proposa com una mesura per avaluar la degradació de la relació senyal a soroll (SNR) en els Mapes de Delay-Doppler (DDM) a causa del impacte de les RFI. La major contribució d'aquesta tesi doctoral és el sistema FENIX. FENIX és un sistema innovador de detecció i mitigació de RFI i inhibidors de freqüència per aplicacions GNSS i GNSS-R. FENIX utilitza la MFT per eliminar la interferència abans del procés de correlació amb el codi GNSS independentment del tipus de RFI. L'algoritme de mitigació de FENIX s'ha avaluat i comparat amb altres tècniques i els principals components de la seva arquitectura patentada es descriuen. Finalment, un simulador anomenat FENIX-Sim permet avaluar el seu rendiment en diferents escenaris d'interferència. El funcionament en temps real del prototip FENIX ha estat provat utilitzant diferents mètodes. En primer lloc, s'ha creat un analitzador de xarxes per a mesurar la funció de transferència del FENIX en presència de diverses RFI representatives. Els resultats mostren com la funció de transferència s'adapta per mitigar el senyal interferent. A més a més, s'han realitzat diferents proves en temps real amb receptors GNSS compatibles amb els senyals GPS L1 C/A, GPS L2C i Galileo E1OS. Els resultats mostren que FENIX proporciona una resistència addicional contra les RFI i els senyals dels inhibidors de freqüència de fins a 30 dB. A més a més, FENIX s'ha provat amb un sistema comercial de temporització basat en GNSS. En condicions nominals, sense RFI, FENIX introdueix un petit error addicional de tan sols 2-4 ns. Per contra, el biaix màxim mesurat en condicions d'alta interferència (-30 dBm) és de 45 ns, el qual és acceptable per als sistemes de temporització actuals que requereixen nivells de precisió d'uns 100 ns. Finalment, es proposa el disseny d'un sistema robust de seguiment, complementari als GNSS, per a aplicacions que requereixen alta fiabilitat contra RFI.Postprint (published version

A Database of Microwave Single Scattering Properties of Ice Hydrometeors

Microwave remote sensing by satellites is important for global observations of ice hydrometeors. Interpretation of the measurements requires sufficiently accurate knowledge of hydrometeors’ interaction with photons, i.e.\ua0 article scattering and absorption. This presents a challenge for several reasons. Liquid hydrometeors can typically be modelled by spheroids, while the shapes of ice hydrometeors are known to be significantly more complex and variable. Also, the shapes can from a remote sensing perspective generally not be known exactly, as they vary from case to case. Finally, calculating the light scattering properties is challenging and computationally costly.This thesis presents work related to recent efforts in improving the representation of light scattering by ice hydrometeors. A new single scattering database is presented, which includes 34 frequencies in between 1 and 874 GHz, and supports both passive and active microwave applications. A total of 34 different particle models were included, ranging from pristine crystals to aggregates. Complete random orientation is assumed throughout, slightly limiting its usefulness with respect to polarimetric measurements. Most aggregates were generated through simulation of aggregation, by letting particles collide randomly. The database can be considered the most extensive of this type to date, and future versions are intended to include oriented and melting particles. The general intention is to aid existing and future satellite retrievals, and satellite data assimilation into weather prediction models, all requiring accurate modelling of measured radiances. Special attention has been given to the upcoming Ice Cloud Imager (ICI), part of Europe’s next generation of weather satellites.Using the aggregation simulation tools developed for the database, a more dedicated case study was performed, which looked at the impact of different aggregate shape parameters on the resulting scattering properties. Both the amount and aspect ratio of the aggregate constituent crystals was found to have a high impact on both extinction (183, 325 and 664 GHz) and back-scattering (13, 36 and 94 GHz). Effective density and aerodynamic area had a high impact as well. Calculated radar triple frequency signatures were seen to clearly depend on the particle shape, consistent with previous studies. Overall, the results indicate that the particle shape should be considered in both passive and active applications above 13 GHz, and future database development will consider this. A potential application is also retrieval of ice particle shape through remote sensing

Millimeter and sub-millimeter wave radiometers for atmospheric remote sensing from CubeSat platforms

2018 Fall.Includes bibliographical references.To view the abstract, please see the full text of the document

Analyse de la modélisation de l'émission multi-fréquences micro-onde des sols et de la neige, incluant les croutes de glace à l'aide du modèle Microwave Emission Model of Layered Snowpacks (MEMLS).

Résumé : L'étude du couvert nival est essentielle afin de mieux comprendre les processus climatiques et hydrologiques. De plus, avec les changements climatiques observés dans l'hémisphère nord, des événements de dégel-regel ou de pluie hivernale sont de plus en plus courants et produisent des croutes de glace dans le couvert nival affectant les moeurs des communautés arctiques en plus de menacer la survie de la faune arctique. La télédétection micro-ondes passives (MOP) démontre un grand potentiel de caractérisation du couvert nival. Toutefois, a fin de bien comprendre les mesures satellitaires, une modélisation adéquate du signal est nécessaire. L'objectif principal de cette thèse est d'analyser le transfert radiatif (TR) MOP des sols, de la neige et de la glace a fin de mieux caractériser les propriétés géophysiques du couvert nival par télédétection. De plus, un indice de détection des croutes de glace par télédétection MOP a été développé. Pour ce faire, le modèle Microwave Emission Model of Layered Snowpacks (MEMLS) a été étudié et calibré afin de minimiser les erreurs des températures de brillance simulées en présences de croutes de glace. La première amélioration faite à la modélisation du TR MOP de la neige a été la caractérisation de la taille des grains de neige. Deux nouveaux instruments, utilisant la réflectance dans le proche infrarouge, ont été développés afin de mesurer la surface spécifique de la neige (SSA). Il a été démontré que la SSA est un paramètre plus précis et plus objectif pour caractériser la taille des grains de neige. Les deux instruments ont démontré une incertitude de 10% sur la mesure de la SSA. De plus, la SSA a été calibré pour la modélisation MOP a n de minimiser l'erreur sur la modélisation de la température de brillance. Il a été démontré qu'un facteur multiplicatif [phi] = 1.3 appliqué au paramètre de taille des grains de neige dans MEMLS, paramètre dérivé de la SSA, est nécessaire afin de minimiser l'erreur des simulations. La deuxième amélioration apportée à la modélisation du TR MOP a été l'estimation de l'émission du sol. Des mesures radiométriques MOP in-situ ainsi que des profils de températures de sols organiques arctiques gelés ont été acquis et caractérisés a fin de simuler l'émission MOP de ces sols. Des constantes diélectriques effectives à 10.7, 19 et 37 GHz ainsi qu'une rugosité de surface effective des sols ont été déterminés pour simuler l'émission des sols. Une erreur quadratique moyenne (RMSE) de 4.65 K entre les simulations et les mesures MOP a été obtenue. Suite à la calibration du TR MOP du sol et de la neige, un module de TR de la glace a été implémenté dans MEMLS. Avec ce nouveau module, il a été possible de démontré que l'approximation de Born améliorée, déjà implémenté dans MEMLS, pouvait être utilisé pour simuler des croutes de glace pure à condition que la couche de glace soit caractérisée par une densité de 917 kg m[indice supérieur _3] et une taille des grains de neige de 0 mm. Il a aussi été démontré que, pour des sites caractérisés par des croutes de glace, les températures de brillances simulées des couverts de neige avec des croutes de glace ayant les propriétés mesurées in-situ (RMSE=11.3 K), avaient une erreur similaire aux températures de brillances simulées des couverts de neige pour des sites n'ayant pas de croutes de glace (RMSE=11.5 K). Avec le modèle MEMLS validé pour la simulation du TR MOP du sol, de la neige et de la glace, un indice de détection des croutes de glace par télédétection MOP a été développé. Il a été démontré que le ratio de polarisation (PR) était très affecté par la présence de croutes de glace dans le couvert de neige. Avec des simulations des PR à 10.7, 19 et 37 GHz sur des sites mesurés à Churchill (Manitoba, Canada), il a été possible de déterminer des seuils entre la moyenne hivernale des PR et les valeurs des PR mesurés indiquant la présence de croutes de glace. Ces seuils ont été appliqués sur une série temporelle de PR de 33 hivers d'un pixel du Nunavik (Québec, Canada) où les conditions de sols étaient similaires à ceux observés à Churchill. Plusieurs croutes de glace ont été détectées depuis 1995 et les mêmes événements entre 2002 et 2009 que (Roy, 2014) ont été détectés. Avec une validation in-situ, il serait possible de confirmer ces événements de croutes de glace mais (Roy, 2014) a démontré que ces événements ne pouvaient être expliqués que par la présence de croutes de glace dans le couvert de neige. Ces mêmes seuils sur les PR ont été appliqués sur un pixel de l'Île Banks (Territoires du Nord-Ouest, Canada). L'événement répertorié par (Grenfell et Putkonen, 2008) a été détecté. Plusieurs autres événements de croutes de glace ont été détectés dans les années 1990 et 2000 avec ces seuils. Tous ces événements ont suivi une période où les températures de l'air étaient près ou supérieures au point de congélation et sont rapidement retombées sous le point de congélation. Les températures de l'air peuvent être utilisées pour confirmer la possibilité de présence de croutes de glace mais seul la validation in-situ peut définitivement confirmer la présence de ces croutes.Abstract : Snow cover studies are essential to better understand climatic and hydrologic processes. With recent climate change observed in the northern hemisphere, more frequent rain-on-snow and meltrefreeze events have been reported, which affect the habits of the northern comunities and the survival of arctique wildlife. Passive microwave remote sensing has proven to be a great tool to characterize the state of snow cover. Nonetheless, proper modeling of the microwave signal is needed in order to understand how the parameters of the snowpack affect the measured signal. The main objective of this study is to analyze the soil, snow and ice radiative transfer in order to better characterize snow cover properties and develop an ice lens detection index with satellite passive microwave brightness temperatures. To do so, the passive microwave radiative transfer modeling of the Microwave Emission Model of Layered Snowpacks (MEMLS) was improved in order to minimize the errors on the brightness temperature simulations in the presence of ice lenses. The first improvement to passive microwave radiative transfer modeling of snow made was the snow grain size parameterization. Two new instruments, based on short wave infrared reflectance to measure the snow specific surface area (SSA) were developed. This parameter was shown to be a more accurate and objective to characterize snow grain size. The instruments showed an uncertainty of 10% to measure the SSA of snow. Also, the SSA of snow was calibrated for passive microwave modeling in order to reduce the errors on the simulated brightness temperatures. It was showed that a correction factor of φ = 1.3 needed to be applied to the grain size parameter of MEMLS, obtain through the SSA measurements, to minimize the simulation error. The second improvement to passive microwave radiative transfer modeling was the estimation of passive microwave soil emission. In-situ microwave measurements and physical temperature profiles of frozen organic arctic soils were acquired and characterized to improve the modeling of the soil emission. Effective permittivities at 10.7, 19 and 37 GHz and effective surface roughness were determined for this type of soil and the soil brightness temperature simulations were obtain with a minimal root mean square error (RMSE) of 4.65K. With the snow grain size and soil contributions to the emitted brightness temperature optimized, it was then possible to implement a passive microwave radiative transfer module of ice into MEMLS. With this module, it was possible to demonstrate that the improved Born approximation already implemented in MEMLS was equivalent to simulating a pure ice lens when the density of the layer was set to 917 kg m−3 and the grain size to 0 mm. This study also showed that by simulating ice lenses within the snow with there measured properties, the RMSE of the simulations (RMSE= 11.3 K) was similar to the RMSE for simulations of snowpacks where no ice lenses were measured (only snow, RMSE= 11.5 K). With the validated MEMLS model for snowpacks with ice lenses, an ice index was created. It is shown here that the polarization ratio (PR) was strongly affected by the presence of ice lenses within the snowpack. With simulations of the PR at 10.7, 19 and 37 GHz from measured snowpack properties in Chucrhill (Manitoba, Canada), thresholds between the measured PR and the mean winter PR were determined to detect the presence of ice within the snowpack. These thresholds were applied to a timeseries of nearly 34 years for a pixel in Nunavik (Quebec, Canada) where the soil surface is similar to that of the Churchill site. Many ice lenses are detected since 1995 with these thresholds and the same events as Roy (2014) were detected. With in-situ validation, it would be possible to confirm the precision of these thresholds but Roy (2014) showed that these events can not be explained by anything else than the presence of an ice layer within the snowpack. The same thresholds were applied to a pixel on Banks island (North-West Territories, Canada). The 2003 event that was reported by Grenfell et Putkonen (2008) was detected by the thresholds. Other events in the years 1990 and 2000’s were detected with these thresholds. These events all follow periods where the air temperature were warm and were followed by a quick drop in air temperature which could be used to validate the presence of ice layer within the snowpack. Nonetheless, without in-situ validation, these events can not be confirmed