Engineering topics
Liang, Shunlin
Publications and source records attributed to Liang, Shunlin.
The Global LAnd Surface Satellite (GLASS) evapotranspiration product Version 5.0: Algorithm development and preliminary validation
An accurate estimation of spatially and temporally continuous global terrestrial evapotranspiration (ET) is essential in the assessment of surface energy, water and carbon cycles. The Global LAnd Surface Satellite (GLASS) ET product Version 4.0 (v4.0) based on the Bayesian model averaging (BMA) method was generated to estimate global terrestrial ET. However, certain uncertainty for the GLASS ET product v4.0 limits its application. In this study, we introduced the deep neural networks (DNN) merging framework to improve terrestrial ET estimation for GLASS ET product Version 5.0 (v5.0) generation by integrating five satellite-derived ET products [Moderate Resolution Imaging Spectroradiometer (MODIS) ET product (MOD16), Shuttleworth–Wallace dual-source ET product (SW), Priestley–Taylor-based ET product (PT-JPL), modified satellite-based Priestley–Taylor ET product (MS-PT) and simple hybrid ET product (SIM)]. We compared the performance of DNN method against other merging methods, including GLASS ET algorithm v4.0 (BMA), the gradient boosting regression tree (GBRT) method and the random forest (RF) method, based on 195 global eddy covariance (EC) flux towers covering observations from 2000 through 2015. Validations indicated that the DNN had the highest accuracy among four merging methods across different land cover types, yielding the highest average determination coefficients (R 2 , 0.62), root-mean-squared-error (RMSE, 24.1 W/m 2 ) and Kling–Gupta efficiency (KGE, 0.77) with a of 99% confidence interval. Compared with GLASS ET algorithm v4.0, the DNN improved on the R 2 by approximately 7% (p < 0.01) and the KGE by 10%. Based on the DNN, we then generated 8-day GLASS ET product v5.0 globally with a 1 km spatial resolution from 2001 to 2015 driven by GLASS vegetation and surface net radiation (R n ) datasets and Modern-Era Retrospective Analysis for Research and Applications, Version 2 (MERRA2) datasets. Finally, this global terrestrial ET product provides a valuable dataset for monitoring regional and global water resources and environmental changes.
Recent Third Pole's rapid warming accompanies cryospheric melt and water cycle intensification and interactions between monsoon and environment: Multidisciplinary approach with observation, modeling and analysis
The Third Pole (TP) is experiencing rapid warming and is currently in its warmest period in the past 2000 years. This paper reviews the latest development in multi-disciplinary TP research associated with this warming. The rapid warming facilitates intense and broad glacier melt over most of TP, although some glaciers in the northwest are advancing. By heating the atmosphere and reducing snow/ice albedo, aerosols also contribute to the glaciers melting. Glacier melt is accompanied by lake expansion and intensification of the water cycle over the TP. Precipitation has increased over the eastern and northwestern TP. Meanwhile, the TP is greening and most regions are experiencing advancing phenological trends although over the southwest there is a spring phenological delay mainly in response to the recent decline in spring precipitation. Atmospheric and terrestrial thermal and dynamical processes over the TP affect the Asian monsoon at different scales. Recent evidence indicates substantial roles that Meso-scale Convective Systems play in TP's precipitation as well as an association between soil moisture anomalies in TP and the Indian monsoon. Moreover, an increase in geo-hazard events has been associated with recent environmental changes, some of which have had catastrophic consequences caused by glacial lake outbursts and landslides. Active debris flows are growing in both frequency of occurrences and spatial scale. Meanwhile, new types of disasters, such as the twin ice avalanches in Ali in 2016, are now appearing in the region. Adaptation and mitigation measures should be taken to help societies' preparation for future environmental challenges. Some key issues for future TP studies are also discussed.
Comprehensive Assessment of Global Surface Net Radiation Products and Uncertainty Analysis
Abstract Earth surface net radiation ( R n ) characterizes the surface radiation budget and plays a critical role in ecological, biogeochemical, and hydrological processes. The R n products from remote sensing and reanalysis have not been validated comprehensively. In this study, four R n products (Clouds and the Earth's Radiant Energy System [CERES], ERA‐Interim, Modern‐Era Retrospective analysis for Research and Applications version 2, and Japanese 55‐year Reanalysis) were validated using global ground measurements on monthly (255 sites) and annual (172 sites) timescales. These products have similar accuracies, with average root‐mean‐square error (RMSE) ranges of 5.35 W m −2 (monthly) and 2.30 W m −2 (annually). However, varying accuracies and systemic biases exist across different climatic zones. The annual land R n intercomparison illustrates that large uncertainty exists over polar regions and deserts. A significantly negative annual anomaly in the CERES product for the 2001–2008 period is identified when examining annual R n anomalies over the global land surface. Detailed uncertainty analysis indicates that the global CERES R n anomaly is mainly due to different versions of input data such as aerosol optical thickness and atmospheric profiles (in 2006 and 2008) and cloud properties (in 2002). This work demonstrates that temporal analysis provides powerful quality control for global time series satellite products when the validation using ground measurements fails to capture potential issues.
Estimation of high-resolution terrestrial evapotranspiration from Landsat data using a simple Taylor skill fusion method
Not Available
An Empirical Orthogonal Function-Based Algorithm for Estimating Terrestrial Latent Heat Flux from Eddy Covariance, Meteorological and Satellite Observations
Accurate estimation of latent heat flux (LE) based on remote sensing data is critical in characterizing terrestrial ecosystems and modeling land surface processes. Many LE products were released during the past few decades, but their quality might not meet the requirements in terms of data consistency and estimation accuracy. Merging multiple algorithms could be an effective way to improve the quality of existing LE products. In this paper, we present a data integration method based on modified empirical orthogonal function (EOF) analysis to integrate the Moderate Resolution Imaging Spectroradiometer (MODIS) LE product (MOD16) and the Priestley-Taylor LE algorithm of Jet Propulsion Laboratory (PT-JPL) estimate. Twenty-two eddy covariance (EC) sites with LE observation were chosen to evaluate our algorithm, showing that the proposed EOF fusion method was capable of integrating the two satellite data sets with improved consistency and reduced uncertainties. Further efforts were needed to evaluate and improve the proposed algorithm at larger spatial scales and time periods, and over different land cover types.
GLASS daytime all-wave net radiation product: Algorithm development and preliminary validation
Mapping surface all-wave net radiation (R n ) is critically needed for various applications. Several existing R n products from numerical models and satellite observations have coarse spatial resolutions and their accuracies may not meet the requirements of land applications. In this study, we develop the Global LAnd Surface Satellite (GLASS) daytime R n product at a 5 km spatial resolution. Its algorithm for converting shortwave radiation to all-wave net radiation using the Multivariate Adaptive Regression Splines (MARS) model is determined after comparison with three other algorithms. The validation of the GLASS R n product based on high-quality in situ measurements in the United States shows a coefficient of determination value of 0.879, an average root mean square error value of 31.61 Wm -2 , and an average bias of 17.59 Wm -2 . Furthermore, we also compare our product/algorithm with another satellite product (CERES-SYN) and two reanalysis products (MERRA and JRA55), and find that the accuracy of the much higher spatial resolution GLASS R n product is satisfactory. The GLASS R n product from 2000 to the present is operational and freely available to the public.
Validation and Spatiotemporal Analysis of CERES Surface Net Radiation Product
The Clouds and the Earth’s Radiant Energy System (CERES) generates one of the few global satellite radiation products. The CERES ARM Validation Experiment (CAVE) has been providing long-term in situ observations for the validation of the CERES products. However, the number of these sites is low and their distribution is globally sparse, and particularly the surface net radiation product has not been rigorously validated yet. Therefore, additional validation efforts are highly required to determine the accuracy of the CERES radiation products. In this study, global land surface measurements were comprehensively collected for use in the validation of the CERES net radiation (R n ) product on a daily (340 sites) and a monthly (260 sites) basis, respectively. The validation results demonstrated that the CERES R n product was, overall, highly accurate. The daily validations had a Mean Bias Error (MBE) of 3.43 W·m -2 , Root Mean Square Error (RMSE) of 33.56 W·m -2 , and R 2 of 0.79, and the monthly validations had an MBE of 3.40 W·m -2 , RMSE of 25.57 W·m -2 , and R 2 of 0.84. The accuracy was slightly lower for the high latitudes. Following the validation, the monthly CERES R n product, from March 2000 to July 2014, was used for a further analysis. The global spatiotemporal variation of the R n , which occurred during the measurement period, was analyzed. In addition, two hot spot regions, the southern Great Plains and south-central Africa, were then selected for use in determining the driving factors or attribution of the R n variation. We determined that R n over the southern Great Plains decreased by -0.33 W·m -2 per year, which was mainly driven by changes in surface green vegetation and precipitation. In south-central Africa, R n decreased at a rate of -0.63 W·m -2 per year, the major driving factor of which was surface green vegetation
Validity of Five Satellite-Based Latent Heat Flux Algorithms for Semi-arid Ecosystems
Accurate estimation of latent heat flux (LE) is critical in characterizing semiarid ecosystems. Many LE algorithms have been developed during the past few decades. However, the algorithms have not been directly compared, particularly over global semiarid ecosystems. In this paper, we evaluated the performance of five LE models over semiarid ecosystems such as grassland, shrub, and savanna using the Fluxnet dataset of 68 eddy covariance (EC) sites during the period 2000–2009. We also used a modern-era retrospective analysis for research and applications (MERRA) dataset, the Normalized Difference Vegetation Index (NDVI) and Fractional Photosynthetically Active Radiation (FPAR) from the moderate resolution imaging spectroradiometer (MODIS) products; the leaf area index (LAI) from the global land surface satellite (GLASS) products; and the digital elevation model (DEM) from shuttle radar topography mission (SRTM30) dataset to generate LE at region scale during the period 2003–2006. The models were the moderate resolution imaging spectroradiometer LE (MOD16) algorithm, revised remote sensing based Penman–Monteith LE algorithm (RRS), the Priestley–Taylor LE algorithm of the Jet Propulsion Laboratory (PT-JPL), the modified satellite-based Priestley–Taylor LE algorithm (MS-PT), and the semi-empirical Penman LE algorithm (UMD). Direct comparison with ground measured LE showed the PT-JPL and MS-PT algorithms had relative high performance over semiarid ecosystems with the coefficient of determination (R2) ranging from 0.6 to 0.8 and root mean squared error (RMSE) of approximately 20 W/m 2 . Empirical parameters in the structure algorithms of MOD16 and RRS, and calibrated coefficients of the UMD algorithm may be the cause of the reduced performance of these LE algorithms with R2 ranging from 0.5 to 0.7 and RMSE ranging from 20 to 35 W/m 2 for MOD16, RRS and UMD. Sensitivity analysis showed that radiation and vegetation terms were the dominating variables affecting LE Fluxes in global semiarid ecosystem.
Estimation of high-resolution land surface net shortwave radiation from AVIRIS data: Algorithm development and preliminary results
Not Available
Surface daytime net radiation estimation using artificial neural networks
Net all-wave surface radiation (R n ) is one of the most important fundamental parameters in various applications. However, conventional R n measurements are difficult to collect because of the high cost and ongoing maintenance of recording instruments. Therefore, various empirical R n estimation models have been developed. This study presents the results of two artificial neural network (ANN) models (general regression neural networks (GRNN) and Neuroet) to estimate R n globally from multi-source data, including remotely sensed products, surface measurements, and meteorological reanalysis products. R n estimates provided by the two ANNs were tested against in-situ radiation measurements obtained from 251 global sites between 1991–2010 both in global mode (all data were used to fit the models) and in conditional mode (the data were divided into four subsets and the models were fitted separately). Based on the results obtained from extensive experiments, it has been proved that the two ANNs were superior to linear-based empirical models in both global and conditional modes and that the GRNN performed better and was more stable than Neuroet. The GRNN estimates had a determination coefficient (R 2 ) of 0.92, a root mean square error (RMSE) of 34.27 W·m –2 , and a bias of –0.61 W·m –2 in global mode based on the validation dataset. In conclusion, ANN methods are a potentially powerful tool for global R n estimation.
Use of In Situ and Airborne Multiangle Data to Assess MODIS and Landsat-Based Estimates of Directional Reflectance and Albedo
The quantification of uncertainty in satellite-derived global surface albedo products is a critical aspect in producing complete, physically consistent, and decadal land property data records for studying ecosystem change. A challenge in validating albedo measurements acquired from space is the ability to overcome the spatial scaling errors that can produce disagreements between satellite and field-measured values. Here, we present the results from an accuracy assessment of MODIS and Landsat-TM albedo retrievals, based on collocated comparisons with tower and airborne Cloud Absorption Radiometer (CAR) measurements collected during the 2007 Cloud and Land Surface Interaction Campaign (CLASIC). The initial focus was on evaluating intersensor consistency through comparisons of intrinsic bidirectional reflectance estimates. Local and regional assessments were then performed to obtain estimates of the resulting scaling uncertainties, and to establish the accuracy of albedo reconstructions during extended periods of precipitation. In general, the satellite derived estimates met the accuracy requirements established for the high-quality MODIS operational albedos at 500 m (the greater of 0.02 units or +/-10% of surface measured values). However, results reveal a high degree of variability in the root-mean-square error (RMSE) and bias of MODIS visible (0.3-0.7 μm) and Landsat-TM shortwave (0.3-5.0 μm) albedos; where, in some cases, retrieval uncertainties were found to be in excess of 15%. Results suggest that an overall improvement in MODIS shortwave albedo retrieval accuracy of 7.8%, based on comparisons between MODIS and CAR albedos, resulted from the removal of sub-grid scale mismatch errors when directly scaling-up the tower measurements to the MODIS satellite footprint.
An analytic radiative transfer model for a coupled atmosphere and leaf canopy
A new analytical radiative transfer model of a leaf canopy is developed that approximates multiple-scattering radiance by a four-stream formulation. The canopy model is coupled to a homogeneous atmospheric model as well as a non-Lambertian lower boundary soil surface. The same four-stream formulation is also used for the calculation of multiple scattering in the atmosphere. Comparisons of radiance derived from the four-stream model with those calculated by an iterative numerical solution of the radiative transfer equation show that the analytic model has a very high accuracy, even with a turbid atmosphere and a very dense canopy in which multiple scattering dominates. Because the coupling of canopy and atmospheric models fully accommodates anisotropic surface reflectance and atmospheric scattering and its effect on directional radiance, the model is especially useful for application to directional radiance and measurements obtained by remote sensing. Retrieval of biophysical parameters using this model is under investigation.
Calculation of the angular radiance distribution for a coupled atmosphere and canopy
The radiative transfer equations for a coupled atmosphere and canopy are solved numerically by an improved Gauss-Seidel iteration algorithm. The radiation field is decomposed into three components: unscattered sunlight, single scattering, and multiple scattering radiance for which the corresponding equations and boundary conditions are set up and their analytical or iterational solutions are explicitly derived. The classic Gauss-Seidel algorithm has been widely applied in atmospheric research. This is its first application for calculating the multiple scattering radiance of a coupled atmosphere and canopy. This algorithm enables us to obtain the internal radiation field as well as radiances at boundaries. Any form of bidirectional reflectance distribution function (BRDF) as a boundary condition can be easily incorporated into the iteration procedure. The hotspot effect of the canopy is accommodated by means of the modification of the extinction coefficients of upward single scattering radiation and unscattered sunlight using the formulation of Nilson and Kuusk. To reduce the computation for the case of large optical thickness, an improved iteration formula is derived to speed convergence. The upwelling radiances have been evaluated for different atmospheric conditions, leaf area index (LAI), leaf angle distribution (LAD), leaf size and so on. The formulation presented in this paper is also well suited to analyze the relative magnitude of multiple scattering radiance and single scattering radiance in both the visible and near infrared regions.