Stress-dependant characteristics of deep marine sediments recovered from the Ulleung Basin, East Sea, Korea

Similar documents
Electrical and geomechanical Properties of Natural Gas Hydratebearing Sediments from Ulleung Basin, East Sea, Korea

Effect of Gas Hydrate Saturation on Hydraulic Conductivity of Marine Sediments

Effect of porosity and permeability reduction on hydrate production in marine sediments

Effect of Frozen-thawed Procedures on Shear Strength and Shear Wave Velocity of Sands

Experimental study on seismic wave velocity of hydrate-bearing fine-grained sediments

Engineering Geology 117 (2011) Contents lists available at ScienceDirect. Engineering Geology

The Seafloor deformation and well bore stability monitoring during gas production in unconsolidated reservoirs

CE 240 Soil Mechanics & Foundations Lecture 3.2. Engineering Classification of Soil (AASHTO and USCS) (Das, Ch. 4)


Geotechnical Properties of Soil

Chapter 1 - Soil Mechanics Review Part A

Solution:Example 1. Example 2. Solution: Example 2. clay. Textural Soil Classification System (USDA) CE353 Soil Mechanics Dr.

Table of Contents Chapter 1 Introduction to Geotechnical Engineering 1.1 Geotechnical Engineering 1.2 The Unique Nature of Soil and Rock Materials

Downloaded 11/20/12 to Redistribution subject to SEG license or copyright; see Terms of Use at

Rock physics and AVO applications in gas hydrate exploration

A. V T = 1 B. Ms = 1 C. Vs = 1 D. Vv = 1

8. ELASTIC PROPERTY CORRECTIONS APPLIED TO LEG 154 SEDIMENT, CEARA RISE 1

Geology and Soil Mechanics /1A ( ) Mark the best answer on the multiple choice answer sheet.

APPENDIX A. Borehole Logs Explanation of Terms and Symbols

Dashed line indicates the approximate upper limit boundary for natural soils. C L o r O L C H o r O H


Tutorial on Methane Hydrate. Presented by Ad Hoc Group on Methane Hydrate Research March 24, 2004

Calculation of 1-D Consolidation Settlement

Geology 229 Engineering Geology. Lecture 8. Elementary Soil Mechanics (West, Ch. 7)

SWELLING POTENTIAL OF THE RESIDUAL SOIL, BASED ON PLASTICITY INDEX VALUE AT MUAROSIJUNJUNG, WEST SUMATERA

GEOTECHNICAL INVESTIGATION REPORT

The Casagrande plasticity chart does it help or hinder the NZGS soil classification process?

4. Soil Consistency (Plasticity) (Das, chapter 4)

Chapter I Basic Characteristics of Soils

Integrated Ground Behavior

Module 1 GEOTECHNICAL PROPERTIES OF SOIL AND OF REINFORCED SOIL (Lectures 1 to 4)

Gas Hydrate R&D from 1996 to Regional Survey. activities at KIGAM. Byong-Jae Ryu


Effect of Lime on the Compressibility Characteristics of a Highly Plastic Clay

Study of Shear Wave Velocity of Macao Marine Clay under Anisotropic Stress Condition

Clay Robinson, PhD, CPSS, PG copyright 2009

H.1 SUMMARY OF SUBSURFACE STRATIGRAPHY AND MATERIAL PROPERTIES (DATA PACKAGE)

Soil Properties - I. Amit Prashant. Indian Institute of Technology Gandhinagar. Short Course on. Geotechnical Aspects of Earthquake Engineering

On the Role of ph in the Cyclic Behavior of Fine-Grained Soils

INFLUENCE OF MICA CONTENT ON DYNAMIC SHEAR MODULUS OF SANDY SOILS

DEPARTMENT OF GEOSCIENCES

SOIL TYPE INVESTIGATION OF THE DOĞANHİSAR CLAYS, CENTRAL ANATOLIA, TURKEY. İ.İnce* and A. Özdemir

A Study on Dynamic Properties of Cement-Stabilized Soils

GEOL 314 Engineering Geology

Introduction to Soil Mechanics Geotechnical Engineering-II

Exploration, Drilling & Production

The Ohio Department of Transportation Office of Research & Development Executive Summary Report

LABORATORY MEASUREMENTS OF STIFFNESS OF SOFT CLAY USING BENDER ELEMENTS

DEFINITION OF EXPANSION POTENTIAL FOR EXPANSIVE SOIL John D. Nelson 1,2, Kuo-Chieh Chao 2 & Daniel D. Overton 2

Cyclic Triaxial Behavior of an Unsaturated Silty Soil Subjected to Suction Changes

Class Principles of Foundation Engineering CEE430/530

Supplementary Problems for Chapter 7

THE GEOTECHNICAL INDEX PROPERTIES OF SOIL IN WARRI, DELTA STATE, NIGERIA

Unloading Test with Remolded Marine Soil Sample and Disturbance Degree Assessment

INTERNATIONAL JOURNAL OF CIVIL AND STRUCTURAL ENGINEERING Volume 1, No 4, 2011

Soil Mechanics Brief Review. Presented by: Gary L. Seider, P.E.

QUANTITATIVE ASSESSMENT OF GAS HYDRATES IN THE MALLIK L-38 WELL, MACKENZIE DELTA, N.W.T., CANADA

SOME GEOTECHNICAL PROPERTIES OF KLANG CLAY

Consolidation Properties of NAPL Contaminated Sediments

Correlation of unified and AASHTO soil classification systems for soils classification

VARIATION OF CONSOLIDATION COEFFICIENT OF EXPANSIVE CLAYS AT HIGH INITIAL WATER CONTENT

G 3. AN ELECTRONIC JOURNAL OF THE EARTH SCIENCES Published by AGU and the Geochemical Society

Soil Behaviour Type from the CPT: an update

Seismic Velocity Dispersion and the Petrophysical Properties of Porous Media

EFFECT OF SOIL TYPE LOCATION ON THE LATERALLY LOADED SINGLE PILE

The Consolidation and Strength Behavior of Mechanically Compressed Fine-Grained Sediments

Oedometer and direct shear tests to the study of sands with various viscosity pore fluids

UNIVERSITY OF CALGARY. The Role of THF Hydrate Veins on the Geomechanical Behaviour of Hydrate-Bearing Fine. Grained Soils. Jiechun Wu A THESIS

Methane hydrate rock physics models for the Blake Outer Ridge

Boreholes. Implementation. Boring. Boreholes may be excavated by one of these methods: 1. Auger Boring 2. Wash Boring 3.

Evaluation of liquefaction resistance of non-plastic silt from mini-cone calibration chamber tests

GEOTECHNICAL SUBSOIL PROPERTIES FOR WASIT GOVERNORATE

The Effects of Different Surcharge Pressures on 3-D Consolidation of Soil

TIME DOMAIN REFLECTOMETRY (TDR) IN MEASURING WATER CONTENTS AND HYDRATE SATURATIONS IN MARINE SEDIMENTS

Advances in Environmental Biology

19 th ICSMGE - Seoul 2017 EGEL. Geo-Education. J. Carlos Santamarina KAUST

Prediction of torsion shear tests based on results from triaxial compression tests

ADVANCES in NATURAL and APPLIED SCIENCES

B-1 BORE LOCATION PLAN. EXHIBIT Drawn By: 115G BROOKS VETERINARY CLINIC CITY BASE LANDING AND GOLIAD ROAD SAN ANTONIO, TEXAS.

Module 3. DYNAMIC SOIL PROPERTIES (Lectures 10 to 16)


Site Characterization & Hydrogeophysics

Appendix J. Geological Investigation

December 5, Junction Gateway, LLC 7551 W. Sunset Boulevard #203 Los Angeles, CA Mr. James Frost P: Dear Mr.

PRINCIPLES OF GEOTECHNICAL ENGINEERING

UNDRAINED FLOW CHARACTERISTICS OF PARTIALLY SATURATED SANDY SOILS IN TRIAXIAL TESTS

Course Scheme -UCE501: SOIL MECHANICS L T P Cr

The relationship between silica diagenesis and physical properties in the East/Japan Sea: ODP Legs 127/128

A study on nonlinear dynamic properties of soils

Characterization of Sand Formation from a Crude Oilfield

Electrical Resistivity of Compacted Kaolin and its Relation with Suction

DATA REPORT GEOTECHNICAL INVESTIGATION GALVESTON CRUISE TERMINAL 2 GALVESTON, TEXAS

THE IMPACT OF HYDRATE SATURATION ON THE MECHANICAL, ELECTRICAL, AND THERMAL PROPERTIES OF HYDRATE-BEARING SAND, SILTS, AND CLAY

Research Article Study on p-wave Attenuation in Hydrate-Bearing Sediments Based on BISQ Model

Mallik 2002 Gas Hydrate Production Research Well Program

POST CYCLIC SHEAR STRENGTH OF FINE GRAINED SOILS IN ADAPAZARI TURKEY DURING 1999 KOCAELI EARTHQUAKE

Engineer. Engineering. Engineering. (in-ja-neer ) A person trained and skilled in any of the various branches of engineering: a civil engineer

SASKATCHEWAN STRATIGRAPHY GLACIAL EXAMPLE BOULDERS IN GLACIAL DEPOSITS

Cone Penetration Testing in Geotechnical Practice

ARCH 1250 APPLIED ENVIRONMENTAL STUDIES

Transcription:

The 2012 World Congress on Advances in Civil, Environmental, and Materials Research (ACEM 12) Seoul, Korea, August 26-30, 2012 Stress-dependant characteristics of deep marine sediments recovered from the Ulleung Basin, East Sea, Korea Young Moon Kim 1), Jong-Sub Lee 2), Joo Yong Lee 3), and *Changho Lee 4) 1), 2), 4) School of Civil, Environmental and Architectural Engineering, Korea University, Seoul 136-713, Korea 3) Marine and Petroleum Resources Division, Korea Institute of Geoscience and Mineral Resources, Daejeon 305-350, Korea 4) shinesky@korea.ac.kr ABSTRACT Although the necessity of exploration in deep sea increases to develop the natural resources, not many geotechnical studies have been performed to explore the geotechnical characterizations of the sediment in Korean deep sea. Comprehensive laboratory tests are conducted to investigate index properties and stress-dependant characteristics of the deep marine sediments recovered from the hydrate occurrence regions during the Ulleung Basin Gas Hydrate Expedition 2 (UBGH2) in the Ulleung Basin, East Sea, Korea. Stress-dependant characteristics are measured by using an instrumented oedometer cell incorporated with shear wave measurements. The index properties of specimens are measured and compared with previous results: the measured water contents range within 25-35%, which are close to a low boundary of previously reported values; the sediments tested are classified as either CL/OH or CH/OH with high silt fraction (~80%). The stress-dependant mechanical properties represent a bi-linear behavior with the vertical effective stress. The hydrate presence significantly affects on the mechanical properties. The experimental results provide important geotechnical properties of deep marine sediments retrieved from the Ulleung Basin and design parameters for construction of the production facility. 1) 2) 3) 4) Graduate Student Associate Professor Research Engineer Research Professor

1. INTRODUCTION A drilling exploration in the deep-sea area is widely processed to develop a new energy resource advance in the world. Many geological studies including the stratigraphy, sedimentary environment, chemical composition, and micropaleontological studies have been carried out with the exploration of the oil resources from the late 1970s in the Ulleung Basin (Chough et al. 1985; Lee et al 1993; Hiller et al. 1996; Bahk et al. 2000; Kim et al. 2003). Since the presence of the gas hydrate had been confirmed during the Ulleung Basin Gas Hydrate Expedition 1 (UBGH1) in 2007, geotechnical properties of the deep marine sediment are important factors to access the safety of gas production facility and productivity from the hydrate-bearing sediment. However, with the exception of some studies (Kwon et al. 2011; Lee et al. 2011), few studies were performed to explore the geotechnical characterization of the deep marine in a domestic sea. In this study, comprehensive laboratory tests are conducted to investigate the stress-dependant characteristics of deep oceanic sediments recovered from the hydrate occurrence regions during the Ulleung Basin Gas Hydrate Expedition 2 (UBGH2) in the Ulleung Basin, East Sea. The index properties of the specimens including the specific gravity, water content, and Atterberg limits are measured and compared with previous results. Stress-dependant characteristics are measured by using an instrumented oedometer cell incorporated with the shear wave velocity measurement. 2. SAMPLE DESCRIPTION AND INDEX PROPERTIES Natural deep marine sediments were recovered using the drilling vessel D/V Fugro Synergy at 10 sites in the Ulleung Basin. 8 core samples from two boreholes (UBGH2-6B and UBGH2-6C) were selected and used to characterize index properties and stress-dependant properties in this study. Sample information including the hole-coresection numbers and corresponding sample depth are summarized in Table 1. Index properties were obtained for all tested samples. Consolidation test were performed for selected samples to determine consolidation characteristics and stress-dependant properties. It is well known that index properties can give a baseline to understand the behavior of the specimen. The mechanical properties can be estimated by using the well-known empirical correlation with index properties. Furthermore, index properties can be used to predict the hydrate occurrence in the sediments because the grain size, its chemical composition, and pore fluid characteristics affects on the morphology, extent, and growth of natural gas hydrate (Winters et al. 1999). 2.1. Water contents and specific gravity Soil media are composed of interacting discrete solid elements. The voids left by the particle are filled by a second phase that often consists of one or more fluids (e.g. air, water, electrolyte, and/or organic compound). The content of second phase in whole media is important because the macroscale behavior of particular material is controlled

by the interactions between particles, affected by the quantity of the fluid phase (Santamarina et al. 2001). Water content ω can be defined in terms of weight or volume. The gravimetric water content is determined by oven drying method (ASTM D2216). Dry temperature is reduced to 80, considering organic materials in the specimens. The measured water contents are in a relatively narrow range of 25-35%, except 6B-8H-3a. These values are close to a low boundary of previously reported values: ω=24.98-161.02% for sediments in the Ulleung Basin (Kwon et al. 2011). The specific gravity of a solid particle G s, which is defined as the unit weight of the solid material normalized with respect to the unit weight of water, depends on the composition material of the sediment. The specific gravity is measured using a pycnometer (ASTM D854). The measured values range G s =2.58-2.65. Previously reported values for the Ulleung Basin are in similar range (Lee et al. 1993; Kwon et al. 2011; Lee et al. 2011). The porosity n can be evaluated from the measured natural water content and specific gravity when the degree of saturation S is assumed as 100%: S n/(1-n)= ω G s. The porosity values are summarized in Table 1. The porosity determined by the laboratory test is in good agreement with the neutron porosity logging data in the field (data not shown). 2.2. Atterberg limits and soil classification Atterberg limits are commonly used in the geotechnical engineering for the classification, identification, and description of fine-grained soils. Selected specimens were tested to measure Atterberg limits (ASTM D4318). The measured values of liquid limit LL=46-68% and plastic limit PL=19-25% are plotted on the plasticity chart with previously reported values for this region in Fig. 1. Measured values for the specimens fall above A-line and the sediments are classified as either CL/OH or CH/OH, which corresponds to a diatomaceous clay reported by Lee et al. (2011). Hole-Core- Section Water depth [m] Depth [mbsf] ω [%] Atterberg limits LL [%] PL [%] PI [%] USCS G s Porosity σ pf [S/m] 6B-8H-3a 2157 54 45.04 46.27 25.48 20.79 CL or OL 2.63 0.54 41.3 6B-9H-1a 2157 62 35.86 52.28 23.41 22.87 CH or OH 2.62 0.48 41.9 6B-14H-3a 2157 102 25.79 48.87 17.29 31.58 CL or OL 2.58 0.40 41.9 6B-16H-3a2 2157 112 30.45 48.31 19.12 29.19 CL or OL 2.61 0.44 42.4 6B-26H-1b 2157 173 35.87 67.37 21.42 45.95 CH or OH 2.61 0.48 42.4 6B-27H-3b 2156 181 36.24 68.56 24.20 44.36 CH or OH 2.60 0.49 42.4 6C-9H-3c 2157 162 32.04 12.13 NP - ML 2.65 0.46 41.9 6C-10H-2c2 2157 169 33.18 18.47 NP - ML 2.64 0.47 42.2 Table. 1. Index properties of the Ulleung Basin Note: ω is the water content. LL, PL, and PI denotes the liquid limit, plastic limit, and plastic index, respectively. USGS is the unified soil classification system. G S is the specific gravity, σ pf is the electrical conductivity of the pore fluid. NP means non-plastic. The samples in the gray rows were recovered from the hydrate occurrence depth.

100 80 6B-8H 6B-14H 6B-26H 6B-9H 6B-16H 6B-27H U-line PI=0.9(LL-8) A-line PI=0.73(LL-20) Plastic index 60 40 20 0 0 30 60 90 120 150 Liquid limit [%] UBGH1-4 (Kwon et al. 2011) UBGH1-9 (Kwon et al. 2011; Lee et al. 2011) UBGH1-10 (Kwon et al. 2011; Lee et al. 2011) Fig. 1. Plasticity chart. 3. STRESS-DEPENDANT MECHANICAL PROPERTIES 3.1. Experiment equipment and process The stress-dependant properties are measured by using an instrumented oedometer cell, which is housed shear wave transducers. To minimize additional disturbance, a core liner was used as the oedometer cell in this study. After the specimen was placed into the oedometer cell, the specimen was saturated for 24 hrs with an electrolyte. The electrical conductivity of electrolyte is identical with that of the pore fluid measured in the field. Consolidation test (ASTM D 2435) was conducted to measure the stress-volume response. After settlement was measured, shear wave was immediately measured. 3.2. Compressibility Fig. 2 shows the stress-volume responses during loading and unloading for four selected specimens (6B-14H-3a, 6B-16H-3a2, 6C-9H-3c, and 6B-26H-1b). A bi-linear behavior is observed with the vertical effective stress, as shown in Fig. 2. Relatively high initial void ratio for 6B-9H and 26H result partly from the hydrate existence. As shown in Fig. 2, the hydrate existence significantly affects on the initial void ratio, compression index C C, and swelling index C S of the specimens. High specific surface from internal pore structures also cause high compressibility in hydrate-bearing sediments (Lee et al. 2010). The internal pore structures of specimens from Ulleung basin are primarily due to the abundance of diatoms. The values are even higher than those from Gulf of Mexico (Lee et al 2008), implying our specimens have significantly high portion of diatoms.

3.3. Shear wave velocity For selected specimen, shear wave velocities calculated from the measured shear wave signal are plotted with the vertical effective stress in Fig. 3. Bender elements are used to transmit and receive the shear waves in this study. Shear wave velocities measured at each loading stage show a power function increase with the vertical effective stress: V s =48.5(σ' V ) 0.218 for 6B-16H and V s =39.4(σ' V ) 0.275 for 6C-9H. Relative high value of the exponent β=0.275 for 6C-9H reflects changes in the interparticle contact during loading, probably due to the disturbance of the specimen by the hydrate dissociation. Similar results are reported in specimens from hydrate-bearing regions in Gulf of Mexico (Lee et al, 2008). Higher specific surface also cause the higher β exponent in hydrate-bearing sediments (Lee et al. 2010). (a) (b) 2.5 6B-14H 5.4 6B-26H Void ratio 2.0 1.5 1.0 C S =0.239 C S =0.381 6B-16H C C =1.233 C C =1.241 Void ratio 5.3 5.2 5.1 5.0 C S =0.029 C S =0.026 6C-9H C C =0.169 C C =0.137 0.5 10 100 1000 10000 Vertical effective stress [kpa] 4.9 10 100 1000 10000 Vertical effective stress [kpa] Fig. 2. Stress-volume response during loading (solid symbols) and unloading (empty symbols): (a) 6B-14H and 16H; (b) 6B-9H and 26H.

400 6B-16H Shear wave velocity [m/s] 300 200 100 6C-9H V S =39.4(σ' V ) 0.275 V S =48.5(σ' V ) 0.218 0 10 100 1000 10000 Vertical effective stress [kpa] Fig. 3. Evolution of shear wave velocities during loading. 4. CONCLUSION In this study, a series of laboratory are conducted to investigate index properties and stress-dependent geotechnical properties of the deep marine sediments retrieved from the hydrate occurrence regions during the Ulleung Basin Gas Hydrate Expedition 2 (UBGH 2) in the Ulleung Basin, East Sea. The measured water contents are in a relatively narrow range of 25-35%, that are close to a low boundary of previously reported values. The sediments tested in this study are classified as either CL/OH or CH/OH with high silt fraction (~80%). Relative high initial void ratio is observed for 6B-9H and 26H resulted from the hydrate existence. The hydrate existence significantly affects on the consolidation characteristics of the specimens. Void ratio, and shear wave velocity show a bi-linear behavior with the effective vertical stress: shear wave velocity increases with increasing in vertical effective stress, and relative high value of the exponent for 6C-9H reflects changes in the interparticle contact during loading, probably due to the disturbance of the specimen by the hydrate dissociation. 5. REFERENCES ASTM D854 (2010). "Standard test methods for specific gravity of soil solids by water pycnometer" Annual Book of ASTM Standards, Vol. 04. 08. ASTM D2216 (2010). "Standard test methods for laboratory determination of water (moisture) content of soil and rock by mass" Annual Book of ASTM Standards, Vol. 04. 08. ASTM D2435/D2435M-11 (2011). "Standard test methods for one-dimensional consolidation properties of soils using incremental loading" Annual Book of ASTM Standards, Vol. 04. 08.

ASTM D4318 (2010). "Standard test methods for liquid limit, plastic limit, and plasticity index of soils" Annual Book of ASTM Standards, Vol. 04. 08. Bahk, J.J., Chough, S.K., Han, S.J. (2000). "Origins and paleoceanographic significance of laminated muds from the Ulleung Basin, East Sea (Sea of Japan)", Marine Geology, Vol. 162, 459-477. Chough, S.K., Jeong, K.S., Honza, E. (1985). "Zone facies of mass-flow deposits in the Ulleung (Tsushima) Basin, East Sea (Sea of Japan)", Marine Geology, Vol. 65, 113-125. Hardin, B.O., Richart, F. E. (1963). Elastic wave velocities in granular soils, Journal of Mechanic and Foundation Division, Vol. 89(1), 33-65. Hillier, S., Son, B.K., Velde, B. (1996). "Effects of hydrothermal activity on clay mineral diagenesis in Miocene shales and sandstones from the Ulleung (Tsushima) backark basin, East Sea (Sea of Japan), Korea", Clay Minerals Vol. 31, 113-126. Kim, I.S., Park, M.H., Lee, Y., Ryu, B.J., Yu, K.M. (2003). "Geological and geochemical studies on the late quaternary sedimentary environment of the southwestern Ulleung Basin, East Sea", Economic and Environmental Geology, Vol. 36(1), 9-15. Kwon, T.H., Lee, K.R., Cho, G.C., Lee, J.Y., (2011). "Geotechnical properties of deep sediments recovered from the hydrate occurrence regions in the Ulleung Basin, East Sea, offshore Korea", Marine and Petroleum Geology, Vol. 28, 1870-1883 Lee, C., Yun, T.S., Lee, J.S., Bahk, J.J., Santamarina, J.C. (2011). "Geotechnical characterization of marine sediments in the Ulleung Basin, East Sea", Engineering Geology, Vol. 117, 151-158. Lee, H. J., Chun, S. S., Yoon, S. H., Kim, S. R. (1993). Slope stability and geotechnical properties of sediment of the southern margin of Ulleung Basin, East Sea (Sea of Japan)", Marine Geology, Vol. 110, 31-45. Lee, J.Y., Santamarina, J.C., Ruppel, C. (2008). Mechanical and electromagnetic properties of northern Gulf of Mexico sediments with and without THF hydrates, Marine and Petroleum Geology, Vol. 25, 884-895 Lee, J.Y., Santamarina, J.C., Ruppel, C. (2010). Volume change associated with formation and dissociation of hydrate in sediment, Geochemistry, Geophysics, Geosystems, Vol. 11, doi:10.1029/2009gc002667 Lee, J.Y., Francisca, F. M., Santamarina, J.C., Ruppel, C. (2010). Parametric study of the physical properties of hydrate-bearing sand, silt, and clay sediments: 2. Small-strain mechanical properties, Journal of geophysical Research, Vol. 115, B11105, doi:10.1029/2009jb006670, Santamarina, J.C., Klein, K.A., Fam, M.A. (2001). Soils and Waves-Particulate Materials Behavior, Characterization and Process Monitoring, John Wiley and Sons, New York. Winters, W.J., Dallimore, S.R., Collett, T.S., Katsube, T.J., Jenner, K.A., Cranston, R.E.,Wright, F.M., Nixon, F.M., Uchida, T. (1999). "Physical properties of sediments from the JAPEX/JNOC/GSC Mallik 2L-38 gas hydrate research well", In: Dallimore, S.R. et al. (Eds.), Scientific Results from JAPEX/JNOC/GSC Mallik 2L-38 Gas Hydrate Research Well, Mackenzie Delta, Northwest Territories, Canada. Bull. Geol. Surv. Can., Vol. 544, 95-100.