Geology topics

Christopher E. Soulard

Publications and source records attributed to Christopher E. Soulard.

At least 19 recordsLinked to original sources

USGS2026

Satellite embeddings for crop type classification: A comparative examination

Embedding datasets encode complex relationships among multiple sources of Earth observation data into a compact format. Here, we evaluated the utility of a 10-m global Satellite Embedding product (SE) for classifying crop types in central California for the year 2020. We compared the classification accuracy of a random forest model based exclusively on the SE layer to an existing random forest model with multiple imagery inputs. Our results showed the SE-based classification had higher agreement with the reference dataset (California Department of Water Resources crop map) than the classification based on Landsat and National Agricultural Imagery Program inputs (94.7% versus. 91.9% overall accuracy, respectively). The performance of individual crop types was consistent across models, ranging from high agreement for rice (98.4% versus 98% accuracy) to lower agreement for pasture, grain, and fallow/young perennial classes (< 65% accuracy in both models). The SE-based workflow used three times less cloud-based computational resources and represented substantial savings of predictor development time. Geospatial embedding products can aid classification efforts by reducing predictor development time and processing demands while maintaining classification accuracy.

California↗

USGS2026

Multi-temporal mapping and analyses of post-wildfire surface soil moisture modeled from Landsat Thermal-IR and Sentinel-1 SAR

Wildfires alter surface soil moisture by consuming vegetation and exposing the soil to increased solar radiation, increased precipitation throughfall, and reduced infiltration from soil hydrophobicity. This can affect the water balance within the watershed and, in turn, influence post-fire nutrient cycling and ecosystem succession. Modeling surface soil moisture over large areas using satellite remote sensing can aid monitoring of soil moisture patterns and dynamics. We trained a random forest regression model using an in situ soil moisture network and examined pre- and post-fire soil moisture dynamics for three wildfires that burned in 2016 and 2017 in coastal California, U.S.A. Model inputs included Sentinel-1A synthetic aperture radar (SAR) and Landsat thermal infrared imagery (TIR), which produced a high resolution (30 m) soil moisture map time series with 18–20 estimates across 2 years. An unbiased root mean square error of 0.065 and an r 2 of 0.624 indicates a good fit between the model and reference data. We found a significant interaction between treatment (wildfire or unburned) and period (before or after wildfire) by examining the modeled surface soil moisture in high-severity burned areas. Before-after control-impact metrics showed that modeled surface soil moisture was lower after each fire within wildfire areas compared to unburned controls, and overall changes in soil moisture were greater in the wildfire perimeters than unburned areas. These results demonstrate the potential of Sentinel-1 SAR and Landsat TIR data for the development of high-resolution soil moisture time series products to support monitoring and assessment of post-fire soil hydrologic conditions.

California↗

USGS2025

Assessing causes and consequences of winter surface water dynamics in California’s Central Valley using satellite remote sensing

California's Central Valley is increasingly vulnerable to winter floods. A comprehensive spatial baseline of flood extents is critical for inundation analyses that can enhance future flood predictions, but cloud cover has prevented the regular observation of surface water extents with optical satellite imagery. In this study, we leveraged the daily resolution of Moderate Resolution Imaging Spectroradiometer (MODIS) satellite data to create a continuous series of monthly Dynamic Surface Water Extent (DSWEmod) images across the Central Valley from January 2003 to January 2023. We used the timeseries to assess the climatic driving forces of winter (Oct–April) surface water variability at sub-basin and pixel scales. At the sub-basin scale, we evaluated the influences of winter precipitation, occurrence of atmospheric rivers, and antecedent soil moisture on monthly surface water extents and found that the greatest correspondence occurs in mid-winter (Dec–Feb); in contrast, non-precipitation drivers such as water management play a stronger role in autumn and spring. The pixel-level analysis identified the probabilities of precipitation-driven surface water occurrences in the Sacramento basin are highest along rivers, conveyance channels, and floodways, with higher probabilities under wetter antecedent soil moisture conditions. Precipitation-driven surface water occurrences are also common in leveed areas and outside flood boundaries designated by state and federal agencies where exposure of structures to inundation was larger in terms of their value. Finally, areas with more frequent precipitation-driven flooding have poor recharge potential but are commonly within 5 km of areas classified as having good potential. This study demonstrates a novel approach for exploring the utility of MODIS for understanding surface water dynamics in mid-winter, a period characterized by peak precipitation, flood risk, and surface water extent. This information can provide valuable insights for (1) assessing flood risks for infrastructure and populations, (2) identifying areas most suited to strategic water management investments to increase recharge, and (3) analyzing precipitation thresholds that trigger flooding to allow proactive water management strategies to minimize damage and maximize recharge.

California↗

USGS2025

Linking fire radiative power to land cover, fire history, and environmental setting in Alaska, 2003–2022

Background Fire radiative power (FRP) shows promise as a diagnostic and predictive indicator of fire behavior and post-fire effects in Alaska, USA. Aims To investigate relationships between FRP, vegetation functional groups, and environmental settings in Alaska (2003–2022) under various fire history conditions. Methods We tested for distinctness of MODIS FRP distributions associated with vegetation classes and fire legacies (frequency and number of previous burns). We used a random forest model to examine relative importance of vegetation class for FRP versus bottom-up biophysical and temporal parameters. Key results FRP distributions are statistically distinct among vegetation functional groups with contrasting fuel biomass, or within functional groups with contrasting burn characteristics. Location and topography, which constrain vegetation class, strongly determine FRP, and fire history is of lesser importance over the 19-year analysis period. Conclusions FRP can be used to identify wildfire consumption in dissimilar vegetation classes but is highly conditioned by geographic location. The complex and evolving vegetation composition of post-fire boreal landscapes precludes a clear association of expected FRP at distinct stages of recovery. Implications These results can inform further study of FRP as an indicator of fire behavior and fuel consumption and for informing dynamics of post-fire recovery across Alaska.

Alaska↗

USGS2025

Land change, fire, and climate weaken carbon sink in the conterminous U.S.

The land carbon sink of the conterminous United States was evaluated using a bottom-up modeling framework and 30-meter land change data from 1985 to 2020. This cross-scale, cross-landscape, and cross-system approach tracked fractional land cover changes and applied regional model calibration. Results show average terrestrial and aquatic carbon sinks of +110 ± 37 and +19 ± 0.5 teragrams of carbon per year, respectively. The terrestrial carbon sink, showing no clear trend, peaked in the 1990s, with more years as a carbon source since 2000, contradicting recent national and global studies. Land change had the largest impact (−70 ± 5.5 teragrams of carbon per year), exceeding impacts of climate (−33 ± 48 teragrams of carbon per year), wildfire (−7.7 ± 2.4 teragrams of carbon per year), and erosion transport (−1.9 ± 0.13 teragrams of carbon per year). The positive CO2 fertilization effect (+69 ± 12 teragrams of carbon per year) was insufficient to maintain the carbon sink strength. Our framework reveals key paths of carbon loss, with implications for carbon budget and energy policies in the United States and beyond.

conterminous United States↗

USGS2024

The feasibility of using national-scale datasets for classifying wetlands in Arizona with machine learning

The advent of machine learning techniques has led to a proliferation of landscape classification products. These approaches can fill gaps in wetland inventories across the United States (U.S.) provided that large reference datasets are available to develop accurate models. In this study, we tested the feasibility of expediting the classification process by sourcing requisite training and testing data from existing national-scale land cover maps instead of customized sample sets. We created a single map of water and wetland presence by intersecting water and wetland classes from available land cover products (National Wetland Inventory, Gap Analysis Project, National Land Cover Database and Dynamic Surface Water Extent) across the U.S. state of Arizona, which has fewer wetland-specific mapping products than other parts of the U.S. We derived classified samples for four wetland classes from the combined map: open water, herbaceous wetlands, wooded wetlands and non-wetland cover. In Google Earth Engine, we developed a random forest model that combined the training data with spatial predictor variables, including vegetation greenness indices, wetness indices, seasonal index variation, topographic parameters and vegetation height metrics. Results show that the final model separates the four classes with an overall accuracy of 86.2%. The accuracy suggests that existing datasets can be effectively used to compile machine learning training samples to map wetlands in arid landscapes in the U.S. These methods hold promise for the generation of wetland inventories at more frequent intervals, which could allow more nuanced investigations of wetland change over time in response to anthropogenic and climatic drivers.

Arizona↗

USGS2024

Crop type classification, trends, and patterns of central California agricultural fields from 2005 to 2020

California produces many key agricultural products in the United States. Current geospatial agricultural datasets are limited in mapping accuracy, spatial context, or observation period. This study uses machine learning and high-resolution imagery to produce a time series of crop maps to assess crop type trends and patterns across central California from 2005 to 2020. National Agriculture Imagery Program and Landsat imagery were used to classify nine crop types that are common in the study region: grain crops, field crops, rice, citrus and subtropical, deciduous fruit and nut, vineyard, berry and vegetable, pasture, and fallow/young perennial crop types. To create labeled data, we sampled 1253 fields and manually identified crop types for each examined year using high-resolution imagery and Landsat normalized difference vegetation index time series. We applied a random forest machine learning algorithm in Google Earth Engine. Results show that the mean overall classification accuracy of the nine-class map was 93.1%, with individual accuracies ranging from 99.3% (rice) to 89.5% (fallow/young perennial). Mann–Kendall trend tests showed significant ( p less than 0.05) declines in field crop and pasture area during the study period, while deciduous fruit and nut, citrus and subtropical, and fallow/young perennial crop types experienced significant increases. At an aggregate level, there was a general shift from annual crop types to perennial crop types. These data provide a 16-year time span of spatially explicit crop type classifications, trends, and patterns in central California that can be used to aid managers and decision makers for resource planning or hazard mitigation.

California↗

USGS2022

Analysis of surface water trends for the conterminous United States using MODIS satellite data, 2003–2019

Satellite imagery is commonly used to map surface water extents over time, but many approaches yield discontinuous records resulting from cloud obstruction or image archive gaps. We applied the Dynamic Surface Water Extent (DSWE) model to downscaled (250-m) daily Moderate Resolution Imaging Spectroradiometer (MODIS) data in Google Earth Engine to generate monthly surface water maps for the conterminous United States (US) from 2003 through 2019. The aggregation of daily observations to monthly maps of maximum water extent produced records with diminished cloud and cloud shadow effects across most of the country. We used the continuous monthly record to analyze spatiotemporal surface water trends stratified within Environmental Protection Agency Ecoregions. Although not all ecoregion trends were significant ( p < 0.05), results indicate that much of the western and eastern US underwent a decline in surface water over the 17-year period, while many ecoregions in the Great Plains had positive trends. Trends were also generated from monthly streamgage discharge records and compared to surface water trends from the same ecoregion. These approaches agreed on the directionality of trend detected for 54 of 85 ecoregions, particularly across the Great Plains and portions of the western US, whereas trends were not congruent in select western deserts, the Great Lakes region, and the southeastern US. By describing the geographic distribution of surface water over time and comparing these records to instrumented discharge data across the conterminous US, our findings demonstrate the efficacy of using satellite imagery to monitor surface water dynamics and supplement traditional instrumented monitoring.

Water Resources Research↗

USGS2022

DSWEmod - The production of high-frequency surface water map composites from daily MODIS images

Optical satellite imagery is commonly used for monitoring surface water dynamics, but clouds and cloud shadows present challenges in assembling complete water time series. To test whether the daily revisit rate of Moderate Resolution Imaging Spectroradiometer (MODIS) satellite imagery can reduce cloud obstruction and improve high-frequency surface water mapping, we compared map results derived from Landsat (30-m) and MODIS (250-m) data across the state of California for 2003–2019. We adapted the Dynamic Surface Water Extent (DSWE) model in Google Earth Engine to generate surface water map composites from MODIS imagery every 5, 10, 15, and 30 days, and compared products to monthly Landsat-based DSWE maps. Results for DSWEmod (DSWE MODIS) in California suggest that more than 5% data loss (cloud obstruction, etc.) was present in only 2% of the 15-day time series, as compared to 32% of the monthly Landsat DSWE time series. The five-day DSWEmod composites averaged 8.4% obscuration in the winter months. Area estimates derived from cloud-filtered MODIS and Landsat monthly products have the highest linear correlations compared to streamgage discharge records, suggesting that monthly scale analyses best explain the relationship between surface water area and general streamflow dynamics. Shorter-interval DSWEmod products have lower correlations but utility for understanding the timing of surface water peaks and past flood events.

Journal of the American Water Resources Associatio↗

USGS2022

Using Landsat and MODIS satellite collections to examine extent, timing, and potential impacts of surface water inundation in California croplands☆

The state of California, United States of America produces many crop products that are both utilized domestically and exported throughout the world. With nearly 39,000 km 2 of croplands, monitoring unintentional and intentional surface water inundation is important for water resource management and flood hazard readiness. We examine inundation dynamics in California croplands from 2003 to 2020 by intersecting monthly surface water maps (n = 216 months) derived using two satellite remote sensing platforms (Landsat and Moderate Resolution Imaging Spectroradiometer [MODIS]) with a high-quality cropland map generated by the California Department of Water Resources. Surface water maps were produced using the Dynamic Surface Water Extent model, in which satellite image pixels are classified into different levels of detection confidence. Our analysis focused on calculating monthly and annual occurrence of “high confidence” water for each satellite collection across eight cropland types and 58 counties. Results indicate that 49.9% (MODIS) to 56.4% (Landsat) of croplands were inundated at least once during the 18-year timespan. Rice crops, due to their unique need of consistent surface water and dominance as a crop type in CA, had the highest proportion of and mean annual inundation area, while citrus crops had the lowest. Mean monthly inundation patterns in most croplands followed California's precipitation patterns with high inundation during the winter and spring rainy season. At the county level, croplands in the southern Central Valley typically had high occurrences of inundation in conjunction with large crop areas. Exposure and sensitivity of inundation for three crop types (citrus, deciduous, and vineyards) that are typically less associated with intentional inundation were geographically variable, but overall were generally highest in counties in the southern Central Valley, California's primary agricultural region. Flood and precipitation related crop insurance claims indicated that rice had the highest mean indemnity payout for any month with damages typically occurring in March and April. Insurance claims were also high in deciduous fruit and nut crops, which had peak damages in February. A comparison between inundation results and insurance claims suggests that the inundation mapped by our process coincides with claim activity. These data elucidate water inundation patterns across the state that can serve to inform farmers, insurers, decision makers, resource managers, and flood mitigation professionals.

California↗

USGS2020

Application of empirical land-cover changes to construct climate change scenarios in federally managed lands

Sagebrush-dominant ecosystems in the western United States are highly vulnerable to climatic variability. To understand how these ecosystems will respond under potential future conditions, we correlated changes in National Land Cover Dataset “Back-in-Time” fractional cover maps from 1985-2018 with Daymet climate data in three federally managed preserves in the sagebrush steppe ecosystem: Beaty Butte Herd Management Area, Hart Mountain National Antelope Refuge, and Sheldon National Wildlife Refuge. Future (2018 to 2050) abundance and distribution of vegetation cover were modeled at a 300-m resolution under a business-as-usual climate (BAU) scenario and a Representative Concentration Pathway (RCP) 8.5 climate change scenario. Spatially explicit map projections suggest that climate influences may make the landscape more homogeneous in the near future. Specifically, projections indicate that pixels with high bare ground cover become less bare ground dominant, pixels with moderate herbaceous cover contain less herbaceous cover, and pixels with low shrub cover contain more shrub cover. General vegetation patterns and composition do not differ dramatically between scenarios despite RCP 8.5 projections of + 1.2 °C mean annual minimum temperatures and +7.6 mm total annual precipitation. Hart Mountain National Antelope Refuge is forecast to undergo the most change, with both models projecting larger declines in bare ground and larger increases in average herbaceous and shrub cover compared to Beaty Butte Herd Management Area and Sheldon National Wildlife Refuge. These scenarios present plausible future outcomes intended to guide federal land managers to identify vegetation cover changes that may affect habitat condition and availability for species of interest.

Remote Sensing↗

USGS2020

Implementation of a surface water extent model in Cambodia using cloud-based remote sensing

Mapping surface water over time provides the spatially explicit information essential for hydroclimatic research focused on droughts and flooding. Hazard risk assessments and water management planning also rely on accurate, long-term measurements describing hydrologic fluctuations. Stream gages are a common measurement tool used to better understand flow and inundation dynamics, but gage networks are incomplete or non-existent in many parts of the world. In such instances, satellite imagery may provide the only data available to monitor surface water changes over time. Here, we describe an effort to extend the applicability of the USGS Dynamic Surface Water Extent (DSWE) model to non-US regions. We leverage the multi-decadal archive of the Landsat satellite in the Google Earth Engine (GEE) cloud-based computing platform to produce and analyze 372 monthly composite maps and 31 annual maps (January 1988–December 2018) in Cambodia, a flood-prone country in Southeast Asia that lacks a comprehensive stream gage network. DSWE relies on a series of spectral water indices and elevation data to classify water into four categories of water inundation. We compared model outputs to existing surface water maps and independently assessed DSWE accuracy at discrete dates across the time series. Despite considerable cloud obstruction and missing imagery across the monthly time series, the overall accuracy exceeded 85% for all annual tests. The DSWE model consistently mapped open water with high accuracy, and areas classified as “high confidence” water correlate well to other available maps at the country scale. Results in Cambodia suggest that extending DSWE globally using a cloud computing framework may benefit scientists, managers, and planners in a wide array of applications across the globe.

Remote Sensing↗

USGS2019

Phenology patterns indicate recovery trajectories of ponderosa pine forests after high-severity fires

Post-fire recovery trajectories in ponderosa pine (Pinus ponderosa Laws.) forests of the US Southwest are increasingly shifting away from pre-burn vegetation communities. This study investigated whether phenological metrics derived from a multi-decade remotely sensed imagery time-series could differentiate among grass, evergreen shrub, deciduous, or conifer-dominated replacement pathways. We focused on 10 fires that burned ponderosa pine forests in Arizona and New Mexico, USA before the year 2000. A total of 29 sites with discernable post-fire recovery signals were selected within high-severity burn areas. At each site, we used Google Earth Engine to derive time-series of normalized difference vegetation index (NDVI) signals from Landsat Thematic Mapper, Enhanced Thematic Mapper+, and Operational Line Imager data from 1984 to 2017. We aggregated values to 8- and 16-day intervals, fit Savitzy-Golay filters to each sequence, and extracted annual phenology metrics of amplitude, base value, peak value, and timing of peak value in the Timesat analysis package. Results show that relative to post-fire conditions, pre-burn ponderosa pine forests exhibit significantly lower mean NDVI amplitude (0.14 vs. 0.21), higher mean base NDVI (0.47 vs. 0.22), higher mean peak NDVI (0.60 vs. 0.43), and later mean peak NDVI (day of year 277 vs. 237). Vegetation succession exhibits distinct phenometric characteristics as early as year five (amplitude) and as late as year 20 (timing of peak NDVI). This study confirms the feasibility of leveraging phenology metrics derived from long-term imagery time series to identify and monitor ecological outcomes. This information may be of benefit to land resource managers who seek indicators of future landscape composition to inform management strategies.

Arizona, New Mexico↗

USGS2019

Landsat time series assessment of invasive annual grasses following energy development

Invasive annual grasses are of concern in many areas of the Western United States because they tolerate resource variability and have high reproductive capacity, with propagules that are readily dispersed in disturbed areas like those created and maintained for energy development. Early-season invasive grasses “green up” earlier than the most native plants, producing a distinct pulse of greenness in the early spring that can be exploited to identify their location using multi-date imagery. To determine if invasive annual grasses increased around energy development areas after the construction phase, we calculated an invasive index using Landsat TM and ETM+ imagery for a 34-year time period (1985–2018) and assessed trends for 1755 wind turbines installed between 1988 and 2013 in the Southern California Desert. The index uses the maximum normalized difference vegetation index (NDVI) for early-season greenness (January–June) and mean NDVI (July–October) for the later dry season. We estimated the relative cover of invasive annuals each year at turbine locations and control sites, and tested for changes before and after each turbine was installed. The time series was also mapped across the region and temporal trends were assessed relative to seasonal precipitation. The results showed an increase in early-season invasives at turbine sites after installation, but also an increase in many of the surrounding control areas. Maps of the invasive index show a region-wide increase starting at around 1998, and a great deal of the increase occurred in areas surrounding wind development sites. These results suggest that invasions around the energy developments occurred within the context of a larger regional invasion, and while the development did not necessarily initiate the invasion, annual grasses were more prevalent around the development areas.

California↗

USGS2019

Spatial patterns of meadow sensitivities to interannual climate variability in the Sierra Nevada

Conservation of montane meadows is a high priority for land and water managers given their critical role in buffering the effects of climate variability and their vulnerability to increasing temperatures and evaporative demands. Recent advances in cloud computing have provided new opportunities to examine ecological responses to climate variability over the past few decades, and at large spatial scales. In this study we characterized the sensitivities (magnitude and direction of the slope) of meadow vegetation responses to interannual variations in climate. We calculated sensitivity as the regression slope between a 35-year (1985-2016) time series of Landsat-derived vegetation indices characterizing late-season vegetation vigor and water balance variables from the Basin Characterization Model. We identified April 1 snowpack as the climate variable the majority of meadows were most sensitive to. We assessed how vegetation sensitivities to snowpack varied with hydrogeomorphic context (e.g., climate, geology, soils, watershed geometry and land cover) across the Sierra Nevada mountain range using factor analysis to reduce the dimensionality of the hydrogeomorphic data, and multiple linear regression to model sensitivity responses. We found that meadow sensitivities to snowpack varied with long-term average meadow climate, indicators of watershed subsurface water storage capacity, and indicators of meadow vegetation composition. Alpine and sub-alpine meadows with high average annual precipitation, but limited catchment subsurface storage exhibited the largest sensitivities. Our results provide a novel regional perspective on spatial patterns of meadow sensitivities to climate variability and the landscape-scale hydrogeomorphic factors that influence late-season water availability in meadow ecosystems in the Sierra Nevada.

Sierra Nevada↗

USGS2019

Integrating stream gage data and Landsat imagery to complete time-series of surface water extents in Central Valley, California

Accurate monitoring of surface water location and extent is critical for the management of diverse water resource phenomena. The multi-decadal archive of Landsat satellite imagery is punctuated by missing data due to cloud cover during acquisition times, hindering the assembly of a continuous time series of inundation dynamics. This study investigated whether streamflow volume measurements could be integrated with satellite data to fill gaps in monthly surface water chronologies for the Central Valley region of California, USA, from 1984 to 2015. We aggregated measurements of maximum monthly water extent within each of the study area’s 50 8-digit hydrologic unit code [HUC] watersheds from two Landsat-derived datasets: the European Commission’s Joint Research Centre (JRC) Monthly Water History and the U.S. Geological Survey Dynamic Surface Water Extent (DSWE). We calculated Spearman rank correlation coefficients between water extent values in each HUC and streamflow discharge data. Linear regression fits of the water extent/streamflow data pairs with the highest correlations served as the basis for interpolation of missing imagery surface water values on a HUC-wise basis. Results show strong (ρ > 0.7) maximum correlations in 11 (22.4%) and 25 (51.0%) HUCs for the DSWE and JRC time series, respectively, when comparisons were restricted to imagery and gages co-located in each HUC. Strong maximum correlations occurred in 39 (79.6%; DSWE) and 42 (85.7%; JRC) HUCs when imagery was paired with discharge data from any study area gage, providing a solid basis for reconstruction of water extent values. We generated continuous time series of 30+ years in 35 HUCs, demonstrating that this technique can provide quantitative estimates of historical surface water extents and elucidate flooding or drought events over the period of data collection. Results of a non-parametric trend analysis of the long-term time series on an annual, seasonal, and monthly basis varied among HUCs, though most trends indicate an increase in surface water over the past 30 years.

California↗

USGS2018

Removing rural roads from the National Land Cover Database to create improved urban maps for the United States, 1992-2011

Quantifying change in urban land provides important information to create empirical models examining the effects of human land use. Maps of developed land from the National Land Cover Database ( NLCD ) of the conterminous United States include rural roads in the developed land class and therefore overestimate the amount of urban land. To better map the urban class and understand how urban lands change over time, we removed rural roads and small patches of rural development from the NLCD developed class and created four wall-to-wall maps (1992, 2001, 2006, and 2011) of urban land. Removing rural roads from the NLCD developed class involved a multi-step filtering process, data fusion using geospatial road and developed land data, and manual editing. Reference data classified as urban or not urban from a stratified random sample was used to assess the accuracy of the 2001 and 2006 urban and NLCD maps. The newly created urban maps had higher overall accuracy (98.7 percent) than the NLCD maps (96.2 percent). More importantly, the urban maps resulted in lower commission error of the urban class (23 percent versus 57 percent for the NLCD in 2006) with the trade-off of slightly inflated omission error (20 percent for the urban map, 16 percent for NLCD in 2006). The removal of approximately 230,000 km 2 of rural roads from the NLCD developed class resulted in maps that better characterize the urban footprint. These urban maps are more suited to modeling applications and policy decisions that rely on quantitative and spatially explicit information regarding urban lands.

Photogrammetric Engineering and Remote Sensing↗

USGS2017

Forest harvest patterns on private lands in the Cascade Mountains, Washington, USA

Forests in Washington State generate substantial economic revenue from commercial timber harvesting on private lands. To investigate the rates, causes, and spatial and temporal patterns of forest harvest on private tracts throughout the Cascade Mountains, we relied on a new generation of annual land-use/land-cover (LULC) products created from the application of the Continuous Change Detection and Classification (CCDC) algorithm to Landsat satellite imagery collected from 1985 to 2014. We calculated metrics of landscape pattern using patches of intact and harvested forest in each annual layer to identify changes throughout the time series. Patch dynamics revealed four distinct eras of logging trends that align with prevailing regulations and economic conditions. We used multiple logistic regression to determine the biophysical and anthropogenic factors that influence fine-scale selection of harvest stands in each time period. Results show that private lands forest cover became significantly reduced and more fragmented from 1985 to 2014. Variables linked to parameters of site conditions, location, climate, and vegetation greenness consistently distinguished harvest selection for each distinct era. This study demonstrates the utility of annual LULC data for investigating the underlying factors that influence land cover change.

Washington↗