Downscaling Aerosol Optical Thickness from Satellite Observations: Physics and Machine Learning Approaches

In recent years, the satellite observation of aerosol properties has been greatly improved. As a result, the derivation of Aerosol Optical Thickness (AOT), one of the most popular atmospheric parameters used in air pollution monitoring, over ocean and continents from satellite observations shows comparable quality to ground-based measurements. Satellite AOT products is often applied for monitoring at global scale because of its coarse spatial resolution. However, monitoring at local scale such as over cities requires more detailed AOT information. The increase spatial resolution to suitable level has potential for applications of air pollution monitoring at global-to-local scale, detecting emission sources, deciding pollution management strategies, localizing aerosol estimation, etc. In this thesis, we investigated, proposed, implemented and validated algorithms to derive AOT maps with spatial resolution increased up to 1×1 km2 from MODerate resolution Imaging Spectrometer (MODIS) observations provided by National Aeronautics and Space Administration (NASA), while MODIS standard aerosol products provide maps at 10×10 km2 of spatial resolution. The solutions are considered on two perspectives: dynamical downscaling by improving the algorithm for remote sensing of tropospheric aerosol from MODIS and statistical downscaling using Support Vector Regression

A Geostatistical Data Fusion Technique for Merging Remote Sensing and Ground-Based Observations of Aerosol Optical Thickness

Particles in the atmosphere reflect incoming sunlight, tending to cool the Earth below. Some particles, such as soot, also absorb sunlight, which tens to warm the ambient atmosphere. Aerosol optical depth (AOD) is a measure of the amount of particulate matter in the atmosphere, and is a key input to computer models that simulate and predict Earth's changing climate. The global AOD products from the Multi-angle Imaging SpectroRadiometer (MISR) and the MODerate resolution Imaging Spectroradiometer (MODIS), both of which fly on the NASA Earth Observing System's Terra satellite, provide complementary views of the particles in the atmosphere. Whereas MODIS offers global coverage about four times as frequent as MISR, the multi-angle data makes it possible to separate the surface and atmospheric contributions to the observed top-of-atmosphere radiances, and also to more effectively discriminate particle type. Surface-based AERONET sun photometers retrieve AOD with smaller uncertainties than the satellite instruments, but only at a few fixed locations. So there are clear reasons to combine these data sets in a way that takes advantage of their respective strengths. This paper represents an effort at combining MISR, MODIS and AERONET AOD products over the continental US, using a common spatial statistical technique called kriging. The technique uses the correlation between the satellite data and the "ground-truth" sun photometer observations to assign uncertainty to the satellite data on a region-by-region basis. The larger fraction of the sun photometer variance that is duplicated by the satellite data, the higher the confidence assigned to the satellite data in that region. In the Western and Central US, MISR AOD correlation with AERONET are significantly higher than those with MODIS, likely due to bright surfaces in these regions, which pose greater challenges for the single-view MODIS retrievals. In the east, MODIS correlations are higher, due to more frequent sampling of the varying AOD. These results demonstrate how the MISR and MODIS aerosol products are complementary. The underlying technique also provides one method for combining these products in such a way that takes advantage of the strengths of each, in the places and times when they are maximal, and in addition, yields an estimate of the associated uncertainties in space and time

Remote sensing, numerical modelling and ground truthing for analysis of lake water quality and temperature

Freshwater accounts for just 2.5% of the earth’s water resources, and its quality and availability are becoming an issue of global concern in the 21st century. Growing human population, over-exploitation of water sources and pressures of global warming mean that both water quantity and quality are affected. In order to effectively manage water quality there is a need for increased monitoring and predictive modelling of freshwater resources. To address these concerns in New Zealand inland waters, an approach which integrates biological and physical sciences is needed. Remote sensing has the potential to allow this integration and vastly increase the temporal and spatial resolution of current monitoring techniques, which typically involve collecting grab-samples. In a complementary way, lake modelling has the potential to enable more effective management of water resources by testing the effectiveness of a range of possible management scenarios prior to implementation. Together, the combination of remote sensing and modelling data allows for improved model initialisation, calibration and validation, which ultimately aid in understanding of complex lake ecosystem processes. This study investigated the use of remote sensing using empirical and semi-analytical algorithms for the retrieval of chlorophyll a (chl a), tripton, suspended minerals (SM), total suspended sediment (SS) and water surface temperature. It demonstrated the use of spatially resolved statistical techniques for comparing satellite estimated and 3-D simulated water quality and temperature. An automated procedure was developed for retrieval of chl a from Landsat Enhanced Thematic Mapper (ETM+) imagery, using 106 satellite images captured from 1999 to 2011. Radiative transfer-based atmospheric correction was applied to images using the Second Simulation of the Satellite in the Solar Spectrum model (6sv). For the estimation of chl a over a time series of images, the use of symbolic regression resulted in a significant improvement in the precision of chl a hindcasts compared with traditional regression equations. Results from this investigation suggest that remote sensing provides a valuable tool to assess temporal and spatial distributions of chl a. Bio-optical models were applied to quantify the physical processes responsible for the relationship between chl a concentrations and subsurface irradiance reflectance used in regression algorithms, allowing the identification of possible sources of error in chl a estimation. While the symbolic regression model was more accurate than traditional empirical models, it was still susceptible to errors in optically complex waters such as Lake Rotorua, due to the effect of variations of SS and CDOM on reflectance. Atmospheric correction of Landsat 7 ETM+ thermal data was carried out for the purpose of retrieval of lake water surface temperature in Rotorua lakes, and Lake Taupo, North Island, New Zealand. Atmospheric correction was repeated using four sources of atmospheric profile data as input to a radiative transfer model, MODerate resolution atmospheric TRANsmission (MODTRAN) v.3.7. The retrieved water temperatures from 14 images between 2007 and 2009 were validated using a high-frequency temperature sensor deployed from a mid-lake monitoring buoy at the water surface of Lake Rotorua. The most accurate temperature estimation for Lake Rotorua was with radiosonde data as an input into MODTRAN, followed by Moderate Resolution Imaging Spectroradiometer (MODIS) Level 2, Atmospheric Infrared Sounder (AIRS) Level 3, and NASA data. Retrieved surface water temperature was used for assessing spatial heterogeneity of surface water temperature simulated with a three-dimensional (3-D) hydrodynamic model (ELCOM) of Lake Rotoehu, located approximately 20 km east of Lake Rotorua. This comparison demonstrated that simulations reproduced the dominant horizontal variations in surface water temperature in the lake. The transport and mixing of a geothermal inflow and basin-scale circulation patterns were inferred from thermal distributions from satellite and model estimations of surface water temperature and a spatially resolved statistical evaluation was used to validate simulations. This study has demonstrated the potential of accurate satellite-based thermal monitoring to validate water surface temperature simulated by 3-D hydrodynamic models. Semi-analytical and empirical algorithms were derived to determine spatial and temporal variations in SS in Lake Ellesmere, South Island, New Zealand, using MODIS band 1. The semi-analytical model and empirical model had a similar level of precision in SS estimation, however, the semi-analytical model has the advantage of being applicable to different satellite sensors, spatial locations, and SS concentration ranges. The estimations of SS concentration (and estimated SM concentration) from the semi-analytical model were used for a spatially resolved validation of simulations of SM derived from ELCOM-CAEDYM. Visual comparisons were compared with spatially-resolved statistical techniques. The spatial statistics derived from the Map Comparison Kit allowed a non-subjective and quantitative method to rank simulation performance on different dates. The visual and statistical comparison between satellite estimated and model simulated SM showed that the model did not perform well in reproducing both basin-scale and fine-scale spatial variation in SM derived from MODIS satellite imagery. Application of the semi-analytical model to estimate SS over the lifetime of the MODIS sensor will greatly extend its spatial and temporal coverage for historical monitoring purposes, and provide a tool to validate SM simulated by 1-D and 3-D models on a daily basis. A bio-optical model was developed to derive chl a, SS concentrations, and coloured dissolved organic matter /detritus absorption at 443 nm, from MODIS Aqua subsurface remote sensing reflectance of Lake Taupo, a large, deep, oligotrophic lake in North Island, New Zealand. The model was optimised using in situ inherent optical properties (IOPs) from the literature. Images were atmospherically corrected using the radiative transfer model 6sv. Application of the bio-optical model using a single chl a-specific absorption spectrum (a*ϕ(λ)) resulted in low correlation between estimated and observed values. Therefore, two different absorption curves were used, based on the seasonal dominance of phytoplankton phyla with differing absorption properties. The application of this model resulted in reasonable agreement between modelled and in situ chl a concentrations. Highest concentrations were observed during winter when Bacillariophytes (diatoms) dominated the phytoplankton assemblage. On 4 and 5 March 2004 an unusually large turbidity current was observed originating from the Tongariro River inflow in the south-east of the lake. In order to resolve fine details of the plume, empirical relationships were developed between MODIS band 1 reflectance (250 m resolution) and SS estimated from MODIS bio-optical features (1 km resolution) were used estimate SS at 250 m resolution. Complex lake circulation patterns were observed including a large clockwise gyre. With the development of this bio-optical model MODIS can potentially be used to remotely sense water quality in near real time, and the relationship developed for B1 SS allows for resolution of fine-scale features such turbidity currents

Retrieval of liquid water cloud properties from POLDER-3 measurements using a neural network ensemble approach

This paper describes a neural network algorithm for the estimation of liquid water cloud optical properties from the Polarization and Directionality of Earth's Reflectances-3 (POLDER-3) instrument aboard the Polarization &amp; Anisotropy of Reflectances for Atmospheric Sciences coupled with Observations from a Lidar (PARASOL) satellite. The algorithm has been trained on synthetic multi-angle, multi-wavelength measurements of reflectance and polarization and has been applied to the processing of 1 year of POLDER-3 data. Comparisons of the retrieved cloud properties with Moderate Resolution Imaging Spectroradiometer (MODIS) products show that the neural network algorithm has a low bias of around 2 in cloud optical thickness (COT) and between 1 and 2&thinsp;µm in the cloud effective radius. Comparisons with existing POLDER-3 datasets suggest that the proposed scheme may have enhanced capabilities for cloud effective radius retrieval, at least over land. An additional feature of the presented algorithm is that it provides COT and effective radius retrievals at the native POLDER-3 Level 1B pixel level.</p

Extending Landsat 8: Retrieval of an Orange contra-Band for Inland Water Quality Applications

The Operational Land Imager (OLI) onboard Landsat 8 has found successful application in inland and coastal water remote sensing. Its radiometric specification and high spatial resolution allows quantification of water-leaving radiance while resolving small water bodies. However, its limited multispectral band set restricts the range of water quality parameters that can be retrieved. Identification of cyanobacteria biomass has been demonstrated for sensors with a band centered near 620 nm, the absorption peak of the diagnostic pigment phycocyanin. While OLI lacks such a band in the orange region, superposition of the available multispectral and panchromatic bands suggests that it can be calculated by a scaled difference. A set of 428 in situ spectra acquired in diverse lakes in Belgium and The Netherlands was used to develop and test an orange contra-band retrieval algorithm, achieving a mean absolute percentage error of 5.39 % and a bias of −0.88 % in the presence of sensor noise. Atmospheric compensation error propagated to the orange contra-band was observed to maintain about the same magnitude (13 % higher) observed for the red band and thus results in minimal additional effects for possible base line subtraction or band ratio algorithms for phycocyanin estimation. Generality of the algorithm for different reflectance shapes was tested against a set of published average coastal and inland Optical Water Types, showing robust retrieval for all but relatively clear water types (Secchi disk depth > 6 m and chlorophyll a < 1.6 mg m−3). The algorithm was further validated with 79 matchups against the Ocean and Land Colour Imager (OLCI) orange band for 10 globally distributed lakes. The retrieved band is shown to convey information independent from the adjacent bands under variable phycocyanin concentrations. An example application using Landsat 8 imagery is provided for a knowncyanobacterialbloominLakeErie,US.ThemethodisdistributedintheACOLITEatmospheric correction code. The contra-band approach is generic and can be applied to other sensors with overlapping bands. Recommendations are also provided for development of future sensors with broad spectral bands with the objective to maximize the accuracy of possible spectral enhancement

Long-Term Monitoring of Global Climate Forcings and Feedbacks

A workshop on Long-Term Monitoring of Global Climate Forcings and Feedbacks was held February 3-4, 1992, at NASA's Goddard Institute for Space Studies to discuss the measurements required to interpret long-term global temperature changes, to critique the proposed contributions of a series of small satellites (Climsat), and to identify needed complementary monitoring. The workshop concluded that long-term (several decades) of continuous monitoring of the major climate forcings and feedbacks is essential for understanding long-term climate change

Feasibility Study for an Aquatic Ecosystem Earth Observing System Version 1.2.

International audienceMany Earth observing sensors have been designed, built and launched with primary objectives of either terrestrial or ocean remote sensing applications. Often the data from these sensors are also used for freshwater, estuarine and coastal water quality observations, bathymetry and benthic mapping. However, such land and ocean specific sensors are not designed for these complex aquatic environments and consequently are not likely to perform as well as a dedicated sensor would. As a CEOS action, CSIRO and DLR have taken the lead on a feasibility assessment to determine the benefits and technological difficulties of designing an Earth observing satellite mission focused on the biogeochemistry of inland, estuarine, deltaic and near coastal waters as well as mapping macrophytes, macro-algae, sea grasses and coral reefs. These environments need higher spatial resolution than current and planned ocean colour sensors offer and need higher spectral resolution than current and planned land Earth observing sensors offer (with the exception of several R&D type imaging spectrometry satellite missions). The results indicate that a dedicated sensor of (non-oceanic) aquatic ecosystems could be a multispectral sensor with ~26 bands in the 380-780 nm wavelength range for retrieving the aquatic ecosystem variables as well as another 15 spectral bands between 360-380 nm and 780-1400 nm for removing atmospheric and air-water interface effects. These requirements are very close to defining an imaging spectrometer with spectral bands between 360 and 1000 nm (suitable for Si based detectors), possibly augmented by a SWIR imaging spectrometer. In that case the spectral bands would ideally have 5 nm spacing and Full Width Half Maximum (FWHM), although it may be necessary to go to 8 nm wide spectral bands (between 380 to 780nm where the fine spectral features occur -mainly due to photosynthetic or accessory pigments) to obtain enough signal to noise. The spatial resolution of such a global mapping mission would be between ~17 and ~33 m enabling imaging of the vast majority of water bodies (lakes, reservoirs, lagoons, estuaries etc.) larger than 0.2 ha and ~25% of river reaches globally (at ~17 m resolution) whilst maintaining sufficient radiometric resolution

Retrieval of Liquid Water Cloud Properties from POLDER-3 Measurements Using a Neural Network Ensemble Approach

This paper describes a neural network algorithm for the estimation of liquid water cloud optical properties from the Polarization and Directionality of Earth's Reflectances-3 (POLDER-3) instrument aboard the Polarization & Anisotropy of Reflectances for Atmospheric Sciences coupled with Observations from a Lidar (PARASOL) satellite. The algorithm has been trained on synthetic multi-angle, multi-wavelength measurements of reflectance and polarization and has been applied to the processing of 1 year of POLDER-3 data. Comparisons of the retrieved cloud properties with Moderate Resolution Imaging Spectroradiometer (MODIS) products show that the neural network algorithm has a low bias of around 2 in cloud optical thickness (COT) and between 1 and 2m in the cloud effective radius. Comparisons with existing POLDER-3 datasets suggest that the proposed scheme may have enhanced capabilities for cloud effective radius retrieval, at least over land. An additional feature of the presented algorithm is that it provides COT and effective radius retrievals at the native POLDER-3 Level 1B pixel level

Data assimilation in atmospheric chemistry models: current status and future prospects for coupled chemistry meteorology models

Abstract. Data assimilation is used in atmospheric chemistry models to improve air quality forecasts, construct re-analyses of three-dimensional chemical (including aerosol) concentrations and perform inverse modeling of input variables or model parameters (e.g., emissions). Coupled chemistry meteorology models (CCMM) are atmospheric chemistry models that simulate meteorological processes and chemical transformations jointly. They offer the possibility to assimilate both meteorological and chemical data; however, because CCMM are fairly recent, data assimilation in CCMM has been limited to date. We review here the current status of data assimilation in atmospheric chemistry models with a particular focus on future prospects for data assimilation in CCMM. We first review the methods available for data assimilation in atmospheric models, including variational methods, ensemble Kalman filters, and hybrid methods. Next, we review past applications that have included chemical data assimilation in chemical transport models (CTM) and in CCMM. Observational data sets available for chemical data assimilation are described, including surface data, surface-based remote sensing, airborne data, and satellite data. Several case studies of chemical data assimilation in CCMM are presented to highlight the benefits obtained by assimilating chemical data in CCMM. A case study of data assimilation to constrain emissions is also presented. There are few examples to date of joint meteorological and chemical data assimilation in CCMM and potential difficulties associated with data assimilation in CCMM are discussed. As the number of variables being assimilated increases, it is essential to characterize correctly the errors; in particular, the specification of error cross-correlations may be problematic. In some cases, offline diagnostics are necessary to ensure that data assimilation can truly improve model performance. However, the main challenge is likely to be the paucity of chemical data available for assimilation in CCMM