Fluid Flow Model for Predicting the Intrusion Rate of Subsurface

Jun 25, 2018 - The technique includes (1) subslab gas extraction with flow and vacuum ... the source-building separation distance → 0) and advection...
2 downloads 0 Views 765KB Size
Subscriber access provided by University of Sussex Library

Environmental Modeling

Fluid flow model for predicting the intrusion rate of subsurface contaminant vapors into buildings Todd Arthur McAlary, John Gallinatti, Gordon Thrupp, William Wertz, Darius Mali, and Helen Dawson Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.8b01106 • Publication Date (Web): 25 Jun 2018 Downloaded from http://pubs.acs.org on June 28, 2018

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 28

Environmental Science & Technology

1

Fluid Flow Model for Predicting the Intrusion Rate

2

of Subsurface Contaminant Vapors into Buildings

3

Todd A. McAlary**1, John Gallinatti2, Gordon Thrupp2, William Wertz3, Darius Mali4 and Helen

4

Dawson5

5

1

Geosyntec Consultants, Inc., Toronto, ON

6

2

Geosyntec Consultants, Inc., Oakland, CA

7

3

Geosyntec Consultants, Inc., Albany, NY

8

4

Geosyntec Consultants, Inc., Guelph, ON

9

5

Geosyntec Consultants, Inc. Washington, D.C.

10

Abstract

11

A new method is presented for calculating a building-specific subslab to indoor air attenuation

12

factor for use in assessing subsurface vapor intrusion to indoor air. The technique includes: 1)

13

subslab gas extraction with flow and vacuum measurements and mathematical modeling to

14

characterize the bulk average vertical gas conductivity of the floor slab, 2) monitoring of the

15

ambient pressure gradient across the floor slab with a micromanometer, 3) calculating the *

Corresponding author – Geosyntec Consultants, Inc., 3250 Bloor Street West, Suite 600, Toronto, ON M8X 2X9; phone (905) 339-7066; e-mail: [email protected]

ACS Paragon Plus Environment

1

Environmental Science & Technology

Page 2 of 28

16

volumetric flow of soil gas into the building (Qsoil), and 4) dividing Qsoil by the building

17

ventilation rate (Qbuilding) to calculate a building-specific attenuation factor.

18

calculations using order statistics from 121 individual tests are comparable to the U.S. EPA

19

empirical attenuation factors for residential buildings and the U.S. Navy empirical attenuation

20

factors for commercial/industrial buildings. A case study of a commercial building shows

21

encouraging agreement between the attenuation factors calculated via this method and via

22

conventional subslab and indoor air sampling.

23

Introduction

24

Risk assessments typically assume that people breathe an average of 20,000 L of air per day and

25

drink about 2 L of water per day (1). As a result, the concentrations of chemicals in air must be

26

much lower than in water to yield the same exposure, so subsurface vapor intrusion to indoor air

27

is often the pathway posing the greatest potential risk at sites contaminated with VOCs(2).

28

Radon exposure is estimated to result in about 20,000 deaths per year in the United States

29

alone(3).

30

Regulatory guidance for assessing VOC vapor intrusion (2) and radon (3) has been developed to

31

protect human health in the United States, and some other countries have similar programs.

32

Nazaroff et al. conducted early research into radon intrusion, including a five-month study of a

33

detached house with a basement (4) which concluded that pressure-driven flow is an important

34

mechanism for radon entry.

35

temperature differences, wind loads, barometric pressure changes, and mechanical fans. Gas

36

entry typically occurs at discontinuities in the foundation that are irregular and difficult to

37

characterize or predict. For VOCs, the source of vapors is often at depth below the building and

Sample

Soil gas entry to buildings is therefore an important concern for human health.

Building pressure gradients can be caused by indoor-outdoor

ACS Paragon Plus Environment

2

Page 3 of 28

Environmental Science & Technology

38

vertical diffusion to the region just below the foundation is also an important mechanism, which

39

is often simulated using some form of the Jury model (5).

40

Several attempts have been made to mathematically model soil gas or vapor intrusion over the

41

past few decades. Radon models (6-9) and VOC models (10) have been developed for decades

42

by researchers around the world, which is a testament to the importance and difficulty of this

43

task (see also additional citations in the Supporting Information). The more detailed models are

44

limited by availability of data for important parameters and the simpler models are limited by

45

their ability to simulate all the processes involved. As a result, building-specific modeling of

46

indoor air concentrations attributable to VOC vapor or radon gas intrusion has been elusive.

47

The most commonly-used model for VOC vapor intrusion has been the Johnson and Ettinger

48

Model (11), a 1-dimensional, steady-state model that provides an algebraic expression for the

49

attenuation factor (AF), defined as the indoor air concentration (Cindoor) divided by the

50

subsurface vapor concentration at a specified depth (Csource):

51

AF=(Cindoor/Csource)

(1)

52

For the special condition of a source directly beneath the foundation (i.e., the source-building

53

separation distance → 0) and advection is the dominant mechanism of transport across the floor

54

slab, AF→Qsoil/Qbuilding, where Qsoil is the volumetric flow rate of soil gas into the building

55

and Qbuilding is the total volumetric flow rate of air through the building or basement (i.e., the

56

ventilation rate, which equals the interior volume multiplied by the air exchange rate). The

57

building volume is easily measured or estimated. The air exchange rate can usually be estimated

58

within a factor of about three for a typical residence: the mean air exchange rate for residences is

59

0.45 /hr and the 10th percentile is 0.18 /hr, according to the U.S. EPA Exposure Factors

ACS Paragon Plus Environment

3

Environmental Science & Technology

Page 4 of 28

60

Handbook (1). Chan et al. (12) analyzed a database of 70,000 residential building leakage tests

61

and showed comparable results. Qbuilding can be measured using blower door tests(13) and

62

tracer tests(14,15). The air exchange rates for commercial buildings tend to be higher (mean =

63

1.5 /hr(1)) and sometimes are available from mechanical engineers responsible for the heating,

64

ventilating and air conditioning (HVAC) systems. Qsoil is more difficult to estimate. Nazaroff

65

(6) developed an equation assuming that flow to the expansion joint between the floor and the

66

wall of a building is similar to a line-sink at the bottom of the slab:

67 68 69

Qsoil = 2 π k ∆P Xcrack µ ln(2Zcrack/Rcrack)

(2)

70

where:

71

k = the shallow soil or granular fill permeability to air flow [L2]

72

∆P = the indoor-outdoor pressure difference [M/L-T2]

73

µ = the viscosity of air [M/L-T]

74

Xcrack = the total length of cracks through which soil gas vapors are flowing [L]

75

Rcrack = the effective crack width [L]

76

Zcrack = the crack opening depth below grade, or thickness of the floor slab [L]

77

∆P can be measured with a pressure transducer/datalogger and Zcrack is seldom much different

78

than 0.1 to 0.15 m.

79

practically impossible to accurately measure Xcrack and Rcrack on a building-specific basis, so the

80

Nazaroff model is not very useful for building-specific calculations of Qsoil.

However, direct measurements of k are seldom available and it may be

Johnson (16)

ACS Paragon Plus Environment

4

Page 5 of 28

Environmental Science & Technology

81

argued that because of challenges estimating Qsoil, it would be preferable to combine

82

Qsoil/Qbuilding, which can be reasonably estimated (from mass balance principles) to be equal

83

to the ratio of the indoor air concentrations divided by subslab concentrations (i.e., the subslab-

84

to-indoor air AF).

85

As an alternative to mathematical modeling, the U.S. EPA (17) compiled and analyzed a

86

database of indoor air and subslab sampling results, primarily for residential structures and

87

chlorinated compounds to better understand the degree to which subslab VOC concentrations are

88

attenuated as soil gas entering the building is diluted by the ventilation of the building (i.e., the

89

subslab-to-indoor air AF).

90

concentrations high enough to minimize bias from background sources, they calculated a 95th

91

percentile value of 0.03 for the subslab-to-indoor air AF.

92

adopted 0.03 as a default AF for development of subslab screening levels (SSSL = Cindoor/AF).

93

Considering this is a 95th percentile AF, it would be expected to be overly protective for 19 of 20

94

residential buildings and Song et al. (18) and Yao et al. (19) provide arguments that the EPA

95

database analysis may result in a high bias in the AF values. Furthermore, higher ceilings and

96

higher air exchange rates for commercial buildings tend to result in lower AFs than residential

97

buildings, but agencies are often reluctant to adopt lower default AFs for commercial/industrial

98

buildings despite available empirical subslab and indoor air data that indicates this would be

99

appropriate (20). Regulators usually recommend additional sampling if subslab concentrations

100

exceed the conservative screening levels and some recommend mitigation if subslab vapor

101

concentrations exceed screening levels by a relatively modest factor (e.g., 10x) (e.g., 21,22).

102

Considering the U.S. EPA empirical database shows almost 5 orders of magnitude range in the

After filtering the database to focus on buildings with subslab

Many regulatory agencies have

ACS Paragon Plus Environment

5

Environmental Science & Technology

Page 6 of 28

103

subslab AFs, using a default of 0.03 for all buildings is overly protective by a wide margin for

104

most buildings.

105

There is a distinct need for better methods to determine the susceptibility of specific buildings to

106

VOC vapor and radon gas intrusion.

107

temporal variability (23,24) and bias from background sources (25), and subslab samples are

108

subject to uncertainties from spatial variability, so empirical AF values derived from the EPA

109

database are unavoidably imperfect. The purpose of this article is to develop and demonstrate a

110

method to calculate building-specific attenuation factors using subslab flow and vacuum

111

measurements, building pressure and building ventilation rates, which will provide an

112

independent line of evidence to augment empirical attenuation factors.

113

Conceptual Model

114

In the field of hydrogeology, aquifer discharge (Qaq) is calculated using (26):

115

Indoor air samples are subject to uncertainties from

Qaq = K i A

(3)

116

where:

117

K = hydraulic conductivity [L/T]

118

i = hydraulic gradient [L/L], and

119

A = cross-sectional area perpendicular to flow (L2)

120

Equation 3 can be applied to calculating Qsoil where A is the footprint of the building, i is the

121

pressure gradient across the floor slab and K is the bulk average vertical gas conductivity of the

122

floor slab. Groundwater models can be applied to soil gas flow as long as the pressure gradients

ACS Paragon Plus Environment

6

Page 7 of 28

Environmental Science & Technology

123

are small enough that the gas behaves similar to an incompressible fluid (27). Subslab venting

124

systems typically operate with applied vacuum