Estimation of high-resolution terrestrial evapotranspiration from Landsat data using a simple Taylor skill fusion method
Not Available
Yao, Yunjun (ORCID:0000000338038170)↗
SEARCH · Engineering Papers
Results for “landsat”
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Mapping annual irrigation from Landsat imagery and environmental variables across the conterminous United States
Not provided.
Environmental Sciences & Ecology↗
Landsat near-infrared (NIR) band and ELM-FATES sensitivity to forest disturbances and regrowth in the Central Amazon
Forest disturbance and regrowth are key processes in forest dynamics, but detailed information on these processes is difficult to obtain in remote forests such as the Amazon. We used chronosequences of Landsat satellite imagery (Landsat 5 Thematic Mapper and Landsat 7 Enhanced Thematic Mapper Plus) to determine the sensitivity of surface reflectance from all spectral bands to windthrow, clear-cut, and clear-cut and burned (cut + burn) and their successional pathways of forest regrowth in the Central Amazon. We also assessed whether the forest demography model Functionally Assembled Terrestrial Ecosystem Simulator (FATES) implemented in the Energy Exascale Earth System Model (E3SM) Land Model (ELM), ELM-FATES, accurately represents the changes for windthrow and clear-cut. The results show that all spectral bands from the Landsat satellites were sensitive to the disturbances but after 3 to 6 years only the near-infrared (NIR) band had significant changes associated with the successional pathways of forest regrowth for all the disturbances considered. In general, the NIR values decreased immediately after disturbance, increased to maximum values with the establishment of pioneers and early successional tree species, and then decreased slowly and almost linearly to pre-disturbance conditions with the dynamics of forest succession. Statistical methods predict that NIR values will return to pre-disturbance values in about 39, 36, and 56 years for windthrow, clear-cut, and cut + burn disturbances, respectively. The NIR band captured the observed, and different, successional pathways of forest regrowth after windthrow, clear-cut, and cut + burn. Consistent with inferences from the NIR observations, ELM-FATES predicted higher peaks of biomass and stem density after clear-cuts than after windthrows. ELM-FATES also predicted recovery of forest structure and canopy coverage back to pre-disturbance conditions in 38 years after windthrows and 41 years after clear-cut. The similarity of ELM-FATES predictions of regrowth patterns after windthrow and clear-cut to those of the NIR results suggests the NIR band can be used to benchmark forest regrowth in ecosystem models. Our results show the potential of Landsat imagery data for mapping forest regrowth from different types of disturbances, benchmarking, and the improvement of forest regrowth models.
54 ENVIRONMENTAL SCIENCES↗
Linear downscaling from MODIS to landsat: connecting landscape composition with ecosystem functions
The open and free access to Landsat and MODIS products have greatly promoted scientific investigations on spatiotemporal change in land mosaics and ecosystem functions at landscape to regional scales. Unfortunately, there is a major mismatch in spatial resolution between MODIS products at coarser resolution (≥ 250 m) and landscape structure based on classified Landsat scenes at finer resolution (30 m). Based on practical needs for downscaling popular MODIS products at 500 m resolution to match classified land cover at Landsat 30 m resolution, we proposed an innovative modelling approach so that landscape structure and ecosystem functions can be directly studied for their interconnections. As a proof-of-concept of our downscaling approach, we selected the watershed of the Kalamazoo River in southwestern Michigan, USA as the testbed. MODIS products for three fundamental variables of ecosystem function are downscaled to ensure the approach can be extrapolated to multiple functional measurements. They are blue-sky albedo (0–1), evapotranspiration (ET, mm), and gross primary production (GPP, Mg C ha -1 year -1 ). An object-oriented classification of Landsat images in 2011 was processed to generate a land cover map for landscape structure. The downscaling model was tested for the five Level IV ecoregions within the watershed. We achieved satisfactory downscaling models for albedo, ET, and GPP for all five ecoregions. The adjusted R2 was > 0.995 for albedo, 0.915–0.997 for ET, and 0.902–0.962 for GPP. The estimated albedo, ET, and GPP values appear different in the region. The estimated albedo was the lowest for water (0.076–0.107) and the highest for cropland (0.166–0.172). Estimated ET was the highest for the built-up cover type (525.6–687.1 mm) and the lowest for forest (209.7–459.7 mm). The estimated GPP was the highest for the build-up cover type (8.65–9.85 Mg C ha -1 year -1 ) and the lowest for forest. Estimated values for albedo, ET, and GPP appear reasonable for their ranges in the Kalamazoo River region and are consistent with values reported in the literature. Despite these promising results, the downscaling approach relies on strong assumptions and can carry substantial uncertainty. Finally, it is only valid at a spatial scale where similar climate, soil, and landforms exist (i.e., values in isolated patches of the same cover type are similar). Plausibly, the uncertainties associated with each estimation, as well as the model residuals, can be explored for other pattern-process relationships within the landscape.
54 ENVIRONMENTAL SCIENCES↗
Performance of Landsat-8 and Sentinel-2 surface reflectance products for river remote sensing retrievals of chlorophyll-a and turbidity
Rivers and other freshwater systems play a crucial role in ecosystems, industry, transportation and agriculture. Despite the > 40 years of inland water observations made possible by optical remote sensing, a standardized reflectance product for inland waters is yet forthcoming. The aim of this work is to compare the standard USGS land surface reflectance product to two Landsat-8 and Sentinel-2 aquatic remote sensing reflectance products over the Amazon, Columbia and Mississippi rivers. Landsat-8 reflectance products from all three routines are then evaluated for their comparative performance in retrieving chlorophyll-a and turbidity in reference to ship- borne, underway in situ validation measurements. The land surface product shows the best agreement (4% Mean Absolute Percent Difference) with field measurements of radiometry collected on the Amazon River and gen- erates 36% higher reflectance values in the visible bands compared to aquatic methods (ACOLITE and SeaDAS) with larger differences between land and aquatic products observed in Sentinel-2 (0.01sr-1) compared to Landsat-8 (0.001 sr-1). Choice of atmospheric correction routine can bias Landsat-8 retrievals of chlorophyll-a and turbidity by as much as 59% and 35% respectively. Using a more restrictive time window for matching in situ and satellite imagery can reduce differences by 5–31% depending on correction technique. This work highlights the challenges of satellite retrievals over rivers and underscores the need for future optical and biogeochemical research aimed at improving our understanding of the absorbing and scattering properties of river water and their relationships to remote sensing reflectance.
Amazon↗
Characterizing spatiotemporal dynamics in phenology of urban ecosystems based on Landsat data
Seasonal phenology of vegetation plays an important role in global carbon cycle and ecosystem productivity. In urban environments, vegetation phenology is also important because of its influence on public health (e.g., allergies), and energy demand (e.g. cooling effects). In this study, we studied the potential use of remotely sensed observations (i.e. Landsat data) to derive some phenology indicators for vegetation embedded within the urban core domains in four distinctly different U.S. regions (Washington, D.C., King County in Washington, Polk County in Iowa, and Baltimore City and County in Maryland) during the past three decades. We used all available Landsat observations (circa 3000 scenes) from 1982 to 2015 and a self-adjusting double logistic model to detect and quantify the annual change of vegetation phenophases, i.e. indicators of seasonal changes in vegetation. The proposed model can capture and quantify not only phenophases of dense vegetation in rural areas, but also those of mixed vegetation in urban core domains. The derived phenology indicators show a good agreement with similar indicators derived from the Moderate Resolution Imaging Spectroradiometer (MODIS) and in situ observations, suggesting that the phenology dynamic depicted by the proposed model is reliable. The vegetation phenology and its seasonal and interannual dynamics demonstrate a distinct spatial pattern in urban domains with an earlier (9–14 days) start-of season (SOS) and a later (13–20 days) end-of season (EOS), resulting in an extended (5–30 days) growing season length (GSL) when compared to the surrounding suburban and rural areas in the four study regions. There is a general long-term trend of decreasing SOS (-0.30 day per year), and increasing EOS and GSL (0.50 and 0.90 day per year, respectively) over past three decades for these study regions. The magnitude of these trends varies among the four urban systems due to their diverse local climate conditions, vegetation types, and different urban-rural settings. The Landsat derived phenology information for urban domains provides more details when compared to the coarse-resolution datasets such as MODIS, thus improves our understanding of human-natural systems interactions (or feedbacks) in urban domains. Such information is very valuable for urban planning in light of rapid urbanization and expansion of major metropolitans at the national and global levels.
60 APPLIED LIFE SCIENCES↗
RiverPIXELS: paired Landsat images and expert-labeled sediment and water pixels for a selection of rivers v1.0
RiverPIXELS contains GeoTIFFs of hand-labeled water and sediment pixels from Landsat images containing rivers. Each of the 104 labeled patches contains 256 x 256 Landsat pixels (30 meter resolution). Our aim in releasing RiverPIXELS is to provide an "off-the-shelf" training and testing dataset for building machine-learned models to automatically identify rivers from multispectral imagery. While a number of trained models and/or surface water products already exist, RiverPIXELS aims for pixel-level accuracy in order to precisely identify river boundaries in particular. Our selection of rivers is heavily Arctic, but we include tropical and temperate rivers as well. Patches are provided for the Colville (7), Indigirka (6), Kolyma (4), Ucayali (54), Waitaki (21), and Yana (12) Rivers. For each patch, all surface water pixels are labeled (1) and all in-channel sediment pixels are labeled (2). Sediments not in-channel are considered part of the land (0) class. RiverPIXELS also includes paired surface water data from the Global Surface Water dataset that may be useful as additional features in machine learning models. Each patch therefore contains four aligned GeoTIFFs: labeled, landsat, gswmo, and gswocc.
54 ENVIRONMENTAL SCIENCES↗
Linking repeat lidar with Landsat products for large scale quantification of fire-induced permafrost thaw settlement in interior Alaska
The permafrost–fire–climate system has been a hotspot in research for decades under a warming climate scenario. Surface vegetation plays a dominant role in protecting permafrost from summer warmth, thus, any alteration of vegetation structure, particularly following severe wildfires, can cause dramatic top–down thaw. A challenge in understanding this is to quantify fire-induced thaw settlement at large scales (>1000 km 2 ). In this study, we explored the potential of using Landsat products for a large-scale estimation of fire-induced thaw settlement across a well-studied area representative of ice-rich lowland permafrost in interior Alaska. Six large fires have affected ~1250 km 2 of the area since 2000. We first identified the linkage of fires, burn severity, and land cover response, and then developed an object-based machine learning ensemble approach to estimate fire-induced thaw settlement by relating airborne repeat lidar data to Landsat products. The model delineated thaw settlement patterns across the six fire scars and explained ~65% of the variance in lidar-detected elevation change. Our results indicate a combined application of airborne repeat lidar and Landsat products is a valuable tool for large scale quantification of fire-induced thaw settlement.
54 ENVIRONMENTAL SCIENCES↗
Reference Shapefiles and Pre-trained Random Forest Classification Models for Detecting Aufeis on the North Slope of Alaska in Landsat Imagery
This dataset provides shapefiles and trained machine learning models used for aufeis detection at four sites on the North Slope of Alaska. It includes reference data for evaluating Landsat-based detection methods, supporting research on remote sensing approaches for identifying aufeis. The ReferenceData folder contains ArcGIS shapefiles of semi-automated land cover classifications for 217 Landsat Collection 2 images, categorizing pixels into six classes: aufeis, snow, ground, none, water, and cloud. The SiteBuffers.zip file includes 10-kilometer buffer shapefiles defining regions of interest around four aufeis fields (Canning21, FH1, Firth, and Kuparuk), used to test three detection techniques. Additionally, the TrainedRFModels folder contains six pre-trained Scikit-Learn Random Forest classifiers (100 trees, max depth = 30) designed to predict aufeis presence in Landsat Collection 2 Surface Reflectance images using Red, Blue, SWIR2, NDVI, and NDWI bands. This dataset supports the development and validation of remote sensing methods for mapping aufeis in Arctic environments.The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research.The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska.Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).
54 ENVIRONMENTAL SCIENCES↗
STAIR 2.0: A Generic and Automatic Algorithm to Fuse Modis, Landsat, and Sentinel-2 to Generate 10 m, Daily, and Cloud-/Gap-Free Surface Reflectance Product
Remote sensing datasets with both high spatial and high temporal resolution are critical for monitoring and modeling the dynamics of land surfaces. However, no current satellite sensor could simultaneously achieve both high spatial resolution and high revisiting frequency. Therefore, the integration of different sources of satellite data to produce a fusion product has become a popular solution to address this challenge. Many methods have been proposed to generate synthetic images with rich spatial details and high temporal frequency by combining two types of satellite datasets—usually frequent coarse-resolution images (e.g., MODIS) and sparse fine-resolution images (e.g., Landsat). In this paper, we introduce STAIR 2.0, a new fusion method that extends the previous STAIR fusion framework, to fuse three types of satellite datasets, including MODIS, Landsat, and Sentinel-2. In STAIR 2.0, input images are first processed to impute missing-value pixels that are due to clouds or sensor mechanical issues using a gap-filling algorithm. The multiple refined time series are then integrated stepwisely, from coarse- to fine- and high-resolution, ultimately providing a synthetic daily, high-resolution surface reflectance observations. We applied STAIR 2.0 to generate a 10-m, daily, cloud-/gap-free time series that covers the 2017 growing season of Saunders County, Nebraska. Moreover, the framework is generic and can be extended to integrate more types of satellite data sources, further improving the quality of the fusion product. View Full-Text
47 OTHER INSTRUMENTATION↗
Automated detection of cloud and cloud-shadow in single-date Landsat imagery using neural networks and spatial post-processing
Use of Landsat data to answer ecological questions is contingent on the effective removal of cloud and cloud shadow from satellite images. We develop a novel algorithm to identify and classify clouds and cloud shadow, \textsc{sparcs}: Spacial Procedures for Automated Removal of Cloud and Shadow. The method uses neural networks to determine cloud, cloud-shadow, water, snow/ice, and clear-sky membership of each pixel in a Landsat scene, and then applies a set of procedures to enforce spatial rules. In a comparison to FMask, a high-quality cloud and cloud-shadow classification algorithm currently available, \textsc{sparcs} performs favorably, with similar omission errors for clouds (0.8% and 0.9%, respectively), substantially lower omission error for cloud-shadow (8.3% and 1.1%), and fewer errors of commission (7.8% and 5.0%). Additionally, textsc{sparcs} provides a measure of uncertainty in its classification that can be exploited by other processes that use the cloud and cloud-shadow detection. To illustrate this, we present an application that constructs obstruction-free composites of images acquired on different dates in support of algorithms detecting vegetation change.
cloud masking↗
A Dynamic Landsat Derived Normalized Difference Vegetation Index (NDVI) Product for the Conterminous United States
Satellite derived vegetation indices (VIs) are broadly used in ecological research, ecosystem modeling, and land surface monitoring. The Normalized Difference Vegetation Index (NDVI), perhaps the most utilized VI, has countless applications across ecology, forestry, agriculture, wildlife, biodiversity, and other disciplines. Calculating satellite derived NDVI is not always straight-forward, however, as satellite remote sensing datasets are inherently noisy due to cloud and atmospheric contamination, data processing failures, and instrument malfunction. Readily available NDVI products that account for these complexities are generally at coarse resolution; high resolution NDVI datasets are not conveniently accessible and developing them often presents numerous technical and methodological challenges. Here, we address this deficiency by producing a Landsat derived, high resolution (30 m), long-term (30+ years) NDVI dataset for the conterminous United States. We use Google Earth Engine, a planetary-scale cloud-based geospatial analysis platform, for processing the Landsat data and distributing the final dataset. We use a climatology driven approach to fill missing data and validate the dataset with established remote sensing products at multiple scales. We provide access to the composites through a simple web application, allowing users to customize key parameters appropriate for their application, question, and region of interest.
42 ENGINEERING↗
Seasonally-managed wetland footprint delineation using Landsat ETM+ satellite imagery
One major challenge in water resource management is the estimation of evapotranspiration losses from seasonally managed wetlands. Quantifying these losses is complicated by the dynamic nature of the wetlands' areal footprint during the periods of flood-up and drawdown. In this paper, we present a data-lean solution to this problem using an example application in the San Joaquin Basin, California. Through analysis of high-resolution Landsat Enhanced Thematic Mapper Plus (ETM+) satellite imagery, we develop a metric to better capture the extent of total flooded wetland area. The procedure is validated using year-long, continuously-logged field datasets for two wetlands within the study area. The proposed classification which uses a Landsat ETM + Band 5 (mid-IR wavelength) to Band 2 (visible green wavelength) ratio improves estimates by 30–50% relative to previous wetland delineation studies. Finally, requiring modest ancillary data, the study results provide a practical and efficient option for wetland management in data-sparse regions or un-gauged watersheds.
54 ENVIRONMENTAL SCIENCES↗
Hurricane Maria Puerto Rico Landsat Analysis
Hurricane Maria made landfall as a strong Category 4 storm in southeast Puerto Rico on September 20th, 2018. The powerful storm traversed the island in a northwesterly direction causing widespread destruction. This study focused on a rapid assessment of Hurricane Marias impact to Puerto Ricos forests. Calibrated and corrected Landsat 8 image composites for the entire island were generated using Google Earth Engine for a comparable pre-Maria and post-Maria time period that accounted for phenology. Spectral mixture analysis (SMA) using image-derived end members was carried out on both composites to calculate the change in the non-photosynthetic vegetation (Delta-NPV) spectral response, a metric that quantifies the increased fraction of exposed wood and surface litter associated with tree mortality and crown damage from the storm. Hurricane simulations were also conducted using the Weather Research and Forecasting (WRF) regional climate model to estimate wind speeds associated with forest disturbance. Dramatic changes in forest structure across the entire island were evident from pre- and post-Maria composited Landsat 8 images. A Delta-NPV map for only the forested pixels illustrated significant spatial variability in disturbance, with patterns that associated with factors such as slope, aspect and elevation. An initial order-of-magnitude impact estimate based on previous work indicated that Hurricane Maria may have caused mortality and severe damage to 23-31 million trees. Additional field work and image analyses are required to further detail the impact of Hurricane Maria to Puerto Rico forests. A minor update to this dataset was posted on April 20, 2018. The previous version is being retired. If you need access to the prior version of the data, email ngee-tropics-archive@lbl.gov.
54 ENVIRONMENTAL SCIENCES↗
Multiscale assessment of land surface phenology from harmonized Landsat 8 and Sentinel-2, PlanetScope, and PhenoCam imagery
As the spatial and temporal resolution of remotely sensed imagery has improved over the last four decades, algorithms for monitoring and mapping seasonal changes in surface properties have evolved rapidly. Most recently, the availability of daily PlanetScope imagery has created new opportunities for monitoring the land surface phenology (LSP) of terrestrial ecosystems at high spatial resolution. However, the quality and value of LSP information from PlanetScope imagery have not been systematically examined. In this paper, we evaluate the character and quality of LSP information derived from PlanetScope by comparing time series of vegetation indices and LSP metrics from PlanetScope to corresponding time series and LSP metrics derived from Harmonized Landsat 8 and Sentinel-2 (HLS) imagery and PhenoCams at six sites that span a diverse range of land cover types and climate. Results show that vegetation index time series from all three data sources show high temporal correlation, and LSP metrics derived from HLS, PlanetScope, and PhenoCam show high agreement with negligible bias. Semi-variograms for phenometrics estimated from PlanetScope imagery indicate that the majority of spatial variance captured in PlanetScope phenometrics occurs well below the spatial resolution HLS imagery. At the same time, LSP metrics from HLS are most strongly correlated with the 50–75% quantiles of 3 m LSP metrics from PlanetScope. This indicates that HLS captures the average phenology at sub-pixel scale captured in PlanetScope imagery. Here, our results represent the first comprehensive comparison of LSP metrics estimated from PlanetScope and publicly available moderate spatial resolution imagery, and provide insights regarding: (1) the quality and character of LSP metrics derived from HLS and PlanetScope; and (2) the relative merits and trade-offs associated with the use of each data source for LSP studies.
54 ENVIRONMENTAL SCIENCES↗
Estimating Fine-Resolution Shortwave Broadband Albedo of Croplands from Harmonized Landsat and Sentinel-2 Data
Altered surface albedo due to land-cover conversions and management is a significant driver of global climate change. Albedo can be directly measured at ground stations, and remote sensing data can be used to scale-up albedo values to regional and global levels. Some previous studies have retrieved fine-resolution (10–30 m) instantaneous albedo and coarse-resolution (500–1000 m) daily mean albedo from remote sensing data, but they all required the input of Moderate Resolution Imaging Spectroradiometer (MODIS) albedo information at 500-m resolution, and none have assembled both instantaneous and daily albedo based exclusively on fine-resolution satellite data. Here, to address this issue, we compiled 387 instantaneous and 346 daily albedo records using field net radiometer measurements from the bioenergy croplands at the W. K. Kellogg Biological Station in southwest Michigan. We then connected these albedo records with a suite of variables derived from harmonized Landsat and Sentinel-2 data through two machine learning algorithms (random forest regression and extreme gradient boosting) to retrieve clear-sky instantaneous and daily shortwave broadband albedo. The performance statistics indicate reasonable accuracy of model results [root-mean-square error (RMSE)] around or below 0.03 except for snow-covered surfaces), suggesting that the retrieval of both instantaneous and daily albedo based exclusively on fine-resolution satellite data is promising. To facilitate the use of fine-resolution albedo products at the global level, future efforts need to include more albedo records of diverse surface cover types, as well as to accurately model daily albedo for cloudy days to address the “clear-sky bias.”
Harmonized Landsat and Sentinel-2↗
Impact of Climate Variabilities and Human Activities on Surface Water Extents in Reservoirs of Yongding River Basin, China, from 1985 to 2016 Based on Landsat Observations and Time Series Analysis
Yongding River is the largest river flowing through Beijing, the capital city of China. In recent years, Yongding River Basin (YDRB) has seen increasing human impacts on water resources, posing serious challenges in hydrological and ecological health. In this study, remote sensing techniques and statistical time series approaches for hydrological studies were combined to characterize the dynamics and driving factors of reservoir water extents in YDRB during 1985–2016. First, 107 Landsat 4, 5, 7 and 8 images were used to extract surface water extents in YDRB during 1985–2016 using a combination of water indices and Otsu threshold algorithm. Significant positive correlation was found between water extents and the annual inflow for the two biggest reservoirs, the downstream Guanting and upstream Cetian reservoirs, proving their representativeness of surface water availability in this basin. Then, statistical time series approaches including trend-free pre-whitening Mann-Kendall trend test, Pettit change-point test and double mass curve method, which are frequently used in hydrological studies, were adopted to quantify the trend of reservoir water extents dynamics and the relative contributions of climate variability and human activities. Results showed that the water extents in both reservoirs exhibited significant downward trend with change point occurring in 2001 and 2005 for Guanting and Cetian, respectively. About 74%~75% of the shrinkage during the post-change period can be attributed to human activities, among which GDP, population, electricity power production, raw coal production, steel and crude iron production, value of agriculture output, and urban area were the major human drivers. Hydrological connectivity between the upstream Cetian and downstream Guanting reservoirs declined during the post-change period. Since 2012, water extents in both reservoirs recovered as a result of various governmental water management policies including the South-to-North Water Diversion Project. The methodology introduced in this study can be used for analyzing the dynamics and driving mechanism of surface water resources, especially for un-gauged or poorly-gauged watersheds.
47 OTHER INSTRUMENTATION↗
Integrating a process-based ecosystem model with Landsat imagery to assess impacts of forest disturbance on terrestrial carbon dynamics: Case studies in Alabama and Mississippi
Forest ecosystems in the southern United States are dramatically altered by three major disturbances: timber harvesting, hurricane, and permanent land conversion. Understanding and quantifying effects of disturbance on forest carbon, nitrogen, and water cycles is critical for sustainable forest management in this region. In this study, we introduced a process-based ecosystem model for simulating forest disturbance impacts on ecosystem carbon, nitrogen, and water cycles. Based on forest mortality data classified from Landsat TM/ETM + images, this model was then applied to estimate changes in carbon storage using Mississippi and Alabama as a case study. Mean annual forest mortality rate for these states was 2.37%. Due to frequent disturbance, over 50% of the forest land in the study region was less than 30 years old. Forest disturbance events caused a large carbon source (138.92 Tg C, 6.04 Tg C yr -1 ; 1 Tg = 10 12 g) for both states during 1984–2007, accounting for 2.89% (4.81% if disregard carbon storage changes in wood products) of the total forest carbon storage in this region. Large decreases and slow recovery of forest biomass were the main causes for carbon release. Forest disturbance could result in a carbon sink in few areas if wood product carbon was considered as a local carbon pool, indicating the importance of accounting for wood product carbon when assessing forest disturbance effects. The legacy effects of forest disturbance on ecosystem carbon storage could last over 50 years. Lastly, this study implies that understanding forest disturbance impacts on carbon dynamics is of critical importance for assessing regional carbon budgets.
58 GEOSCIENCES↗