Coupled land surface and radiative transfer models for the analysis of passive microwave satellite observations

Soil moisture is one of the key variables controlling the water and energy exchanges between Earth’s surface and the atmosphere. Therefore, remote sensing based soil moisture information has potential applications in many disciplines. Besides numerical weather forecasting and climate research these include agriculture and hydrologic applications like flood and drought forecasting. The first satellite specifically designed to deliver operational soil moisture products, SMOS (Soil Moisture and Ocean Salinity), was launched 2009 by the European Space Agency (ESA). SMOS is a passive microwave radiometer working in the L-band of the microwave domain, corresponding to a frequency of roughly 1.4 GHz and relies on a new concept. The microwave radiation emitted by the Earth’s surface is measured as brightness temperatures in several look angles. A radiative transfer model is used in an inversion algorithm to retrieve soil moisture and vegetation optical depth, a measure for the vegetation attenuation of the soil’s microwave emission. For the application of passive microwave remote sensing products a proper validation and uncertainty assessment is essential. As these sensors have typical spatial resolutions in the order of 40 – 50 km, a validation that relies solely on ground measurements is costly and labour intensive. Here, environmental modelling can make a valuable contribution. Therefore the present thesis concentrates on the question which contribution coupled land surface and radiative transfer models can make to the validation and analysis of passive microwave remote sensing products. The objective is to study whether it is possible to explain known problems in the SMOS soil moisture products and to identify potential approaches to improve the data quality. The land surface model PROMET (PRocesses Of Mass and Energy Transfer) and the radiative transfer model L-MEB (L-band microwave emission of the Biosphere) are coupled to simulate land surface states, e.g. temperatures and soil moisture, and the resulting microwave emission. L-MEB is also used in the SMOS soil moisture processor to retrieve soil moisture and vegetation optical depth simultaneously from the measured microwave emission. The study area of this work is the Upper Danube Catchment, located mostly in Southern Germany. Since model validation is essential if model data are to be used as reference, both models are validated on different spatial scales with measurements. The uncertainties of the models are quantified. The root mean squared error between modelled and measured soil moisture at several measuring stations on the point scale is 0.065 m3/m3. On the SMOS scale it is 0.039 m3/m3. The correlation coefficient on the point scale is 0.84. As it is essential for the soil moisture retrieval from passive microwave data that the radiative transfer modelling works under local conditions, the coupled models are used to assess the radiative transfer modelling with L-MEB on the local and SMOS scales in the Upper Danube Catchment. In doing so, the emission characteristics of rape are described for the first time and the soil moisture retrieval abilities of L-MEB are assessed with a newly developed LMEB parameterization. The results show that the radiative transfer modelling works well under most conditions in the study area. The root mean squared error between modelled and airborne measured brightness temperatures on the SMOS scale is less than 6 – 9 K for the different look angles. The coupled models are used to analyse SMOS brightness temperatures and vegetation optical depth data in the Upper Danube Catchment in Southern Germany. Since the SMOS soil moisture products are degraded in Southern Germany and in different other parts of the world these analyses are used to narrow down possible reasons for this. The thorough analysis of SMOS brightness temperatures for the year 2011 reveals that the quality of the measurements is degraded like in the SMOS soil moisture product. This points towards radio frequency interference problems (RFI), that are known, but have not yet been studied thoroughly. This is consistent with the characteristics of the problems observed in the SMOS soil moisture products. In addition to that it is observed that the brightness temperatures in the lower look angles are less reliable. This finding could be used to improve the brightness temperature filtering before the soil moisture retrieval. An analysis of SMOS optical depth data in 2011 reveals that this parameter does not contain valuable information about vegetation. Instead, an unexpected correlation with SMOS soil moisture is found. This points towards problems with the SMOS soil moisture retrieval, possibly under the influence of RFI. The present thesis demonstrates that coupled land surface and radiative transfer models can make a valuable contribution to the validation and analysis of passive microwave remote sensing products. The unique approach of this work incorporates modelling with a high spatial and temporal resolution on different scales. This makes detailed process studies on the local scale as well as analyses of satellite data on the SMOS scale possible. This could be exploited for the validation of future satellite missions, e.g. SMAP (Soil Moisture Active and Passive) which is currently being prepared by NASA (National Aeronautics and Space Administration). Since RFI seems to have a considerable influence on the SMOS data due to the gained insights and the quality of the SMOS products is very good in other parts of the world, the RFI containment and mitigation efforts carried out since the launch of SMOS should be continued

Coupled land surface and radiative transfer models for the analysis of passive microwave satellite observations

Soil moisture is one of the key variables controlling the water and energy exchanges between Earth’s surface and the atmosphere. Therefore, remote sensing based soil moisture information has potential applications in many disciplines. Besides numerical weather forecasting and climate research these include agriculture and hydrologic applications like flood and drought forecasting. The first satellite specifically designed to deliver operational soil moisture products, SMOS (Soil Moisture and Ocean Salinity), was launched 2009 by the European Space Agency (ESA). SMOS is a passive microwave radiometer working in the L-band of the microwave domain, corresponding to a frequency of roughly 1.4 GHz and relies on a new concept. The microwave radiation emitted by the Earth’s surface is measured as brightness temperatures in several look angles. A radiative transfer model is used in an inversion algorithm to retrieve soil moisture and vegetation optical depth, a measure for the vegetation attenuation of the soil’s microwave emission. For the application of passive microwave remote sensing products a proper validation and uncertainty assessment is essential. As these sensors have typical spatial resolutions in the order of 40 – 50 km, a validation that relies solely on ground measurements is costly and labour intensive. Here, environmental modelling can make a valuable contribution. Therefore the present thesis concentrates on the question which contribution coupled land surface and radiative transfer models can make to the validation and analysis of passive microwave remote sensing products. The objective is to study whether it is possible to explain known problems in the SMOS soil moisture products and to identify potential approaches to improve the data quality. The land surface model PROMET (PRocesses Of Mass and Energy Transfer) and the radiative transfer model L-MEB (L-band microwave emission of the Biosphere) are coupled to simulate land surface states, e.g. temperatures and soil moisture, and the resulting microwave emission. L-MEB is also used in the SMOS soil moisture processor to retrieve soil moisture and vegetation optical depth simultaneously from the measured microwave emission. The study area of this work is the Upper Danube Catchment, located mostly in Southern Germany. Since model validation is essential if model data are to be used as reference, both models are validated on different spatial scales with measurements. The uncertainties of the models are quantified. The root mean squared error between modelled and measured soil moisture at several measuring stations on the point scale is 0.065 m3/m3. On the SMOS scale it is 0.039 m3/m3. The correlation coefficient on the point scale is 0.84. As it is essential for the soil moisture retrieval from passive microwave data that the radiative transfer modelling works under local conditions, the coupled models are used to assess the radiative transfer modelling with L-MEB on the local and SMOS scales in the Upper Danube Catchment. In doing so, the emission characteristics of rape are described for the first time and the soil moisture retrieval abilities of L-MEB are assessed with a newly developed LMEB parameterization. The results show that the radiative transfer modelling works well under most conditions in the study area. The root mean squared error between modelled and airborne measured brightness temperatures on the SMOS scale is less than 6 – 9 K for the different look angles. The coupled models are used to analyse SMOS brightness temperatures and vegetation optical depth data in the Upper Danube Catchment in Southern Germany. Since the SMOS soil moisture products are degraded in Southern Germany and in different other parts of the world these analyses are used to narrow down possible reasons for this. The thorough analysis of SMOS brightness temperatures for the year 2011 reveals that the quality of the measurements is degraded like in the SMOS soil moisture product. This points towards radio frequency interference problems (RFI), that are known, but have not yet been studied thoroughly. This is consistent with the characteristics of the problems observed in the SMOS soil moisture products. In addition to that it is observed that the brightness temperatures in the lower look angles are less reliable. This finding could be used to improve the brightness temperature filtering before the soil moisture retrieval. An analysis of SMOS optical depth data in 2011 reveals that this parameter does not contain valuable information about vegetation. Instead, an unexpected correlation with SMOS soil moisture is found. This points towards problems with the SMOS soil moisture retrieval, possibly under the influence of RFI. The present thesis demonstrates that coupled land surface and radiative transfer models can make a valuable contribution to the validation and analysis of passive microwave remote sensing products. The unique approach of this work incorporates modelling with a high spatial and temporal resolution on different scales. This makes detailed process studies on the local scale as well as analyses of satellite data on the SMOS scale possible. This could be exploited for the validation of future satellite missions, e.g. SMAP (Soil Moisture Active and Passive) which is currently being prepared by NASA (National Aeronautics and Space Administration). Since RFI seems to have a considerable influence on the SMOS data due to the gained insights and the quality of the SMOS products is very good in other parts of the world, the RFI containment and mitigation efforts carried out since the launch of SMOS should be continued

Synergistic optical and microwave remote sensing approaches for soil moisture mapping at high resolution

Aplicat embargament des de la data de defensa fins al dia 1 d'octubre de 2022Soil moisture is an essential climate variable that plays a crucial role linking the Earth’s water, energy, and carbon cycles. It is responsible for the water exchange between the Earth’s surface and the atmosphere, and provides key information about soil evaporation, plant transpiration, and the allocation of precipitation into runoff, surface flow and infiltration. Therefore, an accurate estimation of soil moisture is needed to enhance our current climate and meteorological forecasting skills, and to improve our current understanding of the hydrological cycle and its extremes (e.g., droughts and floods). L-band Microwave passive and active sensors have been used during the last decades to estimate soil moisture, since there is a strong relationship between this variable and the soil dielectric properties. Currently, there are two operational L-band missions specifically devoted to globally measure soil moisture: the ESA’s Soil Moisture and the Ocean Salinity (SMOS), launched in November 2009; and the NASA’s Soil Moisture Active Passive (SMAP), launched in January 2015. The spatial resolution of the SMOS and SMAP radiometers, in the order of tens of kilometers (~40 km), is adequate for global applications. However, to fulfill the needs of a growing number of applications at local or regional scale, higher spatial detail (< 1 km) is required. To bridge this gap and improve the spatial resolution of the soil moisture maps, a variety of spatial enhancement or spatial (sub-pixel) disaggregation approaches have been proposed. This Ph.D. Thesis focuses on the study of the Earth’s surface soil moisture from remotely sensed observations. This work includes the implementation of several soil moisture retrieval techniques and the development, implementation, validation and comparison of different spatial enhancement or downscaling techniques, applied at local, regional, and continental scale. To meet these objectives, synergies between several active/passive microwave sensors (SMOS, SMAP and Sentinel-1) and optical/thermal sensors (MODIS) have been explored. The results are presented as follows: - Spatially consistent downscaling approach for SMOS using an adaptive moving window A passive microwave/optical downscaling algorithm for SMOS is proposed to obtain fine-scale soil moisture maps (1 km) from the native resolution (~40 km) of the instrument. This algorithm introduces the concept of a shape-adaptive window as a central improvement of the disaggregation technique presented by Piles et al. (2014), allowing its application at continental scales. - Assessment of multi-scale SMOS and SMAP soil moisture products across the Iberian Peninsula The temporal and spatial characteristics of SMOS and SMAP soil moisture products at coarse- and fine-scales are assessed in order to learn about their distinct features and the rationale behind them, tracing back to the physical assumptions they are based upon. - Impact of incidence angle diversity on soil moisture retrievals at coarse and fine scales An incidence angle (32.5°, 42.5° and 52.5°)-adaptive calibration of radiative transfer effective parameters single scattering albedo and soil roughness has been carried out, highlighting the importance of such parameterization to accurately estimate soil moisture at coarse-resolution. Then, these parameterizations are used to examine the potential application of a physically-based active-passive downscaling approach to upcoming microwave missions, namely CIMR, ROSE-L and Sentinel-1 Next Generation. Soil moisture maps obtained for the Iberian Peninsula at the three different angles, and at coarse and fine scales are inter-compared using in situ measurements and model data as benchmarks.La humedad del suelo es una variable climática esencial que juega un papel crucial en la relación de los ciclos del agua, la energía y el carbono de la Tierra. Es responsable del intercambio de agua entre la superficie de la Tierra y la atmósfera, y proporciona información crucial sobre la evaporación del suelo, la transpiración de las plantas y la distribución de la precipitación en escorrentía, flujo superficial e infiltración. Por lo tanto, es necesaria una estimación precisa de la humedad del suelo para mejorar las predicciones climáticas y meteorológicas, y comprender mejor el ciclo hidrológico y sus extremos (v.g., sequías e inundaciones). Los sensores pasivos y activos en banda L se han usado durante las últimas décadas para estimar la humedad del suelo debido a la relación directa que existe entre esta variable y las propiedades dieléctricas del suelo. Actualmente, hay dos misiones operativas en banda L específicamente dedicadas a medir la humedad del suelo a escala global: la misión Soil Moisture and Ocean Salinity (SMOS) de la ESA, lanzada en noviembre de 2009; y la misión Soil Moisture Active Passive (SMAP) de la NASA, lanzada en enero de 2015. La resolución espacial de los radiómetros SMOS y SMAP, del orden de unas decenas de kilómetros (~40 km), es adecuada para aplicaciones a escala global. Sin embargo, para satisfacer las necesidades de un número creciente de aplicaciones a escala local o regional, se requiere más detalle espacial (<1 km). Para solventar esta limitación y mejorar la resolución espacial de los mapas de humedad, se han propuesto diferentes técnicas de mejora o desagregación espacial. Esta Tesis se centra en el estudio de la humedad de la superficie terrestre a partir de datos obtenidos a través de teledetección. Este trabajo incluye la implementación de distintos algoritmos de recuperación de la humedad del suelo y el desarrollo, implementación, validación y comparación de distintas técnicas de desagregación, aplicadas a escala local, regional y continental. Para cumplir estos objetivos, se han explorado sinergias entre diferentes sensores de microondas activos/pasivos (SMOS, SMAP y Sentinel-1) y sensores ópticos/térmicos. Los resultados se presentan de la siguiente manera: - Técnica de desagregación espacialmente consistente, basada en una ventana móvil adaptativa, aplicada a los datos SMOS Se propone un algoritmo de desagregación del píxel basado en datos obtenidos de medidas radiométricas de microondas en banda L y datos ópticos, para mejorar la resolución espacial de los mapas de humedad del suelo desde la resolución nativa del instrumento (~40 km) hasta resoluciones de 1 km. El algoritmo introduce el concepto de una ventana de contorno adaptativo, como mejora principal sobre la técnica de desagregación presentada en Piles et al. (2014), permitiendo su implementación a escala continental. - Análisis multiescalar de productos de humedad del suelo SMAP y SMOS sobre la Península Ibérica Se han evaluado las características temporales y espaciales de distintos productos de humedad del suelo SMOS y SMAP, a baja y a alta resolución, para conocer sus características distintivas y comprender las razones de sus diferencias. Para ello, ha sido necesario rastrear los supuestos físicos en los que se basan. - Impacto del ángulo de incidencia en la recuperación de la humedad del suelo a baja y a alta resolución Se ha llevado a cabo una calibración adaptada al ángulo de incidencia (32.5°, 42.5° y 52.5°) de los parámetros efectivos, albedo de dispersión simple y rugosidad del suelo, descritos en el modelo de transferencia radiativa � − �, incidiendo en la importancia de esta parametrización para estimar la humedad del suelo de forma precisa a baja resolución. El resultado de las mismas se ha utilizado para estudiar la potencial aplicación de un algoritmo activo/pasivo de desagregación basado en la física para las próximas misiones de microondas, llamadas CIMR, ROSE-L y Sentinel-1 Next Generation. Los mapas de humedad recuperados a los tres ángulos de incidencia, tanto a baja como a alta resolución, se han obtenido para la Península Ibérica y se han comparado entre ellos usando como referencia mediciones de humedad in situ.Postprint (published version

Assessing the relationship between microwave vegetation optical depth and gross primary production

At the global scale, the uptake of atmospheric carbon dioxide by terrestrial ecosystems through photosynthesis is commonly estimated through vegetation indices or biophysical properties derived from optical remote sensing data. Microwave observations of vegetated areas are sensitive to different components of the vegetation layer than observations in the optical domain and may therefore provide complementary information on the vegetation state, which may be used in the estimation of Gross Primary Production (GPP). However, the relation between GPP and Vegetation Optical Depth (VOD), a biophysical quantity derived from microwave observations, is not yet known. This study aims to explore the relationship between VOD and GPP. VOD data were taken from different frequencies (L-, C-, and X-band) and from both active and passive microwave sensors, including the Advanced Scatterometer (ASCAT), the Soil Moisture Ocean Salinity (SMOS) mission, the Advanced Microwave Scanning Radiometer for Earth Observation System (AMSR-E) and a merged VOD data set from various passive microwave sensors. VOD data were compared against FLUXCOM GPP and Solar-Induced chlorophyll Fluorescence (SIF) from the Global Ozone Monitoring Experiment-2 (GOME-2). FLUXCOM GPP estimates are based on the upscaling of flux tower GPP observations using optical satellite data, while SIF observations present a measure of photosynthetic activity and are often used as a proxy for GPP. For relating VOD to GPP, three variables were analyzed: original VOD time series, temporal changes in VOD (ΔVOD), and positive changes in VOD (ΔVOD≥0). Results show widespread positive correlations between VOD and GPP with some negative correlations mainly occurring in dry and wet regions for active and passive VOD, respectively. Correlations between VOD and GPP were similar or higher than between VOD and SIF. When comparing the three variables for relating VOD to GPP, correlations with GPP were higher for the original VOD time series than for ΔVOD or ΔVOD≥0 in case of sparsely to moderately vegetated areas and evergreen forests, while the opposite was true for deciduous forests. Results suggest that original VOD time series should be used jointly with changes in VOD for the estimation of GPP across biomes, which may further benefit from combining active and passive VOD data

Statistical analysis and combination of active and passive microwave remote sensing methods for soil moisture retrieval

Knowledge about soil moisture and its spatio-temporal dynamics is essential for the improvement of climate and hydrological modeling, including drought and flood monitoring and forecasting, as well as weather forecasting models. In recent years, several soil moisture products from active and passive microwave remote sensing have become available with high temporal resolution and global coverage. Thus, the validation and evaluation of spatial and temporal soil moisture patterns are of great interest, for improving soil moisture products as well as for their proper use in models or other applications. This thesis analyzes the different accuracy levels of global soil moisture products and identifies the major influencing factors on this accuracy based on a small catchment example. Furthermore, on global scale, structural differences betweenthe soil moisture products were investigated. This includes in particular the representation of spatial and temporal patterns, as well as a general scaling law of soil moisture variability with extent scale. The results of the catchment scale as well as the global scale analyses identified vegetation to have a high impact on the accuracy of remotely sensed soil moisture products. Therefore, an improved method to consider vegetation characteristics in pasive soil moisture retrieval from active radar satellite data was developed and tested. The knowledge gained by this thesis will contribute to improve soil moisture retrieval of current and future microwave remote sensors (e.g. SMOS or SMAP)

The European heat wave 2003: early indicators from multisensoral microwave remote sensing?

An extreme heat wave affected large parts of Europe in 2003 with severe socioeconomic impacts. The extreme warm weather conditions lasted over a couple of months with positive temperature anomalies of 5°C for large parts of Europe. Simulations of the event using regional climate models revealed that a pronounced precipitation deficit in the beginning of the year, together with an early onset of the vegetation, resulted in a severe deficit of the soil water content. This amplified the course of the heat wave due to an increasing sensible heat flux from the land surface. The monitoring of temporal and spatial dynamics of soil water content can be accomplished using remote-sensing-based techniques. The present paper addresses the question whether there have been early indicators for the low soil water content using either physically based land surface modeling or remote-sensing-based monitoring techniques. The course of the spring surface soil moisture evolution is investigated using observations from two different microwave remote sensing sensors. An intercomparison of the high-resolution data from the European ENVISAT satellite and coarse resolution data from the AMSR-E mission is made. Remote-sensing-derived soil moisture products are compared against the results from a deterministic land surface model. The model enables to relate the year 2003 anomalies to a long-term (30 years) climatology. The year 2003 remote sensing derived soil moisture dynamics is compared against a multiyear climatology. The results reveal a negative surface soil moisture anomaly in 2003. The results indicate that there was in general potential to monitor the spatial and temporal dimensions of the low surface soil water content early in 2003 using remote sensing techniques. Both remote sensing data sets indicate a consistent soil moisture decrease in early 2003. A good agreement between the observed surface soil moisture and soil moisture simulations from a land surface process model was found. An outlook to the use of remote-sensing-based soil moisture estimates for large-scale monitoring of surface soil moisture trends is given. Copyright 2009 by the American Geophysical Union

Satellite and in situ observations for advancing global Earth surface modelling: a review

In this paper, we review the use of satellite-based remote sensing in combination with in situ data to inform Earth surface modelling. This involves verification and optimization methods that can handle both random and systematic errors and result in effective model improvement for both surface monitoring and prediction applications. The reasons for diverse remote sensing data and products include (i) their complementary areal and temporal coverage, (ii) their diverse and covariant information content, and (iii) their ability to complement in situ observations, which are often sparse and only locally representative. To improve our understanding of the complex behavior of the Earth system at the surface and sub-surface, we need large volumes of data from high-resolution modelling and remote sensing, since the Earth surface exhibits a high degree of heterogeneity and discontinuities in space and time. The spatial and temporal variability of the biosphere, hydrosphere, cryosphere and anthroposphere calls for an increased use of Earth observation (EO) data attaining volumes previously considered prohibitive. We review data availability and discuss recent examples where satellite remote sensing is used to infer observable surface quantities directly or indirectly, with particular emphasis on key parameters necessary for weather and climate prediction. Coordinated high-resolution remote-sensing and modelling/assimilation capabilities for the Earth surface are required to support an international application-focused effort

Assessing Global Surface Water Inundation Dynamics Using Combined Satellite Information from SMAP, AMSR2 and Landsat

A method to assess global land surface water (fw) inundation dynamics was developed by exploiting the enhanced fw sensitivity of L-band (1.4 GHz) passive microwave observations from the Soil Moisture Active Passive (SMAP) mission. The L-band fw (fw(sub LBand)) retrievals were derived using SMAP H-polarization brightness temperature (Tb) observations and predefined L-band reference microwave emissivities for water and land endmembers. Potential soil moisture and vegetation contributions to the microwave signal were represented from overlapping higher frequency (Tb) observations from AMSR2. The resulting (fw(sub LBand)) global record has high temporal sampling (1-3 days) and 36-km spatial resolution. The (fw(sub LBand)) annual averages corresponded favourably (R=0.84, p<0.001) with a 250-m resolution static global water map (MOD44W) aggregated at the same spatial scale, while capturing significant inundation variations worldwide. The monthly (fw(sub LBand)) averages also showed seasonal inundation changes consistent with river discharge records within six major US river basins. An uncertainty analysis indicated generally reliable (fw(sub LBand)) performance for major land cover areas and under low to moderate vegetation cover, but with lower accuracy for detecting water bodies covered by dense vegetation. Finer resolution (30-m) (fw(sub LBand)) results were obtained for three sub-regions in North America using an empirical downscaling approach and ancillary global Water Occurrence Dataset (WOD) derived from the historical Landsat record. The resulting 30-m (fw(sub LBand)) retrievals showed favourable spatial accuracy for water (70.71%) and land (98.99%) classifications and seasonal wet and dry periods when compared to independent water maps derived from Landsat-8 imagery. The new (fw(sub LBand)) algorithms and continuing SMAP and AMSR2 operations provide for near real-time, multi-scale monitoring of global surface water inundation dynamics and potential flood risk

Proxy Indicators for Mapping the End of the Vegetation Active Period in Boreal Forests Inferred from Satellite-Observed Soil Freeze and ERA-Interim Reanalysis Air Temperature

Triggered by decreases in photoperiod and temperature, evergreen needle-leaved trees in the boreal region downregulate photosynthetic activity and enter dormancy in autumn. Accompanying changes in canopy structure and chlorophyll content are small and precede the cessation of photosynthetic activity. Low solar elevation and cloud cover during this period pose additional challenges for the use of optical satellite instruments. Alternatively, environmental variables that correlate with photosynthesis, such as soil freeze, can be detected from satellite microwave observations independent of weather and illumination conditions. We tested for the first time the usability of satellite-observed soil freeze from the Soil Moisture and Ocean Salinity (SMOS) instrument as a proxy indicator for the end of vegetation active period (VAPend) at six eddy covariance sites in Finland and Canada. The time when soil freeze commenced over the large SMOS pixel can be employed to estimate VAPend (R-2=0.84, RMSE=7.5days), defined as the time when the photosynthetic capacity of the forest drops below 10% of the growing season maximum. In comparison to satellite-based soil freeze timing, an air temperature-based proxy from ERA-Interim reanalysis data showed better performance (R-2=0.92, RMSE=5.2days). VAPend was mapped in the boreal forest zone in Finland and Canada from both indicators based on linear regression models.Peer reviewe