USGS ScienceSearch

SEARCH · USGS Science

Results for “Algorithms”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

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

At least 109 records · Page 6Linked to original sources

Seasonal cultivated and fallow cropland mapping using MODIS-based automated cropland classification algorithm

Increasing drought occurrences and growing populations demand accurate, routine, and consistent cultivated and fallow cropland products to enable water and food security analysis. The overarching goal of this research was to develop and test automated cropland classification algorithm (ACCA) that provide accurate, consistent, and repeatable information on seasonal cultivated as well as seasonal fallow cropland extents and areas based on the Moderate Resolution Imaging Spectroradiometer remote sensing data. Seasonal ACCA development process involves writing series of iterative decision tree codes to separate cultivated and fallow croplands from noncroplands, aiming to accurately mirror reliable reference data sources. A pixel-by-pixel accuracy assessment when compared with the U.S. Department of Agriculture (USDA) cropland data showed, on average, a producer’s accuracy of 93% and a user’s accuracy of 85% across all months. Further, ACCA-derived cropland maps agreed well with the USDA Farm Service Agency crop acreage-reported data for both cultivated and fallow croplands with R-square values over 0.7 and field surveys with an accuracy of ≥95% for cultivated croplands and ≥76% for fallow croplands. Our results demonstrated the ability of ACCA to generate cropland products, such as cultivated and fallow cropland extents and areas, accurately, automatically, and repeatedly throughout the growing season.

California

Permissible Home Range Estimation (PHRE) in restricted habitats: A new algorithm and an evaluation for sea otters

Parametric and nonparametric kernel methods dominate studies of animal home ranges and space use. Most existing methods are unable to incorporate information about the underlying physical environment, leading to poor performance in excluding areas that are not used. Using radio-telemetry data from sea otters, we developed and evaluated a new algorithm for estimating home ranges (hereafter Permissible Home Range Estimation, or “PHRE”) that reflects habitat suitability. We began by transforming sighting locations into relevant landscape features (for sea otters, coastal position and distance from shore). Then, we generated a bivariate kernel probability density function in landscape space and back-transformed this to geographic space in order to define a permissible home range. Compared to two commonly used home range estimation methods, kernel densities and local convex hulls, PHRE better excluded unused areas and required a smaller sample size. Our PHRE method is applicable to species whose ranges are restricted by complex physical boundaries or environmental gradients and will improve understanding of habitat-use requirements and, ultimately, aid in conservation efforts.

PLoS ONE

Automated Cropland Classification Algorithm (ACCA) for California using multi-sensor remote sensing

Increasing pressure to feed the growing population with scarce water resources requires accurate and routine cropland mapping. This paper develops and implements a rule-based automated cropland classification algorithm (ACCA) using multi-sensor remote sensing data. Pixel-by-pixel accuracy assessments showed that ACCA produced an overall accuracy of 96 percent (K hat = 0.8) when tested using independent data layers. Furthermore, ACCA-generated county cropland areas showed high agreement (R-square values 0.94) when compared with three independent data sources: (a) US Department of Agriculture (USDA) cropland data layer derived cropland areas, (b) county specific crop acreage data from the Farm Service Agency, and (c) the Census of Agriculture data for the 58 counties in California. Our results demonstrate the ability of ACCA to generate cropland extent and areas over space and time, in an automated fashion with high degree of accuracies year after year, greatly contributing to food and water security analysis and decision making.

California

An algorithm for locating candidate labeling boxes within a polygon

Vector-based geographic information systems usually require annotation, such as a polygon number or attribute data, in a suitable location within a polygon. Traditional methods usually compute the polygon centroid, test the centroid for inclusion or exclusion, and select some alternative point when the centroid falls outside the polygon. Two problems are associated with this approach: (1) the text can be centered on the point, but may be placed in a visually awkward place, and (2) part of the text may fall outside the polygon and may overlap other polygon boundaries or other text labels. An algorithm is presented that circumvents both of these problems, by computing a number of horizontal candidate labeling rectangles (boxes) within a polygon from which a suitable selection can be made or from which one may conclude that the text label does not fit the polygon.

American Cartographer

Effect of fixing earthquake depth in ShakeAlert algorithms on performance for intraslab earthquakes

We investigate whether assuming a fixed shallow depth in the ShakeAlert network‐based earthquake early warning system is sufficient to produce accurate ground‐motion based alerts for intraslab earthquakes. ShakeAlert currently uses a fixed focal depth of 8 km to estimate earthquake location and magnitude. This is an appropriate way to reduce computational costs without compromising alert accuracy in California, where earthquakes typically occur on shallow crustal faults. In the Pacific Northwest (PNW), however, the most common moderate‐magnitude events occur within the subducting Juan de Fuca slab at depths between ∼35 and 65 km. Using a dataset of seismic recordings from 37 M w 4.5+ intraslab earthquakes from the PNW and Chile, we replay events through the Earthquake Point‐Source Integrated Code and eqInfo2GM algorithms to estimate source parameters and compute modified Mercalli intensity (MMI) alert threshold contours. Each event is replayed twice—once using a fixed 8 km depth and a second time using the actual catalog earthquake depth. For each depth scenario, we analyze MMI III and IV contours using various performance metrics to determine the number of correctly alerted sites and measure warning times. We determine that shallow depth replays are more likely to produce errors in location estimates of greater than 50 km if the event is located outside of a seismic network. When located within a seismic network, shallow and catalog depth replays have similar epicenter estimates. Results show that applying catalog earthquake depth does not improve the accuracy of magnitude estimates or MMI alert threshold contours, or increase warning times. We conclude that using a fixed shallow earthquake depth for intraslab earthquakes will not significantly impact alert accuracy in the PNW.

California, Oregon, Washington

Evaluation of the expected moments algorithm and a multiple low-outlier test for flood frequency analysis at streamgaging stations in Arizona

Flooding is among the costliest natural disasters in terms of loss of life and property in Arizona, which is why the accurate estimation of flood frequency and magnitude is crucial for proper structural design and accurate floodplain mapping. Current guidelines for flood frequency analysis in the United States are described in Bulletin 17B (B17B), yet since B17B’s publication in 1982 (Interagency Advisory Committee on Water Data, 1982), several improvements have been proposed as updates for future guidelines. Two proposed updates are the Expected Moments Algorithm (EMA) to accommodate historical and censored data, and a generalized multiple Grubbs-Beck (MGB) low-outlier test. The current guidelines use a standard Grubbs-Beck (GB) method to identify low outliers, changing the determination of the moment estimators because B17B uses a conditional probability adjustment to handle low outliers while EMA censors the low outliers. B17B and EMA estimates are identical if no historical information or censored or low outliers are present in the peak-flow data. EMA with MGB (EMA-MGB) test was compared to the standard B17B (B17B-GB) method for flood frequency analysis at 328 streamgaging stations in Arizona. The methods were compared using the relative percent difference (RPD) between annual exceedance probabilities (AEPs), goodness-of-fit assessments, random resampling procedures, and Monte Carlo simulations. The AEPs were calculated and compared using both station skew and weighted skew. Streamgaging stations were classified by U.S. Geological Survey (USGS) National Water Information System (NWIS) qualification codes, used to denote historical and censored peak-flow data, to better understand the effect that nonstandard flood information has on the flood frequency analysis for each method. Streamgaging stations were also grouped according to geographic flood regions and analyzed separately to better understand regional differences caused by physiography and climate. The B17B-GB and EMA-MGB RPD-boxplot results showed that the median RPDs across all streamgaging stations for the 10-, 1-, and 0.2-percent AEPs, computed using station skew, were approximately zero. As the AEP flow estimates decreased (that is, from 10 to 0.2 percent AEP) the variability in the RPDs increased, indicating that the AEP flow estimate was greater for EMA-MGB when compared to B17B-GB. There was only one RPD greater than 100 percent for the 10- and 1-percent AEP estimates, whereas 19 RPDs exceeded 100 percent for the 0.2-percent AEP. At streamgaging stations with low-outlier data, historical peak-flow data, or both, RPDs ranged from −84 to 262 percent for the 0.2-percent AEP flow estimate. When streamgaging stations were separated by the presence of historical peak-flow data (that is, no low outliers or censored peaks) or by low outlier peak-flow data (no historical data), the results showed that RPD variability was greatest for the 0.2-AEP flow estimates, indicating that the treatment of historical and (or) low-outlier data was different between methods and that method differences were most influential when estimating the less probable AEP flows (1, 0.5, and 0.2 percent). When regional skew information was weighted with the station skew, B17B-GB estimates were generally higher than the EMA-MGB estimates for any given AEP. This was related to the different regional skews and mean square error used in the weighting procedure for each flood frequency analysis. The B17B-GB weighted skew analysis used a more positive regional skew determined in USGS Water Supply Paper 2433 (Thomas and others, 1997), while the EMA-MGB analysis used a more negative regional skew with a lower mean square error determined from a Bayesian generalized least squares analysis. Regional groupings of streamgaging stations reflected differences in physiographic and climatic characteristics. Potentially influential low flows (PILFs) were more prevalent in arid regions of the State, and generally AEP flows were larger with EMA-MGB than with B17B-GB for gaging stations with PILFs. In most cases EMA-MGB curves would fit the largest floods more accurately than B17B-GB. In areas of the State with more baseflow, such as along the Mogollon Rim and the White Mountains, streamgaging stations generally had fewer PILFs and more positive skews, causing estimated AEP flows to be larger with B17B-GB than with EMA-MGB. The effect of including regional skew was similar for all regions, and the observed pattern was increasingly greater B17B-GB flows (more negative RPDs) with each decreasing AEP quantile. A variation on a goodness-of-fit test statistic was used to describe each method’s ability to fit the largest floods. The mean absolute percent difference between the measured peak flows and the log-Pearson Type 3 (LP3)-estimated flows, for each method, was averaged over the 90th, 75th, and 50th percentiles of peak-flow data at each site. In most percentile subsets, EMA-MGB on average had smaller differences (1 to 3 percent) between the observed and fitted value, suggesting that the EMA-MGB-LP3 distribution is fitting the observed peak-flow data more precisely than B17B-GB. The smallest EMA-MGB percent differences occurred for the greatest 10 percent (90th percentile) of the peak-flow data. When stations were analyzed by USGS NWIS peak flow qualification code groups, the stations with historical peak flows and no low outliers had average percent differences as high as 11 percent greater for B17B-GB, indicating that EMA-MGB utilized the historical information to fit the largest observed floods more accurately. A resampling procedure was used in which 1,000 random subsamples were drawn, each comprising one-half of the observed data. An LP3 distribution was fit to each subsample using B17B-GB and EMA-MGB methods, and the predicted 1-percent AEP flows were compared to those generated from distributions fit to the entire dataset. With station skew, the two methods were similar in the median percent difference, but with weighted skew EMA-MGB estimates were generally better. At two gages where B17B-GB appeared to perform better, a large number of peak flows were deemed to be PILFs by the MGB test, although they did not appear to depart significantly from the trend of the data (step or dogleg appearance). At two gages where EMA-MGB performed better, the MGB identified several PILFs that were affecting the fitted distribution of the B17B-GB method. Monte Carlo simulations were run for the LP3 distribution using different skews and with different assumptions about the expected number of historical peaks. The primary benefit of running Monte Carlo simulations is that the underlying distribution statistics are known, meaning that the true 1-percent AEP is known. The results showed that EMA-MGB performed as well or better in situations where the LP3 distribution had a zero or positive skew and historical information. When the skew for the LP3 distribution was negative, EMA-MGB performed significantly better than B17B-GB and EMA-MGB estimates were less biased by more closely estimating the true 1-percent AEP for 1, 2, and 10 historical flood scenarios.

Arizona

An empirical algorithm for estimating agricultural and riparian evapotranspiration using MODIS Enhanced Vegetation Index and ground measurements of ET. II. Application to the lower Colorado River, U.S.

Large quantities of water are consumed by irrigated crops and riparian vegetation in western U.S. irrigation districts. Remote sensing methods for estimating evaporative water losses by soil and vegetation (evapotranspiration, ET) over wide river stretches are needed to allocate water for agricultural and environmental needs. We used the Enhanced Vegetation Index (EVI) from MODIS sensors on the Terra satellite to scale ET over agricultural and riparian areas along the Lower Colorado River in the southwestern U.S., using a linear regression equation between ET of riparian plants and alfalfa measured on the ground, and meteorological and remote sensing data, with an error or uncertainty of about 20%. The algorithm was applied to irrigation districts and riparian areas from Lake Mead to the U.S./Mexico border. The results for agricultural crops were similar to results produced by crop coefficients developed for the irrigation districts along the river. However, riparian ET was only half as great as crop coefficient estimates set by expert opinion, equal to about 40% of reference crop evapotranspiration. Based on reported acreages in 2007, agricultural crops (146,473 ha) consumed 2.2 × 10 9 m 3 yr −1 of water. All riparian shrubs and trees (47,014 ha) consumed 3.8 × 10 8 m 3 yr −1 , of which saltcedar, the dominant riparian shrub (25,044 ha), consumed 1.8 × 10 8 m 3 yr −1 , about 1% of the annual flow of the river. This method could supplement existing protocols for estimating ET by providing an estimate based on the actual state of the canopy as determined by frequent-return satellite data.

Arizona, California, Nevada

SSEBop evapotranspiration estimates using synthetically derived Landsat data from the continuous change detection and classification algorithm

The operational Simplified Surface Energy Balance (SSEBop) model has been utilized to generate gridded evapotranspiration data from Landsat images. These estimates are primarily driven by two sources of information: reference evapotranspiration and Landsat land surface temperature (LST) values. Hence, SSEBop is limited by the availability of Landsat data. Here, in this proof-of-concept paper, we utilize the Continuous Change Detection and Classification (CCDC) algorithm to generate synthetic Landsat data, which are then used as input for SSEBop to generate evapotranspiration estimates for six target areas in the continental United States, representing forests, shrublands, and irrigated agriculture. These synthetic land cover data are then used to generate the LST data required for SSEBop evapotranspiration estimates. The synthetic LST, evaporative fractions, and evapotranspiration data from CCDC closely mirror the phenological cycles in the observed Landsat data. Across the six sites, the median correlation in seasonal LST was 0.79, and the median correlation in seasonal evapotranspiration was 0.8. The median root mean squared error (RMSE) values were 2.82 °C for LST and 0.50 mm/day for actual evapotranspiration. CCDC predictions typically underestimate the average evapotranspiration by less than 1 mm/day. The average performance of the CCDC evaporative fractions, and corresponding evapotranspiration estimates, were much better than the initial LST estimates and, therefore, promising. Future work could include bias correction to improve CCDC’s ability to accurately reproduce synthetic Landsat data during the summer, allowing for more accurate evapotranspiration estimates, and determining the ability of SSEBop to predict regional evapotranspiration at seasonal timescales based on projected land cover change from CCDC.

Remote Sensing

Evaluation of the initial thematic output from a continuous change-detection algorithm for use in automated operational land-change mapping by the U.S. Geological Survey

The U.S. Geological Survey (USGS) has begun the development of operational, 30-m resolution annual thematic land cover data to meet the needs of a variety of land cover data users. The Continuous Change Detection and Classification (CCDC) algorithm is being evaluated as the likely methodology following early trials. Data for training and testing of CCDC thematic maps have been provided by the USGS Land Cover Trends (LC Trends) project, which offers sample-based, manually classified thematic land cover data at 2755 probabilistically located sample blocks across the conterminous United States. These samples represent a high quality, well distributed source of data to train the Random Forest classifier invoked by CCDC. We evaluated the suitability of LC Trends data to train the classifier by assessing the agreement of annual land cover maps output from CCDC with output from the LC Trends project within 14 Landsat path/row locations across the conterminous United States. We used a small subset of circa 2000 data from the LC Trends project to train the classifier, reserving the remaining Trends data from 2000, and incorporating LC Trends data from 1992, to evaluate measures of agreement across time, space, and thematic classes, and to characterize disagreement. Overall agreement ranged from 75% to 98% across the path/rows, and results were largely consistent across time. Land cover types that were well represented in the training data tended to have higher rates of agreement between LC Trends and CCDC outputs. Characteristics of disagreement are being used to improve the use of LC Trends data as a continued source of training information for operational production of annual land cover maps.

Remote Sensing

Efficient implementation of the Metropolis-Hastings algorithm, with application to the Cormack?Jolly?Seber model

Judicious choice of candidate generating distributions improves efficiency of the Metropolis-Hastings algorithm. In Bayesian applications, it is sometimes possible to identify an approximation to the target posterior distribution; this approximate posterior distribution is a good choice for candidate generation. These observations are applied to analysis of the Cormack?Jolly?Seber model and its extensions.

Environmental and Ecological Statistics

An algorithm for locating candidate labelling boxes within a polygon

Vector-based geographic information systems usually require annotation, such as a polygon number or attribute data, in a suitable location within a polygon. Traditional methods usually compute the polygon centroid, test the centroid for inclusion or exclusion, and select some alternative point when the centroid falls outside the polygon. Two problems are associated with this approach: (1) the text can be centered on the point, but may be placed in a visually awkward place, and (2) part of the text may fall outside the polygon and may overlap other polygon boundaries or other text labels. An algorithm is presented that circumvents both of these problems, by computing a number of horizontal candidate labelling rectangles (boxes) within a polygon from which a suitable selection can be made or from which one may conclude that the text label does not fit the polygon.

Conference Paper

Estimating multivariate ecological variables at high spatial resolution using a cost-effective matching algorithm

Simulation models are valuable tools for estimating ecosystem response to environmental conditions and are particularly relevant for investigating climate change impacts. However, because of high computational requirements, models are often applied over a coarse grid of points or for representative locations. Spatial interpolation of model output can be necessary to guide decision-making, yet interpolation is not straightforward because the interpolated values must maintain the covariance structure among variables. We present methods for two key steps for utilizing limited simulations to generate detailed maps of multivariate and time series output. First, we present a method to select an optimal set of simulation sites that maximize the area represented for a given number of sites. Then, we introduce a multivariate matching approach to interpolate simulation results to detailed maps for the represented area. This approach links simulation output to environmentally analogous matched sites according to user-defined criteria. We demonstrate the methods with case studies using output from (1) an individual-based plant simulation model to illustrate site selection, and (2) an ecosystem water balance simulation model to illustrate interpolation. For the site selection case study, we identified 200 simulation sites that represented 96% of a large study area (1.12 × 10 6 km 2 ) at a ~1-km resolution. For the interpolation case study, we generated ~1-km resolution maps across 4.38 × 10 6 km 2 of drylands in North America from a 10 × 10 km grid of simulated sites. Estimates of interpolation errors using cross validation were low (less than 10% of the range of each variable). Our point selection and interpolation methods, which are available as an easy-to-use R package, provide a means of cost-effectively generating detailed maps of expensive, complex simulation output (e.g., multivariate and time series) at scales relevant for local conservation planning. Our methods are flexible and allow the user to identify relevant matching criteria to balance interpolation uncertainty with areal coverage to enhance inference and decision-making at management-relevant scales across large areas.

Ecosphere

Algorithm and data improvements for version 2.1 of the Climate Hazards center’s InfraRed Precipitation with Stations Data Set

To support global drought early warning, the Climate Hazards Center (CHC) at the University of California, Santa Barbara developed the Climate Hazards center InfraRed Precipitation with Stations (CHIRPS) dataset, in collaboration with the US Geological Survey and NASA SERVIR. Specifically designed to support early warning applications, CHIRPS has high a spatial resolution (0.05°), a long period of record (1981 to the near present), and relatively low latencies. Here we will describe a brief formal analysis of distributional bias in CHIRPS2.0. This analysis reveals, as expected, that CHIRPS2.0 means are very similar to observed station data. However, a closer look suggests that low precipitation values are underestimated and high values are over-estimated in the CHIRPS2.0. We describe a potential correction for this below.

Book chapter

Evaluation of daily stream temperature predictions (1979-2021) across the contiguous United States using a spatiotemporal aware machine learning algorithm

Stream temperature controls a variety of physical and biological processes that affect ecosystems, human health, and economic activities. We used 42 years (1979–2021) of data to predict daily summary statistics of stream temperature across >50,000 stream reaches in the contiguous United States using a recurrent graph convolution network. We comprehensively documented the performance – both across all reaches and by stream type (e.g., reservoir or groundwater influence) – as a baseline for future improvement. The model showed reach-level RMSE of <2 °C with 90 % prediction intervals that contain 90.7 % of observations. We also assessed how the model captured variability in ecologically relevant metrics (e.g., R 2 for annual 7-day maximum = 0.76; R 2 for days exceeding 25 °C = 0.75). This model does not outperform state-of-the-art machine learning efforts (e.g., RMSE ≤1.5 °C) due to a limited input set but does provide the most spatially complete modeling to date to support water availability assessments.

contiguous United States

Automated classifications of topography from DEMs by an unsupervised nested-means algorithm and a three-part geometric signature

An iterative procedure that implements the classification of continuous topography as a problem in digital image-processing automatically divides an area into categories of surface form; three taxonomic criteria–slope gradient, local convexity, and surface texture–are calculated from a square-grid digital elevation model (DEM). The sequence of programmed operations combines twofold-partitioned maps of the three variables converted to greyscale images, using the mean of each variable as the dividing threshold. To subdivide increasingly subtle topography, grid cells sloping at less than mean gradient of the input DEM are classified by designating mean values of successively lower-sloping subsets of the study area (nested means) as taxonomic thresholds, thereby increasing the number of output categories from the minimum 8 to 12 or 16. Program output is exemplified by 16 topographic types for the world at 1-km spatial resolution (SRTM30 data), the Japanese Islands at 270 m, and part of Hokkaido at 55 m. Because the procedure is unsupervised and reflects frequency distributions of the input variables rather than pre-set criteria, the resulting classes are undefined and must be calibrated empirically by subsequent analysis. Maps of the example classifications reflect physiographic regions, geological structure, and landform as well as slope materials and processes; fine-textured terrain categories tend to correlate with erosional topography or older surfaces, coarse-textured classes with areas of little dissection. In Japan the resulting classes approximate landform types mapped from airphoto analysis, while in the Americas they create map patterns resembling Hammond's terrain types or surface-form classes; SRTM30 output for the United States compares favorably with Fenneman's physical divisions. Experiments are suggested for further developing the method; the Arc/Info AML and the map of terrain classes for the world are available as online downloads.

Geomorphology

Identification of spectrally similar materials using the USGS Tetracorder algorithm: The calcite-epidote-chlorite problem

A scheme to discriminate and identify materials having overlapping spectral absorption features has been developed and tested based on the U.S. Geological Survey (USGS) Tetracorder system. The scheme has been applied to remotely sensed imaging spectroscopy data acquired by the Airborne Visible and Infrared Imaging Spectrometer (AVIRIS) instrument. This approach was used to identify the minerals calcite, epidote, and chlorite in the upper Animas River watershed, Colorado. The study was motivated by the need to characterize the distribution of calcite in the watershed and assess its acid-neutralizing potential with regard to acidic mine drainage. Identification of these three minerals is difficult because their diagnostic spectral features are all centered at 2.3 ??m, and have similar shapes and widths. Previous studies overestimated calcite abundance as a result of these spectral overlaps. The use of a reference library containing synthetic mixtures of the three minerals in varying proportions was found to simplify the task of identifying these minerals when used in conjunction with a rule-based expert system. Some inaccuracies in the mineral distribution maps remain, however, due to the influence of a fourth spectral component, sericite, which exhibits spectral absorption features at 2.2 and 2.4 ??m that overlap the 2.3-??m absorption features of the other three minerals. Whereas the endmember minerals calcite, epidote, chlorite, and sericite can be identified by the method presented here, discrepancies occur in areas where all four occur together as intimate mixtures. It is expected that future work will be able to reduce these discrepancies by including reference mixtures containing sericite. ?? 2004 Elsevier Inc. All rights reserved.

Remote Sensing of Environment

A simplified algorithm for temperature conversions

The author is surprised that those who teach temperature conversion formulas have not introduced the easier-to-retain relationship based on the fact that –40 F° is the same temperature as –40 C°.

Journal of Chemical Education

Terrestrial remote sensing science and algorithms planned for EOS/MODIS

The Moderate Resolution Imaging Spectroradiometer (MODIS) will be the primary daily global monitoring sensor on the NASA Earth Observing System (EOS) satellites, scheduled for launch on the EOS-AM platform in June 1998 and the EOS-PM platform in December 2000. MODIS is a 36 channel radiometer covering 0·415-14·235 μm wavelengths, with spatial resolution from 250 m to 1 km at nadir. MODIS will be the primary EOS sensor for providing data on terrestrial biospheric dynamics and process activity. This paper presents the suite of global land products currently planned for EOSDIS implementation, to be developed by the authors of this paper, the MODIS land team (MODLAND). These include spectral albedo, land cover, spectral vegetation indices, snow and ice cover, surface temperature and fire, and a number of biophysical variables that will allow computation of global carbon cycles, hydrologic balances and biogeochemistry of critical greenhouse gases. Additionally, the regular global coverage of these variables will allow accurate surface change detection, a fundamental determinant of global change.

International Journal of Remote Sensing