Skip to main content

Dr. Lusheng Liang

Lusheng LiangDr. Lusheng Liang leads the CERES Angular Distribution Model working group and is responsible for developing the empirical angular distribution models, producing and validating the instantaneous CERES fluxes.

The CERES Angular Distribution Model (ADM) working group is responsible to quantify the relationship between radiance and flux over different scene types under different sun-Earth-satellite geometry. The CERES empirical ADMs are developed using CERES Rotating Azimuth Plane (RAP) scan, in which the CERES instruments rotate in azimuth as they scan in elevation, therefore maximize the angular coverage. Scene types are defined based upon many variables, such as surface type, snow and ice fraction, cloud fraction, cloud optical depth, cloud phase, wind speed, etc. These ADMs are used to convert the CERES measured radiances to fluxes.

Contact Information

NASA Langley Research Center
Mail Stop 420, Hampton, VA 23681-2199

Phone: 757-864-1188

Fax: 757-864-7996

Email: lusheng.liang@nasa.gov

Education

Research Interest

Publications

2024

Liang, Lusheng; Su, Wenying; Sejas, Sergio; Eitzen, Zachary; Loeb, Norman G.Liang, L., W. Su, S. Sejas, Z. Eitzen, N. G. Loeb, 2024: Next-generation radiance unfiltering process for the Clouds and the Earth's Radiant Energy System instrument. Atmospheric Measurement Techniques, 17(7), 2147-2163. doi: 10.5194/amt-17-2147-2024. The filtered radiances measured by the Clouds and the Earth's Radiant Energy System (CERES) instruments are converted to shortwave (SW), longwave (LW), and window unfiltered radiances based on regressions developed from theoretical radiative transfer simulations to relate filtered and unfiltered radiances. This paper describes an update to the existing Edition 4 CERES unfiltering algorithm (Loeb et al., 2001), incorporating the most recent developments in radiative transfer modeling, ancillary input datasets, and increased computational and storage capabilities during the past 20 years. Simulations are performed with the updated Moderate Resolution Atmospheric Transmission (MODTRAN) 5.4 version. Over land and snow, the surface bidirectional reflectance distribution function (BRDF) is characterized by a kernel-based representation in the simulations, instead of the Lambertian surface used in the Edition 4 unfiltering process. Radiance unfiltering is explicitly separated into four seasonally dependent land surface groups based on the spectral radiation similarities of different surface types (defined by the International Geosphere-Biosphere Programme); over snow, it is separated into fresh snow, permanent snow, and sea ice. This differs from the Edition 4 unfiltering process where only one set of regressions was used for land and snow, respectively. The instantaneous unfiltering errors are estimated with independent test cases generated from radiative transfer simulations in which the “true” unfiltered radiances from radiative transfer simulations are compared with the unfiltered radiances calculated from the regressions. Overall, the relative errors are mostly within ±0.5 % for SW, within ±0.2 % for daytime LW, and within ±0.1 % for nighttime LW for both CERES Terra Flight Model 1 (FM1) and Aqua FM3 instruments. The unfiltered radiances are converted to fluxes and compared to CERES Edition 4 fluxes. The global mean instantaneous fluxes for Aqua FM3 are reduced by 0.34 to 0.45 W m−2 for SW and increased by 0.25 to 0.46 W m−2 for daytime LW; for Terra FM1, they are reduced by 0.24 to 0.34 W m−2 for SW and increased by 0.08 to 0.28 W m−2 for daytime LW. Nighttime LW flux differences are negligible for both instruments.

2023

Liang, Lusheng; Su, Wenying; Sejas, Sergio; Eitzen, Zachary A.; Loeb, Norman G.Liang, L., W. Su, S. Sejas, Z. A. Eitzen, N. G. Loeb, 2023: Next-generation radiance unfiltering process for the Clouds and Earth’s Radiant Energy System instrument. EGUsphere, 1-25. doi: 10.5194/egusphere-2023-1670. Abstract. The filtered radiances measured by the Clouds and the Earth’s Radiant Energy System (CERES) instruments are converted to shortwave (SW), longwave (LW), and window unfiltered radiances based on regressions developed from theoretical radiative transfer simulations to relate filtered and unfiltered radiances. This paper describes an update to the existing Edition 4 CERES unfiltering algorithm (Loeb et al., 2001), incorporating the most recent developments in radiative transfer modeling, ancillary input datasets, and increased computational and storage capabilities during the past 20 years. Simulations are performed with MODTRAN 5.4. Over land and snow, the surface Bidirectional Reflectance Distribution Function (BRDF) is characterized by a kernel-based representation in the simulations, instead of the Lambertian surface used in the Edition 4 unfiltering process. Radiance unfiltering is explicitly separated into 4 seasonally dependent land surface groups based on the spectral radiation similarities of different surface types (defined by International Geosphere-Biosphere Programme); over snow, it is separated into fresh snow, permanent snow, and sea ice. It contrasts to the Edition 4 unfiltering process that one set of regressions for land and snow, respectively. The instantaneous unfiltering errors are estimated with independent test cases generated from radiative transfer simulations in which the ‘true’ unfiltered radiances from radiative transfer simulations are compared with the unfiltered radiances calculated from the regressions. Overall, the relative errors are mostly within ±0.5 % for SW, within ±0.2 % for daytime LW, and within ±0.1 % for nighttime LW for both CERES Terra Flight Model 1 (FM1) and Aqua FM3 instruments. The unfiltered radiances are converted to fluxes and compared to CERES Edition 4 fluxes. The global mean instantaneous fluxes for Aqua FM3 are reduced by less than 0.42 Wm-2 for SW and increased by less than 0.47 Wm-2 for daytime LW; for Terra FM1, they are reduced by less than 0.31 Wm-2 for SW and increased by less than 0.29 Wm-2 for daytime LW, though regional differences can be as large as 2.0 Wm-2. Nighttime LW flux differences are nearly negligible for both instruments.

2021

Su, Wenying; Liang, Lusheng; Myhre, Gunnar; Thorsen, Tyler J.; Loeb, Norman G.; Schuster, Gregory L.; Ginoux, Paul; Paulot, Fabien; Neubauer, David; Checa-Garcia, Ramiro; Matsui, Hitoshi; Tsigaridis, Kostas; Skeie, Ragnhild B.; Takemura, Toshihiko; Bauer, Susanne E.; Schulz, MichaelSu, W., L. Liang, G. Myhre, T. J. Thorsen, N. G. Loeb, G. L. Schuster, P. Ginoux, F. Paulot, D. Neubauer, R. Checa-Garcia, H. Matsui, K. Tsigaridis, R. B. Skeie, T. Takemura, S. E. Bauer, M. Schulz, 2021: Understanding Top-of-Atmosphere Flux Bias in the AeroCom Phase III Models: A Clear-Sky Perspective. Journal of Advances in Modeling Earth Systems, 13(9), e2021MS002584. doi: 10.1029/2021MS002584. Biases in aerosol optical depths (AOD) and land surface albedos in the AeroCom models are manifested in the top-of-atmosphere (TOA) clear-sky reflected shortwave (SW) fluxes. Biases in the SW fluxes from AeroCom models are quantitatively related to biases in AOD and land surface albedo by using their radiative kernels. Over ocean, AOD contributes about 25% to the S–N mean SW flux bias for the multi-model mean (MMM) result. Over land, AOD and land surface albedo contribute about 40% and 30%, respectively, to the S–N mean SW flux bias for the MMM result. Furthermore, the spatial patterns of the SW flux biases derived from the radiative kernels are very similar to those between models and CERES observation, with the correlation coefficient of 0.6 over ocean and 0.76 over land for MMM using data of 2010. Satellite data used in this evaluation are derived independently from each other, consistencies in their bias patterns when compared with model simulations suggest that these patterns are robust. This highlights the importance of evaluating related variables in a synergistic manner to provide an unambiguous assessment of the models, as results from single parameter assessments are often confounded by measurement uncertainty. Model biases in land surface albedos can and must be corrected to accurately calculate TOA flux. We also compare the AOD trend from three models with the observation-based counterpart. These models reproduce all notable trends in AOD except the decreasing trend over eastern China and the adjacent oceanic regions due to limitations in the emission data set. aerosols; radiative flux; surface albedo

2020

Su, Wenying; Liang, Lusheng; Wang, Hailan; Eitzen, Zachary A.Su, W., L. Liang, H. Wang, Z. A. Eitzen, 2020: Uncertainties in CERES Top-of-Atmosphere Fluxes Caused by Changes in Accompanying Imager. Remote Sensing, 12(12), 2040. doi: 10.3390/rs12122040. The Clouds and the Earth’s Radiant Energy System (CERES) project provides observations of Earth’s radiation budget using measurements from CERES instruments on board the Terra, Aqua, Suomi National Polar-orbiting Partnership (S-NPP), and NOAA-20 satellites. The CERES top-of-atmosphere (TOA) fluxes are produced by converting radiance measurements using empirical angular distribution models, which are functions of cloud properties that are retrieved from imagers flying with the CERES instruments. As the objective is to create a long-term climate data record, not only calibration consistency of the six CERES instruments needs to be maintained for the entire time period, it is also important to maintain the consistency of other input data sets used to produce this climate data record. In this paper, we address aspects that could potentially affect the CERES TOA flux data quality. Discontinuities in imager calibration can affect cloud retrieval which can lead to erroneous flux trends. When imposing an artificial 0.6 per decade decreasing trend to cloud optical depth, which is similar to the trend difference between CERES Edition 2 and Edition 4 cloud retrievals, the decadal SW flux trend changed from − 0.3 5 ± 0.18 Wm − 2 to 0.61 ± 0.18 Wm − 2 . This indicates that a 13% change in cloud optical depth results in about 1% change in the SW flux. Furthermore, different CERES instruments provide valid fluxes at different viewing zenith angle ranges, and including fluxes derived at the most oblique angels unique to S-NPP (>66 ∘ ) can lead to differences of 0.8 Wm − 2 and 0.3 Wm − 2 in global monthly mean instantaneous SW flux and LW flux. To ensure continuity, the viewing zenith angle ranges common to all CERES instruments (<66 ∘ ) are used to produce the long-term Earth’s radiation budget climate data record. The consistency of cloud properties retrieved from different imagers also needs to be maintained to ensure the TOA flux consistency. cloud properties; angular distribution model; climate data record; Earth’s radiation budget
Su, Wenying; Minnis, Patrick; Liang, Lusheng; Duda, David P.; Khlopenkov, Konstantin; Thieman, Mandana M.; Yu, Yinan; Smith, Allan; Lorentz, Steven; Feldman, Daniel; Valero, Francisco P. J.Su, W., P. Minnis, L. Liang, D. P. Duda, K. Khlopenkov, M. M. Thieman, Y. Yu, A. Smith, S. Lorentz, D. Feldman, F. P. J. Valero, 2020: Determining the daytime Earth radiative flux from National Institute of Standards and Technology Advanced Radiometer (NISTAR) measurements. Atmospheric Measurement Techniques, 13(2), 429-443. doi: https://doi.org/10.5194/amt-13-429-2020. Abstract. The National Institute of Standards and Technology Advanced Radiometer (NISTAR) onboard the Deep Space Climate Observatory (DSCOVR) provides continuous full-disk global broadband irradiance measurements over most of the sunlit side of the Earth. The three active cavity radiometers measure the total radiant energy from the sunlit side of the Earth in shortwave (SW; 0.2–4 µm), total (0.4–100 µm), and near-infrared (NIR; 0.7–4 µm) channels. The Level 1 NISTAR dataset provides the filtered radiances (the ratio between irradiance and solid angle). To determine the daytime top-of-atmosphere (TOA) shortwave and longwave radiative fluxes, the NISTAR-measured shortwave radiances must be unfiltered first. An unfiltering algorithm was developed for the NISTAR SW and NIR channels using a spectral radiance database calculated for typical Earth scenes. The resulting unfiltered NISTAR radiances are then converted to full-disk daytime SW and LW flux by accounting for the anisotropic characteristics of the Earth-reflected and emitted radiances. The anisotropy factors are determined using scene identifications determined from multiple low-Earth orbit and geostationary satellites as well as the angular distribution models (ADMs) developed using data collected by the Clouds and the Earth's Radiant Energy System (CERES). Global annual daytime mean SW fluxes from NISTAR are about 6 % greater than those from CERES, and both show strong diurnal variations with daily maximum–minimum differences as great as 20 Wm−2 depending on the conditions of the sunlit portion of the Earth. They are also highly correlated, having correlation coefficients of 0.89, indicating that they both capture the diurnal variation. Global annual daytime mean LW fluxes from NISTAR are 3 % greater than those from CERES, but the correlation between them is only about 0.38.

2018

Loeb, Norman G.; Doelling, David R.; Wang, Hailan; Su, Wenying; Nguyen, Cathy; Corbett, Joseph G.; Liang, Lusheng; Mitrescu, Cristian; Rose, Fred G.; Kato, SeijiLoeb, N. G., D. R. Doelling, H. Wang, W. Su, C. Nguyen, J. G. Corbett, L. Liang, C. Mitrescu, F. G. Rose, S. Kato, 2018: Clouds and the Earth’s Radiant Energy System (CERES) Energy Balanced and Filled (EBAF) Top-of-Atmosphere (TOA) Edition 4.0 Data Product. J. Climate, 31(2), 895–918. doi: 10.1175/JCLI-D-17-0208.1. The Clouds and the Earth’s Radiant Energy System (CERES) Energy Balanced and Filled (EBAF) top-of-atmosphere (TOA) Ed4.0 data product is described. EBAF Ed4.0 is an update to EBAF Ed2.8, incorporating all of the Ed4.0 suite of CERES data product algorithm improvements and consistent input datasets throughout the record. A one-time adjustment to SW and LW TOA fluxes is made to ensure that global mean net TOA flux for July 2005-June 2015 is consistent with the in-situ value of 0.71 W m–2. While global mean all-sky TOA flux differences between Ed4.0 and Ed2.8 are within 0.5 Wm-2, appreciable SW regional differences occur over marine stratocumulus and snow/sea-ice regions. Marked regional differences in SW clear-sky TOA flux occur in polar regions and dust areas over ocean. Clear-sky LW TOA fluxes in EBAF Ed4.0 exceed Ed2.8 in regions of persistent high cloud cover. Owing to substantial differences in global mean clear-sky TOA fluxes, the net cloud radiative effect in EBAF Ed4.0 is -18 Wm-2 compared to -21 Wm-2 in EBAF Ed2.8. We estimate the overall uncertainty in 1°x1° latitude-longitude regional monthly all-sky TOA flux to be 3 Wm-2 (1σ) for the Terra-only period and 2.5 Wm-2 for the Terra-Aqua period both for SW and LW. The SW clear-sky regional monthly uncertainty is estimated to be 6 Wm-2 for the Terra-only period and 5 Wm-2 for the Terra-Aqua period. The LW clear-sky regional monthly uncertainty is 5 Wm-2 for Terra-only and 4.5 Wm-2 for Terra-Aqua.
Su, Wenying; Liang, Lusheng; Doelling, David R.; Minnis, Patrick; Duda, David P.; Khlopenkov, Konstantin V.; Thieman, Mandana M.; Loeb, Norman G.; Kato, Seiji; Valero, Francisco P. J.; Wang, Hailan; Rose, Fred G.Su, W., L. Liang, D. R. Doelling, P. Minnis, D. P. Duda, K. V. Khlopenkov, M. M. Thieman, N. G. Loeb, S. Kato, F. P. J. Valero, H. Wang, F. G. Rose, 2018: Determining the Shortwave Radiative Flux from Earth Polychromatic Imaging Camera. Journal of Geophysical Research: Atmospheres, 123(20), 11,479-11,491. doi: 10.1029/2018JD029390. The Earth Polychromatic Imaging Camera (EPIC) onboard Deep Space Climate Observatory (DSCOVR) provides 10 narrowband spectral images of the sunlit side of the Earth. The blue (443 nm), green (551 nm), and red (680 nm) channels are used to derive EPIC broadband radiances based upon narrowband-to-broadband regressions developed using collocated MODIS equivalent channels and CERES broadband measurements. The pixel-level EPIC broadband radiances are averaged to provide global daytime means at all applicable EPIC times. They are converted to global daytime mean shortwave (SW) fluxes by accounting for the anisotropy characteristics using a cloud property composite based on lower Earth orbiting satellite imager retrievals and the CERES angular distribution models (ADMs). Global daytime mean SW fluxes show strong diurnal variations with daily maximum-minimum differences as great as 20 Wm−2 depending on the conditions of the sunlit portion of the Earth. The EPIC SW fluxes are compared against the CERES SYN1deg hourly SW fluxes. The global monthly mean differences (EPIC-SYN) between them range from 0.1 Wm−2 in July to -4.1 Wm−2 in January, and the RMS errors range from 3.2 Wm−2 to 5.2 Wm−2. Daily mean EPIC and SYN fluxes calculated using concurrent hours agree with each other to within 2% and both show a strong annual cycle. The SW flux agreement is within the calibration and algorithm uncertainties, which indicates that the method developed to calculate the global anisotropic factors from the CERES ADMs is robust and that the CERES ADMs accurately account for the Earth's anisotropy in the near-backscatter direction. CERES; angular distribution model; radiation; DSCOVR; EPIC; Lagrange-1 point

2017

Loeb, Norman G.; Wang, Hailan; Liang, Lusheng; Kato, Seiji; Rose, Fred G.Loeb, N. G., H. Wang, L. Liang, S. Kato, F. G. Rose, 2017: Surface energy budget changes over Central Australia during the early 21st century drought. International Journal of Climatology, 37(1), 159–168. doi: 10.1002/joc.4694. Satellite observations are used to investigate surface energy budget variability over central Australia during the early 21st century drought. Over a large expanse of open shrubland and savanna, surface albedo exhibits a multiyear increase of 0.06 during the drought followed by a sharp decline of 0.08 after heavy rainfall in 2010 broke the drought. The surface albedo variations are associated with increased normalized difference vegetation index (NDVI) during wet years before and after the drought and decreased NDVI during drought years. During the worst drought years (2002–2009), the surface albedo increase is most pronounced in the shortwave infrared region (wavelengths between 1 and 3 µm), implying soil moisture content variability is the likely cause of the albedo changes. At interannual timescales, surface albedo variability is associated with near-surface soil moisture, controlled by episodic precipitation events, whereas the multiyear increase in surface albedo is more closely linked with decreases in soil moisture in deeper surface layers. In addition to a higher surface albedo and lower soil moisture content during the drought, the observations show less evaporation, enhanced reflected shortwave radiation, increased upward emission of thermal infrared radiation, lower downwelling longwave (LW) radiation, reduced net total downward radiation, and higher sensible heating compared with the rainy period following the drought. Upward emission of thermal infrared radiation decreases sharply after the drought with increased surface evaporation. However, the surface energy budget changes during the worst drought years show a stronger relationship between upward emission of thermal radiation and reflected shortwave flux. During this period, evaporative fraction is extremely low and surface albedo is steadily increasing. In such extreme conditions, the surface albedo appears to modulate surface upward LW radiation, preventing it from getting too high. The change in upward LW radiation thus represents a negative feedback as it offsets further decreases in surface net radiation. albedo; radiation; energy budget; Precipitation; drought; Latent heat; Sensible heat
Su, W.; Liang, L.; Miller, W. F.; Sothcott, V. E.Su, W., L. Liang, W. F. Miller, V. E. Sothcott, 2017: The effects of different footprint sizes and cloud algorithms on the top-of-atmosphere radiative flux calculation from the Clouds and Earth's Radiant Energy System (CERES) instrument on Suomi National Polar-orbiting Partnership (NPP). Atmos. Meas. Tech., 10(10), 4001-4011. doi: 10.5194/amt-10-4001-2017. Only one Clouds and Earth's Radiant Energy System (CERES) instrument is onboard the Suomi National Polar-orbiting Partnership (NPP) and it has been placed in cross-track mode since launch; it is thus not possible to construct a set of angular distribution models (ADMs) specific for CERES on NPP. Edition 4 Aqua ADMs are used for flux inversions for NPP CERES measurements. However, the footprint size of NPP CERES is greater than that of Aqua CERES, as the altitude of the NPP orbit is higher than that of the Aqua orbit. Furthermore, cloud retrievals from the Visible Infrared Imaging Radiometer Suite (VIIRS) and the Moderate Resolution Imaging Spectroradiometer (MODIS), which are the imagers sharing the spacecraft with NPP CERES and Aqua CERES, are also different. To quantify the flux uncertainties due to the footprint size difference between Aqua CERES and NPP CERES, and due to both the footprint size difference and cloud property difference, a simulation is designed using the MODIS pixel-level data, which are convolved with the Aqua CERES and NPP CERES point spread functions (PSFs) into their respective footprints. The simulation is designed to isolate the effects of footprint size and cloud property differences on flux uncertainty from calibration and orbital differences between NPP CERES and Aqua CERES. The footprint size difference between Aqua CERES and NPP CERES introduces instantaneous flux uncertainties in monthly gridded NPP CERES measurements of less than 4.0 W m−2 for SW (shortwave) and less than 1.0 W m−2 for both daytime and nighttime LW (longwave). The global monthly mean instantaneous SW flux from simulated NPP CERES has a low bias of 0.4 W m−2 when compared to simulated Aqua CERES, and the root-mean-square (RMS) error is 2.2 W m−2 between them; the biases of daytime and nighttime LW flux are close to zero with RMS errors of 0.8 and 0.2 W m−2. These uncertainties are within the uncertainties of CERES ADMs. When both footprint size and cloud property (cloud fraction and optical depth) differences are considered, the uncertainties of monthly gridded NPP CERES SW flux can be up to 20 W m−2 in the Arctic regions where cloud optical depth retrievals from VIIRS differ significantly from MODIS. The global monthly mean instantaneous SW flux from simulated NPP CERES has a high bias of 1.1 W m−2 and the RMS error increases to 5.2 W m−2. LW flux shows less sensitivity to cloud property differences than SW flux, with uncertainties of about 2 W m−2 in the monthly gridded LW flux, and the RMS errors of global monthly mean daytime and nighttime fluxes increase only slightly. These results highlight the importance of consistent cloud retrieval algorithms to maintain the accuracy and stability of the CERES climate data record.
Su, Wenying; Loeb, Norman G.; Liang, Lusheng; Liu, Nana; Liu, ChuntaoSu, W., N. G. Loeb, L. Liang, N. Liu, C. Liu, 2017: The El Niño-Southern Oscillation Effect on Tropical Outgoing Longwave Radiation: A Daytime Versus Nighttime Perspective. Journal of Geophysical Research: Atmospheres, 122(15), 7820–7833. doi: 10.1002/2017JD027002. Trends of tropical (30° N-30° S) mean daytime and nighttime outgoing longwave radiation (OLR) from CERES and AIRS are analyzed using data from 2003 to 2013. Both the daytime and nighttime OLR from these instruments show decreasing trends because of El Niño conditions early in the period and La Niña conditions at the end. However, the daytime and nighttime OLR decrease at different rates with the OLR decreasing faster during daytime than nighttime. The daytime-nighttime OLR trend is consistent across CERES Terra, Aqua observations, and computed OLR based upon AIRS and MODIS retrievals. To understand the cause of the differing decreasing rates of daytime and nighttime OLR, high cloud fraction and effective temperature are examined using cloud retrievals from MODIS and AIRS. Unlike the very consistent OLR trends between CERES and AIRS, the trends in cloud properties are not as consistent, which is likely due to the different cloud retrieval methods used. When MODIS and AIRS cloud properties are used to compute OLR, the daytime and nighttime OLR trends based upon MODIS cloud properties are approximately half as large as the trends from AIRS cloud properties, but their daytime-nighttime OLR trends are in agreement. This demonstrates that though the current cloud retrieval algorithms lack the accuracy to pinpoint the changes of daytime and nighttime clouds in the tropics, they do provide a radiatively-consistent view for daytime and nighttime OLR changes. The causes for the larger decreasing daytime OLR trend than that for nighttime OLR are not clear and further studies are needed. clouds; 0321 Cloud/radiation interaction; ENSO; outgoing longwave radiation

2015

Su, W.; Corbett, J.; Eitzen, Z.; Liang, L.Su, W., J. Corbett, Z. Eitzen, L. Liang, 2015: Next-generation angular distribution models for top-of-atmosphere radiative flux calculation from CERES instruments: validation. Atmos. Meas. Tech., 8(8), 3297-3313. doi: 10.5194/amt-8-3297-2015. Radiative fluxes at the top of the atmosphere (TOA) from the Clouds and the Earth's Radiant Energy System (CERES) instrument are fundamental variables for understanding the Earth's energy balance and how it changes with time. TOA radiative fluxes are derived from the CERES radiance measurements using empirical angular distribution models (ADMs). This paper evaluates the accuracy of CERES TOA fluxes using direct integration and flux consistency tests. Direct integration tests show that the overall bias in regional monthly mean TOA shortwave (SW) flux is less than 0.2 Wm−2 and the RMSE is less than 1.1 Wm−2. The bias and RMSE are very similar between Terra and Aqua. The bias in regional monthly mean TOA LW fluxes is less than 0.5 Wm−2 and the RMSE is less than 0.8 Wm−2 for both Terra and Aqua. The accuracy of the TOA instantaneous flux is assessed by performing tests using fluxes inverted from nadir- and oblique-viewing angles using CERES along-track observations and temporally and spatially matched MODIS observations, and using fluxes inverted from multi-angle MISR observations. The averaged TOA instantaneous SW flux uncertainties from these two tests are about 2.3 % (1.9 Wm−2) over clear ocean, 1.6 % (4.5 Wm−2) over clear land, and 2.0 % (6.0 Wm−2) over clear snow/ice; and are about 3.3 % (9.0 Wm−2), 2.7 % (8.4 Wm−2), and 3.7 % (9.9 Wm−2) over ocean, land, and snow/ice under all-sky conditions. The TOA SW flux uncertainties are generally larger for thin broken clouds than for moderate and thick overcast clouds. The TOA instantaneous daytime LW flux uncertainties derived from the CERES-MODIS test are 0.5 % (1.5 Wm−2), 0.8 % (2.4 Wm−2), and 0.7 % (1.3 Wm−2) over clear ocean, land, and snow/ice; and are about 1.5 % (3.5 Wm−2), 1.0 % (2.9 Wm−2), and 1.1 % (2.1 Wm−2) over ocean, land, and snow/ice under all-sky conditions. The TOA instantaneous nighttime LW flux uncertainties are about 0.5–1 % (< 2.0 Wm−2) for all surface types. Flux uncertainties caused by errors in scene identification are also assessed by using the collocated CALIPSO, CloudSat, CERES and MODIS data product. Errors in scene identification tend to underestimate TOA SW flux by about 0.6 Wm−2 and overestimate TOA daytime (nighttime) LW flux by 0.4 (0.2) Wm−2 when all CERES viewing angles are considered.
Su, W.; Corbett, J.; Eitzen, Z.; Liang, L.Su, W., J. Corbett, Z. Eitzen, L. Liang, 2015: Next-generation angular distribution models for top-of-atmosphere radiative flux calculation from CERES instruments: methodology. Atmos. Meas. Tech., 8(2), 611-632. doi: 10.5194/amt-8-611-2015. The top-of-atmosphere (TOA) radiative fluxes are critical components to advancing our understanding of the Earth's radiative energy balance, radiative effects of clouds and aerosols, and climate feedback. The Clouds and the Earth's Radiant Energy System (CERES) instruments provide broadband shortwave and longwave radiance measurements. These radiances are converted to fluxes by using scene-type-dependent angular distribution models (ADMs). This paper describes the next-generation ADMs that are developed for Terra and Aqua using all available CERES rotating azimuth plane radiance measurements. Coincident cloud and aerosol retrievals, and radiance measurements from the Moderate Resolution Imaging Spectroradiometer (MODIS), and meteorological parameters from Goddard Earth Observing System (GEOS) data assimilation version 5.4.1 are used to define scene type. CERES radiance measurements are stratified by scene type and by other parameters that are important for determining the anisotropy of the given scene type. Anisotropic factors are then defined either for discrete intervals of relevant parameters or as a continuous functions of combined parameters, depending on the scene type. Significant differences between the ADMs described in this paper and the existing ADMs are over clear-sky scene types and polar scene types. Over clear ocean, we developed a set of shortwave (SW) ADMs that explicitly account for aerosols. Over clear land, the SW ADMs are developed for every 1° latitude × 1° longitude region for every calendar month using a kernel-based bidirectional reflectance model. Over clear Antarctic scenes, SW ADMs are developed by accounting the effects of sastrugi on anisotropy. Over sea ice, a sea-ice brightness index is used to classify the scene type. Under cloudy conditions over all surface types, the longwave (LW) and window (WN) ADMs are developed by combining surface and cloud-top temperature, surface and cloud emissivity, cloud fraction, and precipitable water. Compared to the existing ADMs, the new ADMs change the monthly mean instantaneous fluxes by up to 5 W m−2 on a regional scale of 1° latitude × 1° longitude, but the flux changes are less than 0.5 W m−2 on a global scale.

2014

Su, W.; Corbett, J.; Eitzen, Z.; Liang, L.Su, W., J. Corbett, Z. Eitzen, L. Liang, 2014: Next-generation angular distribution models for top-of-atmosphere radiative flux calculation from the CERES instruments: methodology. Atmos. Meas. Tech. Discuss., 7(8), 8817-8880. doi: 10.5194/amtd-7-8817-2014. The top-of-atmosphere (TOA) radiative fluxes are critical components to advancing our understanding of the Earth's radiative energy balance, radiative effects of clouds and aerosols, and climate feedback. The Clouds and Earth's Radiant Energy System (CERES) instruments provide broadband shortwave and longwave radiance measurements. These radiances are converted to fluxes by using scene type dependent Angular Distribution Models (ADMs). This paper describes the next-generation ADMs that are developed for Terra and Aqua using all available CERES rotating azimuth plane radiance measurements. Coincident cloud and aerosol retrievals, and radiance measurements from the Moderate Resolution Imaging Spectroradiometer (MODIS), and meteorological parameters from Goddard Earth Observing System (GEOS) data assimilation version 5.4.1 are used to define scene type. CERES radiance measurements are stratified by scene type and by other parameters that are important for determining the anisotropy of the given scene type. Anisotropic factors are then defined either for discrete intervals of relevant parameters or as a continuous functions of combined parameters, depending on the scene type. Compared to the existing ADMs, the new ADMs change the monthly mean instantaneous fluxes by up to 5 W m−2 on a regional scale of 1° latitude × 1° longitude, but the flux changes are less than 0.5 W m−2 on a global scale.

Additional Publications

2020

Doelling, D., Wrenn, F., & Liang, L. (2020). Radiative Transfer Models for Deriving Geostationary Broadband Shortwave Radiances Directly from Visible Channels for the Ceres SYN1deg Product. In IGARSS 2020-2020 IEEE International Geoscience and Remote Sensing Symposium (pp. 6690-6693). IEEE.

2019

Fu, D., Di Girolamo, L., Liang, L., & Zhao, G. (2019). Regional biases in MODIS marine liquid water cloud drop effective radius deduced through fusion with MISR. Journal of Geophysical Research: Atmospheres, 124(23), 13182-13196.

2015

Liang, L., L. Di Girolamo and W. Sun (2015), Bias in MODIS cloud drop effective radius for oceanic water clouds as deduced from optical thickness variability across scattering angles. Journal of Geophysical Research: Atmospheres, 120(15), 7661-7681.
Su, W., J. Corbett, Z. Eitzen, and L. Liang (2015): Corrigendum to “Next-generation angular distribution models for top-of-atmosphere radiative flux calculation from CERES instruments: methodology” published in Atmos. Meas. Tech., 8, 611–632, 2015, Atmos. Meas. Tech., 8,1467-1467, doi:10.5194/amt-8-1467-2015, 2015.

2013

Liang, L., and L. Di Girolamo (2013), A global analysis on the view-angle dependence of plane-parallel oceanic liquid water cloud optical thickness using data synergy from MISR and MODIS. Journal of Geophysical Research: Atmospheres, 118(5), 2389-2403.

2010

Di Girolamo, L., L. Liang, and S. Platnick (2010), A global view of one‐dimensional solar radiative transfer through oceanic water clouds, Geophys. Res. Lett., 37, L18809, doi:10.1029/2010GL044094.

2009

Liang, L., L. Di Girolamo and S. Platnick (2009), View-angle consistency in reflectance, optical thickness and spherical albedo of marine water-clouds over the northeastern Pacific through MISR-MODIS fusion, Geophys. Res. Lett., 36, L09811, doi:10.1029/2008GL037124.