Search NASA⌕ Search

SEARCH · Search NASA

Results for “Radiative transfer simulations”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRS2012

Assessment of a Bidirectional Reflectance Distribution Correction of Above-Water and Satellite Water-Leaving Radiance in Coastal Waters

Water-leaving radiances, retrieved from in situ or satellite measurements, need to be corrected for the bidirectional properties of the measured light in order to standardize the data and make them comparable with each other. The current operational algorithm for the correction of bidirectional effects from the satellite ocean color data is optimized for typical oceanic waters. However, versions of bidirectional reflectance correction algorithms specifically tuned for typical coastal waters and other case 2 conditions are particularly needed to improve the overall quality of those data. In order to analyze the bidirectional reflectance distribution function (BRDF) of case 2 waters, a dataset of typical remote sensing reflectances was generated through radiative transfer simulations for a large range of viewing and illumination geometries. Based on this simulated dataset, a case 2 water focused remote sensing reflectance model is proposed to correct above-water and satellite water-leaving radiance data for bidirectional effects. The proposed model is first validated with a one year time series of in situ above-water measurements acquired by collocated multispectral and hyperspectral radiometers, which have different viewing geometries installed at the Long Island Sound Coastal Observatory (LISCO). Match-ups and intercomparisons performed on these concurrent measurements show that the proposed algorithm outperforms the algorithm currently in use at all wavelengths, with average improvement of 2.4% over the spectral range. LISCO's time series data have also been used to evaluate improvements in match-up comparisons of Moderate Resolution Imaging Spectroradiometer satellite data when the proposed BRDF correction is used in lieu of the current algorithm. It is shown that the discrepancies between coincident in-situ sea-based and satellite data decreased by 3.15% with the use of the proposed algorithm.

LISCO(LONG ISLAND SOUND COASTAL OBSERVATORY)↗
NASA NTRS1999

Cloud Thickness from Diffusion of Lidar Pulses in Clouds

Measurements of the distribution of reflected light from a laser beam incident on an aqueous suspension of particles or "cloud" with known thickness and particle size distribution are reported. The distribution is referred to as the "cloud radiative Green's function", G. In the diffusion domain, G is sensitive to cloud thickness, allowing that important quantity to be retrieved. The goal of the laboratory simulation is to provide preliminary estimates of sensitivity of G to cloud thickness,for use in the optimal design of an offbeam Lidar instrument for remote sensing of cloud thickness (THOR, Thickness from Offbeam Returns). These clouds of polystyrene microspheres suspended in water are analogous to real clouds of water droplets suspended in air. The microsphere size distribution is roughly lognormal, from 0.5 microns to 25 microns, similar to real clouds. Density of suspended spheres is adjusted so mean-free-path of visible photons is about 10 cm, approximately 1000 times smaller than in real clouds. The light source is a ND:YAG laser at 530 nm. Detectors are flux and photon-counting Photomultiplier Tube (PMTS), with a glass probe for precise positioning. A Labview 5 VI controls positioning, and data acquisition, via an NI Motion Control board connected to a stepper motor driving an Edmund linear slider, and a 16-channel 16-bit NI-DAQ board. The stepper motor is accurate to 10 microns, and step size is selectable from the VI software. Far from the incident beam, the rate of exponential increase as the direction of the incident beam is approached scales as expected from diffusion theory, linearly with the cloud thickness, and inversely as the square root of the reduced optical thickness, and is independent of particle size. Near the beam the signal begins to increase faster than exponential, due to single and low-order scattering near the backward direction, and here the distribution depends on particle size. Results are being used to verify 3D Monte Carlo radiative transfer simulations, used to estimate signal-to-noise ratios for remotely sensed off beam returns, for both homogeneous and inhomogeneous clouds. Signal-to-noise estimates show that unfiltered observations are straight forward at night, while narrow band pass filters are being studied for day.

Cahalan, Robert F.↗
NASA NTRS2021

Three-dimensional estimation of deciduous forest canopy structure and leaf area using multi-directional, leaf-on and leaf-off airborne lidar data

Airborne laser scanning (ALS) has been widely used to map gap probability and leaf area index (LAI) distribution at plot and landscape scales. As an indirect measurement, most ALS methods to estimate LAI combine waveform or point density information with supporting field measurements such as the leaf angle distribution, gap probability, or direct LAI measures. The development of a more independent estimation approach would facilitate more widespread use of existing ALS data to investigate patterns of forest structure and build realistic 3-D vegetation scenes to simulate remote sensing imagery and energy balance. Here, we develop a data processing workflow (named PVlad) using ALS point cloud apparent reflectance to estimate LAI and voxel-based leaf area density (LAD), aiming to reduce the need for associated field measurements such as the gap probability. The adaptation of the path volume (PV) concept derived from apparent reflectance integrates information from multi-directional ALS pulses, and quantifies the percentage exploration of each voxel for classification and occlusion correction, such that rigorous volumetric sampling approaches can be developed to derive LAI and LAD. The PVlad workflow was applied to discrete-return lidar data (Riegl VQ480i) acquired by NASA Goddard's LiDAR, Hyperspectral and Thermal Imager (G-LiHT) Airborne Imager during leaf-on (summer) and leaf-off (spring) conditions at the Smithsonian Environmental Research Center (SERC). The estimates of LAI and LAD captured structural differences between mature, logged, and intermediate-aged stands over eight deciduous forest plots. The derived LAI values were compared to field litter collection measurements, and the derived LAD vertical distribution was compared to the output of the VoxLAD model using terrestrial laser scan (TLS) field survey data. Using voxel sizes ranging from 0.5 m to 5 m, overall LAI estimation showed linear fitting coefficient bias and for 1 and 2 m voxel sizes, and vertical LAD distribution showed strong correlation with and for 0.5 and 1m voxel sizes. For every forest stand, upper-canopy LAD had a low variance for voxel sizes of ≤ . Application of PVlad to the G-LiHT and other similar ALS data archives enables the development of fine-resolution LAI map products, including voxelization of LAD for ecosystem science and radiative transfer simulations of remote sensing imagery or surface energy balance.

Tiangang Yin↗
NASA NTRS2016

Shortwave Direct Radiative Effects of Above-Cloud Aerosols Over Global Oceans Derived From 8 Years of CALIOP and MODIS Observations

In this paper, we studied the frequency of occurrence and shortwave direct radiative effects (DREs) of above-cloud aerosols (ACAs) over global oceans using 8 years (2007-2014) of collocated CALIOP and MODIS observations. Similar to previous work, we found high ACA occurrence in four regions: southeastern (SE) Atlantic region, where ACAs are mostly light-absorbing aerosols, i.e., smoke and polluted dust according to CALIOP classification, originating from biomass burning over the African Savanna; tropical northeastern (TNE) Atlantic and the Arabian Sea, where ACAs are predominantly windblown dust from the Sahara and Arabian deserts, respectively; and the northwestern (NW) Pacific, where ACAs are mostly transported smoke and polluted dusts from Asia. From radiative transfer simulations based on CALIOP-MODIS observations and a set of the preselected aerosol optical models, we found the DREs of ACAs at the top of atmosphere (TOA) to be positive (i.e., warming) in the SE Atlantic and NW Pacific regions, but negative (i.e., cooling) in the TNE Atlantic Ocean and the Arabian Sea. The cancellation of positive and negative regional DREs results in a global ocean annual mean diurnally averaged cloudy-sky DRE of 0.015 W m(exp. -2) [range of -0.03 to 0.06 W m (exp. -2)] at TOA. The DREs at surface and within the atmosphere are -0.015 W m(exp. -2) [range of -0.09 to -0.21 W m(exp. -2)], and 0.17 W m(exp. -2) [range of 0.11 to 0.24 W m(exp. -2)], respectively. The regional and seasonal mean DREs are much stronger. For example, in the SE Atlantic region, the JJA (July-August) seasonal mean cloudy-sky DRE is about 0.7 W m(exp. -2) [range of 0.2 to 1.2 W m(exp. -2)] at TOA. All our DRE computations are publicly available. The uncertainty in our DRE computations is mainly caused by the uncertainties in the aerosol optical properties, in particular aerosol absorption, the uncertainties in the CALIOP operational aerosol optical thickness retrieval, and the ignorance of cloud and potential aerosol diurnal cycle. In situ and remotely sensed measurements of ACA from future field campaigns and satellite missions and improved lidar retrieval algorithm, in particular vertical feature masking, would help reduce the uncertainty.

direct radiative effects↗
NASA NTRS1979

Radiometric resolution for monitoring vegetation: How many bits are needed?

The significance of the various number of radiometric quantizing levels required for satellite monitoring of vegetation resources was evaluated by using in situ collected spectral reflectance data, an atmospheric radiative transfer simulation model, and a satellite sensor simulation model. Reflectance data were converted to radiance data; passed through a model atmosphere to an altitude of 706 km; and subsequently quantized at 16, 32, 64, 128, 256, and 512 digital count levels for Thematic Mapper bands TM3(0.63 - 0.69 microns) and TM4(0.76 - 0.90 microns). The simulated digital count data were regressed against the in situ biological data to quantify the relationship between quantizing levels.

Tucker, C. J.↗
NASA NTRS1980

Radiometric resolution for monitoring vegetation - how many bits are needed

The number of radiometric quantizing levels required for satellite monitoring of vegetation resources was evaluated by using in situ collected spectral reflectance data, an atmospheric radiative transfer simulation model, and a satellite sensor simulation model. Reflectance data were converted to radiance data, passed through a model atmosphere to an altitude of 706 km, and subsequently quantized at 16, 32, 64, 128, 256 and 512 digital count levels for Thematic Mapper bands TM3 (0.63-0.69 microns) and TM4 (0.76-0.90 microns). The simulated digital count data were regressed against in situ biological data to quantify the relationship between quantizing levels. Results of the analysis demonstrated that solar zenith angle has an effect on the quantization equivalent change in reflectance, that 256 quantizing levels gave a 1-3% improvement per channel over 64 quantizing levels, and that 256 quantizing levels gave a 1% improvement per channel over 128 quantizing levels. No improvements were found for 256 versus 512 quantizing levels.

Tucker, C. J.↗
NASA NTRS1984

Sensitivity of airborne fluorosensor measurements to linear vertical gradients in chlorophyll concentration

A semianalytic Monte Carlo radiative transfer simulation model for airborne laser fluorosensors has been extended to investigate the effects of inhomogeneities in the vertical distribution of phytoplankton concentrations in clear seawater. Simulation results for linearly varying step concentrations of chlorophyll are presented. The results indicate that statistically significant differences can be seen under certain conditions in the water Raman-normalized fluorescence signals between nonhomogeneous and homogeneous cases. A statistical test has been used to establish ranges of surface concentrations and/or verticl gradients in which calibration by surface samples would by inappropriate, and the results are discussed.

Venable, D. D.↗
NASA NTRS1986

Monte Carlo study of the atmospheric spread function

Monte Carlo radiative transfer simulations are used to study the atmospheric spread function appropriate to satellite-based sensing of the earth's surface. The parameters which are explored include the nadir angle of view, the size distribution of the atmospheric aerosol, and the aerosol vertical profile.

Pearce, W. A.↗
NASA NTRS1991

Inferring chromospheric flare heating from hydrogen-line wings

Strong flare heating that penetrates deeply into the solar chromosphere is examined, and a simple model that is supported by a numerical radiative transfer simulation is applied to study how elevated free electron densities produce markedly increased Stark emission in the wings of hydrogen Lyman and Balmer lines. The proportionality that exists between the observed Stark wing intensity and the total energy deposition in the partially ionized layer is calculated for various Lyman and Balmer lines. It is concluded that when the effects of background opacity sources on the wing enhancements are not important, the magnitude of the nonthermal heating that penetrates the chromosphere may be inferred from the observed line wings. It is predicted that Ly-alpha extended wings will exhibit enhancement of Stark shape, even when inherent damping wing opacity exceeds Stark wing opacity over the entire absorption profile.

Gayley, K. G.↗
NASA NTRS1992

Effects of temperature-dependent molecular absorption coefficients on the thermal infrared remote sensing of the earth surface

The effect of temperature-dependent molecular absorption coefficients on thermal infrared spectral signatures measured from satellite sensors is investigated by comparing results from the atmospheric transmission and radiance codes LOWTRAN and MODTRAN and the accurate multiple scattering radiative transfer model ATRAD for different atmospheric profiles. The sensors considered include the operational NOAA AVHRR and two research instruments planned for NASA's Earth Observing System (EOS): MODIS-N (Moderate Resolution Imaging Spectrometer-Nadir-Mode) and ASTER (Advanced Spaceborne Thermal Emission and Reflection Radiometer). The difference in band transmittance is as large as 6 percent for some thermal bands within atmospheric windows and more than 30 percent near the edges of these atmospheric windows. The effect of temperature-dependent molecular absorption coefficients on satellite measurements of sea-surface temperature can exceed 0.6 K. Quantitative comparison and factor analysis indicate that more accurate measurements of molecular absorption coefficients and better radiative transfer simulation methods are needed to achieve SST accuracy of 0.3 K, as required for global numerical models of climate, and to develop land-surface temperature algorithms at the 1-K accuracy level.

Wan, Zhengming↗
NASA NTRS1994

Development of practical multiband algorithms for estimating land-surface temperature from EOS/MODIS data

A practical multiband, hierarchical algorithm for estimating land-surface temperature from NASA's future Earth Observing System (EOS) instruments Moderate Resolution Imaging Spectroradiometer (MODIS) and Advance Spaceborne Thermal Emission and Reflection Radiometer (ASTER) is developed through comprehensive, accurate, radiative transfer simulations at moderate spectral steps of 1-5/cm for wide ranges of atmospheric and surface conditions. The algorithm will accept empirical or estimated information about the surface emissivity and reflectivity and the atmospheric temperature and water-vapor profiles. Ground-based and aircraft measurements are necessary to validate and improve the algorithm and to establish its quality. Its accuracy depends on the calibration accuracy of thermal infrared data, uncertainties in surface heterogeneity, and temperature-dependent atmospheric absorption coefficients. Better knowledge of land-surface spectral emissivities and more accurate coefficients for atmospheric molecular band absorption and water vapor continuum absorption are needed to develop global land-surface temperature algorithms accurate to 1-2 K.

Dozier, J.↗
NASA NTRS2000

Aerosol Radiative Forcing in Asian Continental Outflow

Aerosols in elevated layers were sampled with FSSP-probes and wire impactors over the Pacific ocean aboard the NASA DC-8 aircraft. Analyses of particle size and morphology identifies two distinctly different aerosol types for cases when the mid-visible extinctions exceed 0.2/km. Smaller sizes (effective radii of 0.2 um) and moderate absorption (mid-visible single scattering albedo of.935) are typical for urban-industrial pollution. Larger sizes (effective radii of 0.7 um) and weak absorption (mid-visible single scattering albedo of 0.985) identify dust. This aerosol classification is in agreement with its origin as determined by airmass back trajectory analysis. Based on lidar vertical profiling, aerosol dominated by dust and urban-industrial pollution above 3km were assigned mid-visible optical depths of 0.50 and 0.27, respectively. Radiative transfer simulations, considering a 50% cloud-cover below the aerosol layers, suggest (on a daily tP C)C> basis) small reductions (-4W/m2) to the energy budget at the top of the atmosphere for both aerosol types. For c' 0 dust, more backscattering of sunlight (weaker solar absorption) is compensated by a stronger greenhouse effect due to larger sizes. Forced reductions to the energy budget at the surface are 12W/m2 for both aerosol types. In contrast, impacts on heating rates within the aerosol layers are quite different: While urban-industrial aerosol warms the layer (at +0.6K/day as solar heating dominates), dust cools (at -0.5K/day as infrared cooling dominates). Sensitivity tests show the dependence of the aerosol climatic impact on the optical depth, particle size, absorptivity, and altitude of the layers, as well as clouds and surface properties. Climatic cooling can be eliminated (1) for the urban-industrial aerosol if absorption is increased to yield a mid-visible single scattering albedo of 0.89, or if the ocean is replaced by a land surface; (2) for the dust aerosol if the effective radius is increased from 0.7 to 1.2 um. The removal of low-level clouds doubles the cooling at the top of the atmosphere to about -8W/m2.

Pueschel, R.↗
NASA NTRS2001

TRMM Precipitation Radar and Microwave Imager Observations of Convective and Stratiform Rain Over Land and Their Theoretical Implications

Observations of brightness temperature, Tb made over land regions by the Tropical Rainfall Measuring Mission (TRMM) Microwave Imager (TMI) radiometer have been analyzed along with the nearly simultaneous measurements of the vertical profiles of reflectivity factor, Z, made by the Precipitation Radar (PR) onboard the TRMM satellite. This analysis is performed to explore the interrelationship between the TMI and PR data in areas that are covered predominantly by convective or stratiform rain. In particular, we have compared on a scale of 20 km, average vertical profiles of Z with the averages of Tbs in the 19, 37 and 85 GHz channels. Generally, we find from these data that as Z increases, Tbs in the three channels decrease due to extinction. In order to explain physically the relationship between the Tb and Z observations, we have performed radiative transfer simulations utilizing vertical profiles of hydrometeors applicable to convective and stratiform rain regions. These profiles are constructed taking guidance from the Z observations of PR and recent LDR and ZDR measurements made by land-based polarimetric radars.

Prabhakara, C.↗
NASA NTRS2003

Regional Variability in Convection and Rain Retrievals from the TRMM Microwave Imager (TMI)

Precipitation Radar (PR) on board the TRMM satellite shows that the average height of 30 dBz in convective rain areas of the tropics varies significantly from one region to the other. When the convection is weak this height is shallow and when convection is strong this height extends deeper into the troposphere. The brightness temperature (Tb) measured by the microwave radiometer by itself does not reflect this nature of convection satisfactorily. Radiative transfer simulations of Tbs reveal that this could be due to the variations in the vertical distribution of optically active water and ice hydrometeors and their density, shape, and size. These variations are not coupled uniquely to the strength of the convective updrafts, and as a result the Tbs do not reflect properly the convective strength indicated by PR. Because of this deficiency in the Tbs the rain rate deduced from them differs from that of PR. For this reason, to improve the estimation of rain rate we have developed an empirical method. In this method a parameter based on the areal extent of the Tbs that exceed a certain magnitude is included along with the Tbs. Rain rate deduced with this approach is better correlated with that of PR when compared to the current Version 5 operational algorithm. The percentage of rain volume as a function of rain rate, for a given region of 5deg lat. X 5deg long. over a period of three months, deduced from this method, is also in better agreement with that of the PR.

Prabhakara, C.↗
NASA NTRS2006

Potential Nighttime Contamination of CERES Clear-sky Field of View by Optically Thin Cirrus during the CRYSTAL-FACE Campaign

We investigate the outgoing broadband longwave (LW, 5 to approx. 200 microns) and window (WIN, 8 to approx. 12 microns) channel radiances at the top of atmosphere (TOA) under clear-sky conditions, using data acquired by the Cloud and the Earth s Radiant Energy System (CERES) and Moderate-Resolution Imaging Spectroradiometer (MODIS) instruments onboard the NASA Terra satellite platform. In this study, detailed analyses are performed on the CERES Single Scanner Footprint TOA/Surface Fluxes and Clouds product to understand the radiative effect of thin cirrus. The data are acquired over the Florida area during the Cirrus Regional Study of Tropical Anvils and Cirrus Layers Florida Area Cirrus Experiment (CRYSTAL-FACE) field program. Of particular interest is the anisotropy associated with the radiation field. Measured CERES broadband radiances are compared to those obtained from rigorous radiative transfer simulations. Analysis of results from this comparison indicates that the simulated radiances tend to be larger than their measured counterparts, with differences ranging from 2.1% to 8.3% for the LW band and from 1.7% to 10.6% for the WIN band. The averaged difference in radiance is approximately 4% for both the LW and WIN channels. A potential cause for the differences could be the presence of thin cirrus (i.e., optically thin ice clouds with visible optical thicknesses smaller than approximately 0.3). The detection and quantitative analysis of these thin cirrus clouds are challenging even with sophisticated multispectral instruments. While large differences in radiance between the CERES observations and the theoretical calculations are found, the corresponding difference in the anisotropic factors is very small (0.2%). Furthermore, sensitivity studies show that the influence due to a 1 K bias of the surface temperature on the errors of the LW and WIN channel radiances is of the same order as that associated with a 2% bias of the surface emissivity. The LW and WIN errors associated with a 5% bias of water vapor amount in the lower atmosphere in conjunction with a 50% bias of water vapor amount in the upper atmosphere is similar to that of a 1 K bias of the vertical temperature profile. Even with the uncertainties considered for these various factors, the simulated LW and WIN radiances are still larger than the observed radiances if thin cirrus clouds are excluded.

Lee, Yong-Keun↗
NASA NTRS2014

Retrieval of Aerosol Optical Properties under Thin Cirrus from MODIS

Retrieval of aerosol optical properties using shortwave bands from passive satellite sensors, such as MODIS, is typically limited to cloud-free areas. However, if the clouds are thin enough (i.e. thin cirrus) such that the satellite-observed reflectance contains signals under the cirrus layer, and if the optical properties of this cirrus layer are known, the TOA reflectance can be corrected for the cirrus layer to be used for retrieving aerosol optical properties. To this end, we first correct the TOA reflectances in the aerosol bands (0.47, 0.55, 0.65, 0.86, 1.24, 1.63, and 2.12 micron for ocean algorithm and 0.412, 0.47, and 0.65 micron for deep blue algorithm) for the effects of thin cirrus using 1.38 micron reflectance and conversion factors that convert cirrus reflectance in 1.38 micron band to those in aerosol bands. It was found that the conversion factors can be calculated by using relationships between reflectances in 1.38 micron band and minimum reflectances in the aerosol bands (Gao et al., 2002). Refer to the example in the figure. Then, the cirrus-corrected reflectance can be calculated by subtracting the cirrus reflectance from the TOA reflectance in the optically thin case. A sensitivity study suggested that cloudy-sky TOA reflectances can be calculated with small errors in the form of simple linear addition of cirrus-only reflectances and clear-sky reflectances. In this study, we correct the cirrus signals up to TOA reflectance at 1.38 micron of 0.05 where the simple linear addition is valid without extensive radiative transfer simulations. When each scene passes the set of tests shown in the flowchart, the scene is corrected for cirrus contamination and passed into aerosol retrieval algorithms.

MODIS↗
NASA NTRS2014

Evaluating Radiometric Sensitivity of LandSat 8 Over Coastal-Inland Waters

The operational Land Imager (OLI) aboard Landsat 8 was launched in February 2013 to continue the Landsat's mission of monitoring earth resources at relatively high spatial resolution. Compared to Landsat heritage sensors, OLI has an additional 443-nm band (termed coastal/aerosol (CA) band), which extends its potential for mapping/monitoring water quality in coastal/inland waters. In addition, OLI's pushbroom design allows for longer integration time and, as a result, higher signal-to-noise ratio (SNR). Using a series of radiative transfer simulations, we provide insights into the radiometric sensitivity of OLI when studying coastal/inland waters. This will address how the changes in water constituents manifest at top-of-atmosphere (TOA) and whether the changes are resolvable at TOA (focal plane) relative to OLI's overall noise.

radiometric sensitivity↗
NASA NTRS2014

An Approach for the Long-Term 30-m Land Surface Snow-Free Albedo Retrieval from Historic Landsat Surface Reflectance and MODIS-based A Priori Anisotropy Knowledge

Land surface albedo has been recognized by the Global Terrestrial Observing System (GTOS) as an essential climate variable crucial for accurate modeling and monitoring of the Earth's radiative budget. While global climate studies can leverage albedo datasets from MODIS, VIIRS, and other coarse-resolution sensors, many applications in heterogeneous environments can benefit from higher-resolution albedo products derived from Landsat. We previously developed a "MODIS-concurrent" approach for the 30-meter albedo estimation which relied on combining post-2000 Landsat data with MODIS Bidirectional Reflectance Distribution Function (BRDF) information. Here we present a "pre-MODIS era" approach to extend 30-m surface albedo generation in time back to the 1980s, through an a priori anisotropy Look-Up Table (LUT) built up from the high quality MCD43A BRDF estimates over representative homogenous regions. Each entry in the LUT reflects a unique combination of land cover, seasonality, terrain information, disturbance age and type, and Landsat optical spectral bands. An initial conceptual LUT was created for the Pacific Northwest (PNW) of the United States and provides BRDF shapes estimated from MODIS observations for undisturbed and disturbed surface types (including recovery trajectories of burned areas and non-fire disturbances). By accepting the assumption of a generally invariant BRDF shape for similar land surface structures as a priori information, spectral white-sky and black-sky albedos are derived through albedo-to-nadir reflectance ratios as a bridge between the Landsat and MODIS scale. A further narrow-to-broadband conversion based on radiative transfer simulations is adopted to produce broadband albedos at visible, near infrared, and shortwave regimes.We evaluate the accuracy of resultant Landsat albedo using available field measurements at forested AmeriFlux stations in the PNW region, and examine the consistency of the surface albedo generated by this approach respectively with that from the "concurrent" approach and the coincident MODIS operational surface albedo products. Using the tower measurements as reference, the derived Landsat 30-m snow-free shortwave broadband albedo yields an absolute accuracy of 0.02 with a root mean square error less than 0.016 and a bias of no more than 0.007. A further cross-comparison over individual scenes shows that the retrieved white sky shortwave albedo from the "pre-MODIS era" LUT approach is highly consistent (R(exp 2) = 0.988, the scene-averaged low RMSE = 0.009 and bias = −0.005) with that generated by the earlier "concurrent" approach. The Landsat albedo also exhibits more detailed landscape texture and a wider dynamic range of albedo values than the coincident 500-m MODIS operational products (MCD43A3), especially in the heterogeneous regions. Collectively, the "pre-MODIS" LUT and "concurrent" approaches provide a practical way to retrieve long-term Landsat albedo from the historic Landsat archives as far back as the 1980s, as well as the current Landsat-8 mission, and thus support investigations into the evolution of the albedo of terrestrial biomes at fine resolution.

Albedo algorithm↗