05.06.2013 Views

PNNL-13501 - Pacific Northwest National Laboratory

PNNL-13501 - Pacific Northwest National Laboratory

PNNL-13501 - Pacific Northwest National Laboratory

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

Extension of Site Soil/Plant Physical Properties to Landscape and Regional Scales Using<br />

Remote Sensing<br />

Study Control Number: PN00044/1451<br />

Chris J. Murray, Eileen M. Perry<br />

Terrestrial carbon sequestration research requires estimates of plant and soil characteristics across a landscape. Field<br />

measurements, however, are generally limited in number and represent individual sites. Remote sensing imagery may<br />

provide continuous data that can be used to interpolate the limited site data across landscapes. A successful approach to<br />

interpolation would have applications in many areas of environmental sciences.<br />

Project Description<br />

The purpose of this project was to develop new<br />

capabilities in scaling site measurements or modeled point<br />

data over an entire landscape. Remote sensing datasets<br />

were used to derive both 2 m and 15 m estimates of plant<br />

biomass over farm and range test sites. The 2 m<br />

resolution datasets were used to evaluate three approaches<br />

to interpolation of the 2 m data over a landscape: linear<br />

regression, ordinary kriging, and Gaussian simulation<br />

with locally varying means. For a 46-ha farm test site,<br />

ordinary kriging yielded the best results, in terms of both<br />

root mean square error and the simulated green biomass<br />

images produced. The kriging results did improve with<br />

the number of 2 m samples used; 100 was the minimum<br />

number of samples required, and significant improvement<br />

was seen with 200 and 400 samples. The investigators<br />

feel that boundary effects may be the reason that the<br />

Gaussian simulation with locally varying means did not<br />

perform as well as ordinary kriging for this test case.<br />

Introduction<br />

Our <strong>Laboratory</strong> is frequently involved in natural resource<br />

research that requires datasets for spatial models, and<br />

measured and modeled physical soil/plant characteristics<br />

at landscape to regional scales. Murray et al. (1999)<br />

developed geostatistical methods to interpolate shrubsteppe<br />

transect data to landscape scale based on<br />

continuous plant cover data. Our project objective was to<br />

develop new capabilities for scaling site measurements or<br />

modeled data to the landscape or regional scales. This<br />

research was directed to 1) define the current state of<br />

knowledge, and 2) refine or develop new methods to<br />

extend measured or modeled results from sites (plant or<br />

soil carbon) to landscape or regional scales.<br />

Approach<br />

Our approach was to 1) review relevant research; 2) select<br />

a suitable data set; 3) perform an initial evaluation of data;<br />

and 4) perform interpolation based on linear regression,<br />

ordinary kriging, and Gaussian simulation with locally<br />

varying means.<br />

Results and Accomplishments<br />

The literature survey revealed a keen interest in multiscale<br />

issues within the ecology, geosciences, and remote<br />

sensing communities. However, no prior research (other<br />

than performed by the investigators) was directly<br />

applicable to this project. A dataset was selected that had<br />

relevance to a current rangeland health program.<br />

Airborne imagery at 2 m spatial resolution and satellite<br />

imagery at 25 m resolution were acquired during 1999<br />

near Grand Junction, Colorado. These image sets were<br />

used to derive green biomass based on the normalized<br />

difference vegetation index. The purpose of the research<br />

was to evaluate the ability of different methods to<br />

interpolate biomass site data (represented by randomly<br />

selected pixels from the 2 m imagery) to landscape scale.<br />

The first method was linear regression based on fitting the<br />

2 m pixels to 15 m normalized difference vegetation<br />

index values derived from the satellite imagery. The<br />

second method was ordinary kriging (Goovaerts 1997),<br />

where only the 2 m imagery was used; randomly selected<br />

pixels were used to fit a variogram, and then kriging was<br />

used to interpolate to the landscape scale. Ordinary<br />

kriging was performed using datasets with 100, 200, and<br />

400 randomly selected 2 m pixels; the results are shown<br />

in Figure 1. The third method was Gaussian simulation<br />

with locally varying means (Goovaerts 1997), which used<br />

the satellite imagery to partition the landscape into three<br />

Earth System Science 221

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!