REGIONAL GRAVITY INVESTIGATION OF THE EASTERN PORTION OF THE NORTHERN PENINSULA OF MICHIGAN By Erdogan Oray A THESIS Submitted to Michigan State University in partial fulfillment of the requirements for the degree of DOCTOR OF PHILOSOPHY Department of Geology 1971 P L E A S E NOTE: S o m e pages m a y ha ve indi sti nc t print. F i l m e d as rece ive d. University Microfilms, A Xerox Education Company ABSTRACT REGIONAL GRAVITY INVESTIGATION OF THE EASTERN PORTION OF THE NORTHERN PENINSULA OF MICHIGAN By Erdogan Oray A regional gravity investigation was conducted to delineate the Precambrian features of the eastern portion of the Northern Peninsula of Michigan in relation to the tectonic framework of the Lake Superior region and the Southern Peninsula of Michigan. During the survey, an area of 12,000 square miles in the Northern Peninsula of Michigan was covered by a total of 1,003 gravity stations. An additional 851 gravity stations previously observed in the Southern Peninsula of Michigan, Beaver Island, northern Lake Huron, northern Lake Michigan and Sault Ste. Marie area of Canada were used to establish the regional gravity pattern. The Bouguer Gravity Anomaly Map with the aid of a Double Fourier Series Residual Gravity Map indicates that the eastern portion of the Northern Peninsula is associated with two major linear positive gra­ vity anomalies. One of these anomalies trends south from Whitefish Point on the south shore of Lake Superior and is associated with ba­ salts upthrown into a horst within the eastern limb of the Lake Sup­ erior syncline. This anomaly is correlated with the Middle Keweenawan Erdogan Oray volcanics outcropping on Mamainse Point, Ontario. The other anomaly trends southeast from Grand Island in Lake Superior and can be traced magnetically to the Middle Keweenawan volcanics of the Keweenaw Penin­ sula. This anomaly represents the edge of the western limb of the Lake Superior syncline. The edge of the eastern limb of the syncline is al­ so defined by a positive gravity anomaly. These positive gravity anom­ alies which are also associated with positive magnetic anomalies merge in the vicinity of Beaver Island and mark the termination of the Lake Superior syncline. South of Beaver Island, the Keweenawan basalts con­ tinue in a narrow belt and are expressed by the "Mid-Michigan gravity high." The Bouguer Anomaly Map indicates two local gravity minimums in the Whitefish Bay area on the south shore of Lake Superior. These are interpreted to result from a thick accumulation of Upper Keweenawan clastic sediments. The interpretive results obtained from two dimensional model stu­ dies suggest that the Lake Superior syncline in the eastern portion of the Northern Peninsula consists of up to 12,000 feet of basaltic flows overlain by Upper Keweenawan clastic rocks. into a horst in the Whitefish Bay area. The basalts are faulted The presence of yet another horst with a graben to the east occurs in the Grand Traverse Bay region. Theoretical curves from two different geological models can be fitted to the observed anomalies in the northern portion of the South­ ern Peninsula of Michigan. The basalts either extend throughout the northern tip of the Southern Peninsula where they are highly faulted into a series of horsts and grabens or they are confined to the Grand Traverse Bay area in which case pre-Keweenawan extrusives and intrusives underlie the northern tip of the Southern Peninsula. ACKNOWLEDGMENTS The author wishes to express his gratitude to Dr. William J. Hinze for his direction, aid, and suggestions and thorough review of all phases of this study. Special thanks are due to Drs. Hugh F. Bennett, James W. Trow, and Chilton E. Prouty for their helpful suggestions. I am also indebted to Dr. Donald W. Merritt for his assistance in developing computer programs, Robert C. Reed for criticizing portions of the manuscript, and William Marsh for his assistance while collecting the data. Mutch appreciation is due to my wife for her help during the preparation and typing of the manuscript. The financial support which made the gravity survey of the eastern portion of the Northern Peninsula of Michigan possible was furnished by the U. S. Army, Topographical Command and the Cleveland-Cliffs Iron Company. ii TABLE OF CONTENTS Page ACKNOWLEDGMENTS............................. . . . o . . ii LIST OF T A B I E S .................. .................................. LIST OF FIGURES.................................................. vi Chapter I. INTRODUCTION .............................. 1 Nature and Objective............... ... Organization and Scope ........................... Physiography of the A r e a .................... 2 Geology of the A r e a .......................... U General ............................. Lower Precambrian Rocks.................. Middle Precambrian Rocks ..................... Upper Precambrian Rocks. ..................... Lower Keweenawan Rocks . . . . . . . . Middle Keweenawan Rocks................. 9 Upper Keweenawan Rocks ..................... Major Precambrian StructuralFeatures. . . . Origin of the Lake Superior Syncline. . . . Paleozoic Rocks .............................. Cambrian Rocks..........................19 Ordovician Rocks .................... Silurian Rocks. .................... Previous Geophysical Studies ............ . . . Seismic ...................................... Magnetic........................ ... Gravity...................................2U II. U 6 6 7 9 10 12 17 19 20 21 22 22 23 FIELD METHODS AND DATA R E D U C T I O N ............. 27 Survey and I n s t r u m e n t a t i o n ................. 27 Accuracy of Horizontal CoordinatesandElevations . Reduction of D a t a ..........................30 Lake Michigan and Lake Huron GravityData. . . . in. 1 2 INTERPRETATION OF GRAVITY D A T A ................. 32 Sources of A n o m a l i e s ...................... 32 iii \ 29 31 Chapter Page General Characteristics of Gravity Anomalies o o 36 Discussion of Bouguer Anomaly Nap • . 0 » 36 Discussion of the Double Fourier Series Residual Gravity Map 0 • . « 0 . 0 . U2 Comparison of Profiles With the Aid of Upward Continuation .........................1*5 Model Studies o . o « . . . . . o o o 5>2 Profile E*E^ o o o o o o < f l * * « i 55 Profile D - D ^ o o i o o o * t * t « « 5^ An Alternative Geological Model to Profile E-E' £8 Profile E'-E". . . . . • . . • < > • • 60 An Alternative GeologicalModel to Profile E'-E" 62 Profile F—F' 0 . . < > < « . 0 . 0 ^2 An Alternative Geological Model to Profile F-F' 66 Summary of the Results. . . . e o . . . 66 Geological Implications ........................ 68 .. IV. CONCLUSIONS BIBLIOGRAPHY . o o o « « o . s . o o . o 75 o o o o « o . . » . o o . o * o . 78 APPENDIX— FREE AIR GRAVITY ANOMALY MAP OF THE EASTERN PORTION OF THE NORTHERN PENINSULA OF MICHIGAN.................85 iv LIST OF TABLES Page Tentative correlations for the Keweenawan of Northern Peninsula of Michigan (after Green, 1971). • • • « 8 Specific gravities of Precambrian rocks in the Lake Superior Region » . . . • • « . .............. 5>U Parameters of induced, remanent and combined polarization components of the Keweenawan basalts . » • • • • Thicknesses of the Keweenawan rock units in the Lake Superior Region $5 »£6 LIST OF FIGURES Figure Page 1. Location map showing place names ..................... 3 2. Outcrop lithology and basement structures of Michigan and the surrounding area ........................ 5 3. Network of Gravity Base Looping in the eastern portion of the Northern Peninsula of Michigan. . . < > . . . lw Bouguer Gravity Anomaly Map of the eastern portion of the Northern Peninsula of Michigan. . . . . . . 37 Total Magnetic Intensity Anomaly Map of the eastern portion of the Northern Peninsula of Michigan . . . 39 5. 6. Double Fourier Series Residual Gravity Anomaly Map of the Northern Peninsula of Michigan. ............... U3 7-a. The Bouguer and regional gravity curves...................1|6 7-b. Residual values for double Fourier series fit. 8. . . . U6 Selected gravity and magnetic profiles of the eastern Lake Superior Region . . . < > .............. U& Geological model and observed and computed gravity and total magnetic intensity anomaly profiles along E-E' 57 Geological and observed and computed gravity and total magnetic intensity anomaly profiles along D-D 1 . . 57 10. Alternative geological model to profile E-E* . . . . $9 11. Geological model and observed and computed gravity and total magnetic intensity anomaly profiles along E'-E” 61 12. Alternative geological model to profile E'-E". . . . 63 13. Geological model and observed and computed gravity and total magnetic intensity anomaly profiles along F-F' 61+ Alternative geological model to profile F-F' . . . . 67 9-a„ 9-b. Ik* vi Page Figure 15. 16. Schematic map showing interpreted Precambrian structures in the eastern portion of the Northern Peninsula and the northern tip of the Southern Peninsula of Michigan 69 Free-air Gravity Anomaly Map of the eastern portion of the Northern Peninsula of Michigan . . . . . . . 86 CHAPTER I INTRODUCTION Nature and Objective The basement complex of the Midwest has been the subject of con­ siderable recent interest for the purpose of determining the Precambrian geologic history. Exploration for mineral resources has provided val­ uable information about the shallow formations, but few drill holes have penetrated the basement rocks. Therefore, basement studies in the Midwest depend to a large degree upon geophysical surveys and the extra­ polation of information available from a few basement wells. A number of geophysical studies in areas adjacent to the eastern portion of the Northern Peninsula of Michigan where Precambrian rocks are exposed have indicated that considerable gravity and magnetic relief can be correlated with lithologic and structural features of these rocks. The aeromagnetic studies of eastern Lake Superior (Hinze, et al., 1966) which includes the eastern portion of the Northern Peninsula of Michigan have indicated similar results. available in this area. However, only limited gravity data are Therefore, a regional gravity investigation was conducted to supplement the existing geological and geophysical studies of Precambrian geology under the Paleozoic cover in the eastern portion of the Northern Peninsula of Michigan. The general objective of this study is to delineate the Precambrian features of the eastern portion of the Northern Peninsula of Michigan 1 in relation to the tectonic framework of the Lake Superior region and Southern Peninsula of Michigan, while the specific objective is to quan­ titatively investigate the extension of the Lake Superior basin into the area of study. Although available magnetic data have been considered in the quantitative structural interpretations, the main source of in­ formation is the gravity data. Organization and Scope An area of 12,000 square miles in the Northern Peninsula of Mich­ igan was covered during the survey (Figure 1). A total of 1,003 gravity stations were established using a LaCoste and Romberg, Model G, Geodetic Gravimeter. The station spacing varied from one mile in areas where detailed coverage was considered necessary for geological interpretation to a regional coverage of four miles. Observations were taken at road intersections where elevations are specified on U.S. Geological Survey maps or at "bench marks;" Elevations at 93 stations were determined with the aid of aneroid altimeters. An additional 8£l gravity stations previously observed in the Southern Peninsula of Michigan, Beaver Island, northern Lake Huron, northern Lake Michigan and Sault St. Marie area of Canada were used to establish the regional gravity pattern. Specific information regarding the collection, reduction, presen­ tation and methods of interpretation of the gravity data will be dis­ cussed in sections treating each of the subjects individually. Physiography of the Area The study area constitutes the eastern lowland topographic province of the Lake Superior region (Hamblin, 1958) which occupies the entire Northern Peninsula of Michigan east of 87° U5'W. Although the greater part of this area is covered by glacial drift, it has the form of a LOCATION MAP SHOWING PLACE NAMES ’ SCo, V SOUTHERN PENINSULA Figure 1.— Location map showing place names. OP MICHIGAN k gently southward dipping cuesta. Au Train, Laughing Whitefish and Rock Rivers are the only breaks through this cuesta which has formed on the resistant Au Train formation, of Upper Cambrian or possibly Ordovician age. A number of waterfalls have been formed by north-flowing streams in Alger County as they cross the cuesta. From Munising eastward to Beaver Lake the face of the cuesta follows the Lake Superior shoreline and forms the famous Pictured Rocks. The surface elevation within the lowland topographic province no­ where rises more than 900 feet above sea level or 300 feet above Lake Superior. Major elevation differences are due to differential erosion of the resistant Au Train formation. In many places bedrock is covered by lake clays.and recessional and ground moraines as much as 1*00 feet thick. The drainage in general is poor and swamps cover large areas in Schoolcraft and Luce Counties. Several drumlins east of Waucedah and Foster City and also south of Chatham, 30 miles southwest of Munising in Alger County form striking topographic features. A more comprehensive description of the physiography of Northern Peninsula of Michigan is presented by Van Hise and Leith (1911). Geology of the Area General The study area represents the northern extremity of the Michigan structural basin which also includes the Southern Peninsula of Michigan, eastern Wisconsin, northeastern Illinois, northern Indiana, northwestern Ohio and southwestern Ontario. Cambrian, Ordovician and Silurian sedi­ ments conceal the Precambrian basement rocks for an east-west distance of 1^0 miles from Marquette to Sault Ste. Marie, Figure 2. Within this 5 tart«lh« -''"-'•jL— tfflclln* .— l i i i a n t CwtH u I M M M tffa*-niu) i N t ^ r s fM lta UmM m KM Ht0k-wf>* ii*uu M U TlKWt f > T 7 Figure 2.--Outcrop lithology and basement structures of Michigan and the surrounding area0 \ 6 area few basement drill tests are available for a direct investigation of the crystalline basement. However, the general nature of Precambrian formations east and west of the study area are similar in gross aspect. Therefore, geologic information from the adjacent Precambrian exposures is presented to establish the regional Precambrian geology of the area. There is no widely accepted reference framework of eras and periods to correlate Precambrian lithologic units and geologic events. There­ fore, in discussing Precambrian geology from different areas the termi­ nology used by the U.S. Geological Survey will be followed (James, 1958), the Precambrian being subdivided into lower (Archean), middle (Animikie) and upper (Keweenawan) divisions and the Keweenawan in turn subdivided into three groups, lower, middle and upper. Lower Precambrian Rocks-Archean Series The oldest rocks found in Michigan occur in the western extremity of the study area. Kitchi and Mona schist in the Marquette area and Norway Lake granite near Norway in Dickinson County are thought to be of Lower Precambrian age. These rocks mainly consist of gneisses, gran­ ites, greenstones and other metavolcanic rocks. One or more periods of metamorphism and deformation occurred before the end of early Precambrian time and the period was finally brought to an end by the Algoman orogeny (2.U b.y.) causing a major unconformity between the Archean rocks and the overlying Middle Precambrian rocks. Middle Precambrian Rocks-Animikie Series The Middle Precambrian rocks of Michigan are correlated with the Animikie Group of northeastern Minnesota and the Thunder Bay area, Ontario (James, 1958). These rocks consist of a thick sequence of metasediments and metavolcanics and scattered bodies of metagabbro and 7 of granitic rocks. This sequence was referred to as the Huronian series until James established its correlation with Animikie series. Middle Precambrian time was a period of extensive sedimentation with minor volcanic activity. Iron formations deposited during this period are an important source of iron ore. They locally produce magne­ tic anomalies in excess of 20,000 gammas. The Penokean orogeny which caused widespread deformation and regio­ nal metamorphism brought to a close Middle Precambrian time (1.6 b.y.). During this orogeny, emplacement of small dikes and sills of basic rocks, such as diabase and gabbro, was followed by intrusions of granite and pegmatite dikes. Crustal compression then folded, faulted and metamor­ phosed the Animikean rocks resulting in east-west trending folds and metamorphic belts. Upper Precambrian Rocks-Keweenawan Series Upper Precambrian rocks are widely and extensively exposed in the Lake Superior region and from drill hole and geophysical data they are believed to subcrop beneath the Paleozoic sediments over extensive areas of the Midwest. These rocks are best known on the Keweenaw Pen­ insula due to intensive prospecting for copper mineralization. There­ fore, this area serves as the type stratigraphic section for most of the Keweenawan sequence. The Keweenawan rocks around Lake Superior can be correlated either directly or indirectly with the Northern Peninsula sequence as shown in Table 1 (Green, 1971). Middle Keweenawan igneous rocks have been dated isotopically at 1.2-0.9 b.y. (Goldich, et al., 1961). Since the dated interval repre­ sents only a part of Keweenawan time, it can be assumed that Keweenawan time extends from at least 1 .k to between 0.9 and 0.6 b.y. (Halls, 1966). TABLE leoo Tentative correlations for the Keweenawan of Northern Peninsula of Michigan (after Green,1971) MINNESOTA Northern Wisconsin Hie Royate Bayfield Gpi Hinckley Set. UPPER Fond du Lac Sst. Oronto Gp Upper M tc h ig o n Jooobsville Ssl Nonesuch Sh. Eastern L. Super. O n to ri o Joeobsville Sst. sandstone Copper Harbor Conqlom. Portage Lake la v a s o lava series P o rtage Lake leva series Michipicoten I. lavas M i DOLE ^[sandstone? 2 / felsits? Mellen- Hurley lavas Sdsth Range la v * series lavas LOW ER Series ANIMIKIE Rove Fm. Thomson Fm. Virginia Fm Rove Fm quartz sandstone T y le r SI. quaMi sandstone T y le r Slate Alono Bay lavas 9 Lower Keweenawan Rocks-The sediments which are younger than the Ani­ mikie series and underlie the earliest Keweenawan lavas are assigned a Lower Keweenawan age. They are well developed on the Sibley Peninsula, east of Thunder Bay on the north shore of Lake Superior, and consist predominantly of sandstones, marls and shales. These sediments thin rapidly to the southwest and are represented by a thin conglomerate near Duluth. Further south, the Barron quartzite of Wisconsin and the Sioux formation of southern Minnesota are correlated and thought to be of Lower Keweenawan age (Goldich, et al., 1961). Middle Keweenawan Rocks-On the Keweenaw Peninsula, the Middle Kew­ eenawan rocks are known as the Portage Lake lava series. They range in composition from basalt to andesite and are interbedded with thin layers of conglomerate and sandstone. Over 300 flows occur in the upper part of this sequence, with an average individual thickness of 1*3 feet (White, 1966a) although much thicker flows have been recognized. This leads to an estimated total thickness of between 13>,000 to 25,000 feet. This enormous thickness gradually decreases away from the axis of the Lake Superior syncline indicating a progressive subsidence of the syncline contemporaneous with the accumulation of lavas. Another characteristic of the Middle Keweenawan is that movement of lava flows appears to have been outward from the center of the lake according to evidence supplied by bent "pipe" amygdules in the lava (Hotchkiss, 1923j Butler and Burbank, 1929) whereas interflow sediments show current flow in the opposing dir­ ection according to current bedding and imbricate structure studies (Hamblin and Horner, 1961). White (1957, I960) explained this by periodic reversals in slope within the basin. First, flows spread out from the central parts of the basin toward the margins over a flat surface. \ Then 10 due to this load, the surface of the crust subsided forming a physiographic basin which accommodated the deposition of conglomerates and sandstones by currents flowing toward the center of the basin. Intrusives within the Middle Keweenawan are widespread in the Lake Superior region, the largest being the Duluth gabbro with a total thick­ ness of about 20,000 feet (Taylor, 1961*) 0 Numerous diabase dikes have been observed intruding the Keweenawan volcanics. Although some unmet­ amorphosed diabase dikes have been found cutting the Middle Keweenawan sequence of Michigan they are most frequently found in pre-Keweenawan rocks (Halls, 1966)0 Some of these dike swarms exhibit strong reverse magnetization. Patches of Middle Keweenawan lavas also occur along the eastern shore of Lake Superior. This indicates the continuity of these rocks in an easterly direction from Keweenaw Peninsula. North of Whitefish Bay at Mamainse Point, there is an apparent thickness of 16,000 feet of amygdaloidal lava flows with interbedded boulder conglomerates. This sequence is correlated with the Portage Lake lavas of Michigan. At Mamainse Point, the formation strikes generally north-south and is inclined lakeward at angles not exceeding 30°° also been found in this area (Thomson, Dikes and felsite intrusions have 1953)« Extensive exposures of Middle Keweenawan rocks are also found on Michipicoten Island in eastern Lake Superior. lavas on the island reach 11,000 feet. The maximum thickness of This Middle Keweenawan sequence is unique because acidic lavas make up at least one-half of the section (Halls, 1966). Upper Keweenawan Rocks-The Upper Keweenawan sequence consists pre­ dominantly of fine grained sandstones and shales, overlying a thick \ 11 conglomerateo The basal unit, the Copper Harbor conglomerate, lies con­ formably over the Middle Keweenawan lavas. rhyolite cobbles. It is mainly composed of Within the Copper Harbor conglomerate are lenses of fine to coarse grained, red arkosic sandstone. Harbor varies considerably. The thickness of Copper However, on the Keweenaw Peninsula it retains a uniform thickness between U ,000 and 5,000 feet. A sharp lithologic change from conglomerate to siltstone marks the boundary between Copper Harbor conglomerate and the Nonesuch shale. Al­ though there is no structural discordance, this contact is thought to be a disconformity because of the significant difference between the paleomagnetic pole positions of the two formations (Du Bois, 1962). The Nonesuch shale retains a constant thickness of about 600 feet along the Keweenaw Peninsula. The Nonesuch is overlain by 5>000 feet of Freda formation. This formation consists of interbedded arkosic sandstone and red micaceous silty shale. The Freda reaches a thickness of 12,000 feet in Wisconsin and is known as the Oronto group0 of heavy minerals (7»b percent). It contains a large average percentage This characteristic coupled with an abundance of feldspar indicates that the sediments were subjected to in­ complete weathering and short distances of transportation. Therefore, nearby Keweenawan volcanics and pre-Keweenawan rocks are the likely source rocks of Freda formation. Overlying the Freda is a sequence of red quartzose sandstone and micaceous shale more than 2,000 feet thick known as Jacobsville sandstone in Michigan. Seismic reflection data (Bacon, 1961*) indicate thickness of this formation can be as great as 12,000 feet. that the The relation­ ship of the Jacobsville with the Freda is still conjectural, however, the \ 12 Jacobsville appears to be younger as pebbles of red siltstone and mud­ stone possibly derived from the Freda are found in it along the coast from L'Anse to Marquette (Hamblin, 1958)* The Jacobsville may therefore be considered as having a late Upper Keweenawan or Lower Cambrian age. Further evidence to support this conclusion is the low dip angular un­ conformity between the Jacobsville and the Munising formation of Late Cambrian age on the east side’ of Grand Island (Hamblin, 1958)« However, it is impossible to estimate the time gap represented by this unconformity and therefore neither Upper Keweenawan nor Cambrian age can be assigned to Jacobsville with any certainty.. The more recent studies of this long existing Jacobsville problem seem to favor its equivalence with Bayfield group in Wisconsin which is of Upper Keweenawan age 0 The Jacobsville sandstone may be traced eastward from the Keweenaw Peninsula to the Sault Ste. Marie area as narrow discontinuous shoreline outcrops (Hamblin, 1958)* At all these occurrences the dip is very gentle and generally does not exceed 1° or 2° in a lakeward direction. In addition to quartzose Jacobsville sandstones, feldspathic sandstones which resemble the Freda formation lithologically are found in the Batchawana Bay and Sault Ste. Marie areas. However, there is doubt about the time-equivalance of these feldspathic rocks to the Freda formation (Hamblin, 1958). Major Precambrian Structural Features The dominant structural feature in the Lake Superior region is the Lake Superior basin known as the Lake Superior syncline. In the western half of the lake, the syncline trends northeasterly and is strongly as­ ymmetrical toward the south. The Keweenawan rocks within the syncline show dips ranging from only a few degrees to vertical. \ The highest dips 13 occur in the vicinity of major thrust faults which strike parallel to the axis of the syncline. The structure contours clearly show closure of the syncline south of Duluth although Keweenawan rocks continue south­ ward in a narrow trough to Kansas according to gravity, aeromagnetic and borehole data (Thiel, 1956; Lyons, 1959; Craddock, et al., 1963; Zietz, 1965)o rather complex. At the western end of the basin the structural pattern is A north-south trending ridge extends from the vicinity of MBllen, Wisconsin to the Minnesota shoreline (White, 1966a). Wold and Ostenso (1966) detected the presence of another ridge which extends from Isle Royale southwestward to the ridge interpreted by White. Halls (1966) suggested that much of the western half of the Lake Superior syncline developed on the site of a former geosyncline, evolved during the Animikean and terminated by Penokean orogeny which had little effect on the north shore of Lake Superior. The Animikean geosyncline cannot be traced eastward because there are no known Animikean rocks around eastern Lake Superior. However, the Middle Precambrian rocks occurring along the north shore of Lake Huron may be equivalent of Animikean (Young, 1966) indica­ ting that they may form the northern margin of an easterly extension of the Animikean geosyncline. In contrast to this former geosyncline, the axis of the Lake Superior syncline parallels the Keweenaw Peninsula in western Lake Superior and then gradually turns and becomes southerly in the eastern part of the lake according to aeromagnetic data (Hinze, et al., 1966)o The gravity evidence of the present study indicates closure of the syncline east of Beaver Island and continuation of Keweenawan rocks southward in a narrow belt into the Southern Peninsula of Michigan. One of the major structural features of the Lake Superior region is the Keweenaw fault. It extends through the center of the Keweenaw lit Peninsula bringing the Middle Keweenawan flows to the north into contact with the Upper Keweenawan sediments. The age, extension and the amount of displacement of the Keweenaw fault has been a matter of conjecture for many years. The amount of displacement necessary to bring the Middle Keweenawan basalts in contact with the Jacobsville is in the order of 15,000 feet (Butler and Burbank, 1929). Halls (1966) believes that the fault is at least 10,000 feet downthrown to the south in the Keweenaw Peninsula, confirming the earlier studies. However, a much smaller displacement is suggested by Hamblin (1958) who argues that the Keweena­ wan series were tilted and eroded before the deposition of Jacobsville sandstone. Therefore, the minimum amount of displacement on the fault is not the thickness of the Upper Keweenawan rocks but only the thickness of the Jacobsville which is in the order of magnitude of 2,000 feet. Hamblin (1958) also believes that the age of the Keweenaw fault is postJacobsville and probably post-Devonian since Devonian and Jacobsville rocks at the Limestone Mountain area are affected by the same folding process which is related to the compressional forces that caused the faulting. However, generally the Keweenaw fault is believed to be related to the origin of the Lake Superior syncline. Therefore, the post-Devonian fault­ ing at the Limestone Mountain area is likely a rejuvenation of movements along the old lines of weakness. The genetic relationship of the Keweenaw fault to other faults in the area and its possible eastward extension is also a controversial matter. The Lake Owen fault of northern Wisconsin which occurs to the southwest of the Keweenaw fault is thought to be an extension of it. the east, O'Hara (1967) interprets from aeromagnetic data a southeast­ erly trending fault extending from Keweenaw Point through Au Train Bay \ To 15 to the Manistique area on the Lake Michigan shore and suggests that this fault may be an eastward extension of the Keweenaw thrust fault. Thwaites (1935) on the basis of bottom topography studies of eastern Lake Superior suggested the existence of such a fault. Further geological evidence for extending the Keweenaw; fault in this manner came from Oetking's (1951) studies which show minor thrust faults in the Jacobsville sand­ stone on the east side of Laughing Fish Point, Au Train Island and in the southeastern corner of Au Train Bay. In addition, Patenaude (19610 suggested the existence of a nearly vertical fault on Au Train Island. Although the geological evidence for extending the Keweenaw fault further southeast from Au Train Bay is weak, small scale faulting near Seul Choix Point on Lake Michigan shore has been observed to substantiate this interpretation (Patenaude, 1961*). Newcombe (1933) traced the Keweenaw fault to Stannard Rock and suggested that the fault possibly extended eastward® He further stated a genetic relationship between the Keweenaw fault and the Murray fault in Ontario. However, the apparent direction of movement on these two faults is opposite and such a relationship does not appear to be feasible on the basis of geophysical evidence. Patenaude (1961*) as a result of his geophysical studies in the eastern extremity of the eastern portion of the Northern Peninsula pro­ poses the existence of a north-south trending fault of major proportions. This fault extends southward into Lake Huron and is upthrown to the west. Patenaude further suggests that this fault may be related to the Keweenaw fault system. The discontinuous configuration of the magnetic anomalies near the north shore of Northern Peninsula has led O'Hara (1967) to suggest the \ 16 existence of nearly east-west trending faults in this area. Patenaude's (1961|) studies also revealed such an east-west striking fault south of Whitefish Bay. Local east-west trending pre-Keweenawan structures predominate in the western extremity of the study area. Two such structures are the Marquette syncline and Gwinn trough which are considered to be areas of complex folding and faulting, probably the result of a period of a major orogeny. Frantti (1956) on the basis of geophysical data interpreted a similar syncline in the vicinity of Trenary which has a major magnetic anomaly associated with it. The association of iron-formations with folded Animikean rocks in this part of the study area is a rule rather than an exception. A metasedimentary belt extends from Waucedah to Es- canaba (Allen, 191U) and is clearly seen in magnetic maps. Another school of thought concerning the nature of the Precambrian basement between Marquette and Sault Ste 0 Marie has also been considered in the geological literature. Hamblin (1958) ignoring the possibility of a basement of Keweenawan volcanics, endorsed the existence of an eastwest trending ridge of Middle Precambrian rocks based on sediment dis­ persal patterns in the Jacobsville sediments. This structurally posi­ tive area which is called the Northern Michigan Highland is believed to have existed in the eastern portion of the Northern Peninsula of Michi­ gan during Jacobsville time. Hamblin (1961) states that this highland was probably formed near the beginning of Keweenawan time, either as a result of uplift by faulting or arching of the Northern Peninsula or by differential sinking of the region beneath the eastern Lake Superior, and persisted as a highland until Middle Cambrian time. As a further evi­ dence, Hamblin (1958) interprets the broad east-west gravity high of 17 Bacon's (19£7) reconnaissance survey in the Northern Peninsula to be caused by a buried PreCambrian ridge extending from the Wisconsin Arch to Canada* Origin of the Lake Superior Syncline Several theories have been developed to explain the origin of the Lake Superior syncline* by Hotchkiss (1923)* An early, commonly accepted theory was put forth On the basis of estimates of the original dips of the Keweenawan flows, together with flow directions as indicated by bent "pipe" amygdules and sedimentary structures at the base of the extrusions, Hotchkiss concluded that the source of the flood basalts was along the axis of the Lake Superior syncline* He further envisions a batholithic intrusion during pre-Keweenawan time located beneath the present trough which had a surface expression as a topographically positive area. Dur­ ing the Middle Keweenawan time, escaping basaltic material produced great thicknesses of extrusives on the surface which led to an equivalent sub­ sidence of the crust. The escape of basaltic material together with contraction of the remainder of the batholith due to cooling and loss of volatiles resulted in complete collapse of the batholithic roof. Hotchkiss' theory on the origin of the Lake Superior syncline agrees with Van Hise's (1911) conclusion that the Lake Superior syncline pro­ bably originated during Middle Keweenawan time. Van Hise further stated that the folding of the basin was practically complete by the end of the Keweenawan time and large scale strike faulting affected the area in post-Cambrian and possibly in post-Cretaceous time. However, Irving's work (1883) shows an earlier Animikean sedimentary basin with essentially the same structural features. Halls (1966) confirms the early studies of Irving and concludes that the western half of the Lake Superior syn­ cline developed on the site of an east-west trending Animikean geosyncline. 18 For many years it has been hypothesized that the Lake Superior and the Michigan basins are interrelated. Robinson (1923) in an unpublished manuscript suggested an early connection between the Lake Superior syn­ cline and the Michigan basin. He explains the later separation of the two basins in the following way: a synclinal area extended from St. Paul, Minnesota through the present Lake Superior region into the Southern Peninsula of Michigan which was filled with Cambrian sandstone; a period of early Paleozoic erosion removed a large part of the sandstone in the Southern Peninsula but did not similarly affect the Lake Superior area; this was followed by deposition of thick sediments in the basin; a sec­ ond period of erosion in possibly Tertiary time affected the Lake Supe­ rior basin and removed the Paleozoic sediments from the area; the re­ sult of this erosion was the formation of a ridge not far north of the present south shore of Lake Superior which separated the two basins. Newcombe (1933) accepting the idea of Robinson, proposed a new hypothesis to explain the separation of the Lake Superior syncline from the Michi­ gan basin. He concluded that the major Keweenawan movements resulted in the separation rather than a complex erosional history. A possible connection between the eastern Lake Superior basin and the Michigan structural basin also was suggested by Thwaites (1935). His east-west cross section in the eastern portion of the Northern Pen­ insula shows Keweenawan flows slightly downwarped in the center and thrust faulted near either end. Thwaites' other cross sections also show the Keweenaw fault entering at Au Train Point and two other northsouth faults near the eastern extremity of the Peninsula. Most geologists and geophysicists presently working on various seg­ ments of the Lake Superior syncline attribute its origin and development \ 19 to one of two hypotheses. The first suggests that the basin was formed primarily by downwarp of the crust either by loading when mafic lavas were extruded in great volume, or by crustal compression (White, 1966b), The second and recently more favored hypothesis is that, crustal rifting has taken place during the accumulation of the flood basalts and sediments presently filling the trough. This second hypothesis explains very well the thick crust beneath the lake by attributing its cause to isostatic adjustment to the high density basic material which appears to have in­ truded the crustal rocks (Smith, et al., 1966). Hinze, et al. (1971), associating the origin of the Lake Superior syncline to continental rifting, have compared the Lake Superior region with the East African rift system. They find strong similarities be­ tween the two areas and attribute the observed dissimilarities to dif­ ferent stages of development of the rifting process. The above theories exemplify the variation of ideas proposed to explain the origin of the Lake Superior syncline. The structural ex­ tent, framework and time of origin of this feature have been the subject of many geological and geophysical studies. Paleozoic Rocks In the eastern portion of the Northern Peninsula of Michigan sedimentary formations of Paleozoic age overlap the Precambrian basement complex and dip gently south into the Michigan Basin. Beds younger than late Sil­ urian are not present in the Northern Peninsula. Cambrian Rocks-The Miinising formation which is given a Late Cambrian age on the basis of its fossil content represents the Cambrian rocks in the area. The contact between the Munising formation and the older rocks is well exposed in numerous outcrops in Alger and Dickinson Counties. 20 In Alger County the Munising rests upon the Jacobsville formation with a low dip angular unconformity (Hamblin, 1958). In Dickinson County the Jacobsville is absent and the Munising lies upon the highly deformed Animikean rocks with a major unconformity. The Pictured Rocks cliffs constitute the principal exposure of the Munising formation which consists of, in ascending order, a basal con­ glomerate, the Chapel Rock member and the Miner's Castle member. The basal conglomerate is an orthoquartzite and attains a maximum thickness of 15 feet. It wedges out southward and, therefore, is restricted to the northern part of the Northern Peninsula. The Chapel Rock member overlies the basal conglomerate and consists of well sorted medium grained sand­ stone characterized by large-scale cross-bedding. Along the Pictured Rocks the Chapel Rock is 1*0 to 60 feet thick and thins gradually to the easto The contact between the Chapel Rock and Miner's Castle members is sharp as indicated by significant changes in sorting, sedimentary struc­ tures and heavy mineral assemblages. Therefore, the source areas of these two members were probably different and a disconformity separates them (Hamblin, 1958). The Miner's Castle member consists of poorly sorted sandstone and is 11*0 feet thick at the type locality in the Pic­ tured Rocks. Ordovician Rocks-The Munising formation throughout the Northern Peninsula is overlain by a sequence of sandy dolomites and dolomitic sands. The age of these beds is a matter of conjecture. Van Hise and Bayley (1900) proposed the term Hermansville for the limestones which overlie the Lake Superior Sandstones in Wisconsin and the Munising for­ mation (Upper Lake Superior Sandstone) in the Menominee district of the Northern Peninsula of Michigan. These limestones were described as 21 being sandy and dolomitic and assigned a Lower Ordovician age (now Prairie Du Chien group)* The Au Train formation described by Hamblin (19E>8) at Au Train Falls near Munising overlies the Munising formation and has the stratigraphic position and general character of the Trempea­ leau formation of Upper Cambrian age in the Southern Peninsula. However, Oetking (195>1) described fossils which are believed to be of Middle Or­ dovician age from the Au Train formation near Miner's Castle, Munising* Therefore, the age of the Au Train is in doubt and its relationship to the Hermansville is questioned, as there is no contact relations observed in the small, isolated outcrops of the Northern Peninsula* A disconfor- mity must exist between the Munising formation and well established Black River limestones of the Middle Ordovician age as the Lower Middle Ordo­ vician (Chazy) has never been recognized in the Northern Peninsula. The time represented by this disconformity depends on the true age of the Hermansville and the Au Train which at present remains problematical. Overlying the Au Train formation in the Northern Peninsula is either the Black River or Trenton, depending on the age of the Au Train* Major outcrops of Trenton occur in the vicinity of Escanaba, Menominee and the northern part of Drummond Island* The Trenton group consists of dolo­ mitic limestone and limestones alternating with cherty layers. The de­ position of these rocks is followed by the deposition of the Late Ordo­ vician rocks which consist of shales and fossiliferous limestones over extensive areas of the upper Midwest indicating uniform depositional conditions. Silurian Rocks-During Silurian time seas covered the area and laid down a variety of sediments, the most common being carbonates formed in clear seas. basin0 Reef formation was common around the margin of the Michigan Silurian time in Michigan was also a period of considerable down- warpingo Early Silurian rocks consist largely of shales and carbonates where­ as the Middle Silurian rocks are predominantly dolomite and occur in an arcuate belt from St. Martin, Poverty, Big and Little Summer Islands at Green Bay to east of Drummond Island. The uppermost formation of Middle Silurian age, the Engadine dolomite is a scarp former, forming the Gar­ den Peninsula southwest of Manistique and the Drummond and the Manitoulin Islands. The single greatest episode of sinking and the accumulation of thick evaporites, carbonates and shales, mark Upper Silurian time in the Mich­ igan basin. Devonian rocks do not occur in the Northern Peninsula. However, the so called Mackinac Breccia found in the vicinity of St. Ignace con­ tains Devonian fragments indicating that the Devonian rocks were once present in the area. Previous Geophysical Studies Seismic A seismic survey to study the deep crustal structure of the North­ ern Peninsula of Michigan and Wisconsin was undertaken by Slichter as early as 191*0 and was reported in 1951 (Slichter, 1951)» The results of this investigation included one shot point at Palmer in the Marquette area and another at Manistique, but none in the eastern portion of the Northern Peninsula showed that: 1. There is a superficial layer with a congressional wave velocity of 1*.16 km/sec and a thickness of about 3 km0 23 2. The superficial layer is underlain by a homogeneous layer of velocity 6.16 km/sec and a thickness varying from 33 to 38 km. 3. The depth to Moho in the Palmer area is 1*0 km. Slichter used a velocity of 8.17 km/sec for mantle material in calculating the depth to the Mohorovicic discontinuity. In connection with the International Upper Mantle Project, during the summers of 1963 and 1961* refraction data were obtained along a line crossing Lake Superior from Duluth, Minnesota to Otter Cove on the east­ ern shore of Lake Superior in Ontario. The results of various interpre­ tations of these data show that the Mbho discontinuity deepens from about 36 km in the western part of Lake Superior to 60 km in the region north­ west of Keweenaw Point and there is no significant change in depth fur­ ther east (Meyer, 1961*; Smith et al0, 1966; Berry and West, 1966; Meyer and Ocola, 1967; O'Brien, 1968). There is an upper refractor at a depth of 7 km and the compressional wave velocity under the upper refractor is unusually high (6.8 km/sec) while that of the upper mantle is normal (80I km/sec). Magnetic Patenaude (1961*) reported seven east-west and six north-south aeromagnetic traverses in the eastern portion of the Northern Peninsula of Michigan. His work indicates that Keweenawan volcanics extend to the eastern limit of the Peninsula. Hinze, et al., (1966) conducted a regional aeromagnetic survey of the eastern Lake Superior and eastern half of the Northern Peninsula of Michigan. The results of the survey support the geological interpretation that the Lake Superior structural basin consists of thick basic volcanics overlain by clastic sediments and extends southward into the Northern Peninsula of Michigan. 2h Case and Gair (1965) in their aeromagnetic study of parts of Mar­ quette, Dickinson, Baraga, Alger and Schoolcraft Counties, Michigan cor­ related the major magnetic anomalies and broad areas that have character­ istic magnetic patterns with the geology as determined from published reports. This study was extended by U.S. Geological Survey (1970) to cover the neighboring areas. The resultant aeromagnetic maps reflect the east-west trending pre-Keweenawan structures. Kellogg (1971) has studied the basement geology of the Southern Peninsula of Michigan using the results of geophysical anomalies from the periphery of the Michigan basin. He correlates a large positive gra­ vity anomaly and the discontinuous magnetic anomaly, the "Mid-Michigan anomaly," with Middle Keweenawan basic volcanic rocks differentially up­ lifted against the gneisses, granites and metasediments of the older basement provinces. Hinze and Merritt (1969) concluded that this large positive feature extends northwards into the Lake Superior region. O'Hara and Hinze (1971) have reported on an aeromagnetic investi­ gation of Lake Michigan, Green Bay, Traverse Bay and the eastern portion of the Northern Peninsula of Michigan immediately north of Lake Michi­ gan. In this study, they trace the east-west striking magnetic anomalies in northern Michigan into generally north-south striking anomalies which they correlate with Keweenawan basalts in the Lake Superior basin. Gravity A reconnaissance gravity survey of the Northern Peninsula of Mich­ igan was conducted by Bacon (1957). The major feature of Bacon's Bouguer Anomaly Map is a broad southeast trending anomaly which was interpreted as caused by a ridge between the Lake Superior and Michigan basins. Bacon suggested that this gravity anomaly may merge with the "Mid-Michigan 25 gravity anomaly" and was similar to the "Mid-Continent gravity high" which extends from the western end of the Lake Superior basin south through Minnesota, Iowa, Nebraska and Kansas. Patenaude (1961*) analyzed four gravity profiles, two of them in Pickford-Sault Ste. Marie area and the other two in Au Train-Munising area. The results of his model studies show a fault of major proportions in the vicinity of Au Train Point which is genetically related to the Keweenaw fault. Patenaude also suggests a high angle fault about 8 miles east of Pickford which also may be related to the Keweenawan tectonics. A Bouguer Anomaly Map covering the midwestern states and southern Ontario was compiled by Rudman, et al., (1965)*. This map shows that the "Mid-Continent gravity high" is genetically related to the linear "MidMichigan gravity anomaly." Gravity surveys were conducted in Lake Superior during the summers of 1963 and 1961*. The Bouguer Anomaly Map compiled by Weber and Good- acre (1966) is characterized by steep gradients which can be best ex­ plained by shallow gologic structures no more than 10 km in depth from the surface. Berkson (1969) conducted a gravity survey in the vicinity of Michipicoten Island, Lake Superior. His studies show a positive belt of gravity values trending southeast flanked by a negative anomaly to the northeast. The linear positive gravity anomaly is interpreted to result from high density Keweenawan lavas of the north limb of the Lake Superior syncline while the negative anomaly is attributed to a low density area of the pre-Keweenawan basement under the lake. The geological and geophysical studies up to the present indicate that Keweenawan rocks extend southward from the Lake Superior basin and the structures associated with the "Mid-Continent" and "Mid-Michigan" \ 26 anomalies are similar<> These conclusions have led to speculation that the Lake Superior syncline had its origin in a rift of the continental crust. Hinze, et al., (1971) have compared the North American Mid- Continent basement with the East African Rift System and have proposed mechanism for the development of continental rifts. Their proposed pro cess explains the geological and geophysical observations of the East African Rift System to its present stage of development and of the late Precambrian paleo-rift system of north-central United States0 \ CHAPTER II FIELD METHODS AND DATA REDUCTION Survey and Instrumentation The gravity observations of the survey were taken with a La Coste and Romberg, Model G, No. 103 Geodetic Gravimeter. It has a range of 7000 milligals and has a reading accuracy of iO.Ol mi H i gal. The gravi­ meter was provided by the U 0S 0 Army, Topographical Command. Seventeen base stations were established during the course of the survey to expedite the field work. These stations were distributed throughout the survey area near accessible roads and away from busy high­ ways. The magnitude of the instrument drift was found by reoccupation of base stations at four hour intervals. Each base station was triple looped as shown in Figure 3, in four closed systems. The closure errors of the systems are as follows: System 1 2 3 h Station Closure Error (mgal) SB 0B MB KB 0.01 0.005 0.015 0.0 The base station at Ozark was also triple looped with North and South Mackinac Bridge base stations. The station at the south end of Mackinac Bridge was tied to the base station WU 17 in front of the PhysicsAstronomy Building at Michigan State University, East Lansing (Behrendt 27 \ 87 47 47° LAKE S U P E R IO R NEWBERRY GREEN GTw PEN SENEY SYSTEM SHINGLETON SYSTEM 2 RUDYARD SYSTEM 3 OZARK SYSTEM H*SBC STALWART ROCK 46' PIKE LAKE, 46° WATSON NAHI ION NORTH MACKINAC BRIDGE LAKE M ICHIG AN SOUTH MACKINAC BRIDGE SPALDINI KB "STEPtClfeoN NETW ORK OF G RAVITY BASE LOOPING SCALE 20 30 40 MILES 45 87 86' 85* 84* Figure 3o— Network of Gravity Base Looping in the eastern portion of the Northern Peninsula of Michigan.: ro oo 29 and Woollard, 1961). Gravimeter drift for this tie was eliminated by the method described by Woollard and Bonini (l955)* which required es­ tablishing the slope of the drift curve by hourly interval observations at stations at the south end of the Mackinac Bridge, Harrison Rest Area* Mount Pleasant and the Physics-Astronomy Building, Michigan State Univer­ sity. The Harrison Rest Area and Mount Pleasant stations are located at roughly equal intervals between the ends of the loop. This tie was conducted four times with a maximum discrepancy of 0.02 milligal when earth tides were eliminated. The station in front of the Physics-Astronomy Building at Michigan State University was then tied by triple looping to airport base station WA 117 at the Lansing Airport (Behrendt and Woollard, 1961). A total of 26 stations were reoccupied and after applying the drift correction, the maximum discrepancy for reoccupied stations was 0.01* milligals with an average of 0.015 milligals. During the course of the survey, the gravimeter sensitivity and level adjustments were checked once a week. The displacement sensitivity was 0.7 and the level setting remained correct throughout the survey. Accuracy of Horizontal Coordinates and Elevations Gravity stations with spacings varying from 1 to 5 miles were plotted from U.S. Geological Survey maps of 7o5 and 15 minute quadrangles and Michigan State Highway Department maps of scale 1 :62,£00. Highway Department maps were used only where U.S. Geological Survey maps were unavailable. Horizontal coordinates (latitude and longitude for each station) were calculated from these maps using a different latitude constant for 2.5 minute intervals for 7*5 minute series, 5 minute intervals for 15 minute series and 10 minute intervals for State Highway maps. Stations 30 plotted from U.S. Geological Survey maps have an error of *1 second and stations plotted from Highway Department maps have an error of -3 seconds in their horizontal coordinates. The accuracy of locating the stations on the maps is of the order of 2£0 feet. seconds. This causes an error of 3 to The total error in the horizontal coordinates is estimated to be —14. seconds for stations plotted from U.S. Geological Survey maps and -6 seconds for those plotted on State Highway maps. From a total of 1003 gravity stations, 210 are located at U.S. Geological Survey or U.S0 Coast and Geodetic Survey elevation "bench marks." A total of 700 stations are at road intersections whose eleva­ tions are given on U.S. Geological Survey maps. Elevations at road inter­ sections obtained from U.S. Geological Survey maps with a contour inter­ val of 20 feet have an accuracy of roughly i2 feet. A small number, approximately 10, of these stations were established at road intersections and the elevations were interpolated from the 20 feet contour interval. The error of these elevations does not exceed *10 feet. The elevations of 93 stations were established by using four Wallace & Tieman Surveying Altimeters, Type FA-112. two observers was employed. A single base method with Exceptionally clear days were chosen and all reasonable precautions were taken for the altimeter survey. The accuracy of these elevations is therefore, estimated to be ±5 feet and in no case exceed *8 feet. A test survey using the survey altimeters under the conditions described above showed a maximum error of i $ feet. Reduction of Data The observed readings were first corrected for the gravimeter drift relative to the bases which were tied by triple looping. calibration table for the La Coste and Romberg, No. 103 Geodetic \ The 31 Gravimeter was then used to convert the drift corrected readings from scale division to milligals and the values were adjusted to the national gravimetric datum. The theoretical gravity for each station was determined by the 1930 International Gravity Formula and the vertical gradient of gravity of 0.09ii06 milligals per foot was used to calculate the Free-air correction. The Bouguer or mass correction which accounts for the mass of material between the point of observation and sea level is 0.01276p per foot where p milligals is the density of the material above sea level which was assumed to be 2.67 g/cc0 Tfie above corrections were applied to the observed gravity to ob­ tain the Free-air and simple Bouguer anomalies in the equations: Free-air anomaly = observed gravity - sea level gravity + Freeair correction and Simple Bouguer anomaly = observed gravity - sea level gravity + Free-air correction - Bouguer correction© No terrain corrections were needed since the relief in the study area was low and the stations were carefully placed away from any local topographic features. A CDC 3600 Computer was utilized in the calculation of anomalies. Contouring of maps was done both by hand and by machine. The principal facts of gravity stations in the eastern portion of the Northern Peninsula of Michigan are available from the U CS0 Army, Topographical Command in Washington, D0C0 Lake Michigan and Lake Huron Gravity Data Bouguer gravity anomaly values for stations located in the northern Lake Michigan and Lake Huron were obtained from N0 W0 O'Hara of the University of Michigan, Great Lakes research Institute. These values were provided to Dr. O'Hara by the U.S. Topographical Command. CHAPTER III INTERPRETATION OF GRAVITY DATA Sources of Anomalies Geological and geophysical studies indicate that the Precambrian rocks of the upper Midwest are structurally complex and are composed of a broad spectrum of lithologies. Gravity and magnetic maps are particularly useful in deciphering this complex pattern where the Precambrian surface is buried beneath Paleozoic sediments because of the paucity of direct geolo­ gical information from deep drill holes. Specific gravity and magnetic susceptibility contrasts leading to gravity and magnetic anomalies are com­ monly associated with basement lithologic variations. Furthermore, Hinze and Merritt (1969) have shown that gravity and magnetic anomalies originating from intra-basement lithologic and structural variations are an order of magnitude greater than anomalies originating from basement topography, sedimentary structure and bedrock topography. Although there is an ob­ vious relationship between major magnetic and gravity anomalies and base­ ment lithology, assigning specific lithologies and configurations to sources of many anomalies is uncertain without supporting geological information. In general, specific gravity and magnetic susceptibility contrasts are un­ certain and little is known about the configuration of the causative bodies except that they reside within the basement rocks. Even if the physical properties are known the configuration of the anomalous mass cannot be de­ termined with certainty because of the fundamental ambiguity of potential 32 33 field data. Conversely, if the source can be delineated by geologic know­ ledge and the characteristics of the anomaly, the calculated physical pro­ perty does not lead to an unique determination of lithology of the source because the specific gravities and magnetic susceptibilities of rocks gen­ erally overlap (Grant and West, 1965)o The problem is further complicated because remanent magnetic polarization may contribute significantly to the magnetic polarization of the rock0 As a result, rocks which should cause magnetic highs on the basis of their susceptibility may cause magnetic lows or the opposite may occur. In spite of the difficulties of assigning lighologies to intra-basement features in the eastern portion of the Northern Peninsula of Michi­ gan, certain generalizations are possible particularly when both gravity and magnetic information are available. These generalizations are the result of geological and geophysical studies conducted over the vast ex­ posures of igneous and metamorphic rocks of the Canadian Shield, Thiel (1956), Craddock, Thiel and Gross (1963), and Sims and Zietz (1967) have found that large positive gravity and magnetic anomalies are associated with great thicknesses of Middle Keweenawan basic extrusives and intrusives in the western part of the Lake Superior syncline. This large pos­ itive anomaly which extends as far south as Kansas was designated the "Mid-Continent gravity high" by Thiel, Thiel showed that the lavas in the center of the syncline have been thrust upward as a horst, in places, juxtaposed against low specific gravity clastic sediments of Upper Kewee­ nawan age. These sediments give rise to strongly negative gravity and mag­ netic anomalies. However, a portion of the minimums flanking the central high is attributed to regional downwarping of the crust-mantle boundary cau­ sed by the great weight of the overlying volcanics (White, 1966bj Cohen and 3k Meyer, 1966). Keweenawan basalt flows are also responsible for the gra­ vity and magnetic highs in the western portion of the Northern Penin­ sula of Michigan (Bacon, 1966j Meshref and Hinze, 1970). Locally, in the Lake Superior area, both Keweenawan mafic intrusives and extrusives may produce negative magnetic anomalies due to remanent magnetization (Case and Gair, 1965; Corbett, et al., 1967; Meshref and Hinze, 1970). Granite intrusives generally cause negative gravity ano­ malies (Vfeaver, 1967; Gibb and McConnell, 1969) and either negative or positive magnetic anomalies depending on the nature of the country rock (MacLaren and Charbonneau, 1968). Dutton and Bradley's (1970) compilation of the Precambrian geology, and geophysical anomalies in Wisconsin indicate that belts containing metasedimentary rocks, mainly quartzites, are regional gravity and mag­ netic lows and massive granites are related to magnetic highs. Some of these belts can be traced eastward into Lake Michigan where magnetic lows are associated with metasediments and low grade felsic rocks (O'Hara and Hinze, 1971). O'Hara and Hinze have also traced some magnetic highs from Wisconsin into Lake Michigan and they relate these magnetic highs to basic volcanic rocks and also to intermediate to basic intru­ sions which have penetrated a primarily granitic terrane. In summary, the following generalizations will be applied in as­ signing sources to relative gravity and magnetic anomalies in the east­ ern portion of the Northern Peninsula of Michigan and the surrounding area: 1. Positive Bouguer gravity anomalies are caused by accumulations of mafic extrusives and intrusives of Keweenawan age; 2. Negative Bouguer gravity anomalies are commonly associated with metasediments and granitic rocks; \ 35 3. Magnetic highs are often accompanied by Keweenawan basic ex­ trusives and intrusives, iron formations, and with some preKeweenawan granitic rocks; U. Magnetic lows are frequently associated with metasediments and low grade gneisses; 5. Areas showing both marked positive gravity and magnetic ano­ malies are associated with accumulations of Keweenawan basic rocks which lie close to the surface; 6. Areas showing both negative gravity and magnetic anomalies are interpreted as indicating either the absence of Keweenawan ba­ sic rocks and/or the presence of thick accumulations of Kewee­ nawan sediments or the presence of relatively non-magnetic and low density masses of pre-Keweenawan rocks which lie near the surface; 7. Areas showing positive gravity anomalies and negative or no magnetic anomaly are associated with deep intrusions into the pre-Keweenawan complex or large remanent magnetization effects of the Keweenawan volcanics which oppose the induced magnetiza­ tion; 8. Areas showing negative gravity but positive magnetic anomalies are interpreted as indicating the presence of relatively thin near surface sequence of Keweenawan basic extrusives associated with a thick sequence of Keweenawan sediments or pre-Keweenawan rocks. In conclusion, it is possible to assign a lithologic type to a gra­ vity or magnetic anomaly only in the most general sense because of the overlapping physical property parameters among rock types and the ambiguity 36 of potential field methods. In addition, the regional relationship must be used cautiously because the source of anomalies, particularly the gra­ vity anomalies, may be deep seated and, therefore, not truly reflect the basement surface lithology. General Characteristics of Gravity Anomalies Discussion of Bouguer Anomaly Map Interpretation of gravity maps generally begins with a qualitative approach which consists of noting trends and inspecting gradients, ampli­ tude and shape of anomalies. The most important anomaly on the Bouguer Gravity Anomaly Map (Figure It) is a positive anomaly which extends from the Southern Peninsula of Michigan to the vicinity of Beaver Island where it bifucates, with one branch trending northwest and the other north. This anomaly is known as the "Mid-Michigan gravity high" in the Southern Peninsula and extends nearly the whole length of the Peninsula (Hinze, 1963). The width of the anomaly is very narrow in the vicinity of Beaver Island, but rather broad elsewhere. Its magnitude is of the order of 30 mgals in the Southern Peninsula and Northern Lake Michigan but decreases to 20 mgals in both northerly and northwesterly directions. The northwestern branch of the anomaly is strongly asymmetric, with steeper gradients on the southwestern side. Patenaude (l?6ij) and O'Hara (1967) have interpreted the steep gradients to indicate a northwestsoutheast trending fault contact between Keweenawan volcanics and the older Precambrian rocks on the southwest. This anomaly extends to the vicinity of Munising and on the basis of magnetic evidence, which will be discussed in the following section and is presented in Figure 8, represents the continuation of the Middle Keweenawan basic volcanics from the Keweenaw Peninsula (Hinze, et al., 1966). The northern branch BOUGUER GRAVITY ANOMALY MAP SCALE LEGEND 84* o Figure L.— Bouguer Gravity Anomaly Map of the eastern portion of the Northern Peninsula of Michigan. 38 of the anomaly extends across the eastern end of Lake Superior and ties * in with Keweenawan basalts outcropping on Mamainse Point. The anomaly also extends to the east with considerable amplitude, finally decaying in amplitude near the Saint Marys River. The eastern branch gravity an­ omaly is associated with a very strong magnetic anomaly and is believed to indicate the extension of the Keweenawan basalts in an easterly dir­ ection to Saint Marys River. These gravity anomalies correlate in general with positive magnetic anomalies (Figure $), Hinze, et al., (1966) and O'Hara and Hinze (1971) have interpreted the magnetic anomalies to represent the limbs of the Lake Superior basin extending into the Northern Peninsula of Michigan. They further suggest that the basin terminates in the vicinity of Beaver Island where the magnetic positives join to form a single magnetic high which becomes strongly negative in the Grand Traverse Bay area and con­ nects to the south with the "Mid-Michigan gravity and magnetic anomaly*" The intense magnetic minimum in the Grand Traverse Bay region is interpre­ ted to be the magnetic expression of reversely polarized magnetic basalts as observed elsewhere in the Lake Superior basin (Books, 1968). The negative gravity anomaly to the west of the northern branch of the major positive gravity anomaly in the vicinity of Whitefish Point exhibits a magnitude of the order of 1$ mgals. The steep gradients on the eastern side of the anomaly may be due to a fault down thrown to the west. The magnetic values over this gravity low are relatively high suggesting the presence of structurally disturbed volcanics in the area under an increased thickness of the Keweenawan sediments. The negative gravity anomaly in the Whitefish Bay area to the east of the northern branch of the positive gravity anomaly is similar to the RESIDUAL TOTAL MAGNETIC INTENSITY ANOMALY MAP LEGEND Figure 5.--Total Magnetic Intensity Anomaly Map of the eastern portion of the Northern Peninsula of Michigan., ho previously described gravity low and probably has a similar source, that is a fault downthrown to the east and thickening of the Keweenawan sed­ iments. The positive gravity anomaly south of the Mackinac Straits is char­ acterized by relatively steep gradients on both the eastern and the west­ ern sides and has a magnitude of 20 mgals. in a north-south direction. It is considerably elongated The magnetic values over this positive ano­ maly are relatively high, but there is no correlative magnetic high. In their magnetic study of Lake Huron, Secor, et al., (1967) have interpre­ ted this area to be associated with a pre-Keweenawan mafic igneous com­ plex extending south to Rogers City, Michigan. However, it is also pos­ sible that the positive gravity and magnetic anomalies are caused by an upfaulted block of Keweenawan volcanics. The area east of 81*° 30' W o n the Bouguer Gravity Anomaly Map (Fig­ ure it) is characterized by broad negative anomalies, with the exception of a positive gravity anomaly centered in the vicinity of Elliot Lake in Ontario, whose magnitude is of the order of 20 mgals • It is asso - ciated with an east-west trending positive gravity anomaly which extends to Sudbury, Ontario as shown on the Bouguer Gravity Anomaly Map of Canada. This positive belt of gravity anomalies is identified with Precambrian metasedimentary and metavolcanic rocks within a structurally complex belt of folds and faults. the region east of 8k° The Total Magnetic Intensity Anomaly Map of 301 W is rather complex and exhibits in general, east-southeast trending anomalies with steep gradients and amplitudes ranging up to 1000 gammas or more. Hinze and Merritt (1969) have ass­ ociated this complex pattern of magnetic anomalies in the eastern part of the northern tip of Southern Peninsula with metavolcanics, metasediments Ui and gneisses of pre-Keweenawan age. Further north, in the vicinity of North Channel, the major east-west trending negative magnetic anomaly is interpreted to indicate a predominantly granitic province (Secor, et al., 1967). The western part of the Bouguer Gravity Anomaly Map is dominated by a very broad negative gravity anomaly which extends southward into Lake Michigan. There are a number of smaller anomalies within this neg­ ative anomaly which are associated with synclinal structures containing Animikean sediments and iron formations. In the vicinity of Escanaba, the continuation of Msnominee iron range is clearly shown on both gravity and magnetic maps. Another gravity and magnetic high in the vicinity of Trenary, Michigan is associated with an east-west trending syncline fil­ led with Animikean sediments (Frantti, 19!?6). Still another local pos­ itive gravity and magnetic anomaly occurs over Paleozoic sediments near Cooks, Michigan which may be associated with an east-west striking Precambrian iron formation (Hinze, et al., 1966). The amplitude of the broad negative gravity anomaly which encloses all the local anomalies, decreases in a southerly direction. In Green Bay and the western part of northern Lake Michigan where there is no gravimetric coverage, an intense regional magnetic minimum striking east-west from the western shore of Green Bay to the south of Beaver Island characterizes the south­ western portion of the study area (O'Hara and Hinze, 1971). This magnetic minimum is on strike with felsic rocks of the Mountain-Amberg area of Wisconsin which is also characterized by magnetic and gravity minimums. In summary, the Bouguer Gravity Anomaly Map of the eastern portion of the Northern Peninsula of Michigan and the surrounding area is dominated U2 by a major positive gravity anomaly of considerable amplitude which is attributed to Keweenawan basalts and the local gravity minimums in the Vlhitefish Bay area, which are associated with thickening of the Kewee­ nawan clastic sediments. This major anomaly transects roughly east-west trending gravity anomalies whose sources are generally believed to be pre-Keweenawan rocks. Although, the peaks of magnetic and gravity ano­ malies do not always correspond, their signs are, in general, the same. The main exceptions are the negative magnetic anomaly in the Grand Traverse Bay region and the relatively positive magnetic anomaly in the vicinity of Vlhitefish Point. Discussion of the Double Fourier Series Residual Gravity Map A Double Fourier Series Residual Gravity Map of the eastern portion of the Northern Peninsula of Michigan and surrounding area was prepared in an attempt to isolate anomalies which are less obvious in the Bouguer Gravity Map. The double Fourier series has recently come into use in geology and geophysics as an alternative approach to the least-squares polynomial method for isolating anomalies. A method of double Fourier series for surface fitting of irregularly spaced data has been described by James (1966). This method was used to calculate the Double Fourier Series Residual Gravity Map as shown in Figure 6. The residual values were ob­ tained by removing from the Bouguer Gravity Map fundamental wavelengths of 1*75 miles in the east-west direction and 575 miles in the north-south direction. Thus, all wavelengths in excess of 95 miles in the east-west direction and 113 miles in the north-south direction have been removed in the Double Fourier Series Residual Map. mental wavelengths was The selection of the funda­ based on the observation of the half-wavelengths DOUBLE FOURIER SERIES RESIDUAL GRAVITY ANOMALY MAP u c\*i 9 Figure 6.--Double Fourier Series Residual Gravity Anomaly Map of the Northern Peninsula of Michigan,, hh of the major gravity anomalies of the Bouguer Anomaly Map. Anomalies trending significantly different than either east-west or north-south may be considerably distorted as a result of the directions of the fundamen­ tal wavelengths (Whitten, 1969). In addition, the removal of the long wavelength components from predominantly short wavelength anomalies also causes distortion of these anomalies. The predominant anomalies of the Double Fourier Series Residual Map (Figure 6) correlate with the previously discussed major positive gravity and magnetic anomalies on Figures U and A major anomaly extends north from the Southern Peninsula of Michigan and bifucates northeast of Beaver Island, with one branch trending north through Vlhitefish Bay and the other trending northwest through Grand Island. The north trending ano­ maly with large amplitude and steep gradients is much more obvious on Figure 6 than the corresponding anomaly shown on the Bouguer Gravity Anomaly Map of Figure li. subdued in amplitude. However, the northwest trending branch is rather The eastern extension of the major gravity anomaly has been isolated as a positive gravity anomaly northeast of the Mackinac Straits, which in general correlates with the positive magnetic anomaly in the Pickford area. However, its position is slightly to the south of the magnetic anomaly. The well defined gravity lows on either side of the north branch of the major gravity anomaly in the vicinity of Whitefish Bay on the Bouguer Gravity Anomaly Map are shown to have more regional extent in Figure 6 as a result of the isolation of the positive anomalies. Also a positive anomaly has been isolated between Grand Island and Whitefish Bay, The positive gravity anomaly south of the Mackinac Straits is shown to be isolated in the Double Fourier Series Residual Gravity Map, but still has a northwest-southeast trend. \ U5' The east-west trend of the pre-Keweenawan structures in the western part of the area is shown well in Figure 6. Most important, however, is the separation of the positive gravity anomaly south of Marquette from the northwest branch of the major gravity anomaly0 This anomaly has a distinct east-west trend and,therefore, is associated with the preKeweenawan rocks« Figure 7 shows the Bouguer, double Fourier residual and the regional gravity curves along profile E-E'-E" at the eastern portion of the Nor­ thern Peninsula. The regional gravity values were obtained by subtracting the contoured double Fourier series residual values from the Bouguer gravity values. The purpose of Figure 7 is to illustrate the character­ istics of the regional values which have been removed from the map of Figure 6. The regional profile shows long wavelength and relatively large amplitude anomalies which result partly from the longer wavelength components of the volcanics and partly from deep crustal specific gravity variations, perhaps associated with intrusions related to Keweenawan igneous activity. In conclusion, the Double Fourier Series Residual Gravity Map has added information to the anomalous trends of the Bouguer Gravity Anomaly Map. Comparison of Profiles With the Aid of Upward Continuation In this study, five gravity and magnetic profiles were constructed from the Bouguer Gravity Anomaly Map (Figure Magnetic Anomaly Map (Figure $) h) and the Total Intensity of the eastern portion of the Northern Peninsula of Michigan and the surrounding area, the Total Magnetic Inten­ sity Anomaly Map of the Eastern Lake Superior (Hinze, et al., 1966) and the \ \ -10 BOUGUER GRAVITY REGIONAL GRAVITY -20 (a) -4 0 20 w 10 - < 0 s (b) -10 HORIZONTAL SCALE I 10 1 MILES Figure 7®— Profile of a) the Bouguer and regional gravity curves, and b) residual values for double Fourier series fit® (Location of this profile is shown in Figure 8) 1*7 Bouguer Gravity Anomaly Map of the Lake Superior Region (Weber and Goodacre, 1966). The selected profiles are presented in Figure 8. The pro­ files are approximately perpendicular to the major anomalies and have been plotted ‘from a common datum. In addition, all profiles have been analytically upward continued using a 17 point upward continuation for­ mula developed by Hinze (i960), to make the level of gravity observa­ tions for all profiles approximate^ an equal distance from the basement. Upward continuation has made little change in these because of the broad wavelengths of the anomalies. The magnetic anomalies on the southwestern side of the profiles show good continuity, therefore, the correlation of gravity anomalies on the same side of the profiles will be discussed first. On Profile A-A1, the magnetic anomaly with a magnitude of about 1000 gammas originates from the Middle Keweenawan basic volcanics outcropping on the Keweenaw Peninsula. The Keweenaw fault which has a near vertical dip in this area extends along the southern margin of the volcanics. The volcanics are overlain to the north by clastic sediments and the Keweenaw fault brings them into contact with the Upper Keweenawan elastics to the south. This positive magnetic anomaly and the flanking lows associated with it are directly related to a gravity anomaly which becomes obscure to the north because of lack of gravity data. However, the positive gravity anomaly which is flanked by a minimum to the south is interpre­ ted to indicate the presence of thick volcanics in the area. The magnetic and gravity anomalies are continuous from profile A-A' to profile B-B'. However, on profile B-B', the peaks of the magnetic and gravity anomalies do not correspond. This could be attributed to the scarcity of gravity data in the area. However, it is also possible that the volcanics are near the surface under the magnetic high but are at a SELECTED GRAVITY AND MAGNETIC PROFILES OF THE EASTERN LAKE SUPERIOR REGION Os' Figure 80--Selected gravity and magnetic profiles of the eastern Lake Superior Region, h9 considerable depth and thicker under the gravity high. In addition to having offset peaks, the magnetic and gravity anomalies differ in ampli­ tude and gradients from those on profile A-A'. This can be explained by either differing thickness or dip of volcanics in the two areas. Profile C-C' exhibits much the same pattern as the previous two profiles. The magnetic anomaly on profile C-C’ corresponding to the Middle Keweenawan volcanics is broader and slightly lower in amplitude than the same magnetic anomaly on profile A-A1. Hinze, et al„, (1966) have interpreted a fault at the southwestern side of this anomaly and suggested that it is the continuation of the Keweenaw faulto The interpretation of magnetic and gravity highs on profile A-A1, B-B' and C-C' indicates continuation of the Middle Keweenawan volcanics from Keweenaw Point in a southeasterly direction. This is substantiated by volcanic rocks found on Stannard Rock, which occurs in Lake Superior roughly midway between Keweenaw Point and Grand Island. Further southeast from profile C-C1, profile E-E' is characterized by a broad gravity high which correlates well with the corresponding positive gravity anomaly on profile C-C’. It has approximately the same width, amplitude and the gradients, thus indicating the extension of the Middle Keweenawan volcanics into the eastern portion of the Northern Pen­ insula of Michigan. However, the amplitude of the associated magnetic anomaly is considerably less indicating that the volcanics in this area are relatively undisturbed having low dips and under a considerable thickness of sedimentary cover0 A strongly asymmetric positive magnetic anomaly with steeper grad­ ients on the northern side, marks the central part of profile A-A'. This anomaly is interpreted as the extension of the Middle Keweenawan basic 5o volcanics on Isle Royale with Isle Royale thrust fault marking the north­ ern edge of the anomaly (Hinze, et al., 1966). Further north along pro­ file A-A' another positive magnetic anomaly is also believed to be asso­ ciated with basic volcanics dipping toward the axis of the Lake Superior syncline. This positive magnetic anomaly is directly related to a broad positive gravity anomaly. The central part of profile C-C' is associated with a gravity mini­ mum and a magnetic high. erable amplitude. Both these anomalies are broad and of consid­ The broad positive magnetic anomaly is believed to result from faulted and relatively thin near surface volcanics. The pos­ itive gravity and magnetic anomalies on the northeastern end of profile C-C1 are correlated with basic Keweenawan volcanic rocks of Michipicoten Island. The complex nature of the magnetic anomalies in this area is believed to reflect the existence of two roughly northwest-southeast trending faults (Hinze, et al., 1966). The positive gravity anomaly situated in Whitefish Bay, on profile E'-E" is very similar in its characteristics to the gravity anomaly west of Michipicoten Island. This anomaly may represent the continuation of the Keweenawan volcanics of Mamainse Point which Green (1971) has placed in both the Lower and Middle Keweenawan or it may be the extension of the Middle Keweenawan volcanics of Michipicoten Island. If future gravity studies in eastern Lake Superior confirm the latter suggestion, then perhaps the Michipicoten thrust fault extends along the eastern edge of the positive gravity anomaly. The gravity minimums on either side of the positive anomaly suggest a thickening of the Keweenawan sed­ iments over the volcanics in the area. The broad positive magnetic ano­ maly which is associated with the gravity minimum to the northwest of 51 the gravity high is interpreted to result from faulted and steeply dipping volcanics. The southeastern end of profile E'-E" exhibits a positive magnetic anomaly whose amplitude reaches 1000 gammas. The gravity ano­ maly which is associated with this magnetic high is very broad and has a relatively small amplitude. These magnetic and gravity anomalies are interpreted to indicate a thinning of the Keweenawan volcanics to the east with near surface extrusives or intrusives causing the high ampli­ tude of the magnetic anomaly. The positive gravity anomalies on profiles E-E' and E'-E" join to­ gether to form a single gravity high on profile G-G'. The peak of this anomaly is inversely related to a negative magnetic anomaly which has been interpreted to be the magnetic expression of reversely polarized magnetic basalts (O'Hara and Hinze, 1971). Profile G-G' was further com­ pared with an upward continued typical Bouguer gravity anomaly profile of the "Mid-Continent gravity higho" This comparison showed that the gravity anomalies differ in amplitude and gradients but the width of the anomalies are similar. A possible explanation for this is that the vol­ canics are thicker and the structure is very complex in areas associated with the "Mid-Continent gravity high." In summary, the comparison of gravity and magnetic profiles strongly indicate that the Middle Keweenawan volcanics of the Keweenaw Peninsula curve southward through Stannard Rock and extend into the eastern portion of the Northern Peninsula of Michigan, and join the Middle Keweenawan volcanics extending south from Mamainse Point, Ontario, in the vicinity of Beaver Island to become associated with the "Mid-Michigan gravity and magnetic anomaly" of the Southern Peninsula of Michigan. 52 Model Studies The purpose of the gravity model studies in this section is to quan­ tify the qualitative interpretation of the previous sections. The con­ figuration and horizontal and vertical extent of principally the Kewee­ nawan rocks, both the clastic sediments and igneous rocks, are studied by the method of model studies and matching observed and calculated gravity curves. Magnetic effects from gravitationaly constructed models were also calculated and compared with the observed magnetic curves. Four profiles of observed Bouguer gravity and total intensity magnetic anomalies were selected for quantitative investigation. The location of these profiles is shown in Figure 8. Gravity and magnetic effects of basement structures were computed using standard two dimensional theory developed by Talwani, et al., (1959). The assumption of two dimensionality is justified because most anomalies are horizontally linear having greater lengths than widths. However, the results of model studies are not unique because of uncertainties in the physical properties of the rock units. In the models, high density basalts and low density sedimentary rocks overlie the pre-Keweenawan rocks. In addition, the observed gravity and magnetic curves can be matched in more than one way by varying the thicknesses of the rock units. Therefore, in constructing the models, geologically reasonable structures based on studies elsewhere in the Lake Superior region were considered and no attempt was made to precisely match calculated and observed anomaly profiles. Furthermore, only the upper basement structures were considered in the models. Cohen and Meyer (1966) in their study of the "Mid-Continent gravity high" in northeastern Minnesota and northwestern Wisconsin 53 combined the effect of a basalt trough with the effect from a downwarped Mohoc Their analysis is based on loading of the crust and its subsequent deformation by a two dimensional body0 The crust behaves like an elast­ ic beam, while conserving isostatic balance0 Kellogg (1971) has studied a similar model across the "Mid-Michigan gravity high" in the eastern part of the Southern Peninsula of Michigan.) He reports that the gravity gra­ dient attributed to the downwarped Mbho is only 0o2 mgal per mile0 There­ fore, the downwarped Moho can account for only a small portion of the negative anomaly associated with the "Mid-MLchigan high," the remainder coming from specific gravity contrasts within the upper crust0 Specific gravities and magnetic properties used in the models are based on published measurements. Table 2 presents a brief summary of the specific gravities of rocks of the Lake Superior area from previous stud­ ies. Magnetic susceptibilities of pre-Keweenawan rocks and other mag­ netic rocks, with the exception of Keweenawan basalts, were modified to match the observed magnetic profiles. Remanent as well as induced mag­ netic polarization was used in the calculation of the magnetic effect of the Keweenawan basalts. Approximate values given by Du Bois (1962) for the remanent magnetic polarization component of basic Keweenawan volcanics were vectorially combined with the induced magnetization component, to calculate the total magnetic polarization vector. The induced susceptibi­ lity value in Table 3 which summarizes the results,was also taken from Du Bois (1962)o Calculation of the combined vector was made assuming no rotation of the remanent magnetization component of the Keweenawan lavas due to deformation. It is realized that this assumption has introduced error in the calculated anomalies because the basalts have been deformed. Furthermore, negative magnetic polarization of basalts as observed at Mamainse Point can cause error in the calculated profiles. TABLE 2 Specific gravities of Precambrian rocks in the Lake Superior Region Source Jacobsville FormationBayfield Group Pronto Group Middle Keweenawan Volcanics - 2.90 - 0.10 Pre-Keweenawan Rocks_____ Theil, 1956 2.30 - 0.06 White, 1966 2.30 2.37 2.UU 2.62 2.90 2.92 2.67 Bacon, 1966 2.25 2.66 (Freda) 2.88 2.70 2.97 2.70 Weber and Goodacre, 1966 2.36 6.12 average of both groups 2.50 Steinhart and Smith, 1966 Steinhart et al., 1968 2.30-2.36 2.30 2 .U 3 - 2 .5 U 2.66 2. 70- 2.80 55 TABLE 3... Parameters of induced, remanent and combined polarization components of the Keweenawan basalts. Component Magnetic Polarization Inclination Declination Induced 0.00091 emu/cc +71° 0° Remanent 0.0035>U emu/cc +i*5° 285° Combined O.OOU25 emu/cc +53° 292° Interpreting gravity anomalies by the method of matching profiles is ambiguous. Therefore, the thicknesses used in the models must be geo­ logically reasonable. The guide for this is published thickness estimates of geological units which are summarized in Table ii. Profile E-E1 Profile E-E1 is oriented southwest-northeast in the central part of the eastern portion of the Northern Peninsula. The positive gravity ano­ maly on this profile is attributed to the Middle Keweeenawan basalts ex­ tending from the Keweenaw Point. was taken to be 2.95* direction. The specific gravity of these basalts The basalts gradually thin out in a southwesterly A reverse fault with limited vertical displacement was in­ cluded in the model on the basis of previous geophysical interpretation. However, the gravity curve could be equally well matched without the fault (Figure 9a). The thicknesses used in this model are well within the lim­ itations given in Table The magnetic profile, which was calculated from the gravitationally constructed model, agrees in general with the observed magnetic profile. There is no well data available in the vicinity of this profile to check \ TABLE H Thicknesses (km) of the Keweenawan rock units in the Lake Superior Region Source Jacobsville FormationBayfield Group Tyler et al., 19UP 1.3 Bacon, 1966 3oP Pronto Group Middle Keweenawan Volcanics U. 2 Southwestern Lake Superior Keweenaw Fault White, 1966b 6.1-9.1 Halls, 1966 7c6 bo6 Berkson, 1969 (after Pcola, 1968) Berkson, 1969 (After Anzoleaga, 1968) Location Lake Superior Basin North Shore Volcanics Mamainse Point 1.0-1.5 l.P-1.5 a»P-5oP Eastern Lake Superior IcO 1.5-3.8 5=8 Western Lake Superior — OBSERVED MAGNETICS — CALCULATED MAGNETICS 200 r 0 .0 0 3 «m u/cc □ PRE-KEWEENAWAN BASEMENT COMPLEX S . 6 . : 2 . 7 0 , K s 0 .0 0 0 5 «mg/ee -4 0 -SO 20 30 40 50 60 90 SL 140 SL 82 o > cc. O o VO LEGEND ill(( HI EDGE OF >000.. STRUCTURAL CONTOURS AT T>€ TOP OF K E W E E N M N VQLCAMCS AS OBTAINED FROM OUMTY MOOEL STUOCS FAULT CONTROL LME 83* Figure 15.--Schematic map showing interpreted Precambrian structures in the eastern portion of the Northern Peninsula and the northern tip of the Southern Peninsula of Michigan* 70 part of the study area is magnetically complex. Along the eastern shore of Lake Superior in Canada, known geology (Halls, 1966) was used to esta­ blish the edge of the Keweenawan volcanics. To the south, the Keweenawan volcanics become areally restricted and are gravitationally expressed as the "Mid-Michigan gravity high" in the Grand Traverse Bay area of the Southern Peninsula of Michigan. The structure contours on top of the Middle Keweeenawan basalts were interpolated and extrapolated primarily from the model studies. The depths shown in Figure 13> are the minimal depths to the top of the volcanics because Freda sandstone was not considered in the geological cross-sections that incorporate only the Jacobsville sandstone. The structure contours indicate that the limbs of the Lake Superior syncline, which have been traced into the eastern portion of the Northern Peninsula on the basis of gravity and magnetic anomalies as discussed in the pre­ vious sections, merge in the vicinity of Beaver Island. South of the latitude of Beaver Island, the subcrop width of the volcanics gradually narrows and is faulted into a series of horsts and grabens which then merge with the source of the "Mid-Michigan gravity high,," Due to the paucity of gravity data in northern Lake Michigan, no structure contours could be established between Beaver Island and the Lake Michigan shore of the Northern Peninsula. However, on the basis of available control, the structure contours outline two local basins in the vicinity of Whitefish Bay, which are associated with relatively large negative gravity anomalies. The presence of the fault running roughly parallel to the western edge of the Keweenawan volcanics is questionable. However, it may be present with only limited vertical displacement at depth and only minor 71 topographic expression as a fault scarp at the basement surface. In fact, Patenaude (196U) interpreted the gently sloping contact suggested by the gravity data to indicate an eroded fault scarp. The northern tip of the Southern Peninsula is interpreted to be highly faulted and is characterized by a series of horsts. A horst flanked by a graben to the east occurs in the Grand Traverse Bay area0 The vertical displacement on the boundary faults of this graben is of the order of 1*500 feet. Another horst which trends north-south occurs in the Whitefish Bay area. Vertical displacement of the bounding faults of this structure is of the order of 1500 feet0 The relationship of the Whitefish Bay horst to the interpreted faults in the northern tip of the Southern Peninsula is unknown, but this horst is on strike with the in­ terpreted horst associated with the "Mid-Michigan gravity high" east of Beaver Island. The eastern end of the Lake Superior syncline which is schematically illustrated by the structure contours in Figure 15, closely resembles the western end of the syncline,. In the western end of the Lake Superior syncline, positive gravity anomalies correlate with the Keweenawan extrusives and intrusives, and the large negative anomalies on the Bayfield Peninsula and Keweenaw Bay reflect thick accumulation of sediments (Thiel, 1956; Weber and Goodacre, 1966). Furthermore, the St. Croix horst of Wisconsin which is bounded by the Douglas-Isle Royale and Keweenaw faults continues into the Lake Superior syncline (Halls and West, 1971). The Keweenawan structural province of Figure 15 constitutes a key segment of a major geologic feature which extends from southern Kansas along the "Mid-Continent gravity high," through the Lake Superior trough, and southeastward into the Southern Peninsula of Michigan. \ The occurrence 72 of an inverse relationship between the Bouguer gravity anomalies and basement elevations along the "Mid-Continent" and "Mid-Michigan" gravity highs has led Lyons (1959) and Hinze (1963) to speculate that the geolo­ gic features responsible for the geophysical anomalies originate from elastic deformation in response to the added mass of basic rocks, emplaced in the basement complex in late Precambrian time. McGinnis (1970) ex­ panding on these ideas, has suggested that the Lake Superior basin re­ sulted from elastic deformation due to basic intrusives and extrusives emplaced along a Keweenawan rift zone. Hinze, et al. (1971) have pro­ posed an alternative theory for the origin of the Lake Superior basin. They suggest that the crust undergoes compression at the final stage of the continental rifting process resulting in the slow development of a broad basin over a paleo-rift zone. In the case of the Lake Superior basin, widespread basalts were ex­ truded within and probably exterior to the rift zone at the earlier stages of rifting. Evidence for the presence of high density intrusions associated with these stages is provided by high compressional wave velo­ cities under the Lake Superior syncline (Smith, et al., 1966). This excess mass of mafic intrusions and extrusions resulted in elastic de­ formation of the crust and development of local basins. This stage cor­ responds to the deposition of the Oronto group in the Lake Superior syncline. During this stage, the Freda sandstone may have been deposited in the eastern end of the syncline in shallow water and exposed repeatedly to subaerial conditions as indicated by sediment dispersal studies (Hamblin, 196l). The direction of the currents which deposited the Freda was northwest in the Keweenaw Peninsula but north, south and west in the area of Mamainse Point. Keweenawan basalts, Animikean iron formations 73 and granite were the dominant rock types in the source area. The deposi­ tion of Freda sandstone was followed by compression of the rocks within the rift. time. Perhaps some of the observed thrust faulting developed at this The development of the Lake Superior basin was largely completed with the depostion of the Jacobsville sandstone. Hinze, et al., (1971) have associated the development of the Jacobsville with subsidence of a former "floater" which cooled rather rapidly due to voluminous igneous activity, leaving behind an abnormally denser residium. The "floater" is defined as a portion of mantle whose specific gravity is less than 1 the surrounding mantle material because of its relatively high temperature. Contemporaneous subsidence of the surface with that of the "floater" led to the deposition of the Jacobsville sandstone which thins out to the south because the outer edge of the basin subsided less than the central part. The high quartz content, good sorting and simple heavy mineral suite indicate that Jacobsville sandstones were subjected to considerably more transportation than the Freda sandstones. The direction of sediment transport during Jacobsville time was to the northwest near the eastern margins of Lake Superior but due north in an area extending from Grand Marais to Marquette on the south shore of Lake Superior. Since the development of the Lake Superior basin at this time was nearly complete, the subsidence of the basin did not keep pace with sedimentation and, therefore, Jacobsville sedimentation spread southward. Following the Jacobsville sedimentation, early Paleozoic seas entered the eastern portion of the Northern Peninsula and led to the deposition of the Chapel Rock sandstone of St. Croixan age. A barrier to the Cam­ brian seas from the south formed in the northern portion of the Southern Peninsula and the eastern portion of the Northern Peninsula following 7b the termination of Jacobsville sedimentation (Hamblin, 196l). Following this period, the transgressive Cambrian seas covered the whole eastern portion of the Northern Peninsula and the area became the northern most portion of the Michigan basin. CHAPTER IV CONCLUSIONS The gravitational study of the eastern portion of the Northern Pen­ insula of Michigan, integrated with the available magnetic data and geological information has confirmed some previous geological concepts and brought to light new ideas on the tectonics of the eastern end of the Lake Superior syncline. However, the geophysical interpretation of the gravity and magnetic data is subject to error because of inherent ambiguity and assumptions made in quantitative analysis. Therefore, the conclusions presented here serve as a first order approximation. The Bouguer Gravity Anomaly Map of the eastern portion of the North­ ern Peninsula of Michigan is dominated by a major positive gravity anomaly. The source of this anomaly is attributed to Keweenawan basalts. The local gravity minimums in the Whitefish Bay area result from a thick accumulation of Upper Keweenawan clastic sediments. The major gravity anomaly transects roughly east-west trending gravity anomalies whose sources are pre-Keweenawan rocks which are exposed around the periphery of the Lake Superior syncline. The Bouguer and Double Fourier Series Residual Gravity Anomaly Maps indicate that the major positive anomaly in the eastern portion of the Northern Peninsula is comprised of essentially two major linear positive gravity anomalies and a more subdued positive anomaly along the eastern margin of the Northern Peninsula. One of the major linear anomalies IS \ r. . 76 trends south from Miitefish Point on the south shore of Lake Superior and the other trends southeast from Grand Island in Lake Superior. Posi­ tive magnetic anomalies are also found to be associated with these gravity anomalies. Comparison of major anomalies both magnetically and gravimetrically with anomalies of the Lake Superior region suggests that the northwestsoutheast trending anomaly can be correlated to the Middle Keweenawan volcanics of the Keweenaw Peninsula. This anomaly represents the edge of the western limb of the Lake Superior syncline. The north-south tren­ ding anomaly which is associated with upfaulted basalts within the east­ ern limb of the Lake Superior syncline can be correlated with the Middle Keweenawan volcanics outcropping on Mamainse Point, Ontario. The eastern positive gravity anomaly defines the edge of the eastern limb of the syncline. These anomalies merge in the vicinity of Beaver Island and mark the termination of the Lake Superior syncline. South of Beaver Island, the Keweenawan basalts merge with the source of the "Mid-Michigan gravity high." The interpretive results obtained from two dimensional model studies suggest that the Lake Superior syncline in the eastern portion of the Northern Peninsula consists of up to 12,000 feet of basaltic flows overlain by Upper Keweenawan clastic rocks. The basalts thin both to the east and west and, in general, have low dips except in the vicinity of faults. Low specific gravity elastics of Jacobsville age, which thin to the south overlie the basalts. The Freda sandstone which occupies a stratigraphic position between the basalts and the Jacobsville sandstones was also included in some models. i 77 The presence of a fault extending from Au Train Bay to the vicinity of Manistique as suggested by previous geophysical studies is not sup­ ported by the results of this investigation, but this does not preclude the existence of the fault. An approximately north-south trending horst is interpreted to occupy the Miitefish Bay area. The vertical displace­ ment of the bounding faults of this structure is of the order of 1^00 feet. The basalts are also highly faulted in the northern tip of the Southern Peninsula. A horst which is flanked by a graben to the east occurs in the Grand Traverse Bay area. The vertical displacement on the boundary faults of the graben is of the order of h%00 feet. An alternative interpretation of the geology of the northern portion of the Southern Peninsula of Michigan suggests that basalts may be con­ fined to the Grand Traverse Bay area. In this case pre-Keweenawan ex- trusives and intrusives underlie the area of the northern tip of the Southern Peninsula. The structure of the eastern end of the Lake Superior syncline re­ sembles that of the western end. It constitutes a key segment of a major geologic feature which extends from southern Kansas along the "Mid-Continent gravity high," through the Lake Superior trough, and southeastward into the Southern Peninsula of Michigan. The origin of this major geologic feature is thought to be in a rift of the continental crust. v BIBLIOGRAPHY 78 \ BIBLIOGRAPHY Allen, Ro C. 1911*. Relative to an Extension of the Menominee Iron Range Eastward from Waucedah to Escanaba, Michigan* Econ, Geol., 9 pp. 236-2380 Bacon, L« 0o 1957. Relationship of Gravity to Geological Structure in Michigan's Upper Peninsula. Inst, on Lake Superior Geology, Houghton, Michigan, pp. 51*-58. Bacon, L. 0o 1961*. Investigation of the Thickness of the Jacobsville Sandstone by Seismic Reflection Methods-a Progress Report (abstract), Inst, on Lake Superior Geology, Ishpeming, Michigan, p 0 57. Bacon, L. 0. 1966. Geologic Structure East and South of the Keweenaw Fault on thebasis of Geophysical Evidence. Amer. Geophy. Un», Geophy. Mon. 10, pp. 1*2-55* Bass, Mo No 1971o Personal Communication. Bayley, R. W. and Muehlberger, W. R. 1968, United States. U 0 S. Geol. Survey. Basement Rock Map of the Behrendt, J 0 C. and Woollard, G. P. 1961. An Evaluation of the Gravity Control Network in North America. Geophysics, 26, pp, 57-76. Berkson, J. M. 1969. A Gravity Survey in the Vicinity of Michipicoten Island, Lake Superior. M.S. Thesis, University of Wisconsin, 57 pp. Berry, M 0 J. and West, G. F. 1966. An Interpretation of the First Arrival Data of the Lake Superior Experiment by the Time-Term Method. Bullo Seismolo Soc. Amer., 56, pp. 11*1-171. Books, K. G. 1968. Magnetization of the Lowermost Keweenawan Lava Flows in the Lake Superior Area, in Geological Survey Research. U. S. Geol. Surv. Prof. Paper 600-D, pp. D2l*8-251*. Butler, B. S. and U. S. Geol. Burbank, W. S. 1929.The Copper Surv. Prof. Paper 11*1*, 283 pp. Depositsof Michigan. Case,J. E. and Gair, J. E. 1965* Aeromagnetic Map of Parts of Marquette, Dickinson, Baraga, Alger and Schoolcraft Counties, Michigan and its Geologic Interpretation. U. S. Geol. Surv. Geophy. Inv. Map GP-U67. 79 \ 80 Cohen, T. J. and Meyer, R. P. 1966. The Midcontinent Gravity High: Gross Crustal Structure. Amer. Geophy. Un., Geophy. Mon. 10, pp. 11*1-165. Corbett, J. D,, Hinze, W. J., and Secor, G. B 0 1967c A Regional Geophysical Study of the Port Coldwell Complex, Ontario (abstract). Inst, on Lake Superior Geology, East Lansing, Michigan, p. 80 Craddock, C0, Thiel, E 0 C., and Gross, B. 1963. A Gravity Investijgation of the Precambrian of Southeastern Minnesota and Western Wisconsin. Jour. Geophy Res., 68, pp. 6015-6032. Du Bois, Po M. 1962. Paleomagnetism and Correlation of Keweenawan Rocks. Bull. Geol. Surv. Canada, 71, 75 PP« Dutton, C. E. and Bradley, R. E. 1970. Lithologic, Geophysical and Mineral Commodity Maps of Precambrian Rocks in Wisconsin. U. S. Geological Service and University of Wisconsin, Map 1-631. Frantti, G. E. 1956. Geophysical Investigations in the Central Portion of Michigan's Upper Peninsula. Mining Engineer, 8, pp. 70-72. Gibb, R. A 0 and McConnell, R. K. 1969. The Gravity Anomaly Field in Northern Manitoba and Northeastern Saskatchewan with Maps. Publ. Dorn. Obs. of Canada, 28 pp. Girdler, R. W., Fairhead, J. D., Searle, R. C., and Sowerbutts, W. T. C. 1969. The Evolution of Rifting in Africa. Nature, 221*, pp. 1178- 1182. Goldich, So S., Nier, A. 0., Baadsgaard, H., Hoffman, J. H., and Krueger, H. W. 1961. The PrecambrianGeology and Geochronology of Minnesota. Minn. Geol. Surv. Bullo, 1*1, 193 pp. Grant, F. S. and West, G. F. 1965. Interpretation Theory in Applied Geophysics. New York, McGraw-Hill, 5bii pp. Green, J 0 C. 1971. Stratigraphy of the North Shore Volcanic Group Northeast of Silver Bay, Minnesota (abstract). Inst, on Lake Superior Geology, Duluth, Minnesota, pp. 20-22. Halls, H. C. 1966. A review of the Keweenawan Geologyof the Lake Superior Region. Amer. Geophy. Un.,Geophy. Mon. 10, pp. 3-270 Halls, Ho C. and West, G. F. 1971. A Seismic Refraction Survey in Lake Superior. Canadian Jour. Earth Sciences, 8, pp. 6IO-63O. Hamblin, W 0 K. 1958. The Cambrian Sandstones of Northern Michigan. Mich. Geol. Surv. Publ. 51, 11*5 PP» Hamblin, W. K. 1961. Paleogeographic Evolution of the Lake Superior Region from Late Cambrian Time. Geol. Soc. Amer. Bull., 72, pp. 1-18. \ 81 Hamblin, W. K. and Horner, V. J. 1961. Sources of the Keweenawan Con­ glomerates of Northern Michigan. Jour. Geol., 69, pp. 20U-211. Hinze, W. J. I960. Downward Continuation of Two-Dimensional Gravity and Magnetic Anomalies. Address Mich. Acad. Sci., Arts and Letters. Hinze, W. J. 1963. Regional Gravity and Magnetic Anomaly Maps of the Southern Peninsula of Michigan. Mich. Geol. Surv., Rept0 Inv. 1, 26 pp. Hinze, W. J., O'Hara, N. W., Trow, J. W,, and Secor, G. B. 1966. Aeromagnetic Studies of Eastern Lake Superior. Amer. Geophy. Un., Geophy. Mon. 10, pp. 95-110. Hinze, W. J. and Merritt, D. W. 1969. Basement Rocks of the Southern Peninsula of Michigan* Mich. Basin Soc. Guidebook, Studies of the Precamb. of the Mich. Basin, pp. 28-59. Hinze, ¥ 0 J., Davidson, D. M., and Roy, R. F. 1971. Continental Rifts (abstract). Inst, on Lake Superior Geology, Duluth, Minnesota, Po 29. Hotchkiss, W. 0. 1923. The Lake Superior Geosyncline. Bull. 3h> pp. 669-6780 Geol. Soc. Amer. Irving, R. D. 1883. The Copper-Bearing Rocks of Lake Superior. Geol. Surv. Mon, 5. U. S. James, H. L. 1958. Stratigraphy of Pre-Keweenawan Rocks in parts of Northern Michigan. U. S. Geol. Surv. Prof. Paper 31U-C, pp. 27-kk° James, W, R. 1966. Fortran IV Program Using Double Fourier Series for Surface Fitting of Irregularly Spaced Data. Computer Contribution 5# State Geol. Surv., University of Kansas, Lawrence, Kansas. Kellogg, R. L. 1971o An Aeromagnetic Investigation of the Southern Penin­ sula of Michigan. Ph. D0 Thesis, Michigan State University, 161 pp. Lyons, P. L. 1959. The Greenleaf Anomaly, a Significant Gravity Feature. Kansas Geol. Surv. Bull., 137, pp. 105-120. MacLaren, A. S. and Charbonneau, B. W. 1968. Characteristics of Magnetic Data Over Major Subdivisions of the Canadian Shield. Geol. Assoc. Canada Proc. 19, pp. 57-65* McGinnis, L. D. 1970. Tectonics and the Gravity Field in the Continental Interior. Jour. Geophy. Res0, 75* ppo 317-331® Meshref, W. M. and Hinze, W. J. 1970. Geologic Interpretation of Aero­ magnetic Data in Western Upper Peninsula of Michigan. Mtch. Geol. Surv., Rept. Inv. 12, 25 pp. 82 Meyer, R. P. 1961*. North American Seismic Group Experiment (abstract). Proceedings of the 36th Annual Meeting of the Eastern Section of the Seismological Society of America. Meyer, Ro P. and Ocola, L. 1967. Structural Dislocations in Lake Super­ ior from Seismic Profile Methods. Trans. Amer. Geophy. Un., i*8, p. 196. Newcombe, R. B. 1933* Surv. Pub. 38. Oil and Gas Fields of Michigan. Mich. Geol. O'Brien, P. N. S. 1968. Lake Superior Crustal Structure-A Reinterpretation of the 1963 Seismic Experiment. Jour. Geophy. Res., 73, pp. 2669-2689. Oetking, P. F 0 1951° The Relation of the Lower Paleozoic to the Older Rocks in the Northern Peninsula of Michigan. Ph. D. Thesis, University of Wisconsin. O'Hara, N. W. 1967* An Aeromagnetic and Geophysical Interpretation of the Precambrian Framework and Tectonic Structure of the Eastern Lake Superior Region. Ph. D. Thesis, Michigan State University, 260 pp. O'Hara, N. W. and Hinze, W. J 0 1971. Lake Michigan Aeromagnetic Survey (abstract). Inst, on Lake Superior Geology, Duluth, Minnesota, P° 56. Patenaude, R. W. 1961*. Results of Regional Aeromagnetic Surveys of Eastern Upper Michigan, Central Lower Michigan, and Southeastern Illinois. The University of Wisconsin Geophysical and Polar Research Center, Research Report Series 6U-2, 5l PP. Robinson, W. I. 1923° The Southern Limits of the Lake Superior Syncline. (Unpub. manscr.), Mich. Geol. Surv. Files. Rudman, A. J., Summerson, C. H., and Hinze, W. J. 1965. Geology of Base­ ment in Midwestern United States. Amer. Assoc. Pet. Geol. Bull., k9 , pp. 89U-9OJ4. Secor, G„ B., Hinze, W. J., O'Hara, N. W. , and Trow, J. W. 1967. An Aeromagnetic Survey of Lake Huron (abstract). Inst, on Lake Superior Geology, East Lansing, Michigan, p. 36. Sims, P. K. and Zietz, I. 1967. Aeromagnetic and Inferred Precambrian Paleogeologic Map of East-Central Minnesota and Part of Wisconsin. U. S. Geol. Surv. Geophy. Inv. Map GP-563. Slichter, L. B. 1951. Crustal Structure in the Wisconsin Area. of Naval Research Reports N9, 0NR-86200. Office 83 Smith, T. Jo, Steinhart, Jc S., and Aldrich, L0 T 0 1966. Lake Superior Crustal Studies! A Comparison of Several Techniques0 Journ« Geophy, Res0, 71, PP. 111*1-1172. Steinhart, J. S, and Smith, To Jo 1966. Time Terms and Structure in Western Lake Superior Region. Amer0 Geophy0 Un., Geophy0 Ifono 10, pp. 198-201*. Steinhart, J. S., Hart, S. R 0, and Smith, T. J. 1968. Heat Flow. An­ nual Report of the Director, Department of Terrestrial Magnetism, Carnegie Institution Year Book 66, pp. 130-11*0. Talwani, M 0, Worzel, J. L., and Landisman, M. 1999. Rapid Gravity Com­ putations for Two Dimensional Bodies with Applications to the Mendocino Submarine Fracture Zone. Joum. Geophy. Res., 61*, pp. 1i 9 - 99 . Taylor, R. B. 1961*. Geology of the Duluth Gabbro Complex. Surv. Bull., 1*1*, 61 pp. Minn. Geol. Thiel, E. C. 1996. Correlation of Gravity Anomalies with the Keweenawan Geology of Wisconsin and Minnesota. Geol. Soc. Amer. Bull., 67, pp0 1079-1100o Thomson, J0 E0 1993* Geology of Mamainse Point Copper Area. Report, Ontario, Dept, of Mines, 62, 29 pp. Annual Thwaites, F. To 1939* Sublacustrine Topographic Map of the Bottom of Lake Superior. Kansas Geol. Soc. Guidebook, Ninth Ann0 Field Confo, pp. 226-228. Trow, Jo Wo 1971* Personal Communication0 Tyler, S. A0, Marsden, Ro W 0, Grout, F» F., and Thiel, G0 A 0 191*0. Studies of the Lake Superior Precambrian by Accessory Mineral Methods. Geol. Soco Amer. Bull., 91, pp. 11*29-1937. Uo S. Geological Survey. 1970. Aeromagnetic Map of the MenomineeNorthland Area, Dickinson, Marquette, and Menominee Counties, Michigan and Marinette County, Wisconsin0 Geophy. Inv0 Map GP-711. Van Hise, C. R. and Bayley, W. S. 1900. Description of the Menominee Quadrangle. U 0 S. Geol. Surv. Geol. Atlas, Menominee Folio, 620 Van Hise, C0 R. and Leith, C. K 0 1911» The Geology of the Lake Superior Region. U. S. Geol. Surv„ Mono0 92, 61*1 pp. Weaver, D0 F. 1967. A Geological Interpretation of the Bouguer Anomaly Field of Newfoundland. Publ. Dorn. Obs. of Canada, 39, pp. 221-291* Weber, J. R. and Goodacre, A0 K. 1966. A Reconnaissance Underwater Gravity Survey of Lake Superior. Amer. Geophy. Un., Geophy. Mon. 10, pp. $6-6$, 8U White, V/0 So 1957. Regional Structural Setting of the Michigan Native Copper District. Inst, on Lake Superior Geology, Houghton, Michi­ gan, pp. 3-I60 White, W. S. I960. The Keweenawan Lavas of Lake Superior, an Example of Flood Basalts. Amer0 Jour, of Sci., 258-A, pp. 367 -371*. White, W. S. 1966a. Tectonics of the Keweenawan Basin, Western Lake Superior Region. U 0 S. Geol. Surv. Profo Paper, 52l*-E, 23pp. White, Wo S. 1966b. Geologic Evidence for Crustal Structure in the Western Lake Superior Basin. Amer. Geophy. Un 0 Geophy. Mon. 10, pp. 28-1*1 . Whitten, E 0 H. and Beckman, W. A 0 1969. Fold Geometry Within Part of Michigan Basin, Michigan. Amer. Assoc. PeteGeol. Bull., £ 3, pp. 101*3-1057. ' Wold, R. J. and Ostenso, N. A. 1966. Aeromagnetic, Gravity and Sub-Bottom Profiling Studies in Western Lake Superior# Amer. Geophy. Un., Geophy. Mon. 10, pp. 66-91*. Woollard, G. P. and Bonini, W. E. 1955= A Review of the Factors Affecting the Accuracy of Long Distance and Intercontinental Gravimeter Mea­ surements. Trans. Amer. Geophy. Un., 36, pp. 575-583. Young, G. M. 19660 Huronian Stratigraphy of the McGregor Bay Area, Ontario: Relevance to the Paleogeography of the Lake Superior Region. Can. Jour. Earth Sci., 3, pp. 203-210. Zletz, I. 1965# Aeromagnetic Study of the Midcontinent Gravity Anomaly (abstract). U. S. Progress Report Intern0 Upper Mantle Project. Zietz, I. and Kirby, J. R. 1971. Aeromagnetic Map of the Western Part of the Northern Peninsula, Michigan and Part of Northern Wisconsin. U. S. Geol. Surv. Geophy. Inv. Map GP-750q APPENDIX FREE AIR GRAVITY ANOMALY MAP OF THE EASTERN PORTION OF THE NORTHERN PENINSULA OF MICHIGAN FREE AIR GRAVITY ANOMALY MAP Figure 16.-Free-air Gravity Anomaly Hap of the eastern portion of the Northern Peninsula of Michigan,,