Improving Satellite Leaf Area Index Estimation Based On Various Integration Methods

Leaf Area Index (LAI) is an important land surface biophysical variable that is used to characterize vegetation amount and activity. Current satellite LAI products, however, do not satisfy the requirements of the modeling community due to their large uncertainties and frequent missing values. Each LAI product is currently generated from only one satellite sensor data. There is an urgent need for advanced methods to integrate multiple LAI products to improve the product's accuracy and integrality for various applications. To meet this need, this study proposes four methods, including the Optimal Interpolation (OI), Bayesian Maximum Entropy (BME), Multi-Resolution Tree (MRT) and Empirical Orthogonal Function (EOF), to integrate multiple LAI products. Three LAI products have been considered in this study: Moderate Resolution Imaging Spectroradiometer (MODIS), Multi-angle Imaging SpectroRadiometer (MISR) and Carbon cYcle and Change in Land Observational Products from an Ensemble of Satellites (CYCLOPES) LAI. As the basis of data integration, this dissertation first validates and intercompares MODIS and CYCLOPES LAI products and also evaluates their geometric accuracies. The CYCLOPES LAI product has smoother temporal profiles and fewer spatial variations, but tends to produce spurious large errors in winter. The Locally Adjusted Cubic-spline Capping algorithm is revised to smooth multiple years' average and variance. Although OI, BME and MRT based methods have been used in other fields, this is the first research to employ them in integrating multiple LAI products. This dissertation also presents a new integration method based on EOF to solve the problem of large data volume and inconsistent temporal resolution of different datasets. High resolution LAI reference maps generated with ground measurements are used to validate these algorithms. Validation results show that all of these four methods can fill data gaps and reduce the errors of the existing LAI products. The data gaps are filled with information from adjacent pixels and background. These algorithms remove the spurious large temporal and spatial variation of the original LAI products. The combination of multiple satellite products significantly reduces bias. OI and BME can reduce the RMSE from 1.0 (MODIS) to 0.7 and reduce the bias from +0.3 (MODIS) and -0.2 (CYCLOPES) to -0.1. MRT can produce similar results with OI but with significantly improved efficiency. EOF also generates the results with the RMSE of 0.7 but zero bias. Limited ground measurement data hardly prove which methods outperform the others. OI and BME theoretically produce statistically optimal results. BME relaxes OI's linear and Gaussian assumption and explicitly considers data error, but bears a much higher computational burden. MRT has improved efficiency but needs strict assumptions on the scale transfer function. EOF requires simpler model identification, while it is more "empirical" than "statistical". The original contributions of this study mainly include: 1) a new application of several different integration methods to incorporate multiple satellite LAI products to reduce uncertainties and improve integrality, 2) an enhancement of the Locally Adjusted Cubic-spline Capping by revising the end condition, 3) a novel comprehensive comparison of MODIS C5 LAI product with other satellite products, 4) the development of a new LAI normalization scheme by assuming the linear relationship between measurement error and LAI natural variance to account for the inconsistency between products, and finally, 5) the creation of a new data integration method based on EOF

Satellite observations of lake surface state to improve weather forecasting in a lake-rich region

The thermal and dynamic properties of water bodies are important factors affecting the structure of the atmospheric boundary layer which stores and transports energy and mass. The storage and heat transfer of lakes play an essential role in energy and water exchanges with the atmosphere. At high latitudes, the effects of lake ice on climate mostly occur at the local/regional scale, with the degree of influence dependent on the magnitude, timing, location, duration of ice cover and the size of the water body. Ground-based lake temperature and ice observations have been used to investigate the role of lakes in the weather and climate, and the response of lakes to climate. However, in the last two to three decades, it has been observed that the number of ground-based observations, lake ice in particular, has been decreasing dramatically in several countries across the northern hemisphere. In this context, remotely sensed earth observations represent a practical tool in support of the scientific and operational modeling communities, permitting to monitor Lake Surface Water Temperature (LSWT) and ice cover. Data assimilation methods have been used widely to solve the initial value problem in numerical weather prediction (NWP) models. There is a variety of users and applications of space-borne observations in NWP systems; however, not much attention has been paid on the assimilation of remotely-sensed LSWT data in pre-operational NWP environments for improvement of the weather forecast using the optimal interpolation method. This thesis aimed to demonstrate how retrieved remotely-sensed LSWT observations can improve the representation of lake-atmosphere interactions in NWP models. More specifically, LSWT observations were used to improve the representation of lake surface state in the High Resolution Limited Area Model (HIRLAM), a three-dimensional numerical weather prediction (NWP) model. To attain this goal, satellite-derived LSWT observations from the Moderate Resolution Imaging Spectroradiometer (MODIS) and the Along-Track Scanning Radiometer (AATSR) sensors onboard the Terra/Aqua and ENVISAT satellites, respectively, were first evaluated against in-situ measurements collected by the Finnish Environment Institute (SYKE) for a selection of large to medium-size lakes during the open-water season. Results show a good agreement between MODIS and in-situ measurements from 22 Finnish lakes, with a mean bias of -1.13 ˚C determined over five open water seasons (2007-2011). Evaluation of MODIS during an overlapping period (2007-2009) with the AATSR-L2 product currently distributed by the European Space Agency (ESA) shows a mean bias error of -0.93 ˚C for MODIS and a warm mean bias of 1.08 ˚C for AATSR-L2. Two additional LSWT retrieval algorithms were applied to produce more accurate AATSR products. The algorithms use ESA’s AATSR-L1B brightness temperature product to generate new L2 products: one based on Key et al. (1997) and the other on Prata (2002) with a finer resolution water mask than used in the creation of the AATSR-L2 product distributed by ESA. The accuracies of LSWT retrievals are improved with the Key and Prata algorithms with biases of 0.78 ˚C and -0.11 ˚C, respectively, compared to the original AATSR-L2 product (3.18 ˚C). The impact of remotely-sensed LSWT observations in the analysis of lake surface state of HIRLAM forecasting system was then investigated. Data assimilation experiments were performed with the HIRLAM model. Selected thermal remote-sensing LSWT observations provided by MODIS and AATSR sensors were included into the assimilation. The domain of the experiments, which focused on two winters (2010-2011 and 2011-2012), covered northern Europe. Validation of the resulting objective analyses against independent observations demonstrated that the description of the lake surface state can be improved by the introduction of space-borne LSWT observations, compared to the result of pure prognostic parameterizations or assimilation of the available limited number of in-situ lake temperature observations. Further development of the data assimilation methods and solving of several practical issues were found to be necessary in order to fully benefit from the space-borne observations of lake surface state for the improvement of the operational weather forecast. Lastly, the lake-specific autocorrelation function based on LSWT remotely sensed observations was approximated in HIRLAM. A new autocorrelation function of lake pairs was approximated and compared against the original function utilized in current version of HIRLAM to investigate potential improvements demonstrated through HIRLAM sensitivity experiments. The autocorrelation function is calculated based on distance and lake depth differences for each lake pairs. Results show that large lakes are more sensitive to the impact of the autocorrelation. These results also suggest that the high concentrated observations can improve the enhanced result; however, ground-based observations of LSWT are barely available for NWP applications. Overall, results from this thesis clearly demonstrate the benefits of assimilating space-borne LSWT observations into a weather forecasting system such as HIRLAM, and that comprehensive assimilation of LSWT observations can improve NWP results

Review of the CALIMAS Team Contributions to European Space Agency's Soil Moisture and Ocean Salinity Mission Calibration and Validation

Camps, Adriano ... et al.-- 38 pages, 22 figuresThis work summarizes the activities carried out by the SMOS (Soil Moisture and Ocean Salinity) Barcelona Expert Center (SMOS-BEC) team in conjunction with the CIALE/Universidad de Salamanca team, within the framework of the European Space Agency (ESA) CALIMAS project in preparation for the SMOS mission and during its first year of operation. Under these activities several studies were performed, ranging from Level 1 (calibration and image reconstruction) to Level 4 (land pixel disaggregation techniques, by means of data fusion with higher resolution data from optical/infrared sensors). Validation of SMOS salinity products by means of surface drifters developed ad-hoc, and soil moisture products over the REMEDHUS site (Zamora, Spain) are also presented. Results of other preparatory activities carried out to improve the performance of eventual SMOS follow-on missions are presented, including GNSS-R to infer the sea state correction needed for improved ocean salinity retrievals and land surface parameters. Results from CALIMAS show a satisfactory performance of the MIRAS instrument, the accuracy and efficiency of the algorithms implemented in the ground data processors, and explore the limits of spatial resolution of soil moisture products using data fusion, as well as the feasibility of GNSS-R techniques for sea state determination and soil moisture monitoringThis work has been performed under research grants TEC2005-06863-C02-01/TCM, ESP2005-06823-C05, ESP2007-65667-C04, AYA2008-05906-C02-01/ESP and AYA2010-22062-C05 from the Spanish Ministry of Science and Innovation, and a EURYI 2004 award from the European Science FoundationPeer Reviewe

Global Modeling and Assimilation Office Annual Report and Research Highlights 2011-2012

Over the last year, the Global Modeling and Assimilation Office (GMAO) has continued to advance our GEOS-5-based systems, updating products for both weather and climate applications. We contributed hindcasts and forecasts to the National Multi-Model Ensemble (NMME) of seasonal forecasts and the suite of decadal predictions to the Coupled Model Intercomparison Project (CMIP5)

Estimation of Surface Moisture Content and Evapotranspiration Using Weightage Approach.

Soil moisture (MC) and evapotranspiration (ET) are considered as the most significant boundary conditions controlling most of the hydrological cycle’s processes. However, monitoring them continuously over large areas using the high temporal-resolution optical satellites is very demanding. Satellites such as the Advanced Very High Resolution Radiometer (AVHRR) and the Moderate Resolution Imaging Spectroradiometer (MODIS), have a coarse spatial resolution in their images. Thus it not only impedes the acquisition of an accurate MC and ET but also represents multispectral reflections from the holistic surface features. This beside their dependence on vegetation and ground coefficient when assessing MC and ET. The study aims to enhance the spatial accuracy by weighting the MC produced from different surface cover classes within the pixel. MC for each pixel is segmented into three (3) different classes namely urban, vegetation and multi surface cover according to their respective MC weightage. Secondly, to generate an improved actual ETa map by overlaying the segmented MC with a rectified ETo. Images from AVHRR and MODIS satellites were selected in order to generate MC and ET maps. Two powerful MC algorithms were used based on land Surface Temperature (Ts), vegetation Indices (VI) and field measurements of MC; which were conducted at variable depths to examine the depth influence on MC and Ts magnitudes

A global long-term (1981–2000) land surface temperature product for NOAA AVHRR

Land surface temperature (LST) plays an important role in the research of climate change and various land surface processes. Before 2000, global LST products with relatively high temporal and spatial resolutions are scarce, despite a variety of operational satellite LST products. In this study, a global 0.05∘×0.05∘ historical LST product is generated from NOAA advanced very-high-resolution radiometer (AVHRR) data (1981–2000), which includes three data layers: (1) instantaneous LST, a product generated by integrating several split-window algorithms with a random forest (RF-SWA); (2) orbital-drift-corrected (ODC) LST, a drift-corrected version of RF-SWA LST; and (3) monthly averages of ODC LST. For an assumed maximum uncertainty in emissivity and column water vapor content of 0.04 and 1.0 g cm−2, respectively, evaluated against the simulation dataset, the RF-SWA method has a mean bias error (MBE) of less than 0.10 K and a standard deviation (SD) of 1.10 K. To compensate for the influence of orbital drift on LST, the retrieved RF-SWA LST was normalized with an improved ODC method. The RF-SWA LST were validated with in situ LST from Surface Radiation Budget (SURFRAD) sites and water temperatures obtained from the National Data Buoy Center (NDBC). Against the in situ LST, the RF-SWA LST has a MBE of 0.03 K with a range of −1.59–2.71 K, and SD is 1.18 K with a range of 0.84–2.76 K. Since water temperature only changes slowly, the validation of ODC LST was limited to SURFRAD sites, for which the MBE is 0.54 K with a range of −1.05 to 3.01 K and SD is 3.57 K with a range of 2.34 to 3.69 K, indicating good product accuracy. As global historical datasets, the new AVHRR LST products are useful for filling the gaps in long-term LST data. Furthermore, the new LST products can be used as input to related land surface models and environmental applications. Furthermore, in support of the scientific research community, the datasets are freely available at https://doi.org/10.5281/zenodo.3934354 for RF-SWA LST (Ma et al., 2020a), https://doi.org/10.5281/zenodo.3936627 for ODC LST (Ma et al., 2020c), and https://doi.org/10.5281/zenodo.3936641 for monthly averaged LST (Ma et al., 2020b)

Preliminary Evaluation of Influence of Aerosols on the Simulation of Brightness Temperature in the NASA's Goddard Earth Observing System Atmospheric Data Assimilation System

This document reports on preliminary results obtained when studying the impact of aerosols on the calculation of brightness temperature (BT) for satellite infrared (IR) instruments that are currently assimilated in a 3DVAR configuration of Goddard Earth Observing System (GEOS)-atmospheric data assimilation system (ADAS). A set of fifteen aerosol species simulated by the Goddard Chemistry Aerosol Radiation and Transport (GOCART) model is used to evaluate the influence of the aerosol fields on the Community Radiative Transfer Model (CRTM) calculations taking place in the observation operators of the Gridpoint Statistical Interpolation (GSI) analysis system of GEOSADAS. Results indicate that taking aerosols into account in the BT calculation improves the fit to observations over regions with significant amounts of dust. The cooling effect obtained with the aerosol-affected BT leads to a slight warming of the analyzed surface temperature (by about 0:5oK) in the tropical Atlantic ocean (off northwest Africa), whereas the effect on the air temperature aloft is negligible. In addition, this study identifies a few technical issues to be addressed in future work if aerosol-affected BT are to be implemented in reanalysis and operational settings. The computational cost of applying CRTM aerosol absorption and scattering options is too high to justify their use, given the size of the benefits obtained. Furthermore, the differentiation between clouds and aerosols in GSI cloud detection procedures needs satisfactory revision

Observation and integrated Earth-system science: a roadmap for 2016–2025

This report is the response to a request by the Committee on Space Research of the International Council for Science to prepare a roadmap on observation and integrated Earth-system science for the coming ten years. Its focus is on the combined use of observations and modelling to address the functioning, predictability and projected evolution of interacting components of the Earth system on timescales out to a century or so. It discusses how observations support integrated Earth-system science and its applications, and identifies planned enhancements to the contributing observing systems and other requirements for observations and their processing. All types of observation are considered, but emphasis is placed on those made from space. The origins and development of the integrated view of the Earth system are outlined, noting the interactions between the main components that lead to requirements for integrated science and modelling, and for the observations that guide and support them. What constitutes an Earth-system model is discussed. Summaries are given of key cycles within the Earth system. The nature of Earth observation and the arrangements for international coordination essential for effective operation of global observing systems are introduced. Instances are given of present types of observation, what is already on the roadmap for 2016–2025 and some of the issues to be faced. Observations that are organised on a systematic basis and observations that are made for process understanding and model development, or other research or demonstration purposes, are covered. Specific accounts are given for many of the variables of the Earth system. The current status and prospects for Earth-system modelling are summarized. The evolution towards applying Earth-system models for environmental monitoring and prediction as well as for climate simulation and projection is outlined. General aspects of the improvement of models, whether through refining the representations of processes that are already incorporated or through adding new processes or components, are discussed. Some important elements of Earth-system models are considered more fully. Data assimilation is discussed not only because it uses observations and models to generate datasets for monitoring the Earth system and for initiating and evaluating predictions, in particular through reanalysis, but also because of the feedback it provides on the quality of both the observations and the models employed. Inverse methods for surface-flux or model-parameter estimation are also covered. Reviews are given of the way observations and the processed datasets based on them are used for evaluating models, and of the combined use of observations and models for monitoring and interpreting the behaviour of the Earth system and for predicting and projecting its future. A set of concluding discussions covers general developmental needs, requirements for continuity of space-based observing systems, further long-term requirements for observations and other data, technological advances and data challenges, and the importance of enhanced international co-operation