Numerical Groundwater Modelling Report
Numerical Groundwater Modelling Report
APPLICATION FOR AN
ENVIRONMENTAL ASSESSMENT CERTIFICATE /
ENVIRONMENTAL IMPACT STATEMENT
ASSESSMENT OF POTENTIAL ENVIRONMENTAL EFFECTS
Appendix 5.3.5A
Numerical Groundwater Monitoring
Report
Section 5
NEW GOLD INC.
BLACKWATER GOLD PROJECT
PREPARED FOR:
New Gold Inc.
Suite 1800, Two Bentall Centre
555 Burrard Street
Vancouver, BC, V7X 1M9
PREPARED BY:
Knight Piésold Ltd.
Suite 1400 – 750 West Pender Street
Vancouver, BC V6C 2T8 Canada
p. +1.604.685.0543 • f. +1.604.685.0147
VA101-457/6-13
Knight Piésold
Rev 0 CONSULTING
January 17, 2014 www.k n i g h t p i e s o l d .com
NEW GOLD INC.
BLACKWATER GOLD PROJECT
EXECUTIVE SUMMARY
Knight Piésold Ltd. (KP) was retained by New Gold Inc. (New Gold) to provide a representation of
baseline groundwater conditions and to evaluate potential effects of the Blackwater Gold Project (the
Project) on hydrogeological conditions. To achieve this objective, a three-dimensional steady-state,
regional-scale numerical groundwater model was developed using MODFLOW-SURFACT to
simulate baseline hydrogeological conditions at the Project site. The baseline model was then
modified to include proposed mine facilities in order to assess hydrogeological conditions during
mine operations, mine closure, and the post-closure period.
The steady-state baseline model was calibrated to average annual hydrogeologic conditions. The
modelled area encompasses the project site and surrounding drainages, including Davidson Creek,
Turtle Creek, Creek 661, Creek 705 and a portion of Chedakuz Creek. The calibrated baseline
model was then modified to create three numerical Mine Effects Models representing key phases of
Project development:
A transient Operations Model
A transient Closure/Post-Closure model, and
A steady state Post-Closure Model.
Baseline Model and Calibration
The baseline model was calibrated to average annual groundwater elevations at 18 on-site
groundwater monitoring wells, 22 vibrating wire piezometers and to estimates of average annual
baseflow at 14 locations within the study area. Baseflow estimates were obtained from the results of
a baseline watershed model developed for the Project (KP 2013f). The baseline model was
successfully calibrated by iteratively adjusting hydraulic conductivity and groundwater recharge
values until a suitable match between observed and simulated conditions was achieved. Recharge
applied to the calibrated baseline model was varied according to the distribution of surficial materials
and elevation. The calibrated model achieved groundwater levels with a normalized root mean
square error of 2%.
The simulated baseline water table generally mimics the surface topography with groundwater
elevations ranging from 1,760 meters above sea level (masl) in the high elevation region west of the
mine site to 920 masl at the downstream extent of the modelled Chedakuz Creek. Within the active
model domain, groundwater recharge occurs along topographic highs and flows to groundwater
discharge zones located within the valleys.
Mine Effects Models and Predicted Effects on Hydrogeological Conditions
Proposed major mine facilities were represented in the mine effects models, including the open pit,
Tailings Storage Facility (TSF), east and west waste rock dumps, a Low-Grade Ore (LGO) stockpile,
a freshwater reservoir and seepage collection measures consisting of an Environmental Control Dam
(ECD), groundwater interception trenches and engineered drainage ditches. The TSF facility was
represented using boundary conditions specified on the top layer of the model and included discrete
representation of the supernatant ponds, tailings beach and TSF embankments.
From model results, reductions in groundwater flow contributing to the Blackwater River catchment
were predicted to be negligible at the end of active dewatering. Average annual baseflow
3
contribution to these Blackwater River tributary streams was estimated to be 20 m /d (0.4%) lower
than baseline conditions at the end of active dewatering (Year 13) and to be the same as baseline
conditions in Post-Closure.
Open Pit and Pit Lake Simulation Results
A transient Operations Model was constructed to assess hydrogeologic conditions associated with
active mine dewatering and to simulate groundwater inflow rates to the open pit. Results of the
Operations Model indicate that simulated groundwater inflow rates to the proposed open pit during
operational dewatering are expected to increase from the start of operations through Year 13 as the
open pit increases in size and depth. Annual average inflows to the open pit were estimated to be
50 L/s with a maximum annual inflow rate of approximately 60 L/s.
During closure and after pit filling, groundwater elevations directly surrounding the Pit Lake are
expected to recover to the elevation of the Pit Lake water surface. Modelled groundwater inflow to
the Pit Lake during Post-Closure when the Pit Lake is at its maximum elevation is predicted to be
4.5 L/s. Seepage from the Pit Lake is predicted to be 1.3 L/s. Model results indicate that seepage
from the Pit Lake is expected to contribute to the Davidson Creek and Creek 661 watersheds and
not to the Blackwater River catchment.
Pit Filling Simulation Results
A transient Closure/Post-Closure Model was developed to estimate the length of time for the Pit Lake
to fill to the open pit spillway elevation, taking into account both the open pit volume and the loss of
incoming water to the groundwater system. The Pit Lake was estimated to fill 21 years after the end
of operational dewatering. This prediction assumes that water is pumped to the Pit Lake from TSF
Site D a rate of 362 L/s to assist groundwater inflows in rapid pit filling.
Seepage Pathway Assessment
A seepage analysis was conducted to assess pathways of potential seepage originating from the
proposed mine facilities including, the TSF, Pit Lake, east and west waste rock dumps and the plant
site. MODPATH particle tracking was implemented to delineate flow directions and estimate
seepage travel times to discharge locations from key mine infrastructure. Results are presented
showing the estimated groundwater seepage pathways from each facility, the discharge location of
seepage pathways, and seepage travel times. Seepage flux rates from facilities to downstream
discharge locations were estimated by generating a water budget for each facility using the
Hydrostratigraphic Unit (HSU) package in Groundwater Vistas combined with the results of the
MODPATH particle tracking.
From model results, all seepage originating from the west waste rock dump and the plant site was
predicted to be collected or discharge to drainages that flow to the TSF. Seepage was predicted to
bypass seepage collection measures and discharge to downstream locations from these facilities:
Approximately 0.4 L/s of seepage originating from TSF Site D was predicted to bypass seepage
collection measures and discharge to Davidson Creek. Approximately 0.2 L/s seepage was
predicted to discharge to Creek 661 and 0.1 L/s was predicted to discharge to the TSF spillway
channel.
Seepage of 0.5 L/s originating from the Pit Lake was predicted to travel along local groundwater
flow paths through the upper bedrock and discharge to Creek 661 with travel times of tens of
years. A small portion of seepage (0.01 L/s) was predicted to discharge to Davidson Creek
following travel paths though deep bedrock with travel times exceeding 400 years.
A seepage amount of 1.7 L/s from the east waste rock dump was predicted to flow within the
overburden and shallow bedrock under the engineered drainage ditches and discharge to
Creek 661.
The results of baseline and mine operations numerical groundwater models were used to inform
environmental effects assessment as part of the EIA submission.
TABLE OF CONTENTS
PAGE
EXECUTIVE SUMMARY........................................................................................................................ I
1 – INTRODUCTION ............................................................................................................................. 1
1.1 PROJECT DESCRIPTION .................................................................................................... 1
1.2 NUMERICAL MODELLING OBJECTIVES ........................................................................... 1
1.3 BASELINE DATA SOURCES ............................................................................................... 1
TABLES
FIGURES
APPENDICES
ABBREVIATIONS
1 – INTRODUCTION
Reconnaissance Terrain and Terrain Stability Mapping – Reconnaissance Terrain and Terrain
Stability Mapping (KP 2013a)
2012 Site Investigations – 2012 Site Investigation Report (KP 2013b)
Open Pit Investigation – Feasibility Open Pit Slope Design Report (KP 2013c)
Hydrology and Meteorology Data – 2013 Hydrometeorology Report (KP 2013d)
2013 Site Investigations – 2013 Site Investigation Report (KP 2013e)
Watershed Modelling – Watershed Modelling Report (KP 2013f)
Geotechnical Characterization – Geotechnical Characterization Report (KP 2013g), and
Open Pit Hydrogeology – Open Pit Water Management Report (KP 2013h).
Average annual streamflows estimated using the watershed model developed for the Project (KP 2013f)
were used as calibration targets in the baseline numerical model. Conceptual hydrogeologic models
developed as part of the watershed modelling study (KP 2013f) and open pit inflow assessment
(KP 2013h) were incorporated into the conceptual hydrogeologic model developed for the numerical
groundwater model in this study.
2.2 CLIMATE
Climate at the Blackwater property is sub-continental and characterized by warm summers and cold
winters. The climate is influenced by cold arctic air and moisture-laden weather systems moving west
along the Kitimat Ranges. Meterological parameters estimated for the Project have been estimated using
data collected at two climate stations in the immediate project area and correlated with data from regional
climate stations, and estimates based on watershed modelling conducted for the Project (KP 2013f).
Mean monthly temperatures range from -7.7°C in January to 12.5°C in July at the project elevation of
1,470 masl (KP 2013d). Watershed modelling results indicate the average annual precipitation calculated
from 1998 through 2012 is 640 mm at the project site elevation (KP 2013f). An equal proportion of
precipitation is estimated to fall as rain and as snow. Watershed model results also estimate the mean
annual potential evapotranspiration (PET) at 470 mm and actual evapotranspiration (AET) at 280 mm.
Average groundwater recharge across the modelled areas was estimated as 11% of total precipitation (an
equivalent area weighted average depth of 70 mm) based on the results of the watershed model
(KP 2013f).
2.3.1 Geomorphology
Surficial deposits and landforms in the project area are primarily associated with the Fraser Glaciation,
the last period of continental ice sheet glaciation in British Columbia. Surficial landforms such as
drumlins, eskers and other streamlined glacial landforms evidence that the localized ice flow direction in
the project area at the peak of glaciation was toward the northeast (KP 2013a). Glacial ice appears to
have stagnated in Davidson Creek valley during late deglaciation producing ice‐stagnation landforms
such as kettles and kames. Esker complexes are present on the north sides of Davidson Creek and the
headwaters of Creek 661. Meltwater channels provide evidence for water flow beneath and from the
margin of the receding ice sheet and major meltwater channels have a northeast trend. The valleys of
both Davidson Creek and Creek 505659 (a headwater tributary to Creek 661) contain a succession of
meltwater channels, expressed by a series of up to six terraces. The terraces provide evidence of
sequential downcutting by meltwater streams. An estimated 80% of the surficial materials in the
Davidson Creek valley is classified as lodgement glacial till (Plouffe
et. al. 2004), with the other 20% made up of ablation till, glaciofluvial, glaciolacustrine, fluvial, and organic
material. Ablation till is uncommon and predominantly found at higher elevations on the valley sides.
Bedrock exposure in the project area is rare and restricted to higher elevations. Soil cover is generally
thick within the Davidson Creek watershed and has an average thickness greater than 60 m (KP 2013g).
Bedrock is deepest along the Davidson Creek valley bottom and east of the deposit, where it was
encountered at depths of over 100 m.
Appendix Figure A.2. A summary of each unit comprising the surficial geology model in the Geotechnical
Characterization Report (KP 2013g) and Clague (2013) is summarized below.
Glacial Till Deposits - Glacial till deposits are the most dominant surficial material in the region and
consist of compact to very dense lodgement till with uncommon loose to compact ablation till. Glacial
till thickness is variable, ranging from a few to tens of metres. The material is predominantly well
graded, stiff to very dense, sandy silt to silty sand with some gravel and trace clay and cobbles.
Lodgement till is dense or stiff and contains a significant percentage of fines (silt and clay) that greatly
lowers the permeability. Ablation till is less dense and may contain less fines. Lodgement till was the
dominate material encountered in the valley basin of the Davidson Creek watershed, and ablation till
is found in a few locations on site at higher elevations on the valley sides.
Glaciolacustrine Deposits - Glaciolacustrine sediments were deposited locally in an ephemeral lake
that formed between the advancing Cordilleran ice sheet and higher ground to the west and south.
Sediment‐laden meltwater flowing along the margin of the ice sheet entered the lake, and silt-sized
particles settled out of suspension onto the lake bed. The lake was overridden by the advancing
Fraser ice sheet, which terminated glaciolacustrine deposition and compacted the deposits.
Glaciolacustrine deposits up to 20 m thick were identified in most of the drill holes in the upper
Davidson Creek valley (KP 2013g). The glaciolacustrine deposits are very dense, massive sandy silt,
and did not exhibit any fine laminated layers or weaker clayey laminations. Glaciolacustrine deposits
consist of massive silts with trace clay, sand, and poorly graded gravel. Where encountered,
glaciolacustrine layers consistently lie below the Fraser glacial till deposits.
2.3.2.5 Reworked Regolith (Reworked Completely Weathered Bedrock) and In-situ Regolith
(Completely Weathered Bedrock)
The older glacial sequence rests on reworked and in-situ regolith horizon (completely weathered
bedrock). The reworked regolith comprises poorly graded sediments containing abundant weathered
bedrock clasts. It is presumed that gravitational processes and recorded landscape instability controlled
the deposition, potentially during the onset of cold climatic conditions during the early
Pleistocene, ca. 2.6 Ma.
The original bedrock texture or fabric was evident in in-situ regolith in the majority of drill holes advanced
within Davidson Creek valley. The boundary between the reworked and in‐situ regolith is difficult to
discern in all drill holes. The reworked and in-situ regolith was found to range in thickness from a few
metres to over 30 m, with an average thickness of approximately 15 m. The regolith is thin or absent in
topographically high areas and thicker within topographic lows, indicating that the Davidson Creek
watershed may have been shielded from glacial erosion. A wide range of gradation is observed,
indicative of the various states of decomposition of the weathered bedrock. The presence of this stratum
is unusual in British Columbia, as it is typically scoured by the process of glaciation. A white to light
brown zone of silt and clay sized sediments near the top of the layer is either a soil horizon or a
weathered tuff (volcanic ash).
2.3.3 Bedrock
Geology in the project area consists of a lower unit of Upper Jurassic volcaniclastic, sedimentary, and
mafic to felsic volcanic rocks of the Bowser Lake Group. The Bowser Lake group are intruded by Late
Cretaceous granitic to granodioritic plutons. Widespread Eocene volcanic arc-related extensional felsic
volcanic rocks and minor sedimentary rocks of the Ootsa Lake Group overlie the Bowser Lake Group and
are themselves overlain on higher ridges by basalt and andesite of the Eocene Endako Group. Intact
bedrock exposure is rare and restricted to higher elevations in the area.
Bedrock geology encountered during site investigations conducted in the Davidson Creek watershed
consists of an andesite from the Cenozoic Ootsa Lake Formation and fragmentals volcanics from
Cretaceous Volcaniclastics and Flows (KP 2013b,e). Bedrock in the west part of the TSF footprint also
belongs to the Ootsa Lake Formation, but comprises rhyolites and felsic volcanic rocks. The bedrock in
the southeast portion of the project area, including the deposit area, is rhyolites and felsic volcanic rocks
of the Entiako Formation, which belongs to the Middle Jurassic Hazelton Group. Bedrock to the west is
basement rock from the Bowser Lake Group (KP 2013b,c,e).
Bedrock geology in the project area is shown on Appendix Figure A.3. The bedrock geology encountered
during site investigations conducted in the Davidson Creek watershed consists of andesite from the
Eocene Ootsa Lake Formation and fragmentals volcanics from Cretaceous Volcaniclastics and Flows
(KP 2013b,e). Bedrock in the west part of the TSF footprint also belongs to the Ootsa Lake Formation
but is comprised of rhyolites and felsic volcanic rocks.
The bedrock in the deposit area consists of rhyolites and felsic volcanic rocks of the Entiako Formation,
which belongs to the Middle Jurassic Hazelton Group (KP 2013c). Bedrock to the west of the deposit is
basement sedimentary rock from the Bowser Lake Group.
Bedrock beneath the regolith is characterized by a weathering profile that is distinguished based on rock
characteristics such as discolouration, intactness, weakness and clay content. The weathered bedrock
transition from highly to moderately to slightly weathered with increasing depth is accompanied by an
increase in intactness and decrease in clay infill (KP 2013b,e).
A bedrock elevation contour map was generated using the 2012 and 2013 site investigation data
(KP 2013g) and New Gold’s 2012 condemnation drilling data as shown on Appendix Figure A.3. Intact
bedrock exposure in the area is restricted to higher elevations. Bedrock was encountered at depths
greater than 100 m below ground surface (mbgs) in drill holes advanced on the eastern slope of
Mt. Davidson.
2.4 HYDROGEOLOGY
Baseline hydrogeology data were obtained from drilling information, in-situ hydraulic conductivity testing
(packer testing and response testing), and recorded groundwater levels at monitoring wells, standpipe
piezometers, and vibrating wire piezometers (VWP). Two pumping tests were conducted in the open pit
area: one test was conducted in an area of higher permeability bedrock within the deposit and the second
test was conducted in lower permeability bedrock south of the proposed open pit footprint (KP 2013h).
Two pumping wells and 12 observation wells with multipoint VWPs were installed to support the pumping
test program.
-11
Results of laboratory permeability testing of till samples reported permeability ranges of 10 to
-7 -8 -6
10 m/s for 27 constant head tests and 10 to 10 m/s for 7 falling head tests. A hydraulic
conductivity estimate using the mean grain size distribution of the tests and the Kozeny-Carman
-8
equation (Carman 1956) is 4x10 m/s, which falls within the range of laboratory tested values. A
-7
hydraulic conductivity value of 1x10 m/s has been adopted as a representative value for a till deposit
based on laboratory testing of the till material and grain size distributions.
Glaciolacustrine and lacustrine deposits: Glaciolacustrine sediments consist of thinly-bedded
sandy silt, gravelly silt, and silt. Glaciolacustrine sediments are present beneath the upper till layer in
Davidson Creek valley. Recent lacustrine deposits are present at the surface beneath a few isolated
lakes, such as Tatelkuz Lake. One response test was conducted in sediments identified as lacustrine
-8
(MW12-06D), which resulted in a hydraulic conductivity estimate of <1x10 m/s. Laboratory testing of
40 glaciolacustrine samples collected from drill holes and test pits for particle size analysis typically
reported 60 to 95% fines (5th to 95th percentile distribution) within the samples.
Bedrock: Bedrock within the model is assumed to be a homogeneous unit even though several types
of bedrock are present at the study site. This approach is considered sufficient for the purpose of this
hydrogeology assessment. Results of hydraulic conductivity tests conducted in bedrock are plotted
with depth on Figures 2.3 and 2.4. Bedrock at the project site is divided into units according to
weathering and intactness:
o Completely weathered bedrock: Completely weathered bedrock (in-situ regolith and reworked
regolith) consists of a silt and clay matrix with abundant weathered bedrock clasts. This unit is
inferred to be present only in Davidson Creek valley. Laboratory testing of
32 in-situ and reworked regolith samples collected from drill holes for particle size analysis
reported 20 to 85% fines (5th to 95th percentile distribution within the samples. No in-situ
hydraulic conductivity tests were conducted exclusively within this unit; completion zones of
monitoring wells screened in the completely weathered bedrock unit also spanned the contact
with the overlying sand and gravel unit. The completely weathered bedrock unit is expected to be
-8
a low permeability unit with a representative hydraulic conductivity of approximately 1x10 m/s.
o Weathered bedrock: The profile of weathering within the bedrock was distinguished based on
characteristics of rock fracture spacing, intactness, and discolouration noted on drill core
(KP 2013b,e). The weathering profile grades from highly weathered to moderately weathered to
slightly weathered with depth. Hydraulic conductivity test results conducted in weathered
bedrock within drill holes located outside the deposit area are provided in Table 2.1. Results only
include tests conducted a distance from the deposit area so the hydraulic conductivity statistics
are not influenced by tests conducted within higher permeability bedrock found in the deposit
area (see discussion below). The maximum depth of tests conducted in weathered bedrock was
generally less than 60 m below the top of the bedrock surface, although testing in one drill hole
extended to a depth of 90 m below top of bedrock. The geometric mean hydraulic conductivity
-8 -8
increases slightly as weathering in the bedrock decreases from 4x10 to 7x10 m/s. This
increase in hydraulic conductivity is attributed to a decrease in clay infill within the weathered
spaces (Deer and Patton 1971). Only response tests could be conducted in the highly weathered
bedrock zone since difficulty seating the packer prohibited packer testing within the zone. The
lower portion of the weathered bedrock is considered to be a permeable pathway and a hydraulic
-7
conductivity of 1x10 m/s is considered representative for this zone.
o Higher permeability bedrock zone in the deposit area: Results of in-situ hydraulic conductivity
testing and pumping tests indicate bedrock within the central portion of the deposit area has a
higher permeability than the surrounding bedrock. The extent of this area is estimated in the
Open Pit Water Management Report (KP 2013h) and is closely, but not exactly, related to the
‘broken zone’ in the Feasibility Open Pit Slope Design Report (KP 2013c). A bulk hydraulic
-6
conductivity of 5x10 m/s was estimated for this higher permeability bedrock zone based on the
results of a pumping test. (KP 2013h). The area of the higher permeability bedrock is inferred to
be approximately 1 km wide and to extend to depths of 500 m below the bedrock surface. The
higher permeability bedrock zone is mainly located within the limits of the proposed open pit and
will be excavated during mining operations.
o Competent bedrock: Competent bedrock is present everywhere beneath the weathered
bedrock. Hydraulic conductivity tests were conducted within deeper bedrock in the deposit area
as part of open pit site investigation (KP 2013c). The geometric mean hydraulic conductivity of
-8
45 packer tests conducted in drill holes surrounding the open pit is 8x10 m/s. A bulk hydraulic
-7
conductivity of 1x10 m/s was estimated for the lower permeability bedrock zone surrounding the
deposit based on the results of a pumping test (KP 2013h). Given the proximity of the packer
and pumping test to the higher permeability bedrock associated with the deposit, these test
results for deeper bedrock are expected to be higher than elsewhere across the project site. Drill
holes for the geotechnical site investigations were advanced to depths until two consecutive
-7
packer tests yielded an estimated hydraulic conductivity value on the order of 10 m/s or less. As
a result, tests from these site investigations are generally not considered to be part of the intact
bedrock unit. Hydraulic conductivity in the competent bedrock is assumed to decrease with
-8
depth. A bulk hydraulic conductivity of 2x10 m/s is considered to be a representative value for
the upper portion of the competent bedrock zone.
Geomechanical
2 - 129 - 5E-09 6E-05 8E-07
Investigations
5E-06
Hydrogeological
3 - - 30 3E-08 2E-05 9E-07
Investigation
Lower Permeability Bedrock Zone in the Deposit Area
Geomechanical
2 - 45 - 1E-09 4E-06 8E-08
Investigations
1E-07
Hydrogeological
3 31 4E-08 8E-06 3E-07
Investigation -
NOTES:
1. SOURCE: BLACKWATER 2012 GEOTECHNICAL SITE INVESTIGATION REPORT (KP 2013c) AND 2013 GEOTECHNICAL
SITE INVESTIGATION REPORT (KP 2013e).
2. SOURCE: FEASIBILITY OPEN PIT SLOPE DESIGN REPORT (2013c).
3. SOURCE: OPEN PIT WATER MANAGEMENT REPORT (KP 2013h).
4. HYDRAULIC CONDUCTIVITY TEST RESULTS REPORTED AS NO TAKE OR BELOW MEASURABLE TESTING LIMIT ARE
INCLUDED IN THE TABLE AS 1x10-8 m/s FOR RESPONSE TESTS AND 1x10-9 m/s FOR PACKER TESTS.
5. MINIMUM, MAXIMUM, AND LOG MEAN VALUES ARE CALCULATED USING RESPONSE, PACKER, AND AIRLIFT TESTING
RESULTS.
Figure 2.3 Hydraulic Conductivity with Depth in Bedrock (Inferred Higher Permeability Zone in
the Deposit Area)
NOTES:
1. HYDRAULIC CONDUCTIVITY PACKER AND AIRLIFT TESTING RESULTS FROM KP 2013b,c,e,h.
2. DEPTHS ADJUSTED TO METERS BELOW GROUND LEVEL FOR INCLINED DRILL HOLES (KPL 2013c).
3. HYDRAULIC CONDUCTIVITY TEST RESULTS REPORTED AS NO TAKE ARE INCLUDED IN THE TABLE AS 1x10-9 m/s.
Figure 2.4 Hydraulic Conductivity with Depth in Bedrock (Inferred Lower Permeability Zone in
the Deposit Area)
NOTES:
1. HYDRAULIC CONDUCTIVITY PACKER AND AIRLIFT TESTING RESULTS FROM KP 2013b,c,e,h.
2. DEPTHS ADJUSTED TO METERS BELOW GROUND LEVEL FOR INCLINED DRILL HOLES (KPL 2013c).
3. HYDRAULIC CONDUCTIVITY TEST RESULTS REPORTED AS NO TAKE ARE INCLUDED IN THE TABLE AS 1x10-9 m/s.
3.1 OVERVIEW
A steady-state, regional-scale numerical groundwater model was developed to simulate baseline
hydrogeological conditions and to provide the basis required to assess potential effects of the Project on
the local groundwater system. The model was developed using the MODFLOW-SURFACT computer
code run in the Groundwater Vistas (version 6.20; ESI, 2011) graphical user interface. MODFLOW-
SURFACT is a three-dimensional finite-difference flow model developed by the U.S. Geological Survey
and HGL Software Systems that has become an industry standard for groundwater modelling applications
(Hydrogeologic Inc., 1996).
The baseline model was calibrated to hydraulic head data collected from on-site groundwater monitoring
wells and vibrating wire piezometers and to average annual baseflows (groundwater contribution to
streamflow) estimated using a watershed model constructed for the Project (KP 2013f).
The baseline model simulates pre-development hydrogeological conditions including groundwater flow
directions, distribution of hydraulic head and groundwater/surface water interaction on a project-site
scale. Baseline model development, calibration and results are discussed in the sections that follow.
370,000
380,000
390,000
5,910,000 5,910,000
1-TC
1-DC
BWL
4-DC
Ta
H6
te
lku
eek
z
Cr
La
rtle ek
ke
Tu re
C H4B
on
ds
H3
vi
a
D
e MW12-09S
ak 700
5,900,000 pL Creek 5,900,000
To
1-661
MW12-07D/S 661
MW12-08D/S Creek
H2 1-505659
MW12-05D/S MW12-12D/S H1
MW12-06S
k
ee
11-DC MW12-02D/S
Cr
4-705
BWH
nie
aw
705 PH13-1-3
ek MW12-13D/S
Cre
F
PH13-1-1 PH12-2-1
PH12-4-3
H7 PH12-3-2
1-705
MW12-04D
MW12-11D/S
MW12-10D
PH12-4-2
5,890,000 5,890,000
SAVED: M:\1\01\00457\06\A\GIS\Figs\Report13_NumericalGroundwaterModellingeport\Fig301_Modflow.mxd; Jan 17, 2014 9:18 AM; cczembor
5,880,000
360,000
370,000
380,000
390,000
2 1 0 2 4 6 8 km
SCALE
Elevation (masl)
Cross-Section Along Row 150
1800
1600
Cross-Section Along Row 150
1400
1200
1000
800
600
400
200
Elevation (masl) 0
0 5000 10000 15000 20000 25000 30000
2000
1800
1600
1400
1200
1000
0 0
5000 5000
Model Y-Coordinate (m)
10000 10000
20000 20000
NOTES: NEW GOLD INC.
Boundary Conditions 1. AS INDICATED BY THE NORTH ARROW, THE MODEL GRID IS
ROTATED FROM TRUE NORTH SUCH THAT COLUMNS ALIGN WITH AN BLACKWATER GOLD PROJECT
AZIMUTH OF 325 DEGREES.
BASELINE MODEL
2. PROPOSED MINE FACILITY OUTLINES ARE SHOWN IN RED FOR NUMERICAL GRID, LAYERS AND BOUNDARY
SPATIAL REFERENCE ONLY AND ARE NOT SIMULATED IN THE CONDITIONS (LAYER 1)
BASELINE MODEL.
P/A NO. REF. NO.
VA101-457/6-13 13
Hydraulic
Specific Specific Effective
Hydrostratigraphic Unit MODFLOW Conductivity
Layer Kx,y Kz Storage Yield Porosity
(m/s) (m/s) (1/m) (-) (-)
Glaciofluvial: Channel 1 1E-04 1E-04 1E-04 0.3 0.15
Glaciofluvial: Kame 1 1E-05 1E-05 1E-04 0.3 0.15
Glacial Till 1,2 2E-07 2E-07 1E-04 5E-02 0.15
Glaciolacustrine 1,2 1E-07 2E-08 1E-04 5E-02 0.15
Kame/Till Deposit 1 9E-06 2E-07 1E-04 5E-02 0.15
Completely Weathered Bedrock 3 1E-08 1E-08 1E-06 1E-03 0.001
Weathered Bedrock 3 or 4 1E-07 1E-07 1E-06 1E-03 0.001
Competent Bedrock 4, 5 2E-08 2E-08 1E-06 1E-04 0.0001
Competent Bedrock 6,7,8 1E-08 1E-08 1E-06 1E-04 0.0001
Competent Bedrock 9,10 1E-09 1E-09 1E-06 1E-04 0.0001
Higher Permeability Bedrock Zone in
Deposit 3 to 8 1E-06 1E-06 1E-05 5E-03 0.005
Bedrock Outcrops 1 1E-07 1E-07 1E-06 1E-03 0.001
Model Layer 1 was subdivided into six hydraulic conductivity zones based on the results of surficial
landform mapping conducted within the immediate project area (KP 2013a) and using surficial geology
maps by the Geologic Survey of Canada (Plouffe and Levson 2001, 2002) where the model extended
beyond that area. The majority of Layer 1 represented a glacial till with a hydraulic conductivity of
-7
2x10 m/s. Other units represented within Layer 1 include slightly weathered bedrock outcrops along
topographic highs and glaciofluvial kame and channel deposits within major surface drainages. The
glaciofluvial channel deposits were assumed to comprise approximately 30% of the width of permeable
sediments in major surface drainage channels in the model where detailed surficial mapping was
unavailable for the Project (KP 2013a). The kame/till deposit was specified in the model along the flanks
of Mt. Davidson in areas where a kame deposit is mapped as part of the surficial landform mapping
(KP 2013a) and where a thin kame unit (approximately 3 m thick) was encountered in test pits
(KP 2013b,e). The kame/till hydraulic conductivity zone was assigned to the model during calibration to
reduce flooding within grid cells at the base of Mt. Davidson, and to simulate a water table approximately
2 to 3 mbgs as observed during the field investigation.
Hydraulic conductivity zones for subsurface layers of the model (Layers 2 through 10) were assigned
based on available drill hole lithology and corresponding with the conceptual hydrogeological model
presented in Section 2.4. Layer 2 is predominantly glacial till everywhere except where a glaciolacustrine
deposit is present within upper Davidson Creek Valley and beneath lakes, and where bedrock is present
beneath outcrops. A plan view of hydraulic conductivity zones assigned to Layer 2 is presented on
Figure B.2 in Appendix B.
Layers 3 through 10 represent bedrock units and generally consist of weathered bedrock in Layer 3 and
competent bedrock in Layers 4 through 10. Exceptions to this include Davidson Creek valley where a
completely weathered bedrock unit was assigned in Layer 3 and weathered bedrock was subsequently
assigned in Layer 4. A higher permeability zone of bedrock is specified in the deposit area in
-6
Layers 3 through 8 with a hydraulic conductivity value of 1x10 m/s. Figures displaying the hydraulic
conductivity zones for Layers 3 through 10 are provided in Appendix B.
3.4.5 Faults
Faults acting as barriers to groundwater flow were represented in the model using the Wall boundary
condition (Horizontal Flow Barrier Package). The wall boundary condition allows thin, low permeability
vertical features to be represented in the model without refining the grid size. Faults were only specified
within the deposit area where the inclusion of a mapped fault helped to improve model calibration. The
best match to observed and simulated data was obtained using three faults located along the southern
and western edges of the deposit. Including faults in the model allowed simulated hydraulic gradients to
more closely represent the steep hydraulic gradients recorded at adjacent VWPs along the south edge of
the pit. Wall boundaries representing faults were specified in model Layers 3 through 10 and were
-9
assigned a horizontal hydraulic conductivity of 1x10 m/s.
run off estimated by the results of watershed modelling (KP 2013f). Recharge to the lower till zone was
subsequently assigned a lower value to reduce flooding within lower elevations of the model.
Initial groundwater recharge values were assigned based on surficial material type and watershed
modelling (KP 2013f) which suggested that approximately 10% of annual precipitation (an areal weighted
average of 70 mm) provides recharge to groundwater. These initial values were varied within an
expected range during calibration of the baseline model. The calibrated groundwater recharge values
assigned to each zone of the model are summarized in Table 3.2.
5000
0.8
(m)
Elevation(m)
Model Y-Coordinate
0.6
10000
0.4
#REF!
0.2
15000
0
0 0.2 0.4 0.6 0.8 1 1.2
Area (ha)
20000
Horizontal Hydraulic Vertical Hydraulic
Zone Color Material Type
Conductivity Kh (m /sec) Conductivity Kv (m /sec)
Zone 6 Glaciofluvial: Channel 1.E-04 1.E-04
Zone 2 Glaciofluvial: Kame 1.E-05 1.E-05 NEW GOLD INC.
Zone 1 Glacial Till 2.E-07 2.E-07
Zone 4 Glaciolacustrine 1.E-07 2.E-08 BLACKWATER GOLD PROJECT
Zone 12 Kame/Till Deposit 9.E-06 2.E-07
Zone 11 Bedrock Outcrops 1.E-07 1.E-07 BASELINE MODEL
NOTES: HYDRAULIC CONDUCTIVITY ZONES (LAYER 1)
1. THE CONTOUR LINES SHOWN ABOVE ILLUSTRATE GROUND SURFACE ELEVATION (masl) AS DEFINED IN MODEL
LAYER 1.
2. THE GLACIOLACUSTRINE AND GLACIOFLUVIAL KAME/TILL DEPOSITS HAVE ANISOTROPY RATIOS (Kh/Kv) OF 5 AND
P/A NO. REF. NO.
VA101-457/6-13 13
45, RESPECTIVELY.
0 16JAN'14 ISSUED WITH REPORT KTD CAS KJB REV
REV DATE DESCRIPTION PREP'D CHK'D APP'D
FIGURE 3.3 0
0 5000 10000
Model15000
X-Coordinate (m) 20000 25000 30000
0
5000
Model Y-Coordinate (m)
10000
15000
20000
Groundwater Recharge
Zone Color Material Type
(mm/year)
Zone 1 Till - Lower 22
Zone 2 Till - Mid 40 NEW GOLD INC.
Zone 3 Till -Upper 80
Zone 4 Lacustrine 22 BLACKWATER GOLD PROJECT
Zone 5 Bedrock Outcrop 80
Zone 6 Glaciofluvial - Channel 190 BASELINE MODEL
Zone 7 Glaciofluvial - Kame 190 GROUNDWATER RECHARGE ZONES (LAYER 1)
NOTES:
1. THE CONTOUR LINES SHOWN ABOVE ARE OF GROUND SURFACE P/A NO. REF. NO.
VA101-457/6-13 13
ELEVATION (masl) AS DEFINED IN MODEL LAYER 1.
0 16JAN'14 ISSUED WITH REPORT KTD CAS KJB REV
REV DATE DESCRIPTION PREP'D CHK'D APP'D FIGURE 3.4 0
NEW GOLD INC.
BLACKWATER GOLD PROJECT
hydraulic head of 12 to 15 m from recorded. The goal of model calibration was not to reproduce the
hydraulic head variation at all observation points within the deposit area, but to obtain a general
representation of hydrogeologic conditions. The match to hydraulic heads within the deposit area is
considered suitable for the purpose of the model.
Measured Simulated
Groundwater Groundwater Residual Head
Well I.D.
Elevation Elevation (m)
(masl) (masl)
MW12-02D 1,401 1,405 -4.1
MW12-02S 1,406 1,405 0.8
MW12-04D 1,556 1,565 -8.9
MW12-05D 1,374 1,378 -4.1
MW12-05S 1,373 1,377 -4.4
MW12-06S 1,258 1,261 -2.9
MW12-07D 1,203 1,208 -4.7
MW12-07S 1,204 1,208 -4.3
MW12-08D 1,160 1,156 4.3
MW12-08S 1,158 1,156 1.8
MW12-09S 1,144 1,142 1.8
MW12-10D 1,657 1,663 -5.6
MW12-11D 1,668 1,656 11.9
MW12-11S 1,668 1,659 9.3
MW12-12D 1,243 1,237 6.3
MW12-12S 1,243 1,237 6.4
MW12-13D 1,368 1,365 3.3
MW12-13S 1,359 1,365 -5.5
PH13-1-1_VWP1 1,524 1,514 9.8
PH13-1-1_VWP2 1,525 1,514 10.5
PH13-1-3_VWP1 1,517 1,513 3.6
PH13-1-3_VWP2 1,517 1,513 4.2
PH13-1-3_VWP3 1,518 1,513 4.7
PH12-2-1_VWP2 1,514 1,514 -0.3
PH12-2-1_VWP3 1,518 1,514 3.5
PH12-3-2_VWP1 1,624 1,622 1.8
PH12-3-2_VWP2 1,629 1,622 7.7
PH12-3-2_VWP3 1,633 1,621 12.3
PH12-3-2_VWP4 1,633 1,621 12.0
PH12-3-2_VWP5 1,634 1,619 15.4
PH12-4-2_VWP1 1,611 1,616 -5.4
PH12-4-2_VWP2 1,612 1,616 -4.1
PH12-4-2_VWP3 1,616 1,619 -2.5
PH12-4-2_VWP4 1,617 1,620 -2.7
PH12-4-2_VWP5 1,619 1,620 -1.4
PH12-4-3_VWP1 1,525 1,520 4.6
PH12-4-3_VWP2 1,524 1,520 3.8
PH12-4-3_VWP3 1,535 1,532 2.5
PH12-4-3_VWP4 1,597 1,559 38.4
PH12-4-3_VWP5 1,610 1,559 51.2
RMSE (m) 11.9
NRMSE 2%
NOTES:
1. THE VALUES LISTED AS MEASURED GROUNDWATER ELEVATIONS AT MONITORING WELLS (MW) ARE AVERAGE
ELEVATIONS OVER THE PERIOD OF AVAILABLE MEASURED DATA AND ARE MEASUREMENTS FROM MARCH 2013 AT
VWPS INSTALLED IN OBSERVATION WELLS (‘PH’ SERIES).
2. RMSE = ROOT MEAN SQUARE ERROR; NRMSE = NORMALIZED ROOT MEAN SQUARE ERROR.
NOTES:
1. WATERSHED MODEL VALUES OBTAINED FROM WATERSHED MODELLING REPORT (KP 2013f).
2. BOLD INDICATES THE DOWNSTREAM MOST STATION ON EACH DRAINAGE COURSE.
3. RPD = RELATIVE PERCENT DIFFERENCE.
1,500
MW12-11S
MW12-12D
MW12-12S
MW12-13D
MW12-13S
1,400 PH13-1-1_VWP1
PH13-1-1_VWP2
PH13-1-3_VWP1
Model Underprediction PH13-1-3_VWP2
PH13-1-3_VWP3
PH12-2-1_VWP2
1,300 PH12-2-1_VWP3
PH12-3-2_VWP1
PH12-3-2_VWP2
PH12-3-2_VWP3
PH12-3-2_VWP4
PH12-3-2_VWP5
1,200 PH12-4-2_VWP1
PH12-4-2_VWP2
PH12-4-2_VWP3
PH12-4-2_VWP4
PH12-4-2_VWP5
PH12-4-3_VWP1
1,100 PH12-4-3_VWP2
PH12-4-3_VWP3
1,100 1,200 1,300 1,400 1,500 1,600 1,700 PH12-4-3_VWP4
Measured Hydraulic Head (masl)
NEW GOLD INC.
BLACKWATER GOLD PROJECT
BASELINE MODEL
OBSERVED VS. SIMULATED HYDRAULIC HEAD
125
SECTION ALONG COL
5000
Model Y-Coordinate (m)
10000
20000
Water Table Elevation (masl)
Boundary Conditions
NEW GOLD INC.
NOTES:
BLACKWATER GOLD PROJECT
1. THE CONTOUR LINES SHOWN ABOVE ARE OF MODEL
SIMULATED WATER TABLE ELEVATION (masl).
BASELINE MODEL
2. PROPOSED MINE FACILITY OUTLINES ARE SHOWN FOR SPATIAL SIMULATED WATER TABLE CONTOUR MAP
REFERENCE ONLY AND ARE NOT SIMULATED IN THE BASELINE
MODEL. P/A NO. REF. NO.
VA101-457/6-13 13
1600
1400
1200
1000
800
600
400
200
0
0 5000 10000 15000 20000
1400
1200
1000 SECTION ALONG ROW
800 150
600
400
200
0
0 5000 10000 15000 20000 25000 30000
Model X-Coordinate (m)
NOTES: NEW GOLD INC.
1. CONTOURS REPRESENT HYDRAULIC HEAD EQUIPOTENTIALS. BLACKWATER GOLD PROJECT
2. ARROWS INDICATE GROUNDWATER FLOW DIRECTION BASED
BASELINE MODEL
ON THE SLOPE OF THE EQUPOTENTIAL LINES.
SIMULATED HYDRAULIC HEAD DISTRIBUTION
(CROSS-SECTION)
P/A NO. REF. NO.
VA101-457/6-13 13
4.1 OVERVIEW
A transient Operations Model was developed to assess potential effects of the Blackwater Project on pre-
development hydrogeological conditions. The Operations Model was developed from the calibrated
baseline groundwater model using MODFLOW-SURFACT and Groundwater Vistas.
The objectives of the operations groundwater modelling were to:
Characterize potential effects of mine facilities on baseline hydrogeology during the operational
period
Estimate groundwater inflow rates to the open pit on a yearly basis, and
Delineate the groundwater capture zone surrounding the open pit at its maximum extent.
The results of the transient operations modelling along with the methodology and assumptions used to
develop the models are presented in the sections that follow.
be placed adjacent to the open pit during operations and will be reclaimed by the end of mine operations
(Year 17).
An Environmental Control Dam (ECD) and groundwater interception trenches are planned approximately
1 km downstream of the Site D Main Dam to recover potential seepage from the TSF. One seepage
interception trench will be constructed on each side of Davidson Creek, which will be excavated through
the surficial sand and gravel terraces to LPS downstream of the Site D Main Dam. Seepage to the
collection trenches will report to the ECD pond. Recovered water will be pumped to TSF Site D and the
collection pond will be maintained in a dewatered condition to the maximum extent practical.
Engineered drainage ditches will be constructed downslope of the east and west waste rock dumps and
LGO stockpile to collect surface runoff and shallow groundwater seepage from facilities. Drainage
ditches will direct water to the TSF.
Fresh water for the project will be sourced from Tatelkuz Lake, which is located approximately 20 km
northeast of the mine site. Freshwater for streamflow mitigation in Davidson Creek will be stored in a
Freshwater Reservoir located approximately 800 m downstream of the ECD.
Elevation (masl)
Cross-Section Along Row 150
1800
1600
Cross-Section Along Row 150
1400
1200
1000
800
600
400
200
Elevation (masl) 0
0 5000 10000 15000 20000 25000 30000
2000
1800
1600
1400
1200
1000
0 0
5000 5000
Model Y-Coordinate (m)
10000 10000
20000 20000
Boundary Conditions NOTES: NEW GOLD INC.
1. AS INDICATED BY THE NORTH ARROW, THE MODEL GRID IS BLACKWATER GOLD PROJECT
ROTATED FROM TRUE NORTH SUCH THAT COLUMNS ALIGN WITH AN
AZIMUTH OF 325 DEGREES. OPERATIONS MODEL
NUMERICAL GRID, LAYERS AND BOUNDARY
2. THE ECD SYSTEM IS REPRESENTED WITH DRAIN BOUNDARY
CONDITIONS (LAYER 1)
CONDITIONS AND IS HIGHLIGHTED IN LIGHT BLUE DOWNSTREAM FROM
THE TSF MAIN EMBANKMENT. P/A NO. REF. NO.
VA101-457/6-13 13
elevation. Conductance of the drain cells was calculated using the hydraulic conductivity value assigned
-5
to the interception trenches (1x10 m/s).
Groundwater drawdown was predicted around the camp water supply wells associated with groundwater
extraction through the operational period. A maximum drawdown of approximately 5 m was expected
immediately adjacent to the TW13-01 and TW13-02 well locations, as shown on Figure 4.2. Groundwater
drawdown associated with the camp wells was predicted to intersect with the 1 m drawdown zone of
influence of the open pit.
From model results, reductions in groundwater flow contributing to the Blackwater River catchment were
predicted to be negligible at the end of active dewatering. Groundwater contribution (baseflow) to
3
tributaries of Blackwater River included in the model were 5,585 m /d in baseline conditions. Reductions
in average annual baseflow contribution to these Blackwater River tributary streams were estimated to be
3
20 m /d (0.25 L/s), which is equivalent to a 0.2% decrease in average annual baseflows to Blackwater
River tributaries in the model. Groundwater flows leaving the model domain toward Blackwater River via
the General Head Boundary remained unchanged from baseline conditions.
Simulated groundwater inflows to the open pit during each year of operations are provided on Figure 4.4.
Groundwater inflow rates are plotted for the results of the Operations Model along with estimates based
on an analytical calculation presented in the Pit Water Management Report (KP 2013h). The maximum
groundwater inflow rate is almost 60 L/s with an average inflow for Year 3 through Year 13 of
approximately 50 L/s. The analytical calculation was used to design the dewatering system and includes
groundwater pumping rates required to achieve slope depressurization during the first few years of pit
dewatering. Predicted groundwater inflows to the open pit using the numerical groundwater model do not
consider requirements for slope depressurization.
5000
0.8
Elevation
Model Y-Coordinate (m)
(m)
0.6
10000
0.4
#REF!
0.2
15000
0
0 0.2 0.4 0.6 0.8 1 1.2
Area (ha)
20000
NOTES:
Boundary Conditions 1. CONTOUR INDICATE MODEL SIMULATED WATER TABLE NEW GOLD INC.
ELEVATION (masl).
BLACKWATER GOLD PROJECT
2.THE ECD SYSTEM (DRAIN BOUNDARY CONDITIONS)
ARE HIGHLIGHTED IN LIGHT BLUE DOWNSTREAM FROM OPERATIONS MODEL
THE MAIN TSF EMBANKMENT. GROUNDWATER ELEVATION MAP
3. EMBANKMENT DRAINS ARE NOT ASSIGNED WITHIN
(YEAR 13)
THE FOOTPRINTS OF THE TSF EMBANKMENTS. P/A NO. REF. NO.
VA101-457/6-13 13
TW13-01
Boundary Conditions
Groundwater Extraction Well NEW GOLD INC.
NOTES: BLACKWATER GOLD PROJECT
1. CONTOURS INDICATE PREDICTED WATER TABLE DRAWDOWN (meters) OPERATIONS MODEL
FROM PRE-DEVELOPMENT CONDITIONS. PREDICTED OPEN PIT DRAWDOWN
(YEAR 13)
P/A NO. REF. NO.
VA101-457/6-13 13
60
50
Groundwater Inflow Rate (L/S)
40
30
20
10
0
‐2 ‐1 1 2 3 4 5 6 7 8 9 10 11 12 13
Analytical Pit Inflow Predictions (L/S) Year of Operations
MODFLOW Simulated Pit Inflows (L/S)
NEW GOLD INC.
NOTES: BLACKWATER GOLD PROJECT
1. PIT INFLOW RATES PREDICTED FOR YEARS ‐2 THROUGH 2 USING THE ANALYTICAL METHOD OPERATIONS MODEL
ARE HIGHER THAN THOSE PREDICTED USING THE NUMERICAL MODEL. THE ANALYTICAL PREDICTED PIT INFLOW RATES: NUMERICAL MODEL
METHOD CONSIDERS THE GROUNDWATER PUMPING REQUIREMENTS FOR PIT DEWATERING AND AND ANALYTICAL CALCULATION
PIT WALL DEPRESSURIZATION. THE NUMERICAL MODEL SIMULATIONS DO NOT INCLUDE PUMPING
REQUIREMENTS TO ACHIEVE DEPRESSURIZATION. P/A NO. REF. NO.
VA101-457/6-13 13
0 16JAN'14 ISSUED WITH REPORT KTD CAS KJB
REV
REV DATE DESCRIPTION PREP'D CHK'D APP'D FIGURE 4.4 0
NEW GOLD INC.
BLACKWATER GOLD PROJECT
5 – CLOSURE/POST-CLOSURE SIMULATION
5.1 OVERVIEW
A transient Closure/Post-Closure Model was developed for the Project in order to simulate filling of the
open pit during the closure phase of mining. The pit void will fill with water contributed by groundwater
inflows and water pumped from the TSF D Pond.
The Closure/Post-Closure Model was developed using MODFLOW-SURFACT and Groundwater Vistas.
The model was created by modifying the transient Operations Model to represent the pit as a void
capable of filling with water. The results of the transient closure/post-closure modelling along with the
methodology and assumptions used to develop the models are presented in the sections that follow.
from the simulated pit via spillway drains. The estimated duration of pit filling corresponds well with an
estimate of 20 years previously completed using GoldSim (KP 2013i).
0.000007
0.000006
5000
0.000005
Elevation
Model Y-Coordinate (m)
(m)
0.000004
10000
0.000003
15000
0.000001
0
0 0.5 1 1.5 2 2.5 3 3.5
Area (ha)
20000
NOTES:
Boundary Conditions 1. THE CONTOUR LINES SHOWN ABOVE ARE OF MODEL
SIMULATED WATER TABLE ELEVATION (masl). NEW GOLD INC.
2.THE ECD SYSTEM (DRAIN BOUNDARY CONDITIONS) BLACKWATER GOLD PROJECT
ARE HIGHLIGHTED IN LIGHT BLUE DOWNSTREAM FROM
THE MAIN TSF EMBANKMENT. CLOSURE/POST-CLOSURE MODEL
3. THE ULTIMATE PIT LAKE ELEVATION ASSIGNED TO THE GROUNDWATER ELEVATION MAP
MODEL WAS 1475 masl. (YEAR 55)
3. EMBANKMENT DRAINS ARE NOT ASSIGNED WITHIN P/A NO. REF. NO.
VA101-457/6-13 13
THE FOOTPRINTS OF THE TSF EMBANKMENTS.
0 16JAN'14 ISSUED WITH REPORT KTD CAS KJB REV
REV DATE DESCRIPTION PREP'D CHK'D APP'D FIGURE 5.1 0
250,000,000
200,000,000
Pit Lake Water Storage Volume (m3)
150,000,000
100,000,000
50,000,000
0
0 5 10 15 20 25 30 35 40
Time (Years After End of Operational Dewatering)
6 – POST-CLOSURE SIMULATION
6.1 OVERVIEW
A steady-state Post-Closure Model was developed to characterize potential seepage pathways from key
mine infrastructure during the post-closure period. The Post-Closure Model was developed by modifying
the Closure/Post-Closure Model to represent the elevation of the Pit Lake using a constant head
boundary condition.
The main objectives of the steady-state post-closure modelling were to:
Estimate total seepage rates from the Pit Lake during post-closure
Delineate potential seepage pathways from key mine facilities, including the Pit Lake, TSF, east and
west waste rock dumps, and the plant site using MODPATH, and
Estimate seepage travel times to downstream discharge locations using MODPATH and Endpoint
analysis.
The results of the steady-state post-closure modelling along with the methodology and assumptions used
to develop the model are presented in the sections that follow.
5000
0.8
(m)
Elevation(m)
Model Y-Coordinate
0.6
10000
0.4
#REF!
0.2
15000
0
0 0.2 0.4 0.6 0.8 1 1.2
Area (ha)
20000
Boundary
Conditions NOTES: NEW GOLD INC.
1. THE CONTOUR LINES SHOWN ABOVE ARE OF MODEL BLACKWATER GOLD PROJECT
SIMULATED WATER TABLE ELEVATION (masl).
2.THE ECD SYSTEM (DRAIN BOUNDARY CONDITIONS) STEADY-STATE POST-CLOSURE
ARE HIGHLIGHTED IN LIGHT BLUE DOWNSTREAM FROM GROUNDWATER ELEVATION MAP
THE MAIN TSF EMBANKMENT.
P/A NO. REF. NO.
VA101-457/6-13 13
0 16JAN'14 ISSUED WITH REPORT KTD CAS KJB REV
REV DATE DESCRIPTION PREP'D CHK'D APP'D FIGURE 6.1 0
NEW GOLD INC.
BLACKWATER GOLD PROJECT
6.4.1 Seepage Flow Directions and Travel Times (MODPATH Particle Tracking)
MODPATH particle tracking was implemented to delineate flow directions and estimate seepage travel
times to discharge locations from key mine infrastructure. The following facilities were included in the
MODPATH analysis:
TSF facility, including Site C and Site D
Pit Lake
East and west dumps, and
Plant site.
For the Pit Lake MODPATH analysis, particles were inserted along the downstream rim of the pit shell in
Layers 1 through 8 and within the footprint of the Pit Lake constant head cells in Layer 8. For MODPATH
analyses at all other facilities, particles were inserted at the top of Layer 1 within the facility footprint.
Particles were forward tracked through the groundwater flow system to downstream discharge locations.
Downstream discharge locations are shown on Figure 6.2.
MODPATH combined with Endpoint Analysis was used to determine the discharge location of potential
seepage from each facility. Endpoint analysis allows the user to identify steady-state flow lines that
terminate (discharge) at a cell along the model boundary. The MODPATH simulation can be used to
calculate approximate groundwater travel times along the seepage pathways by taking into consideration
an assumed effective porosity. Effective porosities assigned to the model for the MODPATH velocity
calculations are shown on Table 3.1 and include 0.1% (0.001) for weathered bedrock, 0.01% (0.0001) for
competent bedrock and 15% (0.15) for overburden material. Travel times are representative of advective
transport and do not include effects from dispersion or diffusion.
MODPATH results are sensitive to specification of the “sink strength” input parameter, which defines the
termination criterion for particle traces flowing through boundary cells. All MODPATH scenarios
presented herein adopt a “stop at 50 percent strength” weak sink option to discontinue particle traces in
boundary cells. Conceptually this means that the model terminates a particle trace in a cell if more than
50% of the water in the cell is removed.
The simulated MODPATH particle traces resulting from each simulation are provided in Appendix D along
with figures showing the HSU zones defining the simulated discharge zone used to assess the seepage
rates for each mine facility. Histograms showing the distribution of particle travel times from each facility
to downstream locations are provided in Appendix E.
Results of the seepage analysis for each facility are described in the sections that follow.
4
12
10
11
5 14
5 Lower East Dump Drainage Ditch
6 Natural Channel Reporting to TSF
7 Pit Lake
8 Springs within East Dump Footprint
7 9 Springs within West Dump Footprint
8 10 TSF Main Embankment Drains
11 TSF Site C Pond
Table 6.1 Results of MODPATH Particle Tracking and Advective Travel Times
NOTES:
1. APPROXIMATE SEEPAGE TRAVEL TIMES FROM THE STOCKPILES TO THE DISCHARGE LOCATIONS WERE CALCULATED USING AN
ASSUMED EFFECTIVE POROSITY OF 0.1% (0.001) FOR WEATHERED BEDROCK, 0.01% (0.0001) FOR UNWEATHERED BEDROCK AND
15% (0.15) FOR ALLUVIAL MATERIALS.
2. TRAVEL TIMES ARE PROVIDED FOR THE STEADY-STATE POST-CLOSURE MODEL
3. RAVEL TIMES ARE BASED ON ADVECTIVE TRAVEL ONLY AND DISREGARD THE EFFECTS OF DISPERSION AND DIFFUSION.
4. PARTICLES WITH LONG TRAVEL TIMES ORIGINATE NEAR THE BOUNDARY OF THE TSF SITE D POND DISCHARGE ZONE (SHOWN ON
FIGURE 3). THE PROXIMITY TO THIS BOUNDARY MAY INFLUENCE THE MODELLED TRAVEL TIMES.
7.1 CONCLUSIONS
Based on the data analysis, model calibration and model simulations, the following conclusions can be
drawn from the results of the numerical groundwater flow models:
Groundwater in the southwestern portion of the deposit area behaves in a compartmentalized
manner. The best fit to modelled data was achieved by specifying barriers to groundwater flow
(faults) along the edges of the deposit and a higher permeability bedrock zone in the central portion of
the deposit. Identification of this compartmentalization was the result of and supported by increased
data density in this area. Groundwater compartmentalization is likely present at most locations on the
site and was included in the model calibration process.
Mine dewatering was simulated using a transient Operations Model simulating 15 years of mine
operations (Year -2 through Year 13). Drain cells were specified within the pit shell. The maximum
groundwater inflow rates predicted using the numerical groundwater model were approximately
60 L/s with an average inflow rate of 50 L/s. Numerical model results compare well with estimates
using an analytical calculation.
At the end of active dewatering (Year 13), water table drawdown of 1 m was predicted to extend an
average distance of approximately 1,200 m from the pit edge. The predicted drawdown zone of
influence is irregularly shaped and is elongated beneath the topographic high at Mt. Davidson. The
1 m drawdown contour extends approximately 600 m from the pit edge in a southeast direction
toward Blackwater River.
The groundwater drawdown zone of influence for the open pit at the end of active dewatering as
defined by 1 m drawdown contours was predicted to merge with the drawdown zone of influence from
the camp wells. The assessment was conducted using a camp water demand predicted for the
construction phase which is approximately three times greater than the predicted water demand
during operations and is therefore considered conservative.
From model results, reductions in groundwater flow contributing to the Blackwater River catchment
were predicted to be negligible at the end of active dewatering. Average annual baseflow contribution
3
to these Blackwater River tributary streams was estimated to be 20 m /d (0.4%) lower than baseline
conditions at the end of active dewatering (Year 13) and to be the same as baseline conditions in
Post-Closure.
The Pit Lake was estimated to take 21 years to fill to its maximum volume at the spillway elevation
following the cessation of dewatering. This analysis assumed that water is pumped to the Pit Lake
from TSF Site D at a rate of 362 L/s to assist in rapid pit filling.
Results of the steady-state Post-Closure seepage assessment using MODPATH particle tracking and
a mass balance analysis of facility footprint area contributing seepage to downstream discharge
locations indicated that:
o Total seepage into the foundation materials below TSF Site D was predicted to be 21 L/s. The
majority of seepage originating from TSF Site D was predicted to discharge to the TSF
embankment drains (15 L/s) and to the seepage collection system at the Environmental Control
Dam (5 L/s). Approximately 0.4 L/s of seepage originating from TSF Site D was predicted to
bypass seepage collection measures and discharge to Davidson Creek. Approximately 0.2 L/s
seepage was predicted to discharge to Creek 661 and 0.1 L/s was predicted to discharge to the
TSF spillway.
o Total seepage from the Pit Lake was estimated to be approximately 1.3 L/s and was predicted to
flow in the directions of the Davidson Creek and Creek 661 catchments. Seepage from the Pit
Lake was not predicted to flow toward the Blackwater River catchment. A large portion of the Pit
Lake seepage was predicted to discharge to drainages that flow to TSF Site D (0.7/s). Seepage
of approximately 0.5 L/s was estimated to travel along local groundwater flow paths through the
upper bedrock and discharge to Creek 661. A trace amount of seepage (0.01 L/s) was predicted
to discharge to Davidson Creek following travel paths though deep bedrock.
o Approximately 1.3 L/s of seepage from the east dump was predicted to discharge to springs and
natural channels directed to engineered drainage ditches that convey runoff and toe discharge to
the TSF. A seepage amount of 1.7 L/s from the east waste rock dump was predicted to flow
within the overburden and shallow bedrock under the engineered drainage ditches and discharge
to Creek 661.
o Seepage originating from the west waste rock dump was predicted to discharge primarily to
drainage channels routed to the TSF. Seepage from this facility was not predicted to reach the
downstream environment.
o Seepage originating from the plant site was predicted to flow to natural and engineered drainages
that are routed to the TSF. Seepage from this facility was not predicted to reach the downstream
environment.
Results of the seepage analysis indicated that seepage flow paths originating from TSF Site D, the Pit
Lake and the east waste rock dump converge beneath the TSF spillway and Creek 661. These local
groundwater flow paths would discharge to the overlying drainages. Seepage flow paths originating
from the Pit Lake were modelled to travel via deeper (regional) groundwater flow paths within the
competent bedrock to Davidson Creek.
7.2 RECOMMENDATIONS
The following monitoring is recommended to be collected to verify model results provided in this report,
address uncertainty in predicted effects based on these results and provide guidance for monitoring and
contingency planning:
Monitoring of groundwater levels in the deposit area to evaluate the effects of pit dewatering.
Monitoring of groundwater levels and collection of groundwater quality samples at locations
downgradient of facilities to assess potential seepage and flow reductions. Model results indicate that
shallow seepage pathways within the overburden and shallow bedrock are topographically controlled.
The monitoring network should be designed so groundwater monitoring locations are installed along
these pathways, particularly downgradient of the proposed location of the TSF south abutment. The
monitoring network should also be designed to assess the efficiency of the seepage interception
works.
Streamflow measurements and water quality sampling downgradient of the facilities at Creek 661 and
Davidson Creek.
7.0 REFERENCES
Harbaugh, A.W., Banta, E.R., Hill, M.C., and McDonald, M.G., 2000. MODFLOW-2000, the U.S.
Geological Survey modular ground-water model -- User guide to modularization concepts and the
Ground-Water Flow Process: U.S. Geological Survey Open-File Report 00-92, 121 pp.
Deere, D. U. and Patton, F. D., 1971: Slope stability in residual soils, Fourth Panam. Conf. SMFE, San
Juan, Puerto Rico, 1, pages 87–170.
Freeze, R.A. and Cherry, J.A., 1979. Groundwater, Prentice-Hall.
Heath, R.C., 1989. Basic Ground-Water Hydrology, U.S. Geological Survey Water-Supply Paper 2220.
Hydrogeologic Inc. (HGL), 1996. MODHMS/MODFLOW-SURFACT: A Comprehensive MODFLOW-
Based Hydrologic Modeling System. Reston, VA.
Knight Piésold Ltd. 2013a., Reconnaissance Terrain and Terrain Stability Mapping.
KP Ref No. VA101-457/4-4, Rev 0. February 2012. .
Knight Piésold Ltd., 2013b. 2012 Site Investigation Report. KP Ref No. VA101-457/6-1, Rev 0.
September 2013.
Knight Piésold Ltd., 2013c. Feasibility Open Pit Slope Design. KP Ref No. VA101-457/6-2, Rev 0,
November 2013.
Knight Piésold Ltd., 2013d. 2013 Hydrometerology Report. KP Ref No. VA101-457/6-12, Rev 0,
November 2013.
Knight Piésold Ltd., 2013e. 2013 Site Investigation Report. KP Ref No. VA101-457/6-4, Rev 0,
November 2013.
Knight Piésold Ltd., 2013f. Watershed Modelling Report. KP Ref No. VA101-457/6-6, Rev 0,
November 2013.
Knight Piésold Ltd., 2013g. Geotechnical Characterization Report. KP Ref No. VA101-457/6-8, Rev 0,
November 2013.
Knight Piésold Ltd., 2013h. Open Pit Water Management Plan. KP Ref No. VA101-457/6-8, Rev 0,
November 2013.
Knight Piésold Ltd., 2013i. Updated FS Water Balance Model Memo. KP Ref No. VA13-02400, Rev 0,
December 2013.
Plouffe, A. and Levson V.M., 2001. Surficial geology, Tatelkuz Lake, British Columbia; Geological Survey
of Canada, Open File 4001, scale 1:100,000.
Plouffe, A. and Levson V.M., 2002. Surficial geology, Entiako Lake, British Columbia; Geological Survey
of Canada, Open File 4157, scale 1:100,000.
Plouffe, A., Levson, V.M., and Mate, D.J., 2004. Surficial Geology, Nechako River, British Columbia.
Geological Survey of Canada, Map 2067A, scale 1:250,000.
Pollock, D.W., 1994. User's Guide for MODPATH/MODPATH-PLOT, Version 3: A particle tracking post-
processing package for MODFLOW, the U.S. Geological Survey finite-difference ground-water
flow model: U.S. Geological Survey Open-File Report 94-464, 6 ch.
Western Water Associates, Ltd., 2013. Memorandum regarding New Gold, Blackwater Construction
Camp, Water Supply Feasibility. August 2013.
APPENDIX A
APPENDIX FIGURES
350,000
360,000
370,000
380,000
390,000
400,000
C
he
da
ku
zC
r
ee
k
1-TC
5,910,000 5,910,000
1-DC
BWL
H6 4-DC
Ta
te
eek
H4B
lku
Cr
z
ek
La
rtle H3 re
Tu
ke
C
on
ds
a
vi
D
e
ak 700
pL Creek 1-661
5,900,000
To 5,900,000
H2
k 661
H1 Cree
1-505659
k
11-DC
ee
4-705
Cr
BWH
nie
6-705
705
aw
ek
1-705 Cre
F
H7
SAVED: M:\1\01\00457\06\A\GIS\Figs\Report13_NumericalGroundwaterModellingeport\StreamCoursesProjectHydrologyStns.mxd; Jan 17, 2014 9:14 AM; cczembor
MOUNT
DAVIDSON
5,890,000 5,890,000
a ter Riv er
kw
B la c
5,880,000 5,880,000
340,000
350,000
360,000
370,000
380,000
390,000
400,000
1.75 0.875 0 1.75 3.5 5.25 7 km
SCALE
372,000
374,000
376,000
378,000
380,000
LEGEND:
òN
>
!
KAME
E
" RIVER/CREEK
CONTOUR (5m)
FLUTING
LODGEMENT TILL
ABLATION TILL
GLACIOLACUSTRINE
COLLUVIUM DEPOSITS
BEDROCK
LAKE
5,896,000 5,896,000
1300
SAVED: M:\1\01\00457\06\A\GIS\Figs\Report13_NumericalGroundwaterModellingeport\SurficialGeologyAndLandformsMap.mxd; Dec 20, 2013 11:15 AM; adinca
5,894,000 5,894,000
NOTES:
1. COORDINATE GRID IS IN METRES.
COORDINATE SYSTEM: NAD 1983 UTM ZONE 10N.
1400
2. THIS FIGURE IS PRODUCED AT A NOMINAL SCALE OF 1:35,000
FOR 11x17 (TABLOID) PAPER. ACTUAL SCALE MAY DIFFER
ACCORDING TO CHANGES IN PRINTER SETTINGS OR
PRINTED PAPER SIZE.
5,892,000 5,892,000
00
17 350 175 0 350 700 1,050 1,400 1,750 m
SCALE
372,000
374,000
376,000
378,000
380,000
1800
VA101-457/6 13
FIGURE A.2
0 20NOV'13 ISSUED WITH REPORT JEH AMD JAS KJB REV
REV DATE DESCRIPTION DESIGNED DRAWN CHK'D APP'D 0
NEW GOLD INC.
BLACKWATER GOLD PROJECT
APPENDIX B
0.000007
0.000006
5000
0.000005
Elevation (m)
Model Y-Coordinate (m)
0.000004
10000
0.000003
15000
0.000001
0
0 0.5 1 1.5 2 2.5 3 3.5
Area (ha)
20000
Horizontal Hydraulic Vertical Hydraulic
Zone Color Material Type
Conductivity Kh (m/sec) Conductivity Kv (m/sec)
Zone 6 Glaciof luvial: Channel 1.E-04 1.E-04
Zone 2 Glaciof luvial: Kame 1.E-05 1.E-05 NEW GOLD INC.
Zone 1 Glacial Till 2.E-07 2.E-07
BLACKWATER GOLD PROJECT
Zone 4 Glaciolacustrine 1.E-07 2.E-08
Zone 12 Kame/Till Deposit 9.E-06 2.E-07
BASELINE MODEL
Zone 11 Bedrock Outcrops 1.E-07 1.E-07
HYDRAULIC CONDUCTIVITY ZONES (LAYER 1)
NOTES:
1. THE CONTOUR LINES SHOWN ABOVE ILLUSTRATE GROUND SURFACE ELEVATION (masl) . P/A NO. REF. NO.
VA101-457/6-13 13
REV
0 16JAN'14 ISSUED WITH REPORT KTD CAS KJB FIGURE B.1 0
REV DATE DESCRIPTION PREP'D CHK'D APP'D
Model X-Coordinate (m)
0 5000 10000 15000 20000 25000 30000
0
5000
Model Y-Coordinate (m)
10000
15000
20000
Horizontal Hydraulic Vertical Hydraulic
Zone Color Material Type
Conductivity Kh (m/sec) Conductivity Kv (m /sec)
NEW GOLD INC.
Zone 1 Glacial Till 2.E-07 2.E-07
Zone 4 Glaciolacustrine 1.E-07 2.E-08 BLACKWATER GOLD PROJECT
Zone 5 Competent Bedrock 1.E-08 1.E-08
BASELINE MODEL
HYDRAULIC CONDUCTIVITY ZONES (LAYER 2)
NOTES:
1. THE CONTOUR LINES SHOWN ABOVE ILLUSTRATE THE TOP ELEVATION P/A NO. REF. NO.
OF LAYER 2 (masl). VA101-457/6-13 13
REV
0 16JAN'14 ISSUED WITH REPORT KTD CAS KJB FIGURE B.2 0
REV DATE DESCRIPTION PREP'D CHK'D APP'D
Model X-Coordinate (m)
0 5000 10000 15000 20000 25000 30000
0
5000
Model Y-Coordinate (m)
10000
15000
20000
REV
0 16JAN'14 ISSUED WITH REPORT KTD CAS KJB FIGURE B.3 0
REV DATE DESCRIPTION PREP'D CHK'D APP'D
Model X-Coordinate (m)
0 5000 10000 15000 20000 25000 30000
0
5000
Model Y-Coordinate (m)
10000
15000
20000
REV
0 16JAN'14 ISSUED WITH REPORT KTD CAS KJB FIGURE B.4 0
REV DATE DESCRIPTION PREP'D CHK'D APP'D
Model X-Coordinate (m)
0 5000 10000 15000 20000 25000 30000
0
5000
Model Y-Coordinate (m)
10000
15000
20000
Horizontal Hydraulic Vertical Hydraulic
Zone Color Material Type
Conductivity Kh (m/sec) Conductivity Kv (m /sec)
Zone 10 Competent Bedrock 2.E-08 2.E-08 NEW GOLD INC.
Zone 7 Higher Permeability Bedrock Zone in Deposit 1.E-06 1.E-06 BLACKWATER GOLD PROJECT
REV
0 16JAN'14 ISSUED WITH REPORT KTD CAS KJB FIGURE B.5 0
REV DATE DESCRIPTION PREP'D CHK'D APP'D
Model X-Coordinate (m)
0 5000 10000 15000 20000 25000 30000
0
5000
Model Y-Coordinate (m)
10000
15000
20000
Horizontal Hydraulic Vertical Hydraulic
Zone Color Material Type
Conductivity Kh (m /s ec) Conductivity Kv (m/sec)
Zone 5 Competent Bedrock 1.E-08 1.E-08 NEW GOLD INC.
Zone 7 Higher Permeability Bedrock Zone in Deposit 1.E-06 1.E-06 BLACKWATER GOLD PROJECT
REV
0 16JAN'14 ISSUED WITH REPORT KTD CAS KJB FIGURE B.6 0
REV DATE DESCRIPTION PREP'D CHK'D APP'D
Model X-Coordinate (m)
0 5000 10000 15000 20000 25000 30000
0
5000
Model Y-Coordinate (m)
10000
15000
20000
Horizontal Hydraulic Ve rtical Hydraulic
Zone Color Material Type
Conductivity Kh (m /s ec) Conductivity Kv (m/s ec)
Zone 8 Competent Bedrock 1.E-09 1.E-09 NEW GOLD INC.
BLACKWATER GOLD PROJECT
REV
0 16JAN'14 ISSUED WITH REPORT KTD CAS KJB FIGURE B.7 0
REV DATE DESCRIPTION PREP'D CHK'D APP'D
NEW GOLD INC.
BLACKWATER GOLD PROJECT
APPENDIX C
GENERAL ARRANGEMENTS
APPENDIX D
6000
5000
4000
Elevation (m)
3000
2000
1000
-1000
Area (ha)
NOTES:
1. PARTICLE TRACES ARE COLOURED ACCORDING TO THE MODEL LAYER IT IS TRAVELLING IN.
2. THE CONTOURS SHOWN ABOVE ARE OF WATER TABLE ELEVATION (masl). NEW GOLD INC.
BLACKWATER GOLD PROJECT
TSF D MODPATH PARTICLE ANALYSIS
PARTICLES BY MODEL LAYER
P/A NO. REF. NO.
VA101-457/6 13
NOTES:
1. PARTICLE TRACES ARE COLOURED ACCORDING TO THE LOCATION WHERE THE MODPATH PARTICLE
DISCHARGES FROM THE MODEL DOMAIN:
PINK = CREEK 661
TEAL = DAVIDSON CREEK NEW GOLD INC.
BROWN = TSF MAIN EMBANKMENT DRAINS
RED = TSF SOUTH ABUTMENT DRAINS BLACKWATER GOLD PROJECT
ORANGE = TSF SPILLWAY
PURPLE = TSF SITE D POND TSF D MODPATH PARTICLE ANALYSIS
BLUE = ECD SYSTEM
BLACK = NATURAL CHANNEL REPORTING TO TSF
PARTICLES BY DISCHARGE LOCATIONS
2. THE CONTOURS SHOWN ABOVE ARE OF WATER TABLE ELEVATION (masl).
P/A NO. REF. NO.
VA101-457/6 13
0 16JAN'14 ISSUED WITH REPORT KTD CAS KJB
REV
REV DATE DESCRIPTION PREP'D CHK'D APP'D
FIGURE D.2 0
\\van11\prj_file\1\01\00457\06\A\Report\13 - Numerical Groundwater Modelling Report\Rev0\Appendices\Appendix D\[Appendix D.xls]D3 Print 1/17/2014 11:49 AM
Color MODPATH Discharge Locations
Creek 661
Davidson Creek
TSF Main Embankment Drains
TSF South Abutment Drains
ECD System
TSF Spillway
TSF Site D Pond
Natural Channel Upstream of TSF
NOTES:
1. THE EAST DUMP FOOTPRINT IS COLOURED ACCORDING TO THE LOCATION WHERE THE MODPATH
PARTICLE DISCHARGES. NEW GOLD INC.
2. THE CONTOURS SHOWN ABOVE ARE OF WATER TABLE ELEVATION (masl).
BLACKWATER GOLD PROJECT
TSF D MODPATH PARTICLE ANALYSIS
FOOTPRINT BY PARTICLE DISCHARGE LOCATION
P/A NO. REF. NO.
VA101-457/6 13
NOTES:
NEW GOLD INC.
1. PARTICLE TRACES ARE COLOURED ACCORDING TO THE MODEL LAYER IT IS TRAVELLING IN.
3. PARTICLES WERE INSERTED INTO LAYERS 1 THROUGH 8. BLACKWATER GOLD PROJECT
2. THE CONTOURS SHOWN ABOVE ARE OF WATER TABLE ELEVATION (masl).
PIT LAKE MODPATH PARTICLE ANALYSIS
PARTICLES BY MODEL LAYER
P/A NO. REF. NO.
VA101-457/6 13
NOTES:
1. PARTICLE TRACES ARE COLOURED ACCORDING TO THE LOCATION WHERE THE MODPATH PARTICLE NEW GOLD INC.
DISCHARGES FROM THE MODEL DOMAIN:
PURPLE = CHANNEL TO EAST DUMP DRAINAGE DITCHES BLACKWATER GOLD PROJECT
PINK = CREEK 661
TEAL = DAVIDSON CREEK PIT LAKE MODPATH PARTICLE ANALYSIS
BLUE = LOWER EAST DUMP DRAINAGE DITCH PARTICLES BY DISCHARGE LOCATIONS
BLACK = NATURAL CHANNEL UPSTREAM OF TSF
BROWN = TSF D EMBANKMENT DRAINS
P/A NO. REF. NO.
2. THE CONTOURS SHOWN ABOVE ARE OF WATER TABLE ELEVATION (masl). VA101-457/6 13
NOTES:
Layer 7 Layer 8 1. THE IMAGES OF MODEL 1 LAYERS THROUGH 8 ABOVE
PRESENT MODPATH PARTICLE TRACES BY THE LAYER IN WHICH
THEY ORIGINATE.
2. PARTICLE TRACES ARE COLOURED ACCORDING TO THE
MODEL LAYER IT IS TRAVELLING IN.
3. CONTOURS INDICATE WATER TABLE ELEVATION (masl)
NOTES:
1. PARTICLE TRACES ARE COLOURED ACCORDING TO THE MODEL LAYER IT IS TRAVELLING IN.
2. THE CONTOURS SHOWN ABOVE ARE OF WATER TABLE ELEVATION (masl). NEW GOLD INC.
BLACKWATER GOLD PROJECT
EAST DUMP MODPATH PARTICLE ANALYSIS
PARTICLES BY MODEL LAYER
P/A NO. REF. NO.
VA101-457/6 13
NOTES:
1. PARTICLE TRACES ARE COLOURED ACCORDING TO THE LOCATION WHERE THE MODPATH PARTICLE
DISCHARGES FROM THE MODEL DOMAIN: NEW GOLD INC.
DARK PURPLE = NATURAL CHANNEL REPORTING TO TSF
BLUE = LOWER EAST DUMP DRAINAGE DITCH
BLACKWATER GOLD PROJECT
GREEN = CHANNEL TO EAST DUMP DRAINAGE DITCHES
PURPLE = SPRINGS WITHIN EAST DUMP FOOTPRINT EAST DUMP MODPATH PARTICLE ANALYSIS
BROWN = TSF SPILLWAY PARTICLES BY DISCHARGE LOCATION
2. THE CONTOURS SHOWN ABOVE ARE OF WATER TABLE ELEVATION (masl).
P/A NO. REF. NO.
VA101-457/6 13
WEST DUMP MODPATH PARTICLE ANALYSIS
PARTICLES BY MODEL LAYER
P/A NO. REF. NO.
VA101-457/6 13
PLANT SITE MODPATH PARTICLE ANALYSIS
PARTICLES BY MODEL LAYER
P/A NO. REF. NO.
VA101-457/6 13
APPENDIX E
HISTOGRAMS OF PARTICLE TRAVEL TIMES
PIT LAKE
P/A NO. REF. NO.
VA101-457/6 13
HISTOGRAMS OF PARTICLE TRAVEL TIMES
EAST WASTE ROCK DUMP
P/A NO. REF. NO.
VA101-457/6 13
NOTES:
1. HISTOGRAM BIN SIZE IS 25 YEARS. NEW GOLD INC.
BLACKWATER GOLD PROJECT
HISTOGRAMS OF PARTICLE TRAVEL TIMES
WEST WASTE ROCK DUMP
P/A NO. REF. NO.
VA101-457/6 13
HISTOGRAMS OF PARTICLE TRAVEL TIMES
PLANT SITE
P/A NO. REF. NO.
VA101-457/6 13