• 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

PAH Contamination in Sediments—Uncertainty Analysis (Isatis)

Contact: Gael Plassart, ENVISOL, mailto:[email protected]

gro_case_study_6_2a

Problem Statement

Sediments have been affected because of the loading and unloading of goods carried on boats arriving at the dock in Gaspé Bay harbor. The contamination is likely to migrate from the dock to Gaspé Bay.

Site Background

Gaspé Bay harbor has operated since the beginning of the twentieth century and has had several uses including military, industrial, and commercial.

Several characterization investigations have been performed at the site since 1998. These studies have revealed PAH contamination in sediments that extends from the dock to the Gaspé Bay. In 2011, an additional investigation was conducted throughout the site and consisted of drilling 129 boreholes located in a regular mesh pattern. A first estimate of contaminated sediments was calculated by performing the Thiessen Polygons method. The total volume of contaminated sediments was estimated to be 26,900 m3.

Project Objectives

The goal of the study was to delineate the extent of PAH contamination in sediments and to analyze the uncertainties of the estimates by using geospatial methods. The results were then compared to the initial estimates calculated with the Thiessen Polygons method.

Data Set

gro_case_study_6_2_b_corrected

This study evaluates the data set obtained during the characterization investigation conducted in 2011. In total, 377 samples, collected at a rate of two to three per borehole, were analyzed in the laboratory. Samples were collected between 30 cm and 90 cm in depth. Minimal sampling was performed in the northwest area, and no sampling was performed in the south and east areas, at depths between 60 cm and 90 cm.

Methods

The steps of the study were conducted as follows:

  1. EDA
  2. Construction of the empirical variogram and fitting
  3. Conditional simulation

EDA

EDA was performed in order to visualize general trends and possible anisotropies. The first results highlighted three important points:

  • Deep sampling was not performed at the South and the East part of the site between 60 cm and 90 cm depth (see Figure 85).

gro_fig_6_2Figure 85. Sample locations (samples collected from 60 cm to 90 cm depth).

  • The variogram map showed a possible anisotropy following the N110° (Figure 86). An important question was: is this anisotropy associated with the behavior of the contamination in the sediments or does it correspond to a preferential direction of sampling?

gro_fig_6_3Figure 86. Variogram map.

  • The distribution of the contaminants showed a skewed histogram (Figure 87).

gro_fig_6_4Figure 87. Histogram.

Construction of the Empirical Variogram and Fitting

According to the observations revealed by the EDA, several empirical variograms and models were constructed in order to verify whether the anisotropy was intrinsically associated with the contamination’s behavior or if it corresponded to a sampling effect. The cross-validation revealed that the best variogram and model corresponded to the isotropic model (see Figure 88).

gro_fig_6_5Figure 88. Final variogram using Gaussian-transformed data.

The final variogram was constructed with the Gaussian-transformed data because the conditional simulation method requires this type of transformation.

Conditional Simulation Method

To achieve mapping of contaminated areas and to estimate the volume of contaminated sediments, the conditional simulation were chosen because:

  • The histogram of original data showed a skewed distribution. In these cases, the mean of the distribution is not representative of the data set. Using the kriging method can introduce some bias because it performs the estimates of concentrations by using local means.
  • The conditional simulation method allows to identify and quantify the uncertain volumes.

Figure 89 illustrates the probabilities of exceeding the remediation threshold of 5 ppm set for the PAHs. This chart represents the sediments showing a probability to be contaminated of 30% or more. This 3D map takes into account the contamination sources (high risk of contamination > 60%) and the uncertain areas (30–60%). Furthermore, the map illustrates the most probable volume predicted by the conditional simulation method (median of the risk curve).

gro_fig_6_6Figure 89. Probabilities of exceeding the remediation threshold of 5 ppm for PAH measurements.

The areas of uncertainty are mainly located at the southeast part of the site, where a deep sampling was not performed.

The estimate volume of contaminated sediments, 43,440 m3, was calculated using the conditional simulation and is shown in Figure 90. The uncertain volume is 14,500 m3.

gro_fig_6_7Figure 90. Estimated volume of contaminated sediments.

Results

The volume of contaminated sediments obtained by the Thiessen polygons method was estimated to be 26,900 m3. Comparing this result to the estimate performed by the conditional simulation, the volumes are similar without taking into account the uncertain volume (see Table 11).

Table 11. Comparison of estimated volumes using Thiessen polygons and conditional simulation

Comparison of estimated volumes by performing different methods
Estimated volume of contaminated sediments by performing geostatistics 43,440 m3
Uncertain volume 14,500 m3
Estimated volume of contaminated sediments by performing polygons 26,900 m3

The areas that have not been sampled have not been taken into account for the estimates of contaminated sediments. The lack of analysis, however, does not imply the absence of contamination.

Outcome

Missing the uncertainties associated with each stage of a characterization project may lead to unreliable estimates of contamination. The quantification and localization of uncertain areas allow accurate assessment of the financial risk of the project and, if needed, definition of a complementary sampling characterization appropriate for each area.

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.