Search NASA⌕ Search

SEARCH · Search NASA

Results for “Spatially distributed model evaluation”

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 163 records · Page 9

Comparison of simulated cloud cover with satellite obsrvations over the Western United States

Satellite imagery datasets and regional climate model results are intercompared for evaluation of model accuracy in the simulation of cloud cover. Both monthly average individual simulation times are analyzed. To provide a consistent comparison, satellite data are first mapped into the model's geographic projection, grid domain, and resolution. It is found that September 1988 monthly average cloud fraction results from the modeled simulations correspond to observations, in both spatial pattern and magnitude, with bias less than +/- 20% cloud fraction over the entire inland West. Agreement in the pattern of cloud fraction also is evident for monthly average cloud fraction in July, but there is no negative bias of 10%-30% cloud fraction in the model diagnosis of cloud cover. Correlations between the spatial distributions of model-derived and observed cloud fractions are found to exceed 0.80 for certain geographic regions of the West, and these correlations are largest over mountainous areas during summer. Case studies of a series of daily cloud cover demonstrate the ability of the model to simulate the effects of frontal passage on cloud distribution. The ability of the RegCM1 to simulate daily cloud fraction and diurnal cloud evolution is somewhat weak for the summer convective season. It is anticipated that a more recent version of the regional climate model may improve the simulation of summer season cloud cover, through changes in cloud parameterization and improvements in model resolution.

Wetzel, Melanie A.↗

Investigating Sources of Variability and Error in Simulations of Carbon Dioxide in an Urban Region

Greenhouse gas (GHG) emissions estimation methods that use atmospheric trace gas observations, including inverse modeling techniques, perform better when carbon dioxide (CO2) fluxes are more accurately transported and dispersed in the atmosphere by a numerical model. In urban areas, transport and dispersion is particularly difficult to simulate using current mesoscale meteorological models due, in part, to added complexity from surface heterogeneity and fine spatial/temporal scales. It is generally assumed that the errors in GHG estimation methods in urban areas are dominated by errors in transport and dispersion. Other significant errors include, but are not limited to, those from assumed emissions magnitude and spatial distribution. To assess the predictability of simulated trace gas mole fractions in urban observing systems using a numerical weather prediction model, we employ an Eulerian model that combines traditional meteorological variables with multiple passive tracers of atmospheric CO2 from anthropogenic inventories and a biospheric model. The predictability of the Eulerian model is assessed by comparing simulated atmospheric CO2 mole fractions to observations from four in situ tower sites (three urban and one rural) in the Washington DC/Baltimore, MD area for February 2016. Four different gridded fossil fuel emissions inventories along with a biospheric flux model are used to create an ensemble of simulated atmospheric CO2 observations within the model. These ensembles help to evaluate whether the modeled observations are impacted more by the underlying emissions or transport. The spread of modeled observations using the four emission fields indicates the model's ability to distinguish between the different inventories under various meteorological conditions. Overall, the Eulerian model performs well; simulated and observed average CO2 mole fractions agree within 1% when averaged at the three urban sites across the month. However, there can be differences greater than 10% at any given hour, which are attributed to complex meteorological conditions rather than differences in the inventories themselves. On average, the mean absolute error of the simulated compared to actual observations is generally twice as large as the standard deviation of the modeled mole fractions across the four emission inventories. This result supports the assumption, in urban domains, that the predicted mole fraction error relative to observations is dominated by errors in model meteorology rather than errors in the underlying fluxes in winter months. As such, minimizing errors associated with atmospheric transport and dispersion may help improve the performance of GHG estimation models more so than improving flux priors in the winter months. We also find that the errors associated with atmospheric transport in urban domains are not restricted to certain times of day. This suggests that atmospheric inversions should use CO2 observations that have been filtered using meteorological observations rather than assuming that meteorological modeling is most accurate at certain times of day (such as using only mid-afternoon observations).

CO2↗

A Heuristic Approach to Global Landslide Susceptibility Mapping

Landslides can have significant and pervasive impacts to life and property around the world. Several attempts have been made to predict the geographic distribution of landslide activity at continental and global scales. These efforts shared common traits such as resolution, modeling approach, and explanatory variables. The lessons learned from prior research have been applied to build a new global susceptibility map from existing and previously unavailable data. Data on slope, faults, geology, forest loss, and road networks were combined using a heuristic fuzzy approach. The map was evaluated with a Global Landslide Catalog developed at the National Aeronautics and Space Administration, as well as several local landslide inventories. Comparisons to similar susceptibility maps suggest that the subjective methods commonly used at this scale are, for the most part, reproducible. However, comparisons of landslide susceptibility across spatial scales must take into account the susceptibility of the local subset relative to the larger study area. The new global landslide susceptibility map is intended for use in disaster planning, situational awareness, and for incorporation into global decision support systems.

mapping↗

Role of Updraft Velocity in Temporal Variability of Global Cloud Hydrometeor Number

Understanding how dynamical and aerosol inputs affect the temporal variability of hydrometeor formation in climate models will help to explain sources of model diversity in cloud forcing, to provide robust comparisons with data, and, ultimately, to reduce the uncertainty in estimates of the aerosol indirect effect. This variability attribution can be done at various spatial and temporal resolutions with metrics derived from online adjoint sensitivities of droplet and crystal number to relevant inputs. Such metrics are defined and calculated from simulations using the NASA Goddard Earth Observing System Model, Version 5 (GEOS-5) and the National Center for Atmospheric Research Community Atmosphere Model Version 5.1 (CAM5.1). Input updraft velocity fluctuations can explain as much as 48% of temporal variability in output ice crystal number and 61% in droplet number in GEOS-5 and up to 89% of temporal variability in output ice crystal number in CAM5.1. In both models, this vertical velocity attribution depends strongly on altitude. Despite its importance for hydrometeor formation, simulated vertical velocity distributions are rarely evaluated against observations due to the sparsity of relevant data. Coordinated effort by the atmospheric community to develop more consistent, observationally based updraft treatments will help to close this knowledge gap.

cloud condensation nuclei↗

International Conference on Mesoscale Precipitation: Hydrologic and Meteorological Aspects, 3rd, College Station, TX, Feb. 27-Mar. 1, 1991, Proceedings

Papers presented in this volume discuss a mechanism causing mesoscale organizations of precipitation in midlatitude cyclones, the spatial variability of summer Florida precipitation and its impact on microwave radiometer rainfall-measurement systems, the relation between the mean areal rainfall and the fractional area where it rains above a given threshold, and a possible explanation for low correlation dimension estimates for the atmosphere. Attention is also given to a physically based simulation of radar rainfall data using a space-time rainfall model, the predictability of mesoscale rainfall in the tropics, and an intercomparison of tree satellite infrared rainfall techniques over Japan and surrounding waters. Other papers are on an evaluation of sampling errors of precipitation from spaceborne and ground sensors, an approach to the estimation of the areal rain-rate distribution from spaceborne radar by the use of multiple thresholds, and a comparison of two satellite-based rainfall algorithms using Pacific atoll raingage data.

Valdes, Juan B.↗

Seismology and space-based geodesy

The potential of space-based geodetic measurement of crustal deformation in the context of seismology is explored. The achievements of seismological source theory and data analyses, mechanical modeling of fault zone behavior, and advances in space-based geodesy are reviewed, with emphasis on realizable contributions of space-based geodetic measurements specifically to seismology. The fundamental relationships between crustal deformation associated with an earthquake and the geodetically observable data are summarized. The response and spatial and temporal resolution of the geodetic data necessary to understand deformation at various phases of the earthquake cycle is stressed. The use of VLBI, SLR, and GPS measurements for studying global geodynamics properties that can be investigated to some extent with seismic data is discussed. The potential contributions of continuously operating strain monitoring networks and globally distributed geodetic observatories to existing worldwide modern digital seismographic networks are evaluated in reference to mutually addressable problems in seismology, geophysics, and tectonics.

Tralli, David M.↗

GERMcode: A Stochastic Model for Space Radiation Risk Assessment

A new computer model, the GCR Event-based Risk Model code (GERMcode), was developed to describe biophysical events from high-energy protons and high charge and energy (HZE) particles that have been studied at the NASA Space Radiation Laboratory (NSRL) for the purpose of simulating space radiation biological effects. In the GERMcode, the biophysical description of the passage of HZE particles in tissue and shielding materials is made with a stochastic approach that includes both particle track structure and nuclear interactions. The GERMcode accounts for the major nuclear interaction processes of importance for describing heavy ion beams, including nuclear fragmentation, elastic scattering, and knockout-cascade processes by using the quantum multiple scattering fragmentation (QMSFRG) model. The QMSFRG model has been shown to be in excellent agreement with available experimental data for nuclear fragmentation cross sections. For NSRL applications, the GERMcode evaluates a set of biophysical properties, such as the Poisson distribution of particles or delta-ray hits for a given cellular area and particle dose, the radial dose on tissue, and the frequency distribution of energy deposition in a DNA volume. By utilizing the ProE/Fishbowl ray-tracing analysis, the GERMcode will be used as a bi-directional radiation transport model for future spacecraft shielding analysis in support of Mars mission risk assessments. Recent radiobiological experiments suggest the need for new approaches to risk assessment that include time-dependent biological events due to the signaling times for activation and relaxation of biological processes in cells and tissue. Thus, the tracking of the temporal and spatial distribution of events in tissue is a major goal of the GERMcode in support of the simulation of biological processes important in GCR risk assessments. In order to validate our approach, basic radiobiological responses such as cell survival curves, mutation, chromosomal aberrations, and representative mouse tumor induction curves are implemented into the GERMcode. Extension of these descriptions to other endpoints related to non-targeted effects and biochemical pathway responses will be discussed.

Kim, Myung-Hee Y.↗

Effects of Fuel Preheat on Soot Formation in Microgravity Laminar Diffusion Flames

Nonbuoyant flames offer themselves as an attractive and promising platform to gain a better understanding of soot mechanisms. The effects of buoyancy can be eliminated temporarily in drop towers which sustain brief intervals of reduced gravity-typically lower than 10(exp -3)g- extending up to several seconds at a time. Microgravity facilities have been employed to show that nonbuoyant flames are longer, wider and sootier than their normal-gravity counterparts. Sunderland et al. recently verified the existence of smoke point in laminar nonbuoyant flames. As reported, microgravity flames operating above their smoke point displayed a blunt tip and much broader soot-containing regimes in comparison to their buoyant counterparts. Mortazavi et al. established that residence times in microgravity laminar jet diffusion flames with Re=0(100) tend to be proportional to burner diameter and inversely proportional to burner exit velocity. This offers the capability to alter residence times in nonbuoyant laminar jet diffusion flames when varying the burner exit diameters and velocities. Megaridis et al. presented a quantitative definition of the soot-field structure within laminar microgravity jet diffusion flames which operated well above their smoke point. The experimental methodology involved a full-field laser-light extinction technique and jet diffusion flames of nitrogen-diluted (50% vol.) acetylene fuel burning in quiescent air at atmospheric pressure. The work was conducted at the 2.2s drop tower of the NASA Lewis Research Center (NASA-LeRC). Parallel work on 1-g flames was also presented in (6) to facilitate comparisons on the effect of gravity on the soot fields. As reported, the soot spatial distributions in 0-g flames did not change in a detectable manner after 1s within a typical 2.2s experiment. During that period, the soot field was shown to sustain a pronounced annular structure throughout the luminous nonbuoyant-flame zone. The maximum soot volume fraction measured at 0-g was nearly a factor of two higher than that at 1-g, thus confirming the enhanced sooting tendency of nonbuoyant flames. Greenberg and Ku presented a similar study and reported trends that matched those of for the 50% (vol.) nitrogen-diluted acetylene fuel. Furthermore, they examined pure acetylene flames and reported similar trends with respect to the influence of gravity on maximum soot volume fractions and flame cross-section-averaged soot loadings. Both studies clearly demonstrated the improved spatial resolution of microgravity flames compared to their normal-gravity counterparts. The current study evaluates the influence of moderate fuel preheat on soot formation within 0-g laminar gas jet diffusion flames. While fuel temperature variations have little influence on residence times in 1-g, they have a much more significant effect in 0-g. The primary objective of this program is to quantify this effect and its consequences on sooting by comparing soot volume fraction distributions under preheated and unpreheated-fuel conditions. Furthermore, the current work aims at expanding the limited soot database available for nonbuoyant flames. Soot fields in such flames can be used to perform additional tests of recently developed soot sub-models which have the potential to become powerful predicting tools in combustion design.

Konsur, Bogdan↗

Changes in the Mass Balance of the Greenland Ice Sheet in a Warming Climate During 2003-2009

Mass changes of the Greenland ice sheet (GIS) derived from ICESat and GRACE data both show that the net mass loss from GIS during 2003-2009 is about 175 Gt/year, which contributes 0.5mm/yr global sea-level rise. The rate of mass loss has increased significantly since the 1990's when the GIS was close to mass balance. Even though the GIS was close to mass balance during the 1990's, it was already showing characteristics of responding to8 warmer climate, specifically thinning at the margins and thickening inland at higher elevations. During 2003-2009, increased ice thinning due to increases in melting and acceleration of outlet glaciers began to strongly exceed the inland thickening from increases in accumulation. Over the entire GIS, the mass loss between the two periods, from increased melting and ice dynamics, increased by about 190 Gt/year while the mass gain, from increased precipitation and accumulation, increased by only about 15Gt/year. These ice changes occurred during a time when the temperature on GIS changed at rate of about 2K/decade. The distribution of elevation and mass changes derived from ICESat have high spatial resolution showing details over outlet glaciers, by drainage systems, and by elevation. However, information on the seasonal cycle of changes from ICESat data is limited, because the ICESat lasers were only operated during two to three campaigns per year of about 35 days duration each. In contrast, the temporal resolution of GRACE data, provided by the continuous data collection, is much better showing details of the seasonal cycle and the inter-annual variability. The differing sensitivity of the ICESat altimetry and the GRACE gravity methods to motion of the underlying bedrock from glacial isostatic adjustment (GIA) is used to evaluate the GIA corrections provided by models. The two data types are also combined to make estimates of the partitioning of the mass gains and losses among accumulation, melting, and ice discharge from outlet glaciers.

Zwally, H. Jay↗

New Approach to Monitor Transboundary Particulate Pollution over Northeast Asia

A new approach to more accurately monitor and evaluate transboundary particulate matter (PM) pollution is introduced based on aerosol optical products from Korea's Geostationary Ocean Color Imager (GOCI). The area studied is Northeast Asia (including eastern parts of China, the Korean peninsula and Japan), where GOCI has been monitoring since June 2010. The hourly multi-spectral aerosol optical data that were retrieved from GOCI sensor onboard geostationary satellite COMS (Communication, Ocean, and Meteorology Satellite) through the Yonsei aerosol retrieval algorithm were first presented and used in this study. The GOCI-retrieved aerosol optical data are integrated with estimated aerosol distributions from US EPA Models-3/CMAQ (Community Multi-scale Air Quality) v4.5.1 model simulations via data assimilation technique, thereby making the aerosol data spatially continuous and available even for cloud contamination cells. The assimilated aerosol optical data are utilized to provide quantitative estimates of transboundary PM pollution from China to the Korean peninsula and Japan. For the period of 1 April to 31 May, 2011 this analysis yields estimates that AOD as a proxy for PM2.5 or PM10 during long-range transport events increased by 117-265% compared to background average AOD (aerosol optical depth) at the four AERONET sites in Korea, and average AOD increases of 121% were found when averaged over the entire Korean peninsula. This paper demonstrates that the use of multi-spectral AOD retrievals from geostationary satellites can improve estimates of transboundary PM pollution. Such data will become more widely available later this decade when new sensors such as the GEMS (Geostationary Environment Monitoring Spectrometer) and GOCI-2 are scheduled to be launched.

Data Assimilation↗

NASA Team 2 Sea Ice Concentration Algorithm Retrieval Uncertainty

Satellite microwave radiometers are widely used to estimate sea ice cover properties (concentration, extent, and area) through the use of sea ice concentration (IC) algorithms. Rare are the algorithms providing associated IC uncertainty estimates. Algorithm uncertainty estimates are needed to assess accurately global and regional trends in IC (and thus extent and area), and to improve sea ice predictions on seasonal to interannual timescales using data assimilation approaches. This paper presents a method to provide relative IC uncertainty estimates using the enhanced NASA Team (NT2) IC algorithm. The proposed approach takes advantage of the NT2 calculations and solely relies on the brightness temperatures (TBs) used as input. NT2 IC and its associated relative uncertainty are obtained for both the Northern and Southern Hemispheres using the Advanced Microwave Scanning Radiometer for the Earth Observing System (AMSR-E) TB. NT2 IC relative uncertainties estimated on a footprint-by-footprint swath-by-swath basis were averaged daily over each 12.5-km grid cell of the polar stereographic grid. For both hemispheres and throughout the year, the NT2 relative uncertainty is less than 5%. In the Southern Hemisphere, it is low in the interior ice pack, and it increases in the marginal ice zone up to 5%. In the Northern Hemisphere, areas with high uncertainties are also found in the high IC area of the Central Arctic. Retrieval uncertainties are greater in areas corresponding to NT2 ice types associated with deep snow and new ice. Seasonal variations in uncertainty show larger values in summer as a result of melt conditions and greater atmospheric contributions. Our analysis also includes an evaluation of the NT2 algorithm sensitivity to AMSR-E sensor noise. There is a 60% probability that the IC does not change (to within the computed retrieval precision of 1%) due to sensor noise, and the cumulated probability shows that there is a 90% chance that the IC varies by less than +/-3%. We also examined the daily IC variability, which is dominated by sea ice drift and ice formation/melt. Daily IC variability is the highest, year round, in the MIZ (often up to 20%, locally 30%). The temporal and spatial distributions of the retrieval uncertainties and the daily IC variability is expected to be useful for algorithm intercomparisons, climate trend assessments, and possibly IC assimilation in models.

Sea Ice↗

The impact of Doppler lidar wind observations on a single-level meteorological analysis

Through the use of observation operators, modern data assimilation systems have the capability to ingest observations of quantities that are not themselves model variables, but are mathematically related to those variables. An example of this are the so-called LOS (line of sight) winds that a Doppler wind Lidar can provide. The model - or data assimilation system - needs information about both components of the horizontal wind vectors, whereas the observations in this case only provide the projection of the wind vector onto a given direction. The analyzed value is then calculated essentially based on a comparison between the observation itself and the model-simulated value of the observed quantity. However, in order to assess the expected impact of such an observing system, it is important to examine the extent to which a meteorological analysis can be constrained by the LOS winds. The answer to this question depends on the fundamental character of the atmospheric flow fields that are analyzed, but more importantly it also depends on the real and assumed error covariance characteristics of these fields. A single-level wind analysis system designed to explore these issues has been built at the NASA Data Assimilation Office. In this system, simulated wind observations can be evaluated in terms of their impact on the analysis quality under various assumptions about their spatial distribution and error characteristics and about the error covariance of the background fields. The basic design of the system will be presented along with experimental results obtained with it. In particular, the value of simultaneously measuring LOS winds along two different directions for a given location will be discussed.

Riishojgaard, L. P.↗

Multi-Parameter Observation and Detection of Pre-Earthquake Signals in Seismically Active Areas

The recent large earthquakes (M9.0 Tohoku, 03/2011; M7.0 Haiti, 01/2010; M6.7 L Aquila, 04/2008; and M7.9 Wenchuan 05/2008) have renewed interest in pre-anomalous seismic signals associated with them. Recent workshops (DEMETER 2006, 2011 and VESTO 2009 ) have shown that there were precursory atmospheric /ionospheric signals observed in space prior to these events. Our initial results indicate that no single pre-earthquake observation (seismic, magnetic field, electric field, thermal infrared [TIR], or GPS/TEC) can provide a consistent and successful global scale early warning. This is most likely due to complexity and chaotic nature of earthquakes and the limitation in existing ground (temporal/spatial) and global satellite observations. In this study we analyze preseismic temporal and spatial variations (gas/radon counting rate, atmospheric temperature and humidity change, long-wave radiation transitions and ionospheric electron density/plasma variations) which we propose occur before the onset of major earthquakes:. We propose an Integrated Space -- Terrestrial Framework (ISTF), as a different approach for revealing pre-earthquake phenomena in seismically active areas. ISTF is a sensor web of a coordinated observation infrastructure employing multiple sensors that are distributed on one or more platforms; data from satellite sensors (Terra, Aqua, POES, DEMETER and others) and ground observations, e.g., Global Positioning System, Total Electron Content (GPS/TEC). As a theoretical guide we use the Lithosphere-Atmosphere-Ionosphere Coupling (LAIC) model to explain the generation of multiple earthquake precursors. Using our methodology, we evaluated retrospectively the signals preceding the most devastated earthquakes during 2005-2011. We observed a correlation between both atmospheric and ionospheric anomalies preceding most of these earthquakes. The second phase of our validation include systematic retrospective analysis for more than 100 major earthquakes (M>5.9) in Taiwan and Japan. We have found anomalous behavior before all of these events with no false negatives. Calculated false alarm ratio for the for the same month over the entire period of analysis (2003-2009) is less than 10% and was d as the earthquakes. The commonalities in detecting atmospheric/ionospheric anomalies show that they may exist over both land and sea in regions of maximum stress (i.e., along plate boundaries) Our results indicate that the ISTF model could provide a capability to observe pre-earthquake atmospheric/ionospheric signals by combining this information into a common framework.

Ouzounov, D.↗

How Well Do Global Climate Models Simulate the Variability of Atlantic Tropical Cyclones Associated with ENSO?

The variability of Atlantic tropical cyclones (TCs) associated with El Nino-Southern Oscillation (ENSO) in model simulations is assessed and compared with observations. The model experiments are 28-yr simulations forced with the observed sea surface temperature from 1982 to 2009. The simulations were coordinated by the U.S. CLIVAR Hurricane Working Group and conducted with five global climate models (GCMs) with a total of 16 ensemble members. The model performance is evaluated based on both individual model ensemble means and multi-model ensemble mean. The latter has the highest anomaly correlation (0.86) for the interannual variability of TCs. Previous observational studies show a strong association between ENSO and Atlantic TC activity, as well as distinctions in the TC activities during eastern Pacific (EP) and central Pacific (CP) El Nino events. The analysis of track density and TC origin indicates that each model has different mean biases. Overall, the GCMs simulate the variability of Atlantic TCs well with weaker activity during EP El Nino and stronger activity during La Nina. For CP El Nino, there is a slight increase in the number of TCs as compared with EP El Nino. However, the spatial distribution of track density and TC origin is less consistent among the models. Particularly, there is no indication of increasing TC activity over the U.S. southeast coastal region as in observations. The difference between the models and observations is likely due to the bias of vertical wind shear in response to the shift of tropical heating associated with CP El Nino, as well as the model bias in the mean circulation.

Tropical Cyclones↗

Titan Polar Landscape Evolution

With the ongoing Cassini-era observations and studies of Titan it is clear that the intensity and distribution of surface processes (particularly fluvial erosion by methane and Aeolian transport) has changed through time. Currently however, alternate hypotheses substantially differ among specific scenarios with respect to the effects of atmospheric evolution, seasonal changes, and endogenic processes. We have studied the evolution of Titan's polar region through a combination of analysis of imaging, elevation data, and geomorphic mapping, spatially explicit simulations of landform evolution, and quantitative comparison of the simulated landscapes with corresponding Titan morphology. We have quantitatively evaluated alternate scenarios for the landform evolution of Titan's polar terrain. The investigations have been guided by recent geomorphic mapping and topographic characterization of the polar regions that are used to frame hypotheses of process interactions, which have been evaluated using simulation modeling. Topographic information about Titan's polar region is be based on SAR-Topography and altimetry archived on PDS, SAR-based stereo radar-grammetry, radar-sounding lake depth measurements, and superposition relationships between geomorphologic map units, which we will use to create a generalized topographic map.

Cassini↗

Boreal canopy surfaces from spaceborne stereogrammetry

Surface elevation estimates from high resolution spaceborne image (HRSI) stereogrammetry are used to examine fine-scaled structure of boreal forest canopies. These data can depict detailed spatial patterns of vertical forest structure at remote sites across the circumpolar domain where these estimates would otherwise be unavailable. This work examines where these estimates are most effective at describing vertical forest structure to explain which canopy surfaces they represent. We evaluated the variation in canopy surface estimates captured from four general types of HRSI digital surface models (DSMs) across the full range of boreal canopy cover. These DSMs, classified into 4 types by grouping them according to the acquisition's (1) sun elevation angle (low or high) and (2) seasonality-driven ground surface condition (snow presence/absence), vary with acquisition characteristics and the details of this variation continues to be studied. We explored some of this variation by comparing the distributions of differences in boreal canopy percentile heights derived from reference small footprint lidar in Tanana Valley, Alaska with canopy surface elevations derived from these 4 types of HRSI DSMs. We examined how canopy surface estimates from HRSI DSMs differ according to acquisition characteristics and canopy cover, and ultimately which canopy surfaces are represented in these DSMs. Our results help clarify which boreal canopy surfaces are representative of those captured with HRSI DSMs. They show that in the Tanana Valley (1) DSMs grouped by sun elevation angle and ground surface condition provide different surface estimates of boreal canopies; (2) the two DSM types that appear to most differently capture boreal forest canopy surfaces are DSMs from snow-free images acquired at sun elevation angles <30° (Low sun elev. & snow-free) and those with snow-cover at sun elevation angles ≥30° (High sun elev. & snow-free); (3) DSMs with snow most often do not capture upper canopy surfaces; (4) the “Low sun elev. & snow-free” DSMs resolve surfaces that are most representative of upper canopy surfaces (dense forests >60% cover, 70th–80th percentile heights); and (5) in the most dense forests (>80% cover) where canopy gaps are least likely to bias downward the average surface estimates, the snow-free DSM types are representative of 70th - 80th percentile heights (“Low sun elev. & snow-free”) and 60th–70th percentile heights (“High sun elev. & snow-free”). The combination of horizontal structure (canopy cover) and acquisition characteristics affect the boreal vertical structure (canopy surface height) estimates from spaceborne stereogrammetry. These effects should be considered when analyzing products derived from HRSI DSMs, and as part of a comprehensive approach to spaceborne remote sensing of circumpolar boreal forests.

forest structure↗

Next-Generation Angular Distribution Models for Top-of-Atmosphere Radiative Flux Calculation from CERES Instruments: Validation

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.2Wm(exp -2) and the RMSE is less than 1.1Wm(exp -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.5Wm(exp -2) and the RMSE is less than 0.8Wm(exp -)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.9Wm(exp -2) over clear ocean, 1.6% (4.5Wm(exp -2) over clear land, and 2.0% (6.0Wm(exp -) over clear snow/ice; and are about 3.3% (9.0Wm(exp -2), 2.7% (8.4Wm(exp -2), and 3.7% (9.9Wm(exp -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 CERESMODIS test are 0.5% (1.5Wm(exp -2), 0.8% (2.4Wm(exp -2), and 0.7% (1.3Wm(exp -2) over clear ocean, land, and snow/ice; and are about 1.5% (3.5Wm(exp -2), 1.0% (2.9Wm(exp -2), and 1.1% (2.1Wm(exp -2) over ocean, land, and snow/ice under all-sky conditions. The TOA instantaneous nighttime LW flux uncertainties are about 0.5-1% (<2.0Wm(exp -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.6Wm(exp -2) and overestimate TOA daytime (nighttime) LW flux by 0.4 (0.2)Wm(exp -2) when all CERES viewing angles are considered.

Su, W.↗

Toward Landsat and Sentinel-2 BRDF Normalization and Albedo Estimation: A Case Study in the Peruvian Amazon Forest

The Amazon forest has been the focus of study by the science community during the last few decades. Remote sensing data analysis is the only way to study such a large geographical extent during an extended period of time. Since the launch in 2015 of Sentinel 2 and its increase in temporal resolution through the combination with Landsat sensors, a strong emphasis has been put on exploiting these data. Though these satellites provide near nadir observations, surface reflectance time series are affected by illumination variability throughout the year. These effects can be corrected using a Bidirectional Reflectance Distribution Function (BRDF) model. Franch et al. (2014a) developed a methodology to derive Landsat surface albedo and BRDF. It is based on the BRDF parameters from the MODerate Resolution Imaging Spectroradiometer (MODIS) which are disaggregated at Landsat spatial resolution (30 m). In this work, we apply the Franch et al. (2014a) method to normalize the surface reflectance for BRDF effects using the NASA's Harmonized Landsat Sentinel-2 (HLS) product. We apply this method to the Tambopata region in Peru from 2013 to 2017 and validate it using ground-based albedometer measurements. The results show that the near infrared reflectance can increase up to 0.06 (20%) for low solar angles while the impact on the red range and the NDVI is minor (<0.01). The evaluation of the surface albedo against field measurements shows an error of 0.01.

Landsat 8↗