Engineering topics
Su, Wenying
Publications and source records attributed to Su, Wenying.
Developing an Aircraft-Based Angular Distribution Model of Solar Reflection from Wildfire Smoke to Aid Satellite-Based Radiative Flux Estimation
This study examines the angular distribution of scattered solar radiation associated with wildfire smoke aerosols observed over boreal forests in Canada during the ARCTAS (Arctic Research of the Composition of the Troposphere from Aircraft and Satellites) campaign. First, it estimates smoke radiative parameters (550 nm optical depth of 3.9 and single scattering albedo of 0.90) using quasi-simultaneous multiangular and multispectral airborne measurements by the Cloud Absorption Radiometer (CAR). Next, the paper estimates the broadband top-of-atmosphere radiances that a satellite instrument such as the Clouds and the Earth’s Radiant Energy System (CERES) could have observed, given the narrowband CAR measurements made from an aircraft circling about a kilometer above the smoke layer. This estimation includes both an atmospheric correction that accounts for the atmosphere above the aircraft and a narrowband-to-broadband conversion. The angular distribution of estimated radiances is found to be substantially different than the angular model used in the operational data processing of CERES observations over the same area. This is because the CERES model is a monthly average model that was constructed using observations taken under smoke-free conditions. Finally, a sensitivity analysis shows that the estimated angular distribution remains accurate for a fairly wide range of smoke and underlying surface parameters. Overall, results from this work suggest that airborne CAR measurements can bring some substantial improvements in the accuracy of satellite-based radiative flux estimates.
Surface Irradiances of Edition 4.0 Clouds and the Earth's Radiant Energy System (CERES) Energy Balanced and Filled (EBAF) Data Product
The algorithm to produce the Clouds and the Earth’s Radiant Energy System (CERES) Edition 4.0 (Ed4) Energy Balanced and Filled (EBAF)-surface data product is explained. The algorithm forces computed topof- atmosphere (TOA) irradiances to match with Ed4 EBAF-TOA irradiances by adjusting surface, cloud, and atmospheric properties. Surface irradiances are subsequently adjusted using radiative kernels. The adjustment process is composed of two parts: bias correction and Lagrange multiplier. The bias in temperature and specific humidity between 200 and 500 hPa used for the irradiance computation is corrected based on observations by Atmospheric Infrared Sounder (AIRS). Similarly, the bias in the cloud fraction is corrected based on observations by Cloud–Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO) and CloudSat. Remaining errors in surface, cloud, and atmospheric properties are corrected in the Lagrange multiplier process. Ed4 global annual mean (January 2005 through December 2014) surface net shortwave (SW) and longwave (LW) irradiances increase by 1.3Wm22 and decrease by 0.2Wm22, respectively, compared to EBAF Edition 2.8 (Ed2.8) counterparts (the previous version), resulting in an increase in net SW 1LWsurface irradiance of 1.1Wm22. The uncertainty in surface irradiances over ocean, land, and polar regions at various spatial scales are estimated. The uncertainties in all-sky global annual mean upward and downward shortwave irradiance are 3 and 4Wm22, respectively, and the uncertainties in upward and downward longwave irradiance are 3 and 6Wm22, respectively. With an assumption of all errors being independent, the uncertainty in the global annual mean surface LW 1 SW net irradiance is 8Wm22.
Development of Multi-Sensor Global Cloud and Radiance Composites for DSCOVR EPIC Imager with Subpixel Definition
The Deep Space Climate Observatory (DSCOVR) enables analysis of the daytime Earth radiation budget via the onboard Earth Polychromatic Imaging Camera (EPIC) and National Institute of Standards and Technology Advanced Radiometer (NISTAR). EPIC delivers adequate spatial resolution imagery but only in shortwave bands (317-780 nm), while NISTAR measures the top-of-atmosphere (TOA) whole-disk radiance in shortwave and longwave broadband windows. Accurate calculation of albedo and outgoing longwave flux requires a high-resolution scene identification such as the radiance observations and cloud properties retrievals from low earth orbit (LEO, including NASA Terra and Aqua MODIS, Suomi-NPP VIIRS, and NOAA AVHRR) and geosynchronous (GEO, including GOES east and west, METEOSAT, INSAT-3D, MTSAT-2, and Himawari-8) satellite imagers. The cloud properties are derived using the Clouds and the Earth's Radiant Energy System (CERES) mission Cloud Subsystem group algorithms. These properties have to be co-located with EPIC pixels to provide the scene identification and to select anisotropic directional models (ADMs), which are then used to adjust the NISTAR-measured radiance and subsequently obtain the global daytime shortwave and longwave fluxes. This work presents an algorithm for optimal merging of selected radiance and cloud property parameters derived from multiple satellite imagers to obtain seamless global hourly composites at 5-km resolution. Selection of satellite data for each 5-km pixel is based on an aggregated rating that incorporates five parameters: nominal satellite resolution, pixel time relative to the EPIC time, viewing zenith angle, distance from day/night terminator, and probability of sun glint. To provide a smoother transition in the merged output, in regions where candidate pixel data from two satellite sources have comparable aggregated rating, the selection decision is defined by the cumulative function of the normal distribution so that abrupt changes in the visual appearance of the composite data are avoided. Higher spatial accuracy in the composite product is achieved by using the inverse mapping with gradient search during reprojection and bicubic interpolation for pixel resampling.
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)
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/sq. m for SW (shortwave) and less than 1.0 W/sq. m 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/sq. m when compared to simulated Aqua CERES, and the root-mean-square (RMS) error is 2.2 W/sq. m between them; the biases of daytime and night- time LW flux are close to zero with RMS errors of 0.8 and 0.2 W/sq. m. 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/sq. m 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/sq. m and the RMS error increases to 5.2 W/sq. m. LW flux shows less sensitivity to cloud property differences than SW flux, with uncertainties of about 2 W/sq. m 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.
Development of Multi-Sensor Global Cloud and Radiance Composites for Earth Radiation Budget Monitoring from DSCOVR
The Deep Space Climate Observatory (DSCOVR) enables analysis of the daytime Earth radiation budget via the onboard Earth Polychromatic Imaging Camera (EPIC) and National Institute of Standards and Technology Advanced Radiometer (NISTAR). Radiance observations and cloud property retrievals from low earth orbit and geostationary satellite imagers have to be co-located with EPIC pixels to provide scene identification in order to select anisotropic directional models needed to calculate shortwave and longwave fluxes. A new algorithm is proposed for optimal merging of selected radiances and cloud properties derived from multiple satellite imagers to obtain seamless global hourly composites at 5-kilometer resolution. An aggregated rating is employed to incorporate several factors and to select the best observation at the time nearest to the EPIC measurement. Spatial accuracy is improved using inverse mapping with gradient search during reprojection and bicubic interpolation for pixel resampling. The composite data are subsequently remapped into EPIC-view domain by convolving composite pixels with the EPIC point spread function (PSF) defined with a half-pixel accuracy. PSF-weighted average radiances and cloud properties are computed separately for each cloud phase. The algorithm has demonstrated contiguous global coverage for any requested time of day with a temporal lag of under 2 hours in over 95 percent of the globe.
The Role of DYNAMO in Situ Observations in Improving NASA Ceres-like Daily Surface and Atmospheric Radiative Flux Estimates
The daily surface and atmospheric radiative fluxes from NASA Clouds and the Earths RadiantEnergy System (CERES) Synoptic 1 degree (SYN1deg) Ed3A are among the most widely used data to studycloud-radiative feedback. The CERES SYN1deg data are based on Fu-Liou radiative transfer computations thatuse specific humidity (Q) and air temperature (T) from NASA Global Modeling and Assimilation Office (GMAO)reanalyses as inputs and are therefore subject to the quality of those fields. This study uses in situ Q and Tobservations collected during the Dynamics of the Madden-Julian Oscillation (DYNAMO) field campaign toaugment the input stream used in the NASA GMAO reanalysis and assess the impact on the CERES dailysurface and atmospheric longwave estimates. The results show that the assimilation of DYNAMOobservations considerably improves the vertical profiles of analyzed Q and T over and near DYNAMO stationsby moistening and warming the lower troposphere and upper troposphere and drying and cooling themid-upper troposphere. As a result of these changes in Q and T, the computed CERES daily surface downwardlongwave flux increases by about 5 W m(exp -2), due mainly to the warming and moistening in the lowertroposphere; the computed daily top-of-atmosphere (TOA) outgoing longwave radiation increases by2-3 W m(exp -2) during dry periods only. Correspondingly, the estimated local atmospheric longwave radiativecooling enhances by about 5 W m(exp -2) (7-8 W m(exp -2)) during wet (dry) periods. These changes reduce the bias inthe CERES SYN1deg-like daily longwave estimates at both the TOA and surface and represent animprovement over the DYNAMO region.
Construction of a Matched Global Cloud and Radiance Product from LEO/GEO and EPIC Observations to Estimate Daytime Earth Radiation Budget from DSCOVR
With the launch of the Deep Space Climate Observatory (DSCOVR), new estimates of the daytime Earth radiation budget can be computed from a combination of measurements from the two Earth-observing sensors onboard the spacecraft, the Earth Polychromatic Imaging Camera (EPIC) and the National Institute of Standards and Technology Advanced Radiometer (NISTAR). Although these instruments can provide accurate top-of-atmosphere (TOA) radiance measurements, they lack sufficient resolution to provide details on small-scale surface and cloud properties. Previous studies have shown that these properties have a strong influence on the anisotropy of the radiation at the TOA, and ignoring such effects can result in large TOA-flux errors. To overcome these effects, high-resolution scene identification is needed for accurate Earth radiation budget estimation. Selected radiance and cloud property data measured and derived from several low earth orbit (LEO, including NASA Terra and Aqua MODIS, NOAA AVHRR) and geosynchronous (GEO, including GOES (east and west), METEOSAT, INSAT-3D, MTSAT-2, and HIMAWARI-8) satellite imagers were collected to create hourly 5-km resolution global composites of data necessary to compute angular distribution models (ADM) for reflected shortwave (SW) and longwave (LW) radiation. The satellite data provide an independent source of radiance measurements and scene identification information necessary to construct ADMs that are used to determine the daytime Earth radiation budget. To optimize spatial matching between EPIC measurements and the high-resolution composite cloud properties, LEO/GEO retrievals within the EPIC fields of view (FOV) are convolved to the EPIC point spread function (PSF) in a similar manner to the Clouds and the Earth's Radiant Energy System (CERES) Single Scanner Footprint TOA/Surface Fluxes and Clouds (SSF) product. Examples of the merged LEO/GEO/EPIC product will be presented, describing the chosen radiance and cloud properties and details of how data from the multi-satellite measurements are selected.
Application of the CERES Flux-by-Cloud Type Simulator to GCM Output
The CERES Flux By CloudType data product produces CERES top-of-atmosphere (TOA) fluxes by region and cloud type. Here, the cloud types are defined by cloud optical depth (t) and cloud top pressure (pc), with bins similar to those used by ISCCP (International Satellite Cloud Climatology Project). This data product has the potential to be a powerful tool for the evaluation of the clouds produced by climate models by helping to identify which physical parameterizations have problems (e.g., boundary-layer parameterizations, convective clouds, processes that affect surface albedo). Also, when the flux-by-cloud type and frequency of cloud types are simultaneously used to evaluate a model, the results can determine whether an unrealistically large or small occurrence of a given cloud type has an important radiative impact for a given region. A simulator of the flux-by-cloud type product has been applied to three-hourly data from the year 2008 from the UK Met Office HadGEM2-A model using the Langley Fu-Lour radiative transfer model to obtain TOA SW and LW fluxes.
Consistency of CERES Radiances and Fluxes from Aqua and Suomi-NPP
The Clouds and Earth's Radiant Energy System (CERES) instruments on board Terra, Aqua, and Suomi-NPP have been providing data products critical to advancing our understanding of the effects of clouds and aerosols on radiative energy within the Earth-atmosphere system. The CERES instrument consists of a threechannel broadband scanning radiometer. The scanning radiometer measures radiances in shortwave (SW, 0.3-5 micron), window (WN, 8-12 micron), and total (0.3-200 micron) channels. The longwave (LW) component is derived as the difference between total and SW channels. These measured radiances at a given sun-Earthsatellite geometry are converted to outgoing reflected solar and emitted thermal TOA radiative fluxes by using CERES scene-type dependent angular distribution models (ADMs). The CERES instruments must remain radiometrically stable and correctly inter-calibrated to accurately capture changes in Earth"s radiation budget from interannual to decadal timescales. This presentation will focus on comparisons between collocated radiance measurements from CERES instruments on Aqua and on Suomi-NPP. As we do not have a set of ADMs that is constructed specifically for the CERES instrument on Suomi-NPP, CERES Aqua ADMs are used to invert fluxes from radiance measurements on Suomi-NPP. But the CERES Aqua footprint size is smaller than the CERES Suomi-NPP footprint size and the scene identifications provided by MODIS and VIIRS can also be different from each other. Will using Aqua ADMs for Suomi-NPP flux inversion increase the flux uncertainty? We will examine the deseasonalized flux anomaly time series using Aqua data alone and using combined Aqua and Suomi-NPP data. We will also present a simulation study to assess the Suomi-NPP flux uncertainty from using Aqua ADMs for the flux inversion.
The ENSO Effects on Tropical Clouds and Top-of-Atmosphere Cloud Radiative Effects in CMIP5 Models
The El Nino-Southern Oscillation (ENSO) effects on tropical clouds and top-of-atmosphere (TOA) cloud radiative effects (CREs) in Coupled Model Intercomparison Project Phase5 (CMIP5) models are evaluated using satellite-based observations and International Satellite Cloud Climatology Project satellite simulator output. Climatologically, most CMIP5 models produce considerably less total cloud amount with higher cloud top and notably larger reflectivity than observations in tropical Indo-Pacific (60 degrees East - 200 degrees East; 10 degrees South - 10 degrees North). During ENSO, most CMIP5 models considerably underestimate TOA CRE and cloud changes over western tropical Pacific. Over central tropical Pacific, while the multi-model mean resembles observations in TOA CRE and cloud amount anomalies, it notably overestimates cloud top pressure (CTP) decreases; there are also substantial inter-model variations. The relative effects of changes in cloud properties, temperature and humidity on TOA CRE anomalies during ENSO in the CMIP5 models are assessed using cloud radiative kernels. The CMIP5 models agree with observations in that their TOA shortwave CRE anomalies are primarily contributed by total cloud amount changes, and their TOA longwave CRE anomalies are mostly contributed by changes in both total cloud amount and CTP. The model biases in TOA CRE anomalies particularly the strong underestimations over western tropical Pacific are, however, mainly explained by model biases in CTP and cloud optical thickness (tau) changes. Despite the distinct model cloud biases particularly in tau regime, the TOA CRE anomalies from cloud amount changes are comparable between the CMIP5 models and observations, because of the strong compensations between model underestimation of TOA CRE anomalies from thin clouds and overestimation from medium and thick clouds.
Dome Degradation Characterization of Wide-Field-of-View Nonscanner Aboard ERBE and Its Reprocessing
Earth Radiation Budget Experiment (ERBE) wide-field-of-view (WFOV) nonscanners aboard ERBS and NOAA- 9/NOAA-10 provided broadband shortwave and longwave irradiances from 1985 to 1999. The previous analysis showed dome degradation in the shortwave nonscanner instruments. The correction was performed with a constant spectral (gray assumption) degradation. We suspect that the gray assumption affected daytime longwave irradiance and led to a day-minus-night longwave flux differences (little change in night time longwave) increase over time. Based on knowledge from the CERES process, we will reprocess entire ERBE nonscanner radiation dataset by characterizing shortwave dome transmissivity with spectral dependent degradation using the solar data observed by these instruments. Once spectral dependent degradation is derived, imager derived cloud fraction and the cloud phase as well as surface type over the FOV of nonscanner instruments will be used to model unfiltering coefficients. This poster primarily explains the reprocessing techniques and includes initial comparison of several months of data processed with existing and our recent methods.
Aerosol and Cloud Interaction Observed From High Spectral Resolution Lidar Data
Recent studies utilizing satellite retrievals have shown a strong correlation between aerosol optical depth (AOD) and cloud cover. However, these retrievals from passive sensors are subject to many limitations, including cloud adjacency (or 3D) effects, possible cloud contamination, uncertainty in the AOD retrieval. Some of these limitations do not exist in High Spectral Resolution Lidar (HSRL) observations; for instance, HSRL observations are not a ected by cloud adjacency effects, are less prone to cloud contamination, and offer accurate aerosol property measurements (backscatter coefficient, extinction coefficient, lidar ratio, backscatter Angstrom exponent,and aerosol optical depth) at a neospatial resolution (less than 100 m) in the vicinity of clouds. Hence, the HSRL provides an important dataset for studying aerosol and cloud interaction. In this study, we statistically analyze aircraft-based HSRL profiles according to their distance from the nearest cloud, assuring that all profile comparisons are subject to the same large-scale meteorological conditions. Our results indicate that AODs from HSRL are about 17% higher in the proximity of clouds (approximately 100 m) than far away from clouds (4.5 km), which is much smaller than the reported cloud 3D effect on AOD retrievals. The backscatter and extinction coefficients also systematically increase in the vicinity of clouds, which can be explained by aerosol swelling in the high relative humidity (RH) environment and/or aerosol growth through in cloud processing (albeit not conclusively). On the other hand, we do not observe a systematic trend in lidar ratio; we hypothesize that this is caused by the opposite effects of aerosol swelling and aerosol in-cloud processing on the lidar ratio. Finally, the observed backscatter Angstrom exponent (BAE) does not show a consistent trend because of the complicated relationship between BAE and RH. We demonstrate that BAE should not be used as a surrogate for Angstrom exponent, especially at high RH.
Sunglint Research in the Lab, and Spatial Micro-Complexity
We are beginning a new investigation into polarized sunglint in a laboratory setting, using a wave tank. The wave tank has the capability of generating almost any combination of gravity and capillary waves under computer control, from simple sine waves to essentially random seas. Our aim is to perfect sunglint models, which have advanced little since the work of Cox and Munk in the 1950s using digitized photographs taken from an airplane. Sungllnt patterns, even when viewed from elevated platforms, but more and more so when viewed from low Earth orbit and beyond, have a kind of coherence and stability which transcends a pixel-by-pixel interpretation. The light forming the pattern is a sum of contributions from several or even many waves which are far from adjacent. Thus glint is a phenomenon that can only be understood in the large and never on a pixel-by-pixel basis.