Search

Most land-surface models require parameterization of vertical wind profiles within the atmospheric boundary layer. For vegetated surfaces, it is common to assume a logarithmic profile in the surface layer, which includes estimates of vegetation roughness length for momentum (z0) and zero-plane displacement height (d0). This study finds that remotely-sensed forest canopy heights improve estimates of aerodynamic roughness length for momentum using a previously-developed representation of the roughness sublayer (Raupach 1992; Jasinski et al. 2005). Resulting roughness products consist of two datasets: 1) 14 years of 8-day snapshots of the global land surface at a nominal spatial resolution of 500-meters for users who wish to retain full temporal resolution and interannual variability; and 2) multiyear averages of the 8-day snapshots, here referred to as "climatologies" of roughness, which retain underlying seasonality. Both products are suitable for use in data assimilation and reanalyses such as the National Climate Assessment Land Data Assimilation System (NCA-LDAS), for which these products were initially developed.

NASA;s Ice, Cloud, and Land Elevation Satellite-2 (ICESat-2) mission, scheduled to launch no later than April 2018 (and currently slated for October 2017), is being developed to continue the multiyear observations of the earths surface elevation, ice, and clouds started by ICESat. To increase the use of the satellite data after launch, the ICESat-2 mission invested in an applications program aimed at innovatively applying the data in a variety of fields. The program provides a framework for building a broad and well-defined user community during the prelaunch period to maximize the use of data products after launch and to provide early insight into the range of potential uses of the mission data. Ideas and research on how altimetry data will be used for decision making arise from the end users; there-fore, the ICESat-2 mission is extending itself through its applications program.

The U.S. Global Change Research Program Sustained Assessment process benefits from long-term investments in Earth science research that enable the scientific community to conduct assessment-relevant science. To this end, NASA initiated several research programs over the past five years to support the Earth observation community in developing indicators, datasets, research products, and tools to support ongoing and future National Climate Assessments. These activities complement NASA's ongoing Earth science research programs. One aspect of the assessment portfolio funds four "enabling tools" projects at NASA research centers. Each tool leverages existing capacity within the center, but has developed tailored applications and products for National Climate Assessments. The four projects build on the capabilities of a global atmospheric reanalysis (MERRA-2), a continental U.S. land surface reanalysis (NCA-LDAS), the NASA Earth Exchange (NEX), and a Regional Climate Model Evaluation System (RCMES). Here, we provide a brief overview of each enabling tool, highlighting the ways in which it has advanced assessment science to date. We also discuss how the assessment community can access and utilize these tools for National Climate Assessments and other sustained assessment activities.

The Advanced Topographic Laser Altimeter System (ATLAS) on the Ice, Cloud, and Land Elevation Satellite (ICESat-2) mission is a six beam, low energy, high repetition rate, 532-nanometer laser transmitter with photon counting detectors. Although designed primarily for detecting height changes in ice caps, sea ice, and vegetation, the polar-orbiting satellite will observe global surface water during its designed three-year life span, including inland waterbodies, coasts, and open oceans. In preparation for the mission, an ICESat-2 prototype, the Multiple Altimeter Beam Experimental Lidar (MABEL), was built and flown on high-altitude aircraft experiments over a range of inland and near-shore targets. The purpose was to test the ATLAS concept and to provide a database for developing an algorithm that detects along track surface water height and light penetration under a range of atmospheric and water conditions. The current analysis examines the data sets of three MABEL transects observed from 20 kilometers above ground of coastal and inland waters conducted in 2012 and 2013. Transects ranged from about 2 to 12 kilometers in length and included the middle Chesapeake Bay, the near-shore Atlantic coast at Virginia Beach, and Lake Mead. Results indicate MABEL's high capability for retrieving surface water height statistics with a mean height precision ofapproximately 5-7 centimeters per 100-meter segment length. Profiles of attenuated subsurface backscatter, characterized using a Signal to Background Ratio written in Log10 base, or LSBR (sub 0), were observed over a range of 1.3 to 9.3 meters, depending on water clarity and atmospheric background. Results indicate that observable penetration depth, although primarily dependent on water properties, was greatest when the solar background rate was low. Near-shore bottom reflectance was detected only at the Lake Mead site down to a maximum of 10 meters under a clear night sky and low turbidity of approximately 1.6 Nephelometric Turbidity Units (NTU). The overall results suggest that the feasibility of retrieving operational surface water height statistics from space-based photon counting systems such as ATLAS is very high for resolutions down to about 100 meters, even in partly cloudy conditions. The capability to observe subsurface backscatterprofiles is achievable but requires much longer transects of several hundreds of meters.

During the melting of a snowpack, snow water equivalent (SWE) can be correlated to snow-covered area (SCA) once snow-free areas appear, which is when SCA begins to decrease below 100%. This amount of SWE is called the threshold SWE. Daily SWE data from snow telemetry stations were related to SCA derived from moderate-resolution imaging spectro radiometer images to produce snow-cover depletion curves. The snow depletion curves were created for an 80,000 sq km domain across southern Wyoming and northern Colorado encompassing 54 snow telemetry stations. Eight yearly snow depletion curves were compared, and it is shown that the slope of each is a function of the amount of snow received. Snow-cover depletion curves were also derived for all the individual stations, for which the threshold SWE could be estimated from peak SWE and the topography around each station. A stations peak SWE was much more important than the main topographic variables that included location, elevation, slope, and modelled clear sky solar radiation. The threshold SWE mostly illustrated inter-annual consistency.

Satellite retrievals of snow depth and water equivalent (SWE) are critical for monitoring watershed scale processes around the world. However, the problem is especially challenging in mountainous regions where complex heterogeneities limit the utility of low resolution satellite sensors. The Geoscience Laser Altimeter Sensor (GLAS) aboard the Ice, Cloud, and land Elevation Satellite (ICESat) collected surface elevation data along near-repeat reference transects over land areas from 2003-2009. Although intended for monitoring ice caps and sea ice, the seven year global GLAS data base has provided unprecedented opportunity to test the capability of satellite lidar technology for estimating snow depth over land. GLAS single track and low repeat frequency does not provide data sufficient for operational estimates. However, its comparatively small footprint size of -65 m and its database of seasonal repeat observations during both snow and no-snow conditions have been sufficient to evaluate the potential of spacebased lidar altimetry for estimating snow depth. Recent analysis of ICESat elevations in the Uinta Mountains in NE Utah provide encouraging results for watershed scale estimates of snow depth. Research reported here focuses on the sensitivity of several versions of an ICESat snow depth algorithm to a range of landscape types defined by vegetation cover, slope and roughness. Results are compared to available SNOTEL data.

A method is presented for estimating subaqueous integrated backscatter from the CALIOP lidar. The algorithm takes into account specular reflection of laser light, laser scattering by wind-generated foam as well as sun glint and solar scattering from the foam Analyses show that the estimated subaqueous integrated backscatter is most sensitive to the estimate of transmittance used in the atmospheric correction, and is very insensitive to the estimate of wind speed used. As a case study, CALIOP data over Tampa Bay were compared to MODIS 645 nm remote sensing reflectance, which previously has been shown to be nearly linearly related to turbidity. The results indicate good correlation on nearly all CALIOP clear-free dates during the period 2006 through 2007, particularly those with relatively high atmospheric transmittance. When data are composited over the entire period the correlation is reduced but still statistically significant, an indication of variability in the biogeochemical composition in the water. Overall, the favorable results show promise for the application of satellite lidar integrated backscatter in providing information about subsurface backscatter properties, which can be extracted using appropriate models

Many earth science modeling applications employ continuous input data fields derived from satellite data. Environmental factors, sensor limitations and algorithmic constraints lead to data products of inherently variable quality. This necessitates interpolation of one form or another in order to produce high quality input fields free of missing data. The present research tests several interpolation techniques as applied to satellite-derived leaf area index, an important quantity in many global climate and ecological models. The study evaluates and applies a variety of interpolation techniques for the Moderate Resolution Imaging Spectroradiometer (MODIS) Leaf-Area Index Product over the time period 2001-2006 for a region containing the conterminous United States. Results indicate that the accuracy of an individual interpolation technique depends upon the underlying land cover. Spatial interpolation provides better results in forested areas, while temporal interpolation performs more effectively over non-forest cover types. Combination of spatial and temporal approaches offers superior interpolative capabilities to any single method, and in fact, generation of continuous data fields requires a hybrid approach such as this.

River discharge and changes in lake, reservoir and wetland water storage are critical terms in the global surface water balance, yet they are poorly observed globally and the prospects for adequate observations from in-situ networks are poor (Alsdorf et al., 2003). The NASA-sponsored Surface Water Working Group has established a framework for advancing satellite observations of river discharge and water storage changes which focuses on obtaining measurements of water surface height (stage), slope, and extent. Satellite laser altimetry, which can achieve centimeter-level elevation precision for single, small laser footprints, provides a method to obtain these inland water parameters and contribute to global water balance monitoring. Since its launch in January, 2003 the Ice, Cloud, and land Elevation Satellite (ICESat), a NASA Earth Observing System mission, has achieved over 540 million laser pulse observations of ice sheet, ocean surface, land topography, and inland water elevations and cloud and aerosol height distributions. By recording the laser backscatter from 80 m diameter footprints spaced 175 m along track, ICESat acquires globally-distributed elevation profiles, using a 1064 nm laser altimeter channel, and cloud and aerosol profiles, using a 532 nm atmospheric lidar channel. The ICESat mission has demonstrated the following laser altimeter capabilities relevant to observations of inland water: (1) elevation measurements with a precision of 2 to 3 cm for flat surfaces, suitable for detecting river surface slopes along long river reaches or between multiple crossings of a meandering river channel, (2) from the laser backscatter waveform, detection of water surface elevations beneath vegetation canopies, suitable for measuring water stage in flooded forests, (3) single pulse absolute elevation accuracy of about 50 cm (1 sigma) for 1 degree sloped surfaces, with calibration work in progress indicating that a final accuracy of about 12 cm (1 sigma) will be achieved for clear atmosphere conditions, suitable for detection of stage changes through time, (4) ability to precisely point the spacecraft so as to position the laser profile on the Earth the surface with a cross-track accuracy of 50 m (1 sigma), enabling small water bodies and specific locations to be targeted and re-observed through time, (5) adequate signal levels from specular water surfaces up to 5 degrees off-nadir, enabling complete global access to any location on the Earth's surface from the ICESat repeat orbit by off-nadir pointing, and (6) day and night operation with successful laser ranging to the Earth's surface through thin to moderate cloud cover, enabling more frequent measurements than can be achieved by passive optical sensors. Here we illustrate these capabilities by showing ICESat observations through time for selected river and lake locations.

There are remote sensing and science objectives. The remote sensing objectives are: To develop and test a theoretical method for estimating local momentum aerodynamic roughness length, z(sub 0m), using satellite multispectral imagery. To adapt the method to the MODIS imagery. To develop a high-resolution (approx. 1km) gridded dataset of local momentum roughness for the continental United States and southern Canada, using MODIS imagery and other MODIS derived products. The science objective is: To determine the sensitivity of improved satellite-derived (MODIS-) estimates of surface roughness on the momentum and scalar fluxes, within the context of 3-D atmospheric modeling.

There has been increasing effort in recent years to employ satellite remotely sensed data to identify and map vector habitat and malaria transmission risk in data sparse environments. In the current investigation, available satellite and other land surface climatology data products are employed in short-term forecasting of infection rates in the Mpumalanga Province of South Africa, using a multivariate autoregressive approach. The climatology variables include precipitation, air temperature and other land surface states computed by the Off-line Land-Surface Global Assimilation System (OLGA) including soil moisture and surface evaporation. Satellite data products include the Normalized Difference Vegetation Index (NDVI) and other forcing data used in the Goddard Earth Observing System (GEOS-1) model. Predictions are compared to long- term monthly records of clinical and microscopic diagnoses. The approach addresses the high degree of short-term autocorrelation in the disease and weather time series. The resulting model is able to predict 11 of the 13 months that were classified as high risk during the validation period, indicating the utility of applying antecedent climatic variables to the prediction of malaria incidence for the Mpumalanga Province.

The influence of soil heterogeneity on the partitioning of mesoscale land surface energy fluxes at diurnal time scales is investigated over a 10(exp 6) sq km domain centered on the Little Washita Basin, Oklahoma, for the period June 10 - 18, 1992. The sensitivity study is carried out using MM5/PLACE, the Penn State/NCAR MM5 model enhanced with the Parameterization for Land-Atmosphere-Cloud Exchange or PLACE. PLACE is a one-dimensional land surface model possessing detailed plant and soil water physics algorithms, multiple soil layers, and the capacity to model subgrid heterogeneity. A series of 12-hour simulations were conducted with identical atmospheric initialization and land surface characterization but with different initial soil moisture and texture. A comparison then was made of the simulated land surface energy flux fields, the partitioning of net radiation into latent and sensible heat, and the soil moisture fields. Results indicate that heterogeneity in both soil moisture and texture affects the spatial distribution and partitioning of mesoscale energy balance. Spatial averaging results in an overprediction of latent heat flux, and an underestimation of sensible heat flux. In addition to the primary focus on the partitioning of the land surface energy, the modeling effort provided an opportunity to examine the issue of initializing the soil moisture fields for coupled three-dimensional models. For the present case, the initial soil moisture and temperature were determined from off-line modeling using PLACE at each grid box, driven with a combination of observed and assimilated data fields.

Parameterizations of the frontal area index and canopy area index of natural or randomly distributed plants are developed, and applied to the estimation of local aerodynamic roughness using satellite imagery. The formulas are expressed in terms of the subpixel fractional vegetation cover and one non-dimensional geometric parameter that characterizes the plant's shape. Geometrically similar plants and Poisson distributed plant centers are assumed. An appropriate averaging technique to extend satellite pixel-scale estimates to larger scales is provided. The parameterization is applied to the estimation of aerodynamic roughness using satellite imagery for a 2.3 sq km coniferous portion of the Landes Forest near Lubbon, France, during the 1986 HAPEX-Mobilhy Experiment. The canopy area index is estimated first for each pixel in the scene based on previous estimates of fractional cover obtained using Landsat Thematic Mapper imagery. Next, the results are incorporated into Raupach's (1992, 1994) analytical formulas for momentum roughness and zero-plane displacement height. The estimates compare reasonably well to reference values determined from measurements taken during the experiment and to published literature values. The approach offers the potential for estimating regionally variable, vegetation aerodynamic roughness lengths over natural regions using satellite imagery when there exists only limited knowledge of the vegetated surface.

Parameterizations of the frontal area index and canopy area index of natural or randomly distributed plants are developed, and applied to the estimation of local aerodynamic roughness using satellite imagery. The formulas are expressed in terms of the subpixel fractional vegetation cover and one non-dimensional geometric parameter that characterizes the plant's shape. Geometrically similar plants and Poisson distributed plant centers are assumed. An appropriate averaging technique to extend satellite pixel-scale estimates to larger scales is provided. ne parameterization is applied to the estimation of aerodynamic roughness using satellite imagery for a 2.3 sq km coniferous portion of the Landes Forest near Lubbon, France, during the 1986 HAPEX-Mobilhy Experiment. The canopy area index is estimated first for each pixel in the scene based on previous estimates of fractional cover obtained using Landsat Thematic Mapper imagery. Next, the results are incorporated into Raupach's (1992, 1994) analytical formulas for momentum roughness and zero-plane displacement height. The estimates compare reasonably well to reference values determined from measurements taken during the experiment and to published literature values. The approach offers the potential for estimating regionally variable, vegetation aerodynamic roughness lengths over natural regions using satellite imagery when there exists only limited knowledge of the vegetated surface.

The bulk properties of discontinuous vegetation canopies are estimated at subpixel scales by applying the method of moments to a linear stochastic geometric model of canopy-soil reflectance and one set of multispectral observations without ground truth. The procedure involves the formulation of conditional moments for subsets of pixels that possess similar properties and can be identified through their common orientation in red-infrared scattergrams. The analysis is facilitated by assuming geometric similarity among the canopy elements and by formulating a sampling scale ratio in terms of the bulk geometric scales of the canopy and the pixel. Three versions of the method are demonstrated using two simulated scenes and an actual forested watershed for which aerial radiometric data and corresponding ground truth were obtained.

The functional relation among subpixel canopy cover, illuminated soil, and shadowed soil, which progressively develops with increasing pixel size, is investigated for Poisson distributed plants using a geometric canopy simulation model. An analytical relation among cover components is shown to be applicable when the scale of the pixel is much larger than the scale of the plant and ground shadow. The analysis is facilitated through the use of a nondimensional solar-geometric similarity parameter, eta, equal to the ratio of the area of one plant canopy to its associated ground shadow area, as viewed from nadir. A sampling scale ratio, defined as the ratio of the area of the pixel to the mean area of a single plant shadow, is tested as a quantitative criterion to evaluate when the functional relation among subpixel components occurs. The results of a remote sensing experiment over a natural conifer landscape provide preliminary confirmation of the theoretical analysis.

An apparatus for acquiring bidirectional reflectance-factor data was constructed and used over four surface types. Data sets were obtained over a headed wheat canopy, bare soil having several different roughness conditions, playa (dry lake bed), and gypsum sand. Results are presented in terms of relative bidirectional reflectance factors (BRFs) as a function of view angle at a number of solar zenith angles, nadir BRFs as a function of solar zenith angles, and, for wheat, vegetation indices as related to view and solar zenith angles. The wheat canopy exhibited the largest BRF changes with view angle. BRFs for the red and the NIR bands measured over wheat did not have the same relationship with view angle. NIR/Red ratios calculated from nadir BRFs changed by nearly a factor of 2 when the solar zenith angle changed from 20 to 50 degs. BRF versus view angle relationships were similar for soils having smooth and intermediate rough surfaces but were considerably different for the roughest surface. Nadir BRF versus solar-zenith angle relationships were distinctly different for the three soil roughness levels. Of the various surfaces, BRFs for gypsum sand changed the least with view angle (10 percent at 30 degs).

An analytical framework is provided for examining the physically based behavior of the normalized difference vegetation index (NDVI) in terms of the variability in bulk subpixel landscape components and with respect to variations in pixel scales, within the context of the stochastic-geometric canopy reflectance model. Analysis focuses on regional scale variability in horizontal plant density and soil background reflectance distribution. Modeling is generalized to different plant geometries and solar angles through the use of the nondimensional solar-geometric similarity parameter. Results demonstrate that, for Poisson-distributed plants and for one deterministic distribution, NDVI increases with increasing subpixel fractional canopy amount, decreasing soil background reflectance, and increasing shadows, at least within the limitations of the geometric reflectance model. The NDVI of a pecan orchard and a juniper landscape is presented and discussed.

A stochastic-geometric landsurface reflectance model is formulated and tested for the parameterization of spatially variable vegetation and soil at subpixel scales using satellite multispectral images without ground truth. Landscapes are conceptualized as 3-D Lambertian reflecting surfaces consisting of plant canopies, represented by solid geometric figures, superposed on a flat soil background. A computer simulation program is developed to investigate image characteristics at various spatial aggregations representative of satellite observational scales, or pixels. The evolution of the shape and structure of the red-infrared space, or scattergram, of typical semivegetated scenes is investigated by sequentially introducing model variables into the simulation. The analytical moments of the total pixel reflectance, including the mean, variance, spatial covariance, and cross-spectral covariance, are derived in terms of the moments of the individual fractional cover and reflectance components. The moments are applied to the solution of the inverse problem: The estimation of subpixel landscape properties on a pixel-by-pixel basis, given only one multispectral image and limited assumptions on the structure of the landscape. The landsurface reflectance model and inversion technique are tested using actual aerial radiometric data collected over regularly spaced pecan trees, and using both aerial and LANDSAT Thematic Mapper data obtained over discontinuous, randomly spaced conifer canopies in a natural forested watershed. Different amounts of solar backscattered diffuse radiation are assumed and the sensitivity of the estimated landsurface parameters to those amounts is examined.

A physically based linear stochastic geometric canopy soil reflectance model is presented for characterizing spatial variability of semivegetated landscapes at subpixel and regional scales. Landscapes are conceptualized as stochastic geometric surfaces, incorporating not only the variability in geometric elements, but also the variability in vegetation and soil background reflectance which can be important in some scenes. The model is used to investigate several possible mechanisms which contribute to the often observed characteristic triangular shape of red-infrared scattergrams of semivegetated landscapes. Scattergrams of simulated semivegetated scenes are analyzed with respect to the scales of the satellite pixel and subpixel components. Analysis of actual aerial radiometric data of a pecan orchard is presented in comparison with ground observations as preliminary confirmation of the theoretical results.

The estimation of the spatially variable surface moisture and heat fluxes of natural, semivegetated landscapes is difficult due to the highly random nature of the vegetation (e.g., plant species, density, and stress) and the soil (e.g., moisture content, and soil hydraulic conductivity). The solution to that problem lies, in part, in the use of satellite remotely sensed data, and in the preparation of those data in terms of the physical properties of the plant and soil. The work was focused on the development and testing of a stochastic geometric canopy-soil reflectance model, which can be applied to the physically-based interpretation of LANDSAT images. The model conceptualizes the landscape as a stochastic surface with bulk plant and soil reflective properties. The model is particularly suited for regional scale investigations where the quantification of the bulk landscape properties, such as fractional vegetation cover, is important on a pixel by pixel basis. A summary of the theoretical analysis and the preliminary testing of the model with actual aerial radiometric data is provided.

A physically based linear stochastic geometric canopy soil reflectance model is presented for characterizing spatial variability of semivegetated landscapes at subpixel and regional scales. Landscapes are conceptualized as stochastic geometric surfaces, incorporating not only the variability in geometric elements, but also the variability in vegetation and soil background reflectance which can be important in some scenes. The model is used to investigate several possible mechanisms which contribute to the often observed characteristic triangular shape of red-infrared scattergrams of semivegetated landscapes. Scattergrams of simulated and semivegetated scenes are analyzed with respect to the scales of the satellite pixel and subpixel components. Analysis of actual aerial radiometric data of a pecan orchard is presented in comparison with ground observations as preliminary confirmation of the theoretical results.

The present report summarizes the various approaches relevant to estimating canopy cover at subpixel resolution. The approaches are based on physical models of radiative transfer in non-homogeneous canopies and on empirical methods. The effects of vegetation shadows and topography are examined. Simple versions of the model are tested, using the Taos, New Mexico Study Area database. Emphasis has been placed on using relatively simple models requiring only one or two bands. Although most methods require some degree of ground truth, a two-band method is investigated whereby the percent cover can be estimated without ground truth by examining the limits of the data space. Future work is proposed which will incorporate additional surface parameters into the canopy cover algorithm, such as topography, leaf area, or shadows. The method involves deriving a probability density function for the percent canopy cover based on the joint probability density function of the observed radiances.