Geology topics

X. Chen

Publications and source records attributed to X. Chen.

16 recordsLinked to original sources

USGS2026

Large streamflow differences between forested and urbanized watersheds in the energy-limited eastern United States: The role of evapotranspiration and impervious surfaces

Urban forests and other green infrastructures have been viewed as part of the “Nature-based Solutions” (NbS) to mitigate emerging urban environmental change. This study focuses on the role of evapotranspiration (ET) in regulating water balances of small watersheds in the eastern United States. We compared streamflow and ET patterns at daily, monthly and annual scales and linked these hydrological variables to the physical properties of 11 paired watersheds dominated by forests (FW) or urban (UW) land covers. The annual precipitation ranged from 1028 mm to 1683 mm and potential ET (PET) from 815 mm to 1450 mm. The mean annual flow/precipitation (Q/P) ratios were 0.26 ± 0.13 and 0.41 ± 0.1 for FW and UW, respectively. Overall, UW had lower annual ET (772 mm in UW vs. 947 mm in FW), but higher mean annual and (∼58% higher), monthly water yield (17%–186% higher), and peakflow rates (up to 100 times higher) than FW. The streamflow differences between FW and UW were most pronounced during the growing season and early winter (June-November). The mean Q/P ratios for 30 large hurricane events (2016–2021) were 0.12 ± 0.11 and 0.38 ± 0.23 for FW and UW, respectively. The flow rates in the dormant season (around December-May) in UW were similar or lower than FW. We developed conceptual models to explain the seasonal and storm event streamflow differences using background climate (PET), ET, and land surface characteristics. Urban NbS designs should factor in strategies that maximize ET while minimizing impervious surfaces enhancing watershed “sponge” and “pump” functions.

eastern United States↗

USGS2024

SCEC/USGS Community Stress Drop Validation Study: How spectral fitting approaches influence measured source parameters

Spectral source parameters used to estimate an earthquake’s stress drop (Δσ) can vary significantly across measurement approaches. The Statewide California Earthquake Center/U.S. Geological Survey Community Stress-Drop Validation Study was initiated to compare source parameter estimates, focusing initially on a dataset from the 2019 Ridgecrest earthquake sequence. As part of that validation effort, here we focus on one potential source of uncertainty: whether spectral fitting approaches alone, applied to a common set of spectra from the 2019 Ridgecrest sequence result in different source parameter estimates. By using a common set of benchmark spectra analyzed across a consistent frequency band of 1–40 Hz, we eliminate many sources of variability. A subgroup of validation study participants volunteered to estimate the low-frequency displacement (Ω0) and corner frequency ( f c ) by fitting a smooth function to benchmark displacement spectra. Participants used linear- or log-sampled spectra, assumed a Brune or Boatwright spectral model, and applied different misfit criteria. We compare 17 approaches used to estimate Ω0, f c , and Δσ for 54 earthquake spectra. Our results reveal that 35% of events have Δσ estimates within a factor of two, whereas others exhibit variations exceeding an order of magnitude. The variability in Ω0 and f c can largely be attributed to whether a spectrum is consistent with the smooth function of an idealized simple crack model. The trade-off between Ω0 and f c may be more pronounced when using linearly sampled spectra, as higher frequency spectral bumps control the fits. As expected, methods that assumed a Boatwright model tended to have lower Ω0 and somewhat higher f c compared to those assuming a Brune model, although resulting Δσ estimates are similar. When compared to the overall validation study results, the fitting approach alone may account for between 5% and 90% (25% on average) of the total variability in spectral Δσ.

California↗

USGS2018

The GFDL global atmosphere and land model AM4.0/LM4.0: 1. Simulation characteristics with prescribed SSTs

In this two‐part paper, a description is provided of a version of the AM4.0/LM4.0 atmosphere/land model that will serve as a base for a new set of climate and Earth system models (CM4 and ESM4) under development at NOAA's Geophysical Fluid Dynamics Laboratory (GFDL). This version, with roughly 100 km horizontal resolution and 33 levels in the vertical, contains an aerosol model that generates aerosol fields from emissions and a “light” chemistry mechanism designed to support the aerosol model but with prescribed ozone. In Part 1, the quality of the simulation in AMIP (Atmospheric Model Intercomparison Project) mode—with prescribed sea surface temperatures (SSTs) and sea‐ice distribution—is described and compared with previous GFDL models and with the CMIP5 archive of AMIP simulations. The model's Cess sensitivity (response in the top‐of‐atmosphere radiative flux to uniform warming of SSTs) and effective radiative forcing are also presented. In Part 2, the model formulation is described more fully and key sensitivities to aspects of the model formulation are discussed, along with the approach to model tuning.

Journal of Advances in Modeling Earth Systems↗

USGS2018

The GFDL global atmosphere and land model AM4.0/LM4.0: 2. Model description, sensitivity studies, and tuning strategies

In Part 2 of this two‐part paper, documentation is provided of key aspects of a version of the AM4.0/LM4.0 atmosphere/land model that will serve as a base for a new set of climate and Earth system models (CM4 and ESM4) under development at NOAA's Geophysical Fluid Dynamics Laboratory (GFDL). The quality of the simulation in AMIP (Atmospheric Model Intercomparison Project) mode has been provided in Part 1. Part 2 provides documentation of key components and some sensitivities to choices of model formulation and values of parameters, highlighting the convection parameterization and orographic gravity wave drag. The approach taken to tune the model's clouds to observations is a particular focal point. Care is taken to describe the extent to which aerosol effective forcing and Cess sensitivity have been tuned through the model development process, both of which are relevant to the ability of the model to simulate the evolution of temperatures over the last century when coupled to an ocean model.

Journal of Advances in Modeling Earth Systems↗

USGS2016

An intercomparison of remote sensing river discharge estimation algorithms from measurements of river height, width, and slope

The Surface Water and Ocean Topography (SWOT) satellite mission planned for launch in 2020 will map river elevations and inundated area globally for rivers >100 m wide. In advance of this launch, we here evaluated the possibility of estimating discharge in ungauged rivers using synthetic, daily ‘‘remote sensing’’ measurements derived from hydraulic models corrupted with minimal observational errors. Five discharge algorithms were evaluated, as well as the median of the five, for 19 rivers spanning a range of hydraulic and geomorphic conditions. Reliance upon a priori information, and thus applicability to truly ungauged reaches, varied among algorithms: one algorithm employed only global limits on velocity and depth, while the other algorithms relied on globally available prior estimates of discharge. We found at least one algorithm able to estimate instantaneous discharge to within 35% relative root-mean-squared error (RRMSE) on 14/16 nonbraided rivers despite out-of-bank flows, multichannel planforms, and backwater effects. Moreover, we found RRMSE was often dominated by bias; the median standard deviation of relative residuals across the 16 nonbraided rivers was only 12.5%. SWOT discharge algorithm progress is therefore encouraging, yet future efforts should consider incorporating ancillary data or multialgorithm synergy to improve results

Water Resources Research↗

USGS2016

Simulated high-latitude soil thermal dynamics during the past four decades

Soil temperature (Ts ) change is a key indicator of the dynamics of permafrost. On seasonal and inter-annual time scales, the variability of Ts determines the active layer depth, which regulates hydrological soil properties and biogeochemical processes. On the multi-decadal scale, increasing T 5 s not only drives permafrost thaw/retreat, but can also trigger and accelerate the decomposition of soil organic carbon. The magnitude of permafrost carbon feedbacks is thus closely linked to the rate of change of soil thermal regimes. In this study, we used nine process-based ecosystem models with permafrost processes, all forced by different observation-based climate forcing during the period 1960–2000, to characterize the warming rate of Ts 10 in permafrost regions. There is a large spread of Ts trends at 20 cm depth across the models, with trend values ranging from 0.010 ± 0.003 to 0.031 ± 0.005 ◦C yr−1 . Most models show smaller increase in Ts with increasing depth. Air temperature (Ta ) and longwave downward radiation (LWDR) are the main drivers of Ts trends, but their relative contributions differ 15 amongst the models. Different trends of LWDR used in the forcing of models can explain 61 % of their differences in Ts trends, while trends of Ta only explain 5 % of the differences in Ts trends. Uncertain climate forcing contributes a larger uncertainty in Ts trends (0.021 ± 0.008 ◦C yr−1 , mean ± SD) than the uncertainty of model structure (0.012 ± 0.001 ◦C yr−1 ), diagnosed from the range of response between different mod- 20 els, normalized to the same forcing. In addition, the loss rate of near-surface permafrost area, defined as total area where the maximum seasonal active layer thickness (ALT) is less than 3 m loss rate is found to be significantly correlated with the magnitude of the trends of Ts at 1 m depth across the models (R = −0.85, P = 0.003), but not with the initial total near-surface permafrost area (R = −0.30, P = 0.438). The sensitivity of the total boreal near-surface permafrost area to T 25 s at 1 m, is estimated to be of −2.80 ± 0.67 million km2 ◦C −1 . Finally, by using two long-term LWDR datasets and relationships between trends of LWDR and Ts across models, we infer an observationconstrained total boreal near-surface permafrost area decrease comprised between 39 ± 14 × 103 and 75 ± 14 × 103 km2 yr−1 from 1960 to 2000. This corresponds to 9– 18 % degradation of the current permafrost area.

Cryosphere Discussions↗

USGS2015

A simplified, data-constrained approach to estimate the permafrost carbon–climate feedback

We present an approach to estimate the feedback from large-scale thawing of permafrost soils using a simplified, data-constrained model that combines three elements: soil carbon (C) maps and profiles to identify the distribution and type of C in permafrost soils; incubation experiments to quantify the rates of C lost after thaw; and models of soil thermal dynamics in response to climate warming. We call the approach the Permafrost Carbon Network Incubation–Panarctic Thermal scaling approach (PInc-PanTher). The approach assumes that C stocks do not decompose at all when frozen, but once thawed follow set decomposition trajectories as a function of soil temperature. The trajectories are determined according to a three-pool decomposition model fitted to incubation data using parameters specific to soil horizon types. We calculate litterfall C inputs required to maintain steady-state C balance for the current climate, and hold those inputs constant. Soil temperatures are taken from the soil thermal modules of ecosystem model simulations forced by a common set of future climate change anomalies under two warming scenarios over the period 2010 to 2100. Under a medium warming scenario (RCP4.5), the approach projects permafrost soil C losses of 12.2–33.4 Pg C; under a high warming scenario (RCP8.5), the approach projects C losses of 27.9–112.6 Pg C. Projected C losses are roughly linearly proportional to global temperature changes across the two scenarios. These results indicate a global sensitivity of frozen soil C to climate change ( γ sensitivity) of −14 to −19 Pg C °C −1 on a 100 year time scale. For CH 4 emissions, our approach assumes a fixed saturated area and that increases in CH 4 emissions are related to increased heterotrophic respiration in anoxic soil, yielding CH 4 emission increases of 7% and 35% for the RCP4.5 and RCP8.5 scenarios, respectively, which add an additional greenhouse gas forcing of approximately 10–18%. The simplified approach presented here neglects many important processes that may amplify or mitigate C release from permafrost soils, but serves as a data-constrained estimate on the forced, large-scale permafrost C response to warming.

Philosophical Transactions of the Royal Society A:↗

USGS2015

Assessment of model estimates of land-atmosphere CO 2 exchange across northern Eurasia

A warming climate is altering land-atmosphere exchanges of carbon, with a potential for increased vegetation productivity as well as the mobilization of permafrost soil carbon stores. Here we investigate land-atmosphere carbon dioxide (CO 2 ) cycling through analysis of net ecosystem productivity (NEP) and its component fluxes of gross primary productivity (GPP) and ecosystem respiration (ER) and soil carbon residence time, simulated by a set of land surface models (LSMs) over a region spanning the drainage basin of Northern Eurasia. The retrospective simulations cover the period 1960–2009 at 0.5° resolution, which is a scale common among many global carbon and climate model simulations. Model performance benchmarks were drawn from comparisons against both observed CO 2 fluxes derived from site-based eddy covariance measurements as well as regional-scale GPP estimates based on satellite remote-sensing data. The site-based comparisons depict a tendency for overestimates in GPP and ER for several of the models, particularly at the two sites to the south. For several models the spatial pattern in GPP explains less than half the variance in the MODIS MOD17 GPP product. Across the models NEP increases by as little as 0.01 to as much as 0.79 g C m −2 yr −2 , equivalent to 3 to 340 % of the respective model means, over the analysis period. For the multimodel average the increase is 135 % of the mean from the first to last 10 years of record (1960–1969 vs. 2000–2009), with a weakening CO 2 sink over the latter decades. Vegetation net primary productivity increased by 8 to 30 % from the first to last 10 years, contributing to soil carbon storage gains. The range in regional mean NEP among the group is twice the multimodel mean, indicative of the uncertainty in CO 2 sink strength. The models simulate that inputs to the soil carbon pool exceeded losses, resulting in a net soil carbon gain amid a decrease in residence time. Our analysis points to improvements in model elements controlling vegetation productivity and soil respiration as being needed for reducing uncertainty in land-atmosphere CO 2 exchange. These advances will require collection of new field data on vegetation and soil dynamics, the development of benchmarking data sets from measurements and remote-sensing observations, and investments in future model development and intercomparison studies.

Biogeosciences↗

USGS2011

Faulted terrace risers place new constraints on the late Quaternary slip rate for the central Altyn Tagh fault, northwest Tibet

The active, left-lateral Altyn Tagh fault defines the northwestern margin of the Tibetan Plateau in western China. To clarify late Quaternary temporal and spatial variations in slip rate along the central portion of this fault system (85°–90°E), we have more than doubled the number of dated offset markers along the central Altyn Tagh fault. In particular, we determined offset-age relations for seven left-laterally faulted terrace risers at three sites (Kelutelage, Yukuang, and Keke Qiapu) spanning a 140-km-long fault reach by integrating surficial geologic mapping, topographic surveys (total station and tripod–light detection and ranging [T-LiDAR]), and geochronology (radiocarbon dating of organic samples, 230 Th/U dating of pedogenic carbonate coatings on buried clasts, and terrestrial cosmogenic radionuclide exposure age dating applied to quartz-rich gravels). At Kelutelage, which is the westernmost site (37.72°N, 86.67°E), two faulted terrace risers are offset 58 ± 3 m and 48 ± 4 m, and formed at 6.2–6.1 ka and 5.9–3.7 ka, respectively. At the Yukuang site (38.00°N, 87.87°E), four faulted terrace risers are offset 92 ± 12 m, 68 ± 6 m, 55 ± 13 m, and 59 ± 9 m and formed at 24.2–9.5 ka, 6.4–5.0 ka, 5.1–3.9 ka, and 24.2–6.4 ka, respectively. At the easternmost site, Keke Qiapu (38.08°N, 88.12°E), a faulted terrace riser is offset 33 ± 6 m and has an age of 17.1–2.2 ka. The displacement-age relationships derived from these markers can be satisfied by an approximately uniform slip rate of 8–12 mm/yr. However, additional analysis is required to test how much temporal variability in slip rate is permitted by this data set.

Central Altyn Tagh fault↗

USGS2011

Structural stability of methane hydrate at high pressures

The structural stability of methane hydrate under pressure at room temperature was examined by both in-situ single-crystal and powder X-ray diffraction techniques on samples with structure types I, II, and H in diamond-anvil cells. The diffraction data for types II (sII) and H (sH) were refined to the known structures with space groups Fd 3 m and P 6 3 / mmc , respectively. Upon compression, sI methane hydrate transforms to the sII phase at 120 MPa, and then to the sH phase at 600 MPa. The sII methane hydrate was found to coexist locally with sI phase up to 500 MPa and with sH phase up to 600 MPa. The pure sH structure was found to be stable between 600 and 900 MPa. Methane hydrate decomposes at pressures above 3 GPa to form methane with the orientationally disordered Fm 3 m structure and ice VII ( Pn 3 m ). The results highlight the role of guest (CH 4 )-host (H 2 O) interactions in the stabilization of the hydrate structures under pressure.

Geoscience Frontiers↗

USGS2010

Landslide-dammed lake at Tangjiashan, Sichuan province, China (triggered by the Wenchuan Earthquake, May 12, 2008): Risk assessment, mitigation strategy, and lessons learned

Landslides and rock avalanches triggered by the 2008 Wenchuan Earthquake produced 257 landslide dams, mainly situated along the eastern boundary of the Qinghai-Tibet Plateau where rivers descend approximately 3,000 m into the Sichuan Basin. The largest of these dams blocked the Tongkou River (a tributary of the Fujiang River) at Tangjiashan. The blockage, consisting of 2.04 × 10 7 m 3 of landslide debris, impounded a lake with a projected maximum volume of 3.15 × 10 8 m 3 , potentially inundating 8.92 km 2 of terrain. Its creation during the rainy season and the possibility of an uncontrolled release posed a serious, impending threat to at least 1.3 million people downstream that could add substantially to the total of 69,200 individuals directly killed by the earthquake. Risk assessment of the blockage indicated that it was unlikely to collapse suddenly, and that eventual overtopping could be mitigated by notching the structure in order to create an engineered breach and achieve safe drainage of the lake. In addition to the installation of monitoring and warning instrumentation, for emergency planning we estimated several outburst scenarios equivalent to 20, 25, 33, and 50% of the dam failing suddenly, creating, respectively, 3.35, 3.84, 4.22, and 4.65 km 2 of flooded area, and overbank water depths of 4.6, 5.1, 5.7, and 6.2 m, respectively, in Mianyang, the second largest city in Sichuan Province, 48 km downstream from the blockage. Based on these scenarios, recommendations and plans for excavating a sluiceway, draining the lake, and downstream evacuation were proposed and later were implemented successfully, with the blockage breached by overtopping on June 10, less than a month after dam emplacement. The peak discharge of the release only slightly exceeded the flood of record at Mianyang City. No lives were lost, and significant property damage was avoided. Post-breaching evaluation reveals how future similar mitigation can be improved. Although initial breach erosion was slow, later erosion was judged uncontrollably rapid; increased slope of the engineered channel and adoption of a compound, trapezoid–triangular cross-section can be considered, as can other measures to control the rate of breach incision. Evacuees from Mianyang City spent an unnecessarily long time (12 days) in temporary settlements; more precise risk management planning can reduce this time in the future.

Sichuan Province↗

USGS2009

Crustal structure across the Three Gorges area of the Yangtze platform, central China, from seismic refraction/wide-angle reflection data

We present active-source seismic data recorded along a 300 km-long profile across the Three Gorges area of the western Yangtze platform, central China. From west to east, the profile crosses the Zigui basin, Huangling dome and Jianghan basin. The derived crustal P -wave velocity structure changes significantly across the Tongchenghe fault that lies at the transition from the Huangling dome to the Jianghan basin. West of the Tongchenghe fault, beneath the Zigui basin and the Huangling dome, we observe a ~ 42 km thick crust of relatively low average velocity (6.3–6.4 km/s). In contrast, east of the Tongchenghe fault, beneath the Jianghan basin, the crust is only 30 km thick and has a high average velocity (6.6–6.7 km/s). A west–east variation in crustal composition along the Tongchenghe fault is also inferred. West of the fault, P -wave velocities suggest a felsic composition with an intermediate layer at the base of the crust, whilst, east of the fault, felsic, intermediate, and mafic crustal layers are apparent. Our results suggest that the crust beneath the Jianghan basin has been thinned by rifting, accompanied by intrusion of the lower crust by mafic dikes and sills. The west-to-east division of the crust in the Three Gorges area coincides with first-order geophysical contrasts in gravity, topography, crustal and lithospheric thickness.

Tectonophysics↗

USGS2007

Monitoring boreal forest leaf area index across a Siberian burn chronosequence: A MODIS validation study

Landscapes containing differing amounts of ecological disturbance provide an excellent opportunity to validate and better understand the emerging Moderate Resolution Imaging Spectrometer (MODIS) vegetation products. Four sites, including 1‐year post‐fire coniferous, 13‐year post‐fire deciduous, 24‐year post‐fire deciduous, and >100 year old post‐fire coniferous forests, were selected to serve as a post‐fire chronosequence in the central Siberian region of Krasnoyarsk (57.3°N, 91.6°E) with which to study the MODIS leaf area index (LAI) and vegetation index (VI) products. The collection 4 MODIS LAI product correctly represented the summer site phenologies, but significantly underestimated the LAI value of the >100 year old coniferous forest during the November to April time period. Landsat 7‐derived enhanced vegetation index (EVI) performed better than normalized difference vegetation index (NDVI) to separate the deciduous and conifer forests, and both indices contained significant correlation with field‐derived LAI values at coniferous forest sites ( r 2 = 0.61 and r 2 = 0.69, respectively). The reduced simple ratio (RSR) markedly improved LAI prediction from satellite measurements ( r 2 = 0.89) relative to NDVI and EVI. LAI estimates derived from ETM+ images were scaled up to evaluate the 1 km resolution MODIS LAI product; from this analysis MODIS LAI overestimated values in the low LAI deciduous forests (where LAI<5) and underestimated values in the high LAI conifer forests (where LAI>6). Our results indicate that further research on the MODIS LAI product is warranted to better understand and improve remote LAI quantification in disturbed forest landscapes over the course of the year.

Krasnoyarsk Kray↗

USGS2006

The Short Wave Aerostat-Mounted Imager (SWAMI): A novel platform for acquiring remotely sensed data from a tethered balloon

We describe a new remote sensing system called the Short Wave Aerostat-Mounted Imager (SWAMI). The SWAMI is designed to acquire co-located video imagery and hyperspectral data to study basic remote sensing questions and to link landscape level trace gas fluxes with spatially and temporally appropriate spectral observations. The SWAMI can fly at altitudes up to 2 km above ground level to bridge the spatial gap between radiometric measurements collected near the surface and those acquired by other aircraft or satellites. The SWAMI platform consists of a dual channel hyperspectral spectroradiometer, video camera, GPS, thermal infrared sensor, and several meteorological and control sensors. All SWAMI functions (e.g. data acquisition and sensor pointing) can be controlled from the ground via wireless transmission. Sample data from the sampling platform are presented, along with several potential scientific applications of SWAMI data.

Remote Sensing of Environment↗

USGS2005

A simple and effective radiometric correction method to improve landscape change detection across sensors and across time

Satellite data offer unrivaled utility in monitoring and quantifying large scale land cover change over time. Radiometric consistency among collocated multi-temporal imagery is difficult to maintain, however, due to variations in sensor characteristics, atmospheric conditions, solar angle, and sensor view angle that can obscure surface change detection. To detect accurate landscape change using multi-temporal images, we developed a variation of the pseudoinvariant feature (PIF) normalization scheme: the temporally invariant cluster (TIC) method. Image data were acquired on June 9, 1990 (Landsat 4), June 20, 2000 (Landsat 7), and August 26, 2001 (Landsat 7) to analyze boreal forests near the Siberian city of Krasnoyarsk using the normalized difference vegetation index (NDVI), enhanced vegetation index (EVI), and reduced simple ratio (RSR). The temporally invariant cluster (TIC) centers were identified via a point density map of collocated pixel VIs from the base image and the target image, and a normalization regression line was created to intersect all TIC centers. Target image VI values were then recalculated using the regression function so that these two images could be compared using the resulting common radiometric scale. We found that EVI was very indicative of vegetation structure because of its sensitivity to shadowing effects and could thus be used to separate conifer forests from deciduous forests and grass/crop lands. Conversely, because NDVI reduced the radiometric influence of shadow, it did not allow for distinctions among these vegetation types. After normalization, correlations of NDVI and EVI with forest leaf area index (LAI) field measurements combined for 2000 and 2001 were significantly improved; the r 2 values in these regressions rose from 0.49 to 0.69 and from 0.46 to 0.61, respectively. An EVI "cancellation effect" where EVI was positively related to understory greenness but negatively related to forest canopy coverage was evident across a post fire chronosequence with normalized data. These findings indicate that the TIC method provides a simple, effective and repeatable method to create radiometrically comparable data sets for remote detection of landscape change. Compared to some previous relative radiometric normalization methods, this new method does not require high level programming and statistical skills, yet remains sensitive to landscape changes occurring over seasonal and inter-annual time scales. In addition, the TIC method maintains sensitivity to subtle changes in vegetation phenology and enables normalization even when invariant features are rare. While this normalization method allowed detection of a range of land use, land cover, and phenological/biophysical changes in the Siberian boreal forest region studied here, it is necessary to further examine images representing a wide variety of ecoregions to thoroughly evaluate the TIC method against other normalization schemes. ?? 2005 Elsevier Inc. All rights reserved.

Remote Sensing of Environment↗

USGS2004

Using lidar and effective LAI data to evaluate IKONOS and Landsat 7 ETM+ vegetation cover estimates in a ponderosa pine forest

Structural and functional analyses of ecosystems benefit when high accuracy vegetation coverages can be derived over large areas. In this study, we utilize IKONOS, Landsat 7 ETM+, and airborne scanning light detection and ranging (lidar) to quantify coniferous forest and understory grass coverages in a ponderosa pine ( Pinus ponderosa ) dominated ecosystem in the Black Hills of South Dakota. Linear spectral mixture analyses of IKONOS and ETM+ data were used to isolate spectral endmembers (bare soil, understory grass, and tree/shade) and calculate their subpixel fractional coverages. We then compared these endmember cover estimates to similar cover estimates derived from lidar data and field measures. The IKONOS-derived tree/shade fraction was significantly correlated with the field-measured canopy effective leaf area index (LAI e ) ( r 2 =0.55, p <0.001) and with the lidar-derived estimate of tree occurrence ( r 2 =0.79, p <0.001). The enhanced vegetation index (EVI) calculated from IKONOS imagery showed a negative correlation with the field measured tree canopy effective LAI and lidar tree cover response ( r 2 =0.30, r =−0.55 and r 2 =0.41, r =−0.64, respectively; p <0.001) and further analyses indicate a strong linear relationship between EVI and the IKONOS-derived grass fraction ( r 2 =0.99, p <0.001). We also found that using EVI resulted in better agreement with the subpixel vegetation fractions in this ecosystem than using normalized difference of vegetation index (NDVI). Coarsening the IKONOS data to 30 m resolution imagery revealed a stronger relationship with lidar tree measures ( r 2 =0.77, p <0.001) than at 4 m resolution ( r 2 =0.58, p <0.001). Unmixed tree/shade fractions derived from 30 m resolution ETM+ imagery also showed a significant correlation with the lidar data ( r 2 =0.66, p <0.001). These results demonstrate the power of using high resolution lidar data to validate spectral unmixing results of satellite imagery, and indicate that IKONOS data and Landsat 7 ETM+ data both can serve to make the important distinction between tree/shade coverage and exposed understory grass coverage during peak summertime greenness in a ponderosa pine forest ecosystem.

Remote Sensing of Environment↗