Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reflection coefficient”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 919 records · Page 51

Turbulent diffusion with memories and intrinsic shear

The first part of the present theory is devoted to the derivation of a Fokker-Planck equation. The eddies smaller than the hydrodynamic scale of the diffusion cloud form a diffusivity, while the inhomogeneous, bigger eddies give rise to a nonuniform migratory drift. This introduces an eddy-induced shear which reflects on the large-scale diffusion. The eddy-induced shear does not require the presence of a permanent wind shear and is intrinsic to the diffusion. Secondly, a transport theory of diffusivity is developed by the method of repeated-cascade and is based upon a relaxation of a chain of memories with decreasing information. The full range of diffusion consists of inertia, composite, and shear subranges, for which variance and eddy diffusivities are predicted. The coefficients are evaluated. Comparison with experiments in the upper atmosphere and oceans is made.

Tchen, C. M.↗

Ice-Soil Mixtures: Visual and Near-Infrared Remote Sensing Techniques

The spectral properties of ice soil mixtures are reviewed. Water ice has an absorption coefficient which varies by several orders of magnitude in the visual and near infrared, and has several prominent overtone absorptions. Different wavelengths can be used to probe to different depths in the surface as well as for different mineral impurity concentrations. Empirical and/or theoretical models might be used to derive characteristics grain size of the ice or of the impurity minerals and to derive abundance of the ice and rock or soil components. Quantitative analysis of remotely obtained reflectance spectra can only be performed by using absorption features in the spectra of ice and soil, and not by broadband response. It is suggested that this might be done by special selection of several narrow band filters in the near infrared which will adequately define the ice absorptions.

Clark, R. N.↗

The baroeffect and an appropriate momentum boundary condition

As disparate molecular weight gases isothermally diffuse between two ends of a capillary tube, they can support a pressure gradient. The magnitude of this pressure gradient depends critically on viscous wall stress and becomes a measure of the boundary condition. Here, a one-dimensional analytical model is proposed that allows for specular gas reflection from the wall. It explains anomalous (4/3) correction factors required previously to match experiment to baroeffect models. It predicts a new physical phenomenon, a surface-driven baroeffect for equal molecular weight gases. Diffusive slip contributions exceed the order of convective diffusion for Peclet number Pe less than 1 and approximately equal convective diffusion for Pe values between 1 and 4. For binary gases, this model further extends the baroeffect experiment as a means to find momentum accommodation coefficients, slip lengths, and external friction coefficients. A corrected binary diffusion coefficient is defined.

Noever, David A.↗

On Utilization of NEXRAD Scan Strategy Information to Infer Discrepancies Associated With Radar and Rain Gauge Surface Volumetric Rainfall Accumulations

To evaluate the Tropical Rainfall Measuring Mission (TRMM) monthly Ground Validation (GV) rain map, 42 quality controlled tipping bucket rain gauge data (1 minute interpolated rain rates) were utilized. We have compared the gauge data to the surface volumetric rainfall accumulation of NEXRAD reflectivity field, (converting to rain rates using a 0.5 dB resolution smooth Z-R table). The comparison was carried out from data collected at Melbourne, Florida during the month of July 98. GV operational level 3 (L3 monthly) accumulation algorithm was used to obtain surface volumetric accumulations for the radar. The gauge records were accumulated using the 1 minute interpolated rain rates while the radar Volume Scan (VOS) intervals remain less than or equal to 75 minutes. The correlation coefficient for the radar and gauge totals for the monthly time-scale remain at 0.93, however, a large difference was noted between the gauge and radar derived rain accumulation when the radar data interval is either 9 minute, or 10 minute. This difference in radar and gauge accumulation is being explained in terms of the radar scan strategy information. The discrepancy in terms of the Volume Coverage Pattern (VCP) of the NEXRAD is being reported where VCP mode is ascertained using the radar tilt angle information. Hourly radar and gauge accumulations have been computed using the present operational L3 method supplemented with a threshold period of +/- 5 minutes (based on a sensitivity analysis). These radar and gauge accumulations are subsequently improved using a radar hourly scan weighting factor (taking ratio of the radar scan frequency within a time bin to the 7436 total radar scans for the month). This GV procedure is further being improved by introducing a spatial smoothing method to yield reasonable bulk radar to gauge ratio for the hourly and daily scales.

Roy, Biswadev↗

Impeller flow field characterization with a laser two-focus velocimeter

Use of Computational Fluid Dynamics (CFD) codes, prevalent in the rocket engine turbomachinery industry, necessitates data of sufficient quality and quantity to benchmark computational codes. Existing data bases for typical rocket engine configurations, in particular impellers, are limited. In addition, traditional data acquisition methods have several limitations: typically transducer uncertainties are 0.5% of transducer full scale and traditional pressure probes are unable to provide flow characteristics in the circumferential (blade-to-blade) direction. Laser velocimetry circumvents these limitations by providing +0.5% uncertainty in flow velocity and +0.5% uncertainty in flow angle. The percent of uncertainty in flow velocity is based on the measured value, not full range capability. The laser electronics multiple partitioning capability allows data acquired between blades as the impeller rotates, to be analyzed separately, thus providing blade-to-blade flow characterization. Unlike some probes, the non-intrusive measurements made with the laser velocimeter does not disturb the flow. To this end,, and under Contract (NAS8-38864) to the National Aeronautics and Space Administration (NASA) at Marshall Space Flight Center (MSFC), an extensive test program was undertaken at Rocketdyne. Impellers from two different generic rocket engine pump configurations were examined. The impellers represent different spectrums of pump design: the Space Shuttle Main Engine (SSME) high pressure fuel turbopump (HPFTP) impeller was designed in the 1 1970's the Consortium for CFD application in Propulsion Technology Pump Stage Technology Team (Pump Consortium) optimized impeller was designed with the aid of modern computing techniques. The tester configuration for each of the impellers consisted of an axial inlet, an inducer, a diffuser, and a crossover discharge. While the tested configurations were carefully chosen to be representative of generic rocket engine pumps, several features of both testers were intentionally atypical. A crossover discharge, downstream of the impeller, rather than a volute discharge was used to minimize asymmetric flow conditions that might be reflected in the impeller discharge flow data. Impeller shroud wear ring radial clearances were purposely close to minimize leakage flow, thus increasing confidence in using the inlet data as an input to CFD programs. The empirical study extensively examined the flow fields of the two impellers via performance of laser two-focus velocimeter surveys in an axial plane upstream of the impellers and in multiple radial planes downstream of the impellers. Both studies were performed at the impeller design flow coefficients. Inlet laser surveys that provide CFD code inlet boundary conditions were performed in one axial plane, with ten radial locations surveyed. Three wall static pressures, positioned circumferentially around the impeller inlet, were used to identify asymmetrical pressure distributions in the inlet survey plane. The impeller discharge flow characterization consisted of three radial planes for the SSME HPFTP impeller and two radial planes for the Pump Consortium optimized impeller. Housing wall static pressures were placed to correspond to the radial locations surveyed with the laser velocimeter. Between five and thirteen axial stations across the discharge channel width were examined in each radial plane during the extensive flow mapping. The largely successful empirical flow characterization of two different impellers resulted in a substantial contribution to the limited existing data base, and yielded accurate data for CFD code benchmarking.

Brozowski, L. A.↗

Investigation, development and application of optimal output feedback theory. Vol. 4: Measures of eigenvalue/eigenvector sensitivity to system parameters and unmodeled dynamics

Some measures of eigenvalue and eigenvector sensitivity applicable to both continuous and discrete linear systems are developed and investigated. An infinite series representation is developed for the eigenvalues and eigenvectors of a system. The coefficients of the series are coupled, but can be obtained recursively using a nonlinear coupled vector difference equation. A new sensitivity measure is developed by considering the effects of unmodeled dynamics. It is shown that the sensitivity is high when any unmodeled eigenvalue is near a modeled eigenvalue. Using a simple example where the sensor dynamics have been neglected, it is shown that high feedback gains produce high eigenvalue/eigenvector sensitivity. The smallest singular value of the return difference is shown not to reflect eigenvalue sensitivity since it increases with the feedback gains. Using an upper bound obtained from the infinite series, a procedure to evaluate whether the sensitivity to parameter variations is within given acceptable bounds is developed and demonstrated by an example.

Halyo, Nesim↗

A 10-Year Comparison of Water Levels Measured with a Geodetic GPS Receiver Versus a Conventional Tide Gauge

A standard geodetic GPS receiver and a conventional Aquatrak tide gauge, collocated at Friday Harbor, Washington, are used to assess the quality of 10 years of water levels estimated from GPS sea surface reflections.The GPS results are improved by accounting for (tidal) motion of the reflecting sea surface and for signal propagation delay by the troposphere. The RMS error of individual GPS water level estimates is about 12 cm. Lower water levels are measured slightly more accurately than higher water levels. Forming daily mean sea levels reduces the RMS difference with the tide gauge data to approximately 2 cm. For monthly means, the RMS difference is 1.3 cm. The GPS elevations, of course, can be automatically placed into a well-defined terrestrial reference frame. Ocean tide coefficients, determined from both the GPS and tide gauge data, are in good agreement, with absolute differences below 1 cm for all constituents save K1 and S1. The latter constituent is especially anomalous, probably owing to daily temperature-induced errors in the Aquatrak tide gauge

Larson, Kristine M.↗

Enhancements to the Open Access Spectral Band Adjustment Factor Online Calculation Tool for Visible Channels

With close to 40 years of satellite observations, from which, cloud, land-use, and aerosol parameters can be measured, inter-consistent calibrations are needed to normalize retrievals across satellite records. Various visible-sensor inter-calibration techniques have been developed that utilize radiometrically stable Earth targets, e.g., deep convective clouds and desert/polar ice pseudo-invariant calibration sites. Other equally effective, direct techniques for intercalibration between satellite imagers are simultaneous nadir overpass comparisons and ray-matched radiance pairs. Combining independent calibration results from such varied techniques yields robust calibration coefficients, and is a form of self-validation. One potential source of significant error when cross-calibrating satellite sensors, however, are the often small but substantial spectral discrepancies between comparable bands, which must be accounted for. As such, visible calibration methods rely on a Spectral Band Adjustment Factor (SBAF) to account for the spectral-response function- induced radiance differences between analogous imagers. The SBAF is unique to each calibration method as it is a function of the Earth-reflected spectra. In recent years, NASA Langley pioneered the use of SCIAMACHY-, GOME-2-, and Hyperion-retrieved Earth spectra to compute SBAFs. By carefully selecting hyperspectral footprints that best represent the conditions inherent to an inter-calibration technique, the uncertainty in the SBAF is greatly reduced. NASA Langley initially provided the Global Space-based Inter-calibration System processing and research centers with online SBAF tools, with which users select conditions to best match their calibration criteria. This article highlights expanded SBAF tool capabilities for visible wavelengths, with emphasis on the use of the spectral range filtering for the purpose of separating scene conditions for the channel that the SBAF is needed based on the reflectance values of other bands. In other words, spectral filtering will enable better scene-type selection for bands where scene determination is difficult without information from other channels, which should prove valuable to users in the calibration community.

Scarino, Benjamin↗

The Inter-Calibration of the DSCOVR EPIC Imager 2 with Aqua-MODIS and NPP-VIIRS

The Deep Space Climate Observatory (DSCOVR) Earth Polychromatic Imaging Camera (EPIC) continuously observes the illuminated disk from the Lagrange-1 point. The EPIC sensor was designed to monitor the diurnal variation of ozone, clouds, aerosols, and vegetation, especially those features that benefit from observation near-backscatter conditions. The EPIC sensor does not contain any onboard calibration systems. This study describes the inter-calibration of EPIC channels 5 (0.44 μm), 6 (0.55 μm), 7 (0.68 μm), and 10 (0.78 μm) with respect to Aqua-MODIS and NPP-VIIRS. The calibration is transferred using coincident ray-matched reflectance pairs over all-sky tropical ocean (ATO) and deep convective cloud (DCC) targets. A robust and automated image-alignment technique based on feature matching was formulated to improve the navigation quality of the EPIC images. The EPIC V02 dataset exhibits improved navigation over V01. Because the visible channels display similar spatial features, a single visible channel can be used to co-register the remaining visible bands. The VIIRS-referenced EPIC ATO and DCC ray-matched calibration coefficients are within 0.3%. The EPIC four-year calibration trends based on VIIRS are within 0.15%/year. The MODIS-based EPIC calibration coefficients were compared against the Geogdzhayev and Marshak 2018 published calibration coefficients and are found to be within 1.6%.

DSCOVR↗

TROPOMI Geometry-dependent Lambertian-Equivalent surface Reflectivity (GLER) Product for Improved Trace-Gas Retrieval

Accurate information about the reflectivity of the Earth's surface is required for most satellite retrievals of atmospheric composition, and this information is generally taken from monthly surface reflectivity climatology that neglects angular dependence. Previously we introduced Geometry-dependent Lambertian-equivalent surface reflectivity (GLER) which captures solar and satellite viewing angle dependence as well as daily and seasonal changes. GLER is calculated from simulations of Rayleigh-only top-of-atmosphere (TOA) radiances over non-Lambertian surfaces. We use NASA's Moderate Resolution Imaging Spectroradiometer (MODIS) bidirectional reflectance distribution function (BRDF) product over land and the wind-dependent Cox–Munk wave-facet slope distribution including water-leaving radiance over water to accounts for surface BRDF. We have developed global GLER product, previously for the Ozone Monitoring Instrument (OMI) and recently for Sentinel-5 Precursor (S5P) TROPOspheric Monitoring Instrument (TROPOMI) with several new improvements and updates. We have implemented the near real time daily V006 MODIS MCD43C1 BRDF data and gap-filled with a daily BRDF coefficient climatology created from 2002-2017 V006 MCD43GF data. The NASA’s Global Modelling Initiative hourly 0.25 x 0.25 deg Replay simulations are used for more accurate determination of pixel specific terrain pressure. To improve detection of seasonal snow/ice scenes, we use the 4-km snow cover product from the Interactive Multi-sensor Snow and Ice Mapping System (IMS). Finally, we use an improved version of the vector linearized discrete ordinate radiative transfer (VLIDORT) for update of the top-of-atmosphere (TOA) radiance look-up-tables (LUTs). We demonstrate how the use of GLER is beneficial to TROPOMI’s high spatial resolution (up to 3.5 km x 3.5 km) measurements to monitor atmospheric trace gas pollutants down to the sub-city scale.

TROPOMI↗

A 70th Degree Lunar Gravity Model (GLGM-2) from Clementine and other tracking data

A spherical harmonic model of the lunar gravity field complete to degree and order 70 has been developed from S band Doppler tracking data from the Clementine mission, as well as historical tracking data from Lunar Orbiters 1-5 and the Apollo 15 and 16 subsatellites. The model combines 361,000 Doppler observations from Clementine with 347,000 historical observations. The historical data consist of mostly 60-s Doppler with a noise of 0.25 to several mm/s. The Clementine data consist of mostly 10-s Doppler data, with a data noise of 0.25 mm/s for the observations from the Deep Space Network, and 2.5 mm/s for the data from a naval tracking station at Pomonkey, Maryland. Observations provided Clementine, provide the strongest satellite constraint on the Moon's low-degree field. In contrast the historical data, collected by spacecraft that had lower periapsis altitudes, provide distributed regions of high-resolution coverage within +/- 29 deg of the nearside lunar equator. To obtain the solution for a high-degree field in the absence of a uniform distribution of observations, we applied an a priori power law constraint of the form 15 x 10(exp -5)/sq l which had the effect of limiting the gravitational power and noise at short wavelengths. Coefficients through degree and order 18 are not significantly affected by the constraint, and so the model permits geophysical analysis of effects of the major basins at degrees 10-12. The GLGM-2 model confirms major features of the lunar gravity field shown in previous gravitational field models but also reveals significantly more detail, particularly at intermediate wavelengths (10(exp 3) km). Free-air gravity anomaly maps derived from the new model show the nearside and farside highlands to be gravitationally smooth, reflecting a state of isostatic compensation. Mascon basins (including Imbrium, Serenitatis, Crisium, Smythii, and Humorum) are denoted by gravity highs first recognized from Lunar Orbiter tracking. All of the major mascons are bounded by annuli of negative anomalies representing significant subsurface mass deficiencies. Mare Orientale appears as a minor mascon surrounded by a horseshoe-shaped gravity low centered on the Inner and Outer Rook rings that is evidence of significant subsurface structural heterogeneity. Although direct tracking is not available over a significant part of the lunar farside, GLGM-2 resolves negative anomalies that correlate with many farside basins, including South Pole-Aitken, Hertzsprung, Korolev, Moscoviense, Tsiolkovsky, and Freundlich-Sharonov.

Lemonie, Frank G. R.↗

The correlation between far-IR and radio continuum emission from spiral galaxies

A sample of 30 galaxies selected for their intense IRAS flux at 60 and 100 micron using the Arecibo telescope at 21 cm to measure the continuum and HI line luminosities were observed. The centimeter wave continuum correlates very well with the far-infrared flux, with a correlation coefficient as high as that found for other samples, and the same ratio between FIR and radio luminosities. Weaker correlations are seen between the FIR and optical luminosity and between the FIR and radio continuum. There is very little correlation between the FIR and the HI mass deduced from the integral of the 21 cm line. The strength of the radio continuum correlation suggests that there is little contribution to either the radio and FIR from physical processes not affecting both. If they each reflect time integrals of the star formation rate then the time constants must be similar, or the star formation rate must change slowly in these galaxies.

Dickey, John M.↗

TPSAS-NF1676L-28148-DND

Currently, ocean carbon biomass is estimated from ocean color measurements primarily through particulate backscatter coefficients (bbp). Small non-algal particles contribute significantly to bbp due to the fact that small particles scatter more isotropically while bigger particles scatter mostly in the forward direction. CALIPSO’s ocean subsurface cross polarization measurements (Figure 1) are not sensitive to small particles far less than 2 microns in size since these particles rarely cause depolarization. Thus, cross polarization backscatter is contributed mostly by larger, depolarizing phytoplankton particles. This is why CALIPSO’s cross polarization has much larger contrasts comparing with MODIS reflectance (CALIPSO cross polarization is very weak in Gyres and much stronger in high latitudes, See Figure 1). CALIPSO cross-polarization measurements in the Southern Oceans also suggest that water cloud droplet number concentration covaries with phytoplankton carbon biomass (Figure 2), but it has less correlation with particulate organic carbon (POC) (Figure 2) and chlorophyll concentration. The study also agrees with water cloud droplet size retrieved from Parasol’s polarization measurements.

Yongxiang Hu↗

First Provisional Land Surface Reflectance Product from Geostationary Satellite Himawari-8 AHI

A provisional surface reflectance (SR) product from the Advanced Himawari Imager (AHI) on-board the new generation geostationary satellite (Himawari-8) covering the period between July 2015 and December 2018 is made available to the scientific community. The Multi-Angle Implementation of Atmospheric Correction (MAIAC) algorithm is used in conjunction with time series Himawari-8 AHI observations to generate 1-km gridded and tiled land SR every 10 minutes during day time. This Himawari-8 AHI SR product includes retrieved atmospheric properties (e.g., aerosol optical depth at 0.47μm and 0.51μm), spectral surface reflectance (AHI bands 1–6), parameters of the RTLS BRDF model, and quality assurance flags. Product evaluation shows that Himawari-8 AHI data on average yielded 35% more cloud-free, valid pixels in a single day when compared to available data from the low earth orbit (LEO) satellites Terra/Aqua with MODIS sensor. Comparisons of Himawari-8 AHI SR against corresponding MODIS SR products (MCD19A1) over a variety of land cover types with the similar viewing geometry show high consistency between them, with correlation coefficients (r) being 0.94 and 0.99 for red and NIR bands, respectively. The high-frequency geostationary data are expected to facilitate studies of ecosystems on daily to diurnal time scales, complementing observations from networks such as the FLUXNET.

Himawari-8 AHI↗

SWIM: A Semi-Analytical Ocean Color Inversion Algorithm for Optically Shallow Waters

In clear shallow waters, light that is transmitted downward through the water column can reflect off the sea floor and thereby influence the water-leaving radiance signal. This effect can confound contemporary ocean color algorithms designed for deep waters where the seafloor has little or no effect on the water-leaving radiance. Thus, inappropriate use of deep water ocean color algorithms in optically shallow regions can lead to inaccurate retrievals of inherent optical properties (IOPs) and therefore have a detrimental impact on IOP-based estimates of marine parameters, including chlorophyll-a and the diffuse attenuation coefficient. In order to improve IOP retrievals in optically shallow regions, a semi-analytical inversion algorithm, the Shallow Water Inversion Model (SWIM), has been developed. Unlike established ocean color algorithms, SWIM considers both the water column depth and the benthic albedo. A radiative transfer study was conducted that demonstrated how SWIM and two contemporary ocean color algorithms, the Generalized Inherent Optical Properties algorithm (GIOP) and Quasi-Analytical Algorithm (QAA), performed in optically deep and shallow scenarios. The results showed that SWIM performed well, whilst both GIOP and QAA showed distinct positive bias in IOP retrievals in optically shallow waters. The SWIM algorithm was also applied to a test region: the Great Barrier Reef, Australia. Using a single test scene and time series data collected by NASA's MODIS-Aqua sensor (2002-2013), a comparison of IOPs retrieved by SWIM, GIOP and QAA was conducted.

shallow waters↗

Evaluation of spatial, radiometric and spectral thematic mapper performance for coastal studies

The main emphasis of the research was to determine what effect different wetland plant canopies would have upon observed reflectance in Thematic Mapper bands. The three major vegetation canopy types (broadleaf, gramineous and leafless) produce unique spectral responses for a similar quantity of live biomass. Biomass estimates computed from spectral data were most similar to biomass estimates determined from harvest data when models developed for a specific canopy were used. In other words, the spectral biomass estimate of a broadleaf canopy was most similar to the harvest biomass estimate when a broadleaf canopy radiance model was used. Work is continuing to more precisely determine regression coefficients for each canopy type and to model the change in the coefficients with various combinations of canopy types. Researchers suspect that textural and spatial considerations can be used to identify canopy types and improve biomass estimates from Thematic Mapper data.

Klemas, V.↗

A Principal Component and Machine Learning Approach to Spatially Gap Fill Hyperspectral Ocean Color Satellite Retrievals

Retrievals of ocean color properties from space are important for monitoring the health of the ocean ecosystem but such retrievals tend to be limited in spatial coverage due to conditions such as clouds, aerosols, and sun glint. Gap filling of ocean color retrievals is typically performed by combining retrievals from multiple satellites or temporally averaging multiple days of retrievals but despite these techniques large gaps still exist posing challenges for near real time monitoring of events like harmful algae blooms. To address these limitations, we propose a spatial gap filling approach using machine learning to learn how to perform an atmospheric correction under challenging retrieval conditions. In this approach a principal component analysis is used to decompose the hyperspectral measurements into spectral features that describe the scattering and absorption of the atmosphere as well as the underlying surface. The coefficients of the principal components are then used to train a neural network to predict ocean color properties derived from a standard ocean color algorithm such as the MODIS atmospheric correction algorithm. This machine learning approach is independent of a priori information and does not rely on any radiative transfer modeling. We apply the approach to two hyperspectral UV/VIS instruments, the Ozone Monitoring Instrument (OMI) and TROPOspheric Monitoring Instrument (TROPOMI) to show that it can be used to estimate ocean color properties such as chlorophyll, remote sensing reflectance, and fluorescence line height. This method could be used as a gap-filling technique for the future Ocean Color Instrument (OCI) which will be onboard NASA's Plankton, Aerosol Cloud, ocean Ecosystem (PACE) ocean color satellite to provide additional information for monitoring the health of our global oceans. Additionally, it could be applied to the geostationary satellite Tropospheric Emissions: Monitoring of Pollution (TEMPO) to better understand diurnal variability in ocean ecology.

MODIS atmospheric correction algorithm↗

Satellite-Based Evaluation of AeroCom Model Bias in Biomass Burning Regions

Global models are widely used to simulate biomass burning aerosol (BBA). Exhaustive evaluations on model representation of aerosol distributions and properties are fundamental to assess health and climate impacts of BBA. Here we conducted a comprehensive comparison of Aerosol Comparisons between Observations and Models (AeroCom) project model simulations with satellite observations. A total of 59 runs by 18 models from three AeroCom Phase-III experiments (i.e., biomass burning emissions, CTRL16, and CTRL19) and 14 satellite products of aerosols were used in the study. Aerosol optical depth (AOD) at 550 nm was investigated during the fire season over three key fire regions reflecting different fire dynamics (i.e., deforestation-dominated Amazon, Southern Hemisphere Africa where savannas are the key source of emissions, and boreal forest burning in boreal North America). The 14 satellite products were first evaluated against AErosol RObotic NETwork (AERONET) observations, with large uncertainties found. But these uncertainties had small impacts on the model evaluation that was dominated by modeling bias. Through a comparison with Polarization and Directionality of the Earth’s Reflectances measurements with the Generalized Retrieval of Aerosol and Surface Properties algorithm (POLDER-GRASP), we found that the modeled AOD values were biased by −93 % to 152 %, with most models showing significant underestimations even for the state-of-the-art aerosol modeling techniques (i.e., CTRL19). By scaling up BBA emissions, the negative biases in modeled AOD were significantly mitigated, although it yielded only negligible improvements in the correlation between models and observations, and the spatial and temporal variations in AOD biases did not change much. For models in CTRL16 and CTRL19, the large diversity in modeled AOD was in almost equal measures caused by diversity in emissions, lifetime, and the mass extinction coefficient (MEC). We found that in the AeroCom ensemble, BBA lifetime correlated significantly with particle deposition (as expected) and in turn correlated strongly with precipitation. Additional analysis based on Cloud-Aerosol LIdar with Orthogonal Polarization (CALIOP) aerosol profiles suggested that the altitude of the aerosol layer in the current models was generally too low, which also contributed to the bias in modeled lifetime. Modeled MECs exhibited significant correlations with the Ångström exponent (AE, an indicator of particle size). Comparisons with the POLDER-GRASP-observed AE suggested that the models tended to overestimate the AE (underestimated particle size), indicating a possible underestimation of MECs in models. The hygroscopic growth in most models generally agreed with observations and might not explain the overall underestimation of modeled AOD. Our results imply that current global models contain biases in important aerosol processes for BBA (e.g., emissions, removal, and optical properties) that remain to be addressed in future research.

biomass burning aerosols↗