Geology topics

John L. Dwyer

Publications and source records attributed to John L. Dwyer.

At least 19 recordsLinked to original sources

USGS2021

Semi-centennial of Landsat observations and pending Landsat 9 launch

The first Landsat was placed in orbit on 23 July 1972, followed by a series of missions that have provided nearly continuous, two-satellite 8-day repeat image coverage of the Earth’s land areas for the last half-century. These observations have substantially enhanced our understanding of the Earth’s terrestrial dynamics, both as a major element of the Earth’s physical system, the primary home of humans, and the major source of natural resources that support them. The history of Landsat is complex, reflective of the human systems that sustain it. Despite the conflicted perspectives surrounding the continuation of the program, Landsat has survived based on worldwide recognition of its critical contributions to understanding land dynamics, management of natural resources and Earth system science. Launch of Landsat 9 is anticipated in Fall 2021, and current planning for the next generation, Landsat Next is well underway. The community of Landsat data users is looking forward to another 50 years of the Landsat program.

Photogrammetric Engineering and Remote Sensing (PE↗

USGS2020

Lessons learned implementing an operational continuous United States national land change monitoring capability: The Land Change Monitoring, Assessment, and Projection (LCMAP) approach

Growing demands for temporally specific information on land surface change are fueling a new generation of maps and statistics that can contribute to understanding geographic and temporal patterns of change across large regions, provide input into a wide range of environmental modeling studies, clarify the drivers of change, and provide more timely information for land managers. To meet these needs, the U.S. Geological Survey has implemented a capability to monitor land surface change called the Land Change Monitoring, Assessment, and Projection (LCMAP) initiative. This paper describes the methodological foundations and lessons learned during development and testing of the LCMAP approach. Testing and evaluation of a suite of 10 annual land cover and land surface change data sets over six diverse study areas across the United States revealed good agreement with other published maps (overall agreement ranged from 73% to 87%) as well as several challenges that needed to be addressed to meet the goals of robust, repeatable, and geographically consistent monitoring results from the Continuous Change Detection and Classification (CCDC) algorithm. First, the high spatial and temporal variability of observational frequency led to differences in the number of changes identified, so CCDC was modified such that change detection is dependent on observational frequency. Second, the CCDC classification methodology was modified to improve its ability to characterize gradual land surface changes. Third, modifications were made to the classification element of CCDC to improve the representativeness of training data, which necessitated replacing the random forest algorithm with a boosted decision tree. Following these modifications, assessment of prototype Version 1 LCMAP results showed improvements in overall agreement (ranging from 85% to 90%).

Remote Sensing of Environment↗

USGS2020

The Landsat Burned Area algorithm and products for the conterminous United States

Complete and accurate burned area map data are needed to document spatial and temporal patterns of fires, to quantify their drivers, and to assess the impacts on human and natural systems. In this study, we developed the Landsat Burned Area (BA) algorithm, an update from the Landsat Burned Area Essential Climate Variable (BAECV) algorithm. Here, we present the BA algorithm and products, changes relative to the BAECV algorithm and products, and updated validation metrics. We also present spatial and temporal patterns of burned area across the conterminous U.S., how burned area varies in relation to the number of operational Landsat sensors, and a comparison with other burned area datasets, including the BAECV, Monitoring Trends in Burn Severity (MTBS), GeoMAC, and Moderate Resolution Imaging Spectroradiometer (MODIS) MCD64A1.006 data. The BA algorithm identifies burned areas in analysis ready data (ARD) time-series of Landsat imagery from 1984 through 2018 using machine learning, thresholding, and image segmentation. Validation with reference data from high-resolution commercial satellite imagery resulted in omission and commission error rates averaging 19% and 41%, respectively. In comparison, validation with Landsat reference data had omission and commission error rates averaging 40% and 28%, respectively when burned areas in cultivated crops and pasture/hay land-cover types were excluded. Both validation tests documented lower commission error rates relative to the BAECV products. The amount of burned area detected varies not only in response to climate but also with the number of operational sensors and scenes collected. The combined amount of burned area detected by multiple sensors was larger than from any individual sensor, but there was no significant difference between individual sensors. Therefore, we used BA products from individual sensors to assess trends over time and all available sensors to compare with other existing BA products. From 1984 through 2018, annual burned area averaged 30,000 km2, ranged between 14,000 km2 in 1991 and 46,500 km2 in 2012, and increased over time at a rate of 356 km2/year. Compared to existing burned area products, the new Landsat BA products identified 29% more burned area than the BAECV products (1984–2015), 183% more than the MTBS/GeoMAC products (1984–2018), and 56% more than the MCD64A1.006 products (2003–2018). The products had similar patterns of year-to-year variability; the R2 values of linear regressions between annual burned area were >0.70 with the BAECV products and the MTBS/GeoMAC products, but somewhat lower for the MCD64A1.006 product (R2 = 0.66). The BA products are routinely produced as new Landsat data are collected and provide a unique data source to monitor and assess the spatial and temporal patterns and the impacts of fire.

Remote Sensing of Environment↗

USGS2020

New operational national satellite burned area product

Introduction Lack of consistent spatial and temporal fire information with relevant spatial resolution hinders land management and broad-scale assessments of fire activity, especially in the eastern United States and the Great Plains where fi re is important ecologically and culturally. Remote sensing can be used to monitor fi re activity, augment existing fi re data, and fill information gaps. In particular, Landsat offers one of the most complete time series of remote sensing data sets as the Landsat satellites with spectral bands useful for mapping fi res and burn severity have been operational since 1984. Furthermore, Landsat satellite imagery collect data at a resolution useful for on-the-ground comparisons and management decisions. Methods A gradient-boosting regression model algorithm was used to predict burn probabilities (BP), indicating the likelihood that a pixel had burned in a fire. Then the algorithm translated the burn probability images to burn classification (BC) images using thresholding and region growing. Burned areas smaller than 5 acres were removed to reduce noise. The BP and BC products were generated for Landsat scenes collected from 1984 to present with <= 80% cloud cover and are available for download from USGSs EarthExplorer website (www.earthexplorer. usgs.gov). Annual summaries were also produced from the scene-level BC and BP products. These include: 1) the maximum burn probability value from all scenes in a year, (2) the burn classification count, or the number of times a pixel was classified as burned in a year, (3) the burn date, or the day of year of the first Landsat scene a pixel was classified as burned, and (4) a burn classification which removes burned areas still visible from previous years with a filter. Results The Landsat Burned Area Products identified 183% more burned area than the Monitoring Trends in Burn Severity and 56% more burned area than the MODIS MCD64A1.006 products2. From 1984 through 2018, annual burned area mapped by the Landsat Burned Area Products averaged 30,000 km2, ranged between 14,000 km2 in 1991 and 46,500 km2 in 2012, and increased over time at a rate of 356 km2/year. The Landsat Burned Area products consistently map burned areas 5 acres or larger over time while capturing unburned islands within fi res. Such characteristics are critical for understanding long-term impacts of fi res on human and natural systems. The products are being routinely produced as new Landsat data are collected and demonstrate the utility remote sensing offers for monitoring landscape change.

Report↗

USGS2018

Analysis ready data: Enabling analysis of the Landsat archive

Data that have been processed to allow analysis with a minimum of additional user effort are often referred to as Analysis Ready Data (ARD). The ability to perform large scale Landsat analysis relies on the ability to access observations that are geometrically and radiometrically consistent, and have had non-target features (clouds) and poor quality observations flagged so that they can be excluded. The United States Geological Survey (USGS) has processed all of the Landsat 4 and 5 Thematic Mapper (TM), Landsat 7 Enhanced Thematic Mapper Plus (ETM+), Landsat 8 Operational Land Imager (OLI) and Thermal Infrared Sensor (TIRS) archive over the conterminous United States (CONUS), Alaska, and Hawaii, into Landsat ARD. The ARD are available to significantly reduce the burden of pre-processing on users of Landsat data. Provision of pre-prepared ARD is intended to make it easier for users to produce Landsat-based maps of land cover and land-cover change and other derived geophysical and biophysical products. The ARD are provided as tiled, georegistered, top of atmosphere and atmospherically corrected products defined in a common equal area projection, accompanied by spatially explicit quality assessment information, and appropriate metadata to enable further processing while retaining traceability of data provenance.

Remote Sensing↗

USGS2018

An initial validation of Landsat 5 and 7 derived surface water temperature for U.S. lakes, reservoirs, and estuaries

The United States Harmful Algal Bloom and Hypoxia Research Control Act of 2014 identified the need for forecasting and monitoring harmful algal blooms (HAB) in lakes, reservoirs, and estuaries across the nation. Temperature is a driver in HAB forecasting models that affects both HAB growth rates and toxin production. Therefore, temperature data derived from the U.S. Geological Survey Landsat 5 Thematic Mapper and Landsat 7 Enhanced Thematic Mapper Plus thermal band products were validated across 35 lakes and reservoirs, and 24 estuaries. In situ data from the Water Quality Portal (WQP) were used for validation. The WQP serves data collected by state, federal, and tribal groups. Discrete in situ temperature data included measurements at 11,910 U.S. lakes and reservoirs from 1980 through 2015. Landsat temperature measurements could include 170,240 lakes and reservoirs once an operational product is achieved. The Landsat-derived temperature mean absolute error was 1.34°C in lake pixels >180 m from land, 4.89°C at the land-water boundary, and 1.11°C in estuaries based on comparison against discrete surface in situ measurements. This is the first study to quantify Landsat resolvable U.S. lakes and reservoirs, and large-scale validation of an operational satellite provisional temperature climate data record algorithm. Due to the high performance of open water pixels, Landsat satellite data may supplement traditional in situ sampling by providing data for most U.S. lakes, reservoirs, and estuaries over consistent seasonal intervals (even with cloud cover) for an extended period of record of more than 35 years.

International Journal of Remote Sensing↗

USGS2017

A land product characterization system for comparative analysis of satellite data and products

A Land Product Characterization System (LPCS) has been developed to provide land data and products to the community of individuals interested in validating space-based land products by comparing them with similar products available from other sensors or surface-based observations. The LPCS facilitates the application of global multi-satellite and in situ data for characterization and validation of higher-level, satellite-derived, land surface products (e.g., surface reflectance, normalized difference vegetation index, and land surface temperature). The LPCS includes data search, inventory, access, and analysis functions that will permit data to be easily identified, retrieved, co-registered, and compared statistically through a single interface. The system currently includes data and products available from Landsat 4 through 8, Moderate Resolution Imaging Spectroradiometer (MODIS) Terra and Aqua, Suomi National Polar-Orbiting Partnership (S-NPP)/Joint Polar Satellite System (JPSS) Visible Infrared Imaging Radiometer Suite (VIIRS), and simulated data for the Geostationary Operational Environmental Satellite (GOES)-16 Advanced Baseline Imager (ABI). In addition to the future inclusion of in situ data, higher-level land products from the European Space Agency (ESA) Sentinel-2 and -3 series of satellites, and other high and medium resolution spatial sensors, will be included as available. When fully implemented, any of the sensor data or products included in the LPCS would be available for comparative analysis.

Remote Sensing↗

USGS2017

Mapping tide-water glacier dynamics in east Greenland using landsat data

Landsat multispectral scanner and thematic mapper images were co-registered For the Kangerdlugssuaq Fjord region in East Greenland and were used to map glacier drainage-basin areas, changes in the positions of tide-water glacier termini and to estimate surface velocities of the larger tide-water glaciers. Statistics were compiled to document distance and area changes to glacier termini. The methodologies developed in this study are broadly applicable to the investigation of tide-water glaciers in other areas. The number of images available for consecutive years and the accuracy with which images are co-registered are key factors that influence the degree to which regional glacier dynamics can be characterized using remotely sensed data. Three domains of glacier state were interpreted: net increase in terminus area in the southern part of the study area, net loss of terminus area for glaciers in upper Kangerdlugssuaq Fjord and a slight loss of glacier terminus area northward from Ryberg Fjord. Local increases in the concentrations of drifting icebergs in the fjords coincide with the observed extension of glacier termini positions Ice-surface velocity estimates were derived for several glaciers using automated image cross-correlation techniques The velocity determined for Kangerdlugssuaq Gletscher is approximately 5.0 km a −1 and that for Kong Christian IV Gletscher is 0.9 km a −1 . The continuous presence of icebergs and brash ice in front of these glaciers indicates sustained rates of ice-front calving.

Journal of Glaciology↗

USGS2017

Cloud detection algorithm comparison and validation for operational Landsat data products

Clouds are a pervasive and unavoidable issue in satellite-borne optical imagery. Accurate, well-documented, and automated cloud detection algorithms are necessary to effectively leverage large collections of remotely sensed data. The Landsat project is uniquely suited for comparative validation of cloud assessment algorithms because the modular architecture of the Landsat ground system allows for quick evaluation of new code, and because Landsat has the most comprehensive manual truth masks of any current satellite data archive. Currently, the Landsat Level-1 Product Generation System (LPGS) uses separate algorithms for determining clouds, cirrus clouds, and snow and/or ice probability on a per-pixel basis. With more bands onboard the Landsat 8 Operational Land Imager (OLI)/Thermal Infrared Sensor (TIRS) satellite, and a greater number of cloud masking algorithms, the U.S. Geological Survey (USGS) is replacing the current cloud masking workflow with a more robust algorithm that is capable of working across multiple Landsat sensors with minimal modification. Because of the inherent error from stray light and intermittent data availability of TIRS, these algorithms need to operate both with and without thermal data. In this study, we created a workflow to evaluate cloud and cloud shadow masking algorithms using cloud validation masks manually derived from both Landsat 7 Enhanced Thematic Mapper Plus (ETM +) and Landsat 8 OLI/TIRS data. We created a new validation dataset consisting of 96 Landsat 8 scenes, representing different biomes and proportions of cloud cover. We evaluated algorithm performance by overall accuracy, omission error, and commission error for both cloud and cloud shadow. We found that CFMask, C code based on the Function of Mask (Fmask) algorithm, and its confidence bands have the best overall accuracy among the many algorithms tested using our validation data. The Artificial Thermal-Automated Cloud Cover Algorithm (AT-ACCA) is the most accurate nonthermal-based algorithm. We give preference to CFMask for operational cloud and cloud shadow detection, as it is derived from a priori knowledge of physical phenomena and is operable without geographic restriction, making it useful for current and future land imaging missions without having to be retrained in a machine-learning environment.

Remote Sensing of Environment↗

USGS2017

Mapping burned areas using dense time-series of Landsat data

Complete and accurate burned area data are needed to document patterns of fires, to quantify relationships between the patterns and drivers of fire occurrence, and to assess the impacts of fires on human and natural systems. Unfortunately, in many areas existing fire occurrence datasets are known to be incomplete. Consequently, the need to systematically collect burned area information has been recognized by the United Nations Framework Convention on Climate Change and the Intergovernmental Panel on Climate Change, which have both called for the production of essential climate variables (ECVs), including information about burned area. In this paper, we present an algorithm that identifies burned areas in dense time-series of Landsat data to produce the Landsat Burned Area Essential Climate Variable (BAECV) products. The algorithm uses gradient boosted regression models to generate burn probability surfaces using band values and spectral indices from individual Landsat scenes, lagged reference conditions, and change metrics between the scene and reference predictors. Burn classifications are generated from the burn probability surfaces using pixel-level thresholding in combination with a region growing process. The algorithm can be applied anywhere Landsat and training data are available. For this study, BAECV products were generated for the conterminous United States from 1984 through 2015. These products consist of pixel-level burn probabilities for each Landsat scene, in addition to, annual composites including: the maximum burn probability and a burn classification. We compared the BAECV burn classification products to the existing Global Fire Emissions Database (GFED; 1997–2015) and Monitoring Trends in Burn Severity (MTBS; 1984–2013) data. We found that the BAECV products mapped 36% more burned area than the GFED and 116% more burned area than MTBS. Differences between the BAECV products and the GFED were especially high in the West and East where the BAECV products mapped 32% and 88% more burned area, respectively. However, the BAECV products found less burned area than the GFED in regions with frequent agricultural fires. Compared to the MTBS data, the BAECV products identified 31% more burned area in the West, 312% more in the Great Plains, and 233% more in the East. Most pixels in the MTBS data were detected by the BAECV, regardless of burn severity. The BAECV products document patterns of fire similar to those in the GFED but also showed patterns of fire that are not well characterized by the existing MTBS data. We anticipate the BAECV products will be useful to studies that seek to understand past patterns of fire occurrence, the drivers that created them, and the impacts fires have on natural and human systems.

continental United States↗

USGS2016

Selection and quality assessment of Landsat data for the North American forest dynamics forest history maps of the US

Using the NASA Earth Exchange platform, the North American Forest Dynamics (NAFD) project mapped forest history wall-to-wall, annually for the contiguous US (1986–2010) using the Vegetation Change Tracker algorithm. As with any effort to identify real changes in remotely sensed time-series, data gaps, shifts in seasonality, misregistration, inconsistent radiometry and cloud contamination can be sources of error. We discuss the NAFD image selection and processing stream (NISPS) that was designed to minimize these sources of error. The NISPS image quality assessments highlighted issues with the Landsat archive and metadata including inadequate georegistration, unreliability of the pre-2009 L5 cloud cover assessments algorithm, missing growing-season imagery and paucity of clear views. Assessment maps of Landsat 5–7 image quantities and qualities are presented that offer novel perspectives on the growing-season archive considered for this study. Over 150,000+ Landsat images were considered for the NAFD project. Optimally, one high quality cloud-free image in each year or a total of 12,152 images would be used. However, to accommodate data gaps and cloud/shadow contamination 23,338 images were needed. In 220 specific path-row image years no acceptable images were found resulting in data gaps in the annual national map products.

International Journal of Digital Earth↗

USGS2016

A note on the temporary misregistration of Landsat-8 Operational Land Imager (OLI) and Sentinel-2 Multi Spectral Instrument (MSI) imagery

The Landsat-8 and Sentinel-2 sensors provide multi-spectral image data with similar spectral and spatial characteristics that together provide improved temporal coverage globally. Both systems are designed to register Level 1 products to a reference image framework, however, the Landsat-8 framework, based upon the Global Land Survey images, contains residual geolocation errors leading to an expected sensor-to-sensor misregistration of 38 m (2σ). These misalignments vary geographically but should be stable for a given area. The Landsat framework will be readjusted for consistency with the Sentinel-2 Global Reference Image, with completion expected in 2018. In the interim, users can measure Landsat-to-Sentinel tie points to quantify the misalignment in their area of interest and if appropriate to reproject the data to better alignment.

Remote Sensing of Environment↗

USGS2012

Landsat: building a strong future

Conceived in the 1960s, the Landsat program has experienced six successful missions that have contributed to an unprecedented 39-year record of Earth Observations that capture global land conditions and dynamics. Incremental improvements in imaging capabilities continue to improve the quality of Landsat science data, while ensuring continuity over the full instrument record. Landsats 5 and 7 are still collecting imagery. The planned launch of the Landsat Data Continuity Mission in December 2012 potentially extends the Landsat record to nearly 50 years. The U.S. Geological Survey (USGS) Landsat archive contains nearly three million Landsat images. All USGS Landsat data are available at no cost via the Internet. The USGS is committed to improving the content of the historical Landsat archive though the consolidation of Landsat data held in international archives. In addition, the USGS is working on a strategy to develop higher-level Landsat geo- and biophysical datasets. Finally, Federal efforts are underway to transition Landsat into a sustained operational program within the Department of the Interior and to authorize the development of the next two satellites — Landsats 9 and 10.

Remote Sensing of Environment↗

USGS2012

Development of the Landsat Data Continuity Mission cloud-cover assessment algorithms

The upcoming launch of the Operational Land Imager (OLI) will start the next era of the Landsat program. However, the Automated Cloud-Cover Assessment (CCA) (ACCA) algorithm used on Landsat 7 requires a thermal band and is thus not suited for OLI. There will be a thermal instrument on the Landsat Data Continuity Mission (LDCM)-the Thermal Infrared Sensor-which may not be available during all OLI collections. This illustrates a need for CCA for LDCM in the absence of thermal data. To research possibilities for full-resolution OLI cloud assessment, a global data set of 207 Landsat 7 scenes with manually generated cloud masks was created. It was used to evaluate the ACCA algorithm, showing that the algorithm correctly classified 79.9% of a standard test subset of 3.95 109 pixels. The data set was also used to develop and validate two successor algorithms for use with OLI data-one derived from an off-the-shelf machine learning package and one based on ACCA but enhanced by a simple neural network. These comprehensive CCA algorithms were shown to correctly classify pixels as cloudy or clear 88.5% and 89.7% of the time, respectively.

IEEE Transactions on Geoscience and Remote Sensing↗

USGS2011

Continuity of Landsat observations: Short term considerations

As of writing in mid-2010, both Landsat-5 and -7 continue to function, with sufficient fuel to enable data collection until the launch of the Landsat Data Continuity Mission (LDCM) scheduled for December of 2012. Failure of one or both of Landsat-5 or -7 may result in a lack of Landsat data for a period of time until the 2012 launch. Although the potential risk of a component failure increases the longer the sensor's design life is exceeded, the possible gap in Landsat data acquisition is reduced with each passing day and the risk of Landsat imagery being unavailable diminishes for all except a handful of applications that are particularly data demanding. Advances in Landsat data compositing and fusion are providing opportunities to address issues associated with Landsat-7 SLC-off imagery and to mitigate a potential acquisition gap through the integration of imagery from different sensors. The latter will likely also provide short-term, regional solutions to application-specific needs for the continuity of Landsat-like observations. Our goal in this communication is not to minimize the community's concerns regarding a gap in Landsat observations, but rather to clarify how the current situation has evolved and provide an up-to-date understanding of the circumstances, implications, and mitigation options related to a potential gap in the Landsat data record.

Remote Sensing of Environment↗

USGS2010

An overview of the Landsat Data Continuity Mission

The Landsat Data Continuity Mission (LDCM) is the follow-on mission to Landsat 7 and will be the eighth mission in the Landsat series. The mission is in development via an interagency partnership between the National Aeronautics and Space Administration (NASA) and the Department of Interior (DOI) / United States Geological Survey (USGS). The LDCM satellite will carry two earth-observing sensors, the Operational Land Imager (OLI) to collect image data for nine spectral bands in the reflective portion of the spectrum and the Thermal Infrared Sensor (TIRS) to collect coincident image data for two thermal spectral bands. The LDCM ground segment will control the satellite and will receive, process, archive, and distribute the science data collected by the OLI and TIRS instruments. The USGS Earth Resources Observation & Science Center (EROS) will distribute LDCM data products at no cost to requestors. The mission objective is to continues the Landsat program's collection, archive, and distribution of multispectral imagery affording global, synoptic, and repetitive coverage of the Earth's land surfaces at a scale where natural and human-induced changes can be detected, differentiated, characterized, and monitored over time. The LDCM launch readiness date is currently December, 2012.

Conference Paper↗

USGS2008

Agreement evaluation of AVHRR and MODIS 16-day composite NDVI data sets

Satellite-derived normalized difference vegetation index (NDVI) data have been used extensively to detect and monitor vegetation conditions at regional and global levels. A combination of NDVI data sets derived from AVHRR and MODIS can be used to construct a long NDVI time series that may also be extended to VIIRS. Comparative analysis of NDVI data derived from AVHRR and MODIS is critical to understanding the data continuity through the time series. In this study, the AVHRR and MODIS 16-day composite NDVI products were compared using regression and agreement analysis methods. The analysis shows a high agreement between the AVHRR-NDVI and MODIS-NDVI observed from 2002 and 2003 for the conterminous United States, but the difference between the two data sets is appreciable. Twenty per cent of the total difference between the two data sets is due to systematic difference, with the remainder due to unsystematic difference. The systematic difference can be eliminated with a linear regression-based transformation between two data sets, and the unsystematic difference can be reduced partially by applying spatial filters to the data. We conclude that the continuity of NDVI time series from AVHRR to MODIS is satisfactory, but a linear transformation between the two sets is recommended.

International Journal of Remote Sensing↗

USGS2006

Remotely sensed data available from the US Geological Survey EROS Data Center

The Center for Earth Resources Observation Systems (EROS) is a field center of the geography discipline within the US geological survey (USGS) of the Department of the Interior. The EROS Data Center (EDC) was established in the early 1970s as the nation’s principal archive of remotely sensed data. Initially the EDC was responsible for the archive, reproduction, and distribution of black-and-white and color-infrared aerial photography acquired under numerous mapping programs conducted by various Federal agencies including the USGS, Department of Agriculture, Environmental Protection Agency, and NASA. The EDC was also designated the central archive for data acquired by the first satellite sensor designed for broad-scale earth observations in support of civilian agency needs for earth resource information. A four-band multispectral scanner (MSS) and a return-beam vidicon (RBV) camera were initially flown on the Earth Resources Technology Satellite-1, subsequently designated Landsat-1. The synoptic coverage, moderate spatial resolution, and multi-spectral view provided by these data stimulated scientists with an unprecedented perspective from which to study the Earth’s surface and to understand the relationships between human activity and natural systems.

Book chapter↗