• Skip to main content
itrc_logo

Geospatial Analysis for Optimization at Environmental Sites

Navigating this Website
Overview
Fact Sheets
Fact Sheets Overview
Fact Sheet 1: Do You Need Geospatial Analysis?
Fact Sheet 2: Are Conditions Suitable for Geospatial Analysis?
Fact Sheet 3: How is Geospatial Analysis Applied?
Fact Sheet 4: What Software is Available to Help?
PM's Tool Box
PM's Tool Box Overview
Review Checklist
Choosing Methods
Common Misapplications
Optimization Questions
Geospatial Analysis Support for Optimization Questions in the Project Life Cycle
Data Requirements
General Considerations
Methods for Optimization
Geospatial Methods for Optimization Questions in the Project Life Cycle Stages
Release Detection
Site Characterization
Remediation
Monitoring
Closure
Documenting Results
Fundamental Concepts
Fundamental Concepts for Geospatial Analysis
Basic Data Concepts for Geospatial Analysis
Interpolation Methods and Model Prediction
Uncertainty in Geospatial Analyses
Characteristics of Interpolation Methods
Work Flow
Work Flow for Conducting Geospatial Analysis
Geospatial Analysis Work Flow Overview
Perform Exploratory Data Analysis
Select Geospatial Method
Build Geospatial Model
Evaluate Geospatial Method Accuracy
Generate Geospatial Analysis Results
Using Results
Using Analysis Results for Optimization
Plume Intensity and Extent
Trend Maps
Estimating Quantities
Hot Spot Detection
Sample Spacing
Estimating Concentrations Based on Proxy Data
Background Estimation
Quantifying Uncertainty
Remedial Action Optimization
Monitoring Program Optimization
Examples
Examples Overview
Example 1
Example 2
Example 3
Example 4
Methods
Methods Overview
Simple Geospatial Methods
More Complex Geospatial Methods
Advanced Methods
Index of Methods
Software
Software Overview
Software Comparison Tables
Software Descriptions
Workshops and Short Courses
Case Studies
Case Studies Overview
Superfund Site Monitoring Optimization (MAROS)
PAH Contamination in Sediments—Uncertainty Analysis (Isatis)
Optimization of Long-Term Monitoring at Former Nebraska Ordnance Plant (GTS; Summit Envirosolutions)
Optimization of Lead-Contaminated Soil Remediation at a Former Lead Smelter (EVS/MVS)
Extent of Radiological Contamination in Soil at Four Sites near the Fukushima Daiichi Power Plant, Japan (ArcGIS)
Optimization of Groundwater Monitoring at a Research Facility in New Jersey (GWSDAT)
Optimization of Sediment Sampling at a Tidally Influenced Site (ArcGIS)
Stringfellow Superfund Site Monitoring Optimization (MAROS)
Lead Contamination in Soil (ArcGIS)
Stakeholder Perspectives
Additional Information
Project Life Cycle Stages
History of Remedial Process Optimization
Additional Resources
Acronyms
Glossary
Index of Methods
Acknowledgments
Team Contacts
Document Feedback

 

Geospatial Analysis for Optimization at Environmental Sites
HOME

Extent of Radiological Contamination in Soil at Four Sites near the Fukushima Daiichi Power Plant, Japan (ArcGIS)

Contact: Ted Parks, AMEC Foster Wheeler, [email protected], Alex Mikszewski, AMEC Foster Wheeler, [email protected]

gro_case_study_6_5_a_10_30_16

Problem Statement

This case study illustrates the use of ArcGIS interpolation and spatial statistics tools to delineate the extent of contamination and required remediation at four sites in response to the 2011 radiological release at the Fukushima Daiichi nuclear power plant. This project was completed as part of a demonstration of Orion ScanPlot technology (ScanPlot) (AMEC Foster Wheeler 2015). ScanPlot technology is a mobile data-collection system for surface and shallow-depth radiological examination and characterization of open land. ArcGIS tools were used to process and analyze the ScanPlot data to produce maps of detected radioactivity, including concentrations of individual radioactive isotopes, to inform site remediation activities.

Site Background

Assessment and remediation activities were completed at four separate areas within a 20 km radius around the Fukushima Daiichi nuclear power plant. These areas were located in schoolyards, roads, and agricultural fields. Due to variable terrain of the areas, ScanPlot was deployed in three different configurations: a “buggy”, back-pack unit, and truck-mounted unit (see Figure 103).

gro_fig_6_20Figure 103. ScanPlotSM “buggy” (left) and back-pack deployment (right).

Project Objectives

The primary objective was to delineate the spatial extent of soil contamination and optimize remedial decision making at the four study areas. To accomplish this objective, data obtained from ScanPlot was mapped and spatially analyzed during the field investigation using ESRI’s ArcGIS platform. Geospatial interpolation of the data was used to compute zonal statistics and identify areas above risk-based criteria. A secondary objective of the geospatial analysis was to evaluate the conceptual model for radiological contaminant transport and deposition. Following remediation, interpolation and data plotting verified that cleanup levels were achieved.

Data Set

gro_case_study_6_5

At each study area, radiological contaminant levels and corresponding GPS locations were measured at high density on an approximate 1 m grid, resulting in thousands of data points measuring cesium-137 exposure levels in microsieverts per hour (uSv/hr). Over 13,000 data points were collected. Data from ScanPlot was imported into ArcGIS and plotted using the coordinates collected by the GPS component of the system. IDW interpolation was used in the Spatial Analyst extension to rapidly and accurately contour and visualize the data.

Methods

Data plots were generated as the data were collected during the field assessments. The interpolated surfaces informed the ScanPlot locations and ultimately the remedial footprint delineation. IDW was a good interpolation method for this task because of the large volume, high-density, and even distribution of the data. IDW in Spatial Analyst can contour up to approximately 45 million data points.

The resulting maps provided visualizations of the distribution of radiological contaminants. The maps showed that the contamination existed across the sites with relative uniformity. This result was consistent with the expected fall-out delivery method of the contamination. Sites with more variable topography or pronounced drainage features exhibited signs of surface water transport of radioactive materials. Higher levels of contamination were found in ditches and low-lying spots, presumably washed there by runoff. Sites higher in elevation and more exposed to wind showed overall higher concentrations than those in low-lying areas that were partially shielded by topographic features.

Example contour maps of cesium-137 exposure levels in microsieverts per hour (uSv/hr) for two ScanPlot surveys are presented on Figure 104 (left and right).

To optimize remedial decision making, AMEC used the ArcGIS Zonal Statistics tool set to calculate the mean of the IDW interpolated data for a 10 m grid covering the sites. This approach was used to identify portions of the site to excavate where average levels exceeded risk-based criteria. In targeting just these areas, responders could have avoided removing soil over the entire site and still achieved acceptable exposure-point contamination levels. Maps created using the Zonal Statistics tool set in the Spatial Analyst Extension are presented in Figure 105 (left and right).

gro_fig_6_21Figure 104. Maps of cesium-137 levels created using ESRI’s Spatial Analyst Extension.

gro_fig_6_22Figure 105. Left: Average combined cesium-134 and 137 levels created using the Zonal Statistics tool set. Right: Identification of 10-meter grid cells (red) with levels exceeding risk-based thresholds (red).

Due to the sensitivity of the receptors (farmland and schools), regulatory officials conservatively decided to excavate the entirety of the sites to a depth of 10 cm. Following remedial activities, the sites were scanned a second time. The results of the postremediation scan verified the effectiveness of the remedial activities in reducing levels to below risk-based criteria. Postremediation contour maps generating using IDW in Spatial Analyst are presented in Figure 106 (left and right).

gro_fig_6_23Figure 106. Maps of postremediation cesium-137 levels created using ESRI’s Spatial Analyst Extension.

More advanced analysis and interpolation of the radiological data can be completed using ESRI’s Geostatistical Analyst extension. This extension has a robust set of data exploration tools, including histogram, trend, and semivariogram analysis, to characterize the data set for stationarity, a required condition for kriging. Data are stationary if there is a constant mean and variance across the data domain, and data correlation depends only on separation distance and direction, rather than absolute location. An implicit assumption commonly accepted by statisticians is that kriging requires data to be normally distributed, or that it can be transformed to approximately follow a normal distribution. The data exploration utility in Geostatistical Analyst enables the user to understand how to transform data to improve kriging estimates and reduce interpolation error.

Figure 107 (left and right) illustrates the histogram utility in Geostatistical Analyst, while Figure 108 illustrates the trend analysis utility.

gro_fig_6_24Figure 107. Left: histogram of raw cesium-137 data for one of the ScanPlot sites. Right: histogram of log-transformed data better resembles a normal distribution.

gro_fig_6_25Figure 108. Trend visualization in Geostatistical Analyst with the Z-axis representing cesium-137 levels.

Strong trends are evident moving from high levels to low levels, implying large-scale, deterministic correlation of the data. This trend is likely related to runoff-mediated transport of contaminants.

The knowledge gained from performing data exploration can be applied when implementing kriging interpolation in Geostatistical Analyst. First, Geostatistical Analyst can log-transform the data prior to variography. Then, the second-order trend surface observed in Figure 107 can be removed from the data. This trend surface is presented on Figure 108. With these two steps, the data better approximate stationarity, with a constant variance (sill) evident on the semivariogram depicted in Figure 109.

gro_fig_6_26Figure 109. Detrending utility in Geostatistical Analyst.

Removing a second order polynomial equation from the data set places kriging emphasis on local variability rather than large-scale boundary effects. Data point locations are presented in blue. Kriging occurs on the residuals remaining from the detrending. In this case a “stable” semivariogram model is used, which is a hybrid between typical Gaussian and Exponential models. Note anisotropy can be incorporated into the model, but in this case is omitted.

gro_fig_6_27Figure 110. Semivariogram fitting utility in Geostatistical Analyst.

The resulting kriged surface is presented on Figure 111 (left), and is very similar to that previously produced using IDW interpolation (right). This similarity is expected given the high density of data collection. Geostatistical Analyst automatically back-transforms data to the final raster surface, so the user does not need to manually adjust for the prior log-transformation and detrending. The IDW surface appears smoother and performs more averaging in areas with higher measured levels.

gro_fig_6_28Figure 111. Left: interpolated Cesium-137 surface created using ordinary kriging with log-transformation and detrending. Right: interpolation surface created using IDW.

Geostatistical Analyst automatically performs cross-validation on spatial interpolation models and summarizes the results for the user. Cross-validation allows assessment of the accuracy of an interpolation model by computing and investigating the gridding errors, also referred to as residuals. The gridding errors are calculated by removing the first observation from the data set, and using the remaining data and the specified model to estimate a value at that location. The difference between the data value and the estimated value is the residual. This process is repeated for each individual data value and summary statistics (for example, mean error, root-mean-square error) are then calculated for the residual data set.

For kriging models, Geostatistical Analyst also provides the prediction standard error, or prediction standard deviation (same units as data), of each prediction location. This enables calculation of three additional metrics during cross-validation: mean standardized error (dimensionless), average standard error (analogous to root mean square error), and root-mean-square standardized error (measuring the assessment of prediction variability).

Beyond visual comparison, the performance of kriging versus that of IDW can be evaluated using cross-validation statistics. Geostatistical Analyst has a compare models utility that summarized cross-validation results for different models and plots them side-by-side. A screen-shot of this utility comparing the two interpolated surfaces presented in Figure 111 is presented in Figure 112.

gro_fig_6_29Figure 112. Cross-validation comparison of kriging versus IDW in Geostatistical Analyst.

Kriging has a lower mean residual error and a slightly reduced bias compared to IDW. However, IDW has a lower root-mean-square error, demonstrating its suitability for rapid interpolation of high density data. Note that standardized statistics are not available to deterministic models such as IDW that are not based on spatial correlation of the data. The root-mean-square standardized error for the kriging model would ideally be closer to one, indicating that the kriging model could be improved with additional analysis and revision.

Results

The Spatial Analyst extension to Esri’s ArcGIS was used to rapidly interpolate between thousands of data points collected during assessment of four sites impacted by the radiological release at the Fukushima Daiichi nuclear power plant. These interpolations were performed while field investigations were ongoing in order to inform additional data collection efforts. Zonal statistics were calculated on the interpolated surfaces in Spatial Analyst to develop risk-based remedial optimization scenarios. Advanced data analysis performed in Geostatistical Analyst confirmed the presence of trend and spatial correlation in the preremediation data set and validated the original IDW interpolations. Postremediation (soil excavation) scanning demonstrated the successful elimination of the areas of greatest concern identified using zonal statistics.

image_pdfPrint this page/section



GRO

web document
glossaryGRO Glossary
referencesGRO References
acronymsGRO Acronyms
ITRC
Contact Us
About ITRC
Visit ITRC
social media iconsClick here to visit ITRC on FacebookClick here to visit ITRC on TwitterClick here to visit ITRC on LinkedInITRC on Social Media
about_itrc
Permission is granted to refer to or quote from this publication with the customary acknowledgment of the source (see suggested citation and disclaimer). This web site is owned by ITRC • 1250 H Street, NW • Suite 850 • Washington, DC 20005 • (202) 266-4933 • Email: [email protected] • Terms of Service, Privacy Policy, and Usage Policy ITRC is sponsored by the Environmental Council of the States.