Evaluation of the Surface Water Distribution in North-Central Namibia Based on MODIS and AMSR Series

Semi-arid North-central Namibia has high potential for rice cultivation because large seasonal wetlands (oshana) form during the rainy season. Evaluating the distribution of surface water would reveal the area potentially suitable for rice cultivation. In this study, we detected the distribution of surface water with high spatial and temporal resolution by using two types of complementary satellite data: MODIS (MODerate-resolution Imaging Spectroradiometer) and AMSR-E (Advanced Microwave Scanning Radiometer–Earth Observing System), using AMSR2 after AMSR-E became unavailable. We combined the modified normalized-difference water index (MNDWI) from the MODIS data with the normalized-difference polarization index (NDPI) from the AMSR-E and AMSR2 data to determine the area of surface water. We developed a simple gap-filling method (“database unmixing”) with the two indices, thereby providing daily 500-m-resolution MNDWI maps of north-central Namibia regardless of whether the sky was clear. Moreover, through receiver-operator characteristics (ROC) analysis, we determined the threshold MNDWI (−0.316) for wetlands. Using ROC analysis, MNDWI had moderate performance (the area under the ROC curve was 0.747), and the recognition error for seasonal wetlands and dry land was 21.2%. The threshold MNDWI let us calculate probability of water presence (PWP) maps for the rainy season and the whole year. The PWP maps revealed the total area potentially suitable for rice cultivation: 1255 km2 (1.6% of the study area)

Analysis and quantification of ENSO-linked changes in the tropical Atlantic cloud vertical distribution using 14 years of MODIS observations

A total of 14 years (September 2002 to September 2016) of Aqua Moderate Resolution Imaging Spectroradiometer (MODIS) monthly mean cloud data are used to quantify possible changes in the cloud vertical distribution over the tropical Atlantic. For the analysis multiple linear regression techniques are used. For the investigated time period significant linear changes were found in the domain-averaged cloud-top height (CTH) (−178 m per decade), the high-cloud fraction (HCF) (−0.0006 per decade), and the low-cloud amount (0.001 per decade). The interannual variability of the time series (especially CTH and HCF) is highly influenced by the El Niño–Southern Oscillation (ENSO). Separating the time series into two phases, we quantified the linear change associated with the transition from more La Niña-like conditions to a phase with El Niño conditions (Phase 2) and vice versa (Phase 1). The transition from negative to positive ENSO conditions was related to a decrease in total cloud fraction (TCF) (−0.018 per decade; not significant) due to a reduction in the high-cloud amount (−0.024 per decade; significant). Observed anomalies in the mean CTH were found to be mainly caused by changes in HCF rather than by anomalies in the height of cloud tops themselves. Using the large-scale vertical motion ω at 500 hPa (from ERA-Interim ECMWF reanalysis data), the observed anomalies were linked to ENSO-induced changes in the atmospheric large-scale dynamics. The most significant and largest changes were found in regions with strong large-scale upward movements near the Equator. Despite the fact that with passive imagers such as MODIS it is not possible to vertically resolve clouds, this study shows the great potential for large-scale analysis of possible changes in the cloud vertical distribution due to the changing climate by using vertically resolved cloud cover and linking those changes to large-scale dynamics using other observations or model data

Monitoring of an Indonesian Tropical Wetland by Machine Learning-Based Data Fusion of Passive and Active Microwave Sensors

In this study, a novel data fusion approach was used to monitor the water-body extent in a tropical wetland (Lake Sentarum, Indonesia). Monitoring is required in the region to support the conservation of water resources and biodiversity. The developed approach, random forest database unmixing (RFDBUX), makes use of pixel-based random forest regression to overcome the limitations of the existing lookup-table-based approach (DBUX). The RFDBUX approach with passive microwave data (AMSR2) and active microwave data (PALSAR-2) was used from 2012 to 2017 in order to obtain PALSAR-2-like images with a 100 m spatial resolution and three-day temporal resolution. In addition, a thresholding approach for the obtained PALSAR-2-like backscatter coefficient images provided water body extent maps. The validation revealed that the spatial patterns of the images predicted by RFDBUX are consistent with the original PALSAR-2 backscatter coefficient images (r = 0.94, RMSE = 1.04 in average), and that the temporal pattern of the predicted water body extent can track the wetland dynamics. The PALSAR-2-like images should be a useful basis for further investigation of the hydrological/climatological features of the site, and the proposed approach appears to have the potential for application in other tropical regions worldwide

Harmonization of remote sensing land surface products : correction of clear-sky bias and characterization of directional effects

Tese de doutoramento, Ciências Geofísicas e da Geoinformação (Deteção Remota), Universidade de Lisboa, Faculdade de Ciências, 2018Land surface temperature (LST) is the mean radiative skin temperature of an area of land resulting from the mean energy balance at the surface. LST is an important climatological variable and a diagnostic parameter of land surface conditions, since it is the primary variable determining the upward thermal radiation and one of the main controllers of sensible and latent heat fluxes between the surface and the atmosphere. The reliable and long-term estimation of LST is therefore highly relevant for a wide range of applications, including, amongst others: (i) land surface model validation and monitoring; (ii) data assimilation; (iii) hydrological applications; and (iv) climate monitoring. Remote sensing constitutes the most effective method to observe LST over large areas and on a regular basis. Satellite LST products generally rely on measurements in the thermal infrared (IR) atmospheric window, i.e., within the 8-13 micrometer range. Beside the relatively weak atmospheric attenuation under clear sky conditions, this band includes the peak of the Earth’s spectral radiance, considering surface temperature of the order of 300K (leading to maximum emission at approximately 9.6 micrometer, according to Wien’s Displacement Law). The estimation of LST from remote sensing instruments operating in the IR is being routinely performed for nearly 3 decades. Nevertheless, there is still a long list of open issues, some of them to be addressed in this PhD thesis. First, the viewing position of the different remote sensing platforms may lead to variability of the retrieved surface temperature that depends on the surface heterogeneity of the pixel – dominant land cover, orography. This effect introduces significant discrepancies among LST estimations from different sensors, overlapping in space and time, that are not related to uncertainties in the methodologies or input data used. Furthermore, these directional effects deviate LST products from an ideally defined LST, which should correspond to the ensemble directional radiometric temperature of all surface elements within the FOV. In this thesis, a geometric model is presented that allows the upscaling of in situ measurements to the any viewing configuration. This model allowed generating a synthetic database of directional LST that was used consistently to evaluate different parametric models of directional LST. Ultimately, a methodology is proposed that allows the operational use of such parametric models to correct angular effects on the retrieved LST. Second, the use of infrared data limits the retrieval of LST to clear sky conditions, since clouds “close” the atmospheric window. This effect introduces a clear-sky bias in IR LST datasets that is difficult to quantify since it varies in space and time. In addition, the cloud clearing requirement severely limits the space-time sampling of IR measurements. Passive microwave (MW) measurements are much less affected by clouds than IR observations. LST estimates can in principle be derived from MW measurements, regardless of the cloud conditions. However, retrieving LST from MW and matching those estimations with IR-derived values is challenging and there have been only a few attempts so far. In this thesis, a methodology is presented to retrieve LST from passive MW observations. The MW LST dataset is examined comprehensively against in situ measurements and multiple IR LST products. Finally, the MW LST data is used to assess the spatial-temporal patterns of the clear-sky bias at global scale.Fundação para a Ciência e a Tecnologia, SFRH/BD/9646

River gauging at global scale using optical and passive microwave remote sensing

Recent discharge observations are lacking for most rivers globally. Discharge can be estimated from remotely sensed floodplain and channel inundation area, but there is currently no method that can be automatically extended to many rivers. We examined whether automated monitoring is feasible by statistically relating inundation estimates from moderate to coarse (>0.05°) resolution remote sensing to monthly station discharge records. Inundation extents were derived from optical MODIS data and passive microwave sensors, and compared to monthly discharge records from over 8000 gauging stations and satellite altimetry observations for 442 reaches of large rivers. An automated statistical method selected grid cells to construct “satellite gauging reaches” (SGRs). MODIS SGRs were generally more accurate than passive microwave SGRs, but there were complementary strengths. The rivers widely varied in size, regime, and morphology. As expected performance was low (R  0.6. The best results (R > 0.9) were obtained for large unregulated lowland rivers, particularly in tropical and boreal regions. Relatively poor results were obtained in arid regions, where flow pulses are few and recede rapidly, and in temperate regions, where many rivers are modified and contained. Provided discharge variations produce clear changes in inundated area and gauge records are available for part of the satellite record, SGRs can retrieve monthly river discharge values back to around 1998 and up to present

Mapping gains and losses in woody vegetation across global tropical drylands

Woody vegetation in global tropical drylands is of significant importance for both the interannual variability of the carbon cycle and local livelihoods. Satellite observations over the past decades provide a unique way to assess the vegetation long-term dynamics across biomes worldwide. Yet, the actual changes in the woody vegetation are always hidden by interannual fluctuations of the leaf density, because the most widely used remote sensing data are primarily related to the photosynthetically active vegetation components. Here, we quantify the temporal trends of the nonphotosynthetic woody components (i.e., stems and branches) in global tropical drylands during 2000–2012 using the vegetation optical depth (VOD), retrieved from passive microwave observations. This is achieved by a novel method focusing on the dry season period to minimize the influence of herbaceous vegetation and using MODerate resolution Imaging Spectroradiometer (MODIS) Normalized Difference Vegetation Index (NDVI) data to remove the interannual fluctuations of the woody leaf component. We revealed significant trends (P < 0.05) in the woody component (VODwood) in 35% of the areas characterized by a nonsignificant trend in the leaf component (VODleaf modeled from NDVI), indicating pronounced gradual growth/decline in woody vegetation not captured by traditional assessments. The method is validated using a unique record of ground measurements from the semiarid Sahel and shows a strong agreement between changes in VODwood and changes in ground observed woody cover (r2 = 0.78). Reliability of the obtained woody component trends is also supported by a review of relevant literatures for eight hot spot regions of change. The proposed approach is expected to contribute to an improved assessment of, for example, changes in dryland carbon pools

Multi-annual flood mapping using multi-sensor satellite data in the Iishana Sub-Basin (Namibia/Angola)

In this study a remote sensing approach based on multi-sensor satellite data for the creation of different hydrological products is presented. It includes different methods for deriving general water masks from optical multi-spectral remote sensing data and sensors using synthetic aperture radar (SAR), as well as their joint processing. Result layers are calculated for the maximum flood extent, the relative frequencies of flood observations and the total flood duration for all three study years (10/2007 - 09/2008, 10/2008 - 09/2009, 10/2016 – 09/2017). Since the study combines already validated methods, quality layers provide information about the confidence of the results. Further analyses such as the aggregation of the annual results allow insights into the complex surface water dynamics in the transnational study area of the Iishana Zone (Namibia/Angola)

Evaluating global trends (1988-2010) in harmonized multi-satellite surface soil moisture

[1] Global trends in a new multi-satellite surface soil moisture dataset were analyzed for the period 1988–2010. 27% of the area covered by the dataset showed significant trends (p = 0.05). Of these, 73% were negative and 27% positive. Subtle drying trends were found in the Southern US, central South America, central Eurasia, northern Africa and the Middle East, Mongolia and northeast China, northern Siberia, and Western Australia. The strongest wetting trends were found in southern Africa and the subarctic region. Intra-annual analysis revealed that most trends are not uniform among seasons. The most prominent trend patterns in remotely sensed surface soil moisture were also found in GLDAS-Noah and ERA Interim modeled surface soil moisture and GPCP precipitation, lending confidence to the obtained results. The relationship with trends in GIMMS-NDVI appeared more complex. In areas of mutual disagreement more research is needed to identify potential deficiencies in models and/or remotely sensed products

Large simulated radiative effects of smoke in the south-east Atlantic

A 1200×1200km2 area of the tropical South Atlantic Ocean near Ascension Island is studied with the HadGEM climate model at convection-permitting and global resolutions for a 10-day case study period in August 2016. During the simulation period, a plume of biomass burning smoke from Africa moves into the area and mixes into the clouds. At Ascension Island, this smoke episode was the strongest of the 2016 fire season. The region of interest is simulated at 4km resolution, with no parameterised convection scheme. The simulations are driven by, and compared to, the global model. For the first time, the UK Chemistry and Aerosol model (UKCA) is included in a regional model with prognostic aerosol number concentrations advecting in from the global model at the boundaries of the region. Fire emissions increase the total aerosol burden by a factor of 3.7 and cloud droplet number concentrations by a factor of 3, which is consistent with MODIS observations. In the regional model, the inversion height is reduced by up to 200m when smoke is included. The smoke also affects precipitation, to an extent which depends on the model microphysics. The microphysical and dynamical changes lead to an increase in liquid water path of 60 g m−2 relative to a simulation without smoke aerosol, when averaged over the polluted period. This increase is uncertain, and smaller in the global model. It is mostly due to radiatively driven dynamical changes rather than precipitation suppression by aerosol. Over the 5-day polluted period, the smoke has substantial direct radiative effects of +11.4 W m−2 in the regional model, a semi-direct effect of −30.5 W m−2 and an indirect effect of −10.1 W m−2. Our results show that the radiative effects are sensitive to the structure of the model (global versus regional) and the parameterization of rain autoconversion. Furthermore, we simulate a liquid water path that is biased high compared to satellite observations by 22% on average, and this leads to high estimates of the domain-averaged aerosol direct effect and the effect of the aerosol on cloud albedo. With these caveats, we simulate a large net cooling across the region, of −27.6 W m−2