Geochemical Characterization of Trace MVT Mineralization in Paleozoic Sedimentary Rocks of Northeastern Wisconsin , USA

Disseminated Mississippi Valley-type (MVT) mineralization occurs throughout northeastern Wisconsin, USA, and is recognized as the source of regionally extensive natural groundwater contamination in the form of dissolved arsenic, nickel, and other related metals. Although considerable attention has been given to arsenic contamination of groundwater in the region, limited attention has been focused on characterizing the bedrock sources of these and other metals. A better understanding of the potential sources of groundwater contamination is needed, especially in areas where groundwater is the dominant source of drinking water. This article describes the regional, stratigraphic, and petrographic distribution of MVT mineralization in Paleozoic rocks of northeastern Wisconsin, with a focus on sulfide minerals. Whole-rock geochemical analysis performed on 310 samples of dolomite, sandstone, and shale show detectable levels of arsenic, nickel, cobalt, copper, lead, zinc, and other metals related to various sulfide mineral phases identified using scanning electron microscopy. MVT minerals include pyrite, marcasite, sphalerite, galena, chalcopyrite, fluorite, celestine, barite, and others. We describe the first nickeland cobalt-bearing sulfide mineral phases known from Paleozoic strata in the region. Arsenic, nickel, and cobalt are sometimes present as isomorphous substitutions in pyrite and marcasite, but discrete mineral phases containing nickel and cobalt elements are also observed, including bravoite and vaesite. Locally abundant stratigraphic zones of sulfide minerals occur across the region, especially in the highly enriched Sulfide Cement Horizon at the top of the Ordovician St. Peter Sandstone. Abundant quantities of sulfides also appear near the contact between the Silurian Mayville Formation and the underlying Maquoketa and Neda formations in certain areas along and east of the Niagara escarpment. This article illustrates how a detailed geochemical and mineralogical investigation can yield a better understanding of groundwater quality problems.


Introduction
Mississippi Valley-Type (MVT) lead-zinc ore deposits are carbonate-hosted accumulations of sulfide and associate minerals that are dominated by sphalerite, galena, iron sulfides, and other associated minerals.They typically occur in Phanerozoic dolostone (with some in limestone or sandstone) located on the flanks of sedimentary basins, orogenic forelands, or foreland thrust belts [1][2][3].MVT ore deposits are located throughout the world, but were named for deposits that occur along the along the Mississippi River Valley in North America.Large, well-studied deposits occur in the United States and Canada [1][2][3][4][5], Europe (especially Ireland and Poland) [3,[5][6][7], and elsewhere.Modern interpretations of the origin of MVT ore deposits generally include large-scale brine migration driven within platform carbonates by large-scale tectonic events [3].
The upper Mississippi Valley region of southwest Wisconsin, eastern Iowa, and northwestern Illinois, USA (Figure 1) contains a historic Mississippi Valley-type (MVT) ore district [8].Although the lead mineralization in the region was discovered during the late 17th century, these deposits did not receive major attention until the early 19th century.It was in this region that the first large "metals rush" in the United States took place in what is known as the Upper Mississippi Valley Zinc-Lead Ore District.Galena, and later sphalerite, were the two principal minerals of interest in small ore bodies located near folds and faults hosted by Paleozoic carbonate rocks [8].The likelihood of a hydrothermal origin for these deposits was recognized by the mid-20th century by Heyl et al. and McLimans [8,9].Research by Sverjensky, Garven, Bethke, Rowan and others during the 1980s and 1990s led to the general acceptance that these deposits were the result of northward basin-scale brine migration from the Illinois basin [1,2,10,11].The economic viability of these deposits faded during the later 20th century, and the last mine in the district closed in 1979.Base map modified after [12,13].
Similar MVT mineralization is also present throughout eastern Wisconsin, which was initially thought to represent outlying deposits of the main ore district [14,15].Limited exploration for MVT mineralization took place in southeastern Wisconsin during the early 1980s through subsurface coring program in eight southeastern Wisconsin counties by Mobil Mineral Resources, Inc. [16].However, extensive Pleistocene glacial overburden, coupled with a limited understanding of the along the Mississippi River Valley in North America.Large, well-studied deposits occur in the United States and Canada [1][2][3][4][5], Europe (especially Ireland and Poland) [3,[5][6][7], and elsewhere.Modern interpretations of the origin of MVT ore deposits generally include large-scale brine migration driven within platform carbonates by large-scale tectonic events [3].
The upper Mississippi Valley region of southwest Wisconsin, eastern Iowa, and northwestern Illinois, USA (Figure 1) contains a historic Mississippi Valley-type (MVT) ore district [8].Although the lead mineralization in the region was discovered during the late 17th century, these deposits did not receive major attention until the early 19th century.It was in this region that the first large "metals rush" in the United States took place in what is known as the Upper Mississippi Valley Zinc-Lead Ore District.Galena, and later sphalerite, were the two principal minerals of interest in small ore bodies located near folds and faults hosted by Paleozoic carbonate rocks [8].The likelihood of a hydrothermal origin for these deposits was recognized by the mid-20th century by Heyl et al. and McLimans [8,9].Research by Sverjensky, Garven, Bethke, Rowan and others during the 1980s and 1990s led to the general acceptance that these deposits were the result of northward basin-scale brine migration from the Illinois basin [1,2,10,11].The economic viability of these deposits faded during the later 20th century, and the last mine in the district closed in 1979.Base map modified after [12,13].
Similar MVT mineralization is also present throughout eastern Wisconsin, which was initially thought to represent outlying deposits of the main ore district [14,15].Limited exploration for MVT mineralization took place in southeastern Wisconsin during the early 1980s through subsurface coring program in eight southeastern Wisconsin counties by Mobil Mineral Resources, Inc. [16].However, extensive Pleistocene glacial overburden, coupled with a limited understanding of the along the Mississippi River Valley in North America.Large, well-studied deposits occur in the United States and Canada [1][2][3][4][5], Europe (especially Ireland and Poland) [3,[5][6][7], and elsewhere.Modern interpretations of the origin of MVT ore deposits generally include large-scale brine migration driven within platform carbonates by large-scale tectonic events [3].
The upper Mississippi Valley region of southwest Wisconsin, eastern Iowa, and northwestern Illinois, USA (Figure 1) contains a historic Mississippi Valley-type (MVT) ore district [8].Although the lead mineralization in the region was discovered during the late 17th century, these deposits did not receive major attention until the early 19th century.It was in this region that the first large "metals rush" in the United States took place in what is known as the Upper Mississippi Valley Zinc-Lead Ore District.Galena, and later sphalerite, were the two principal minerals of interest in small ore bodies located near folds and faults hosted by Paleozoic carbonate rocks [8].The likelihood of a hydrothermal origin for these deposits was recognized by the mid-20th century by Heyl et al. and McLimans [8,9].Research by Sverjensky, Garven, Bethke, Rowan and others during the 1980s and 1990s led to the general acceptance that these deposits were the result of northward basin-scale brine migration from the Illinois basin [1,2,10,11].The economic viability of these deposits faded during the later 20th century, and the last mine in the district closed in 1979.Base map modified after [12,13].
Similar MVT mineralization is also present throughout eastern Wisconsin, which was initially thought to represent outlying deposits of the main ore district [14,15].Limited exploration for MVT mineralization took place in southeastern Wisconsin during the early 1980s through subsurface coring program in eight southeastern Wisconsin counties by Mobil Mineral Resources, Inc. [16].However, extensive Pleistocene glacial overburden, coupled with a limited understanding of the Similar MVT mineralization is also present throughout eastern Wisconsin, which was initially thought to represent outlying deposits of the main ore district [14,15].Limited exploration for MVT mineralization took place in southeastern Wisconsin during the early 1980s through subsurface coring program in eight southeastern Wisconsin counties by Mobil Mineral Resources, Inc. [16].However, extensive Pleistocene glacial overburden, coupled with a limited understanding of the distribution and character of geologic structures and the absence of economic petroleum deposits, has hampered efforts to fully characterize the region's subsurface geology.This article focuses on a broad region of Paleozoic-hosted MVT mineralization in northeastern Wisconsin located about 250 km (150 miles) northeast of the Upper Mississippi Valley Zinc-Lead ore district and the Illinois basin (Figure 1).
During the 1970s and 1980s, Paleozoic rocks in Wisconsin received further attention from carbonate sedimentologists who focused on Ordovician rocks as a potential ancient analog for the formation of low-temperature dolomite.This model, known either as the "Dorag" model or the mixing-zone model, involved interaction between fresh water and seawater along a coastal mixing zone in southern Wisconsin [17].Despite an early call for a hydrothermal origin for the dolomite [18], the low-temperature interpretation continued to be favored by most sedimentologists who worked in the region.In fact, the Dorag model became a paradigm for dolomitization models elsewhere, and remained so throughout much of the 1980s and 1990s (e.g., [19][20][21]).
The geochemical interpretations for the origin of the dolomite and MVT mineralization in eastern Wisconsin were in stark contrast with one another until conclusive data were presented that supported a hydrothermal origin for both [22,23].This hydrothermal system was interpreted to be a distinct, unrelated flow system from the one that formed the Upper Mississippi Valley ore district.Subeconomic MVT mineralization extends throughout the entire region of eastern Wisconsin into the Upper Peninsula of Michigan.Luczaj [22,23] used petrographic, isotopic, and fluid-inclusion data to propose a common hydrothermal origin for the MVT mineralization, pervasive epigenetic dolomitization, and K-silicate mineralization in the region.His research mainly focused on an evaluation of the mixing-zone versus hydrothermal epigenetic dolomitization models for the formation of pervasively dolomitized strata in the region.

Geologic Setting of the Study Area
The study area in eastern Wisconsin region lies on the western edge of the ancestral Michigan basin and is located to the east of the Wisconsin arch and to the south of the Superior Craton (Figure 1).As much as 500 to 800 m of Lower and Middle Paleozoic quartz sandstone, dolostone, and shale are present (Figures 2 and 3), and the strata thicken and dip gently toward the Michigan basin (Figure 4) where they are overlain by younger Paleozoic and Mesozoic sedimentary rocks.
Numerous gentle folds, faults, and fractures are observed in Paleozoic rocks in the region, despite a geographic position near the middle of the North American craton away from orogenic belts.Dip-slip and strike-slip faults are present within the study area, based upon water well records and quarry exposures, but few are directly observed because they are concealed by extensive Pleistocene glacial drift [23][24][25].The carbonate rocks throughout the study area are completely dolomitized over all stratigraphic intervals.Faults and joints were certainly preferred conduits for mineralizing groundwater because sulfide minerals are preferentially precipitated along these surfaces throughout most of the study area.Stratigraphic reconstructions and organic maturity data suggest that the thickness of eroded sediments in eastern Wisconsin was probably much less than 1-1.5 km at any time during the burial history of these rocks [23].

Origin of the MVT Mineralization in Eastern Wisconsin
Luczaj [22,23] proposed the existence of a regional hydrothermal system in eastern Wisconsin that was active during the Late Devonian-Mississippian periods.The precise timing is not well defined, but could be earlier than the Late Paleozoic system responsible for mineralization of the Upper Mississippi Valley Zinc-Lead ore district [26].The rocks in eastern Wisconsin preserve a pervasive hydrothermal signature and contain abundant epigenetic dolomite, subeconomic MVT mineralization, and minor authigenic K-silicate mineralization.Fluid-inclusion, stable isotopic, and petrographic evidence suggests that the MVT mineralization, K-silicate mineralization, and dolomitization in eastern Wisconsin are genetically related to the same regional hydrothermal system that operated at temperatures between about 60 and 120 ˝C, with sulfide mineralization likely occurring between about 80 and 110 ˝C [22,23].This ancient system was responsible for precipitation and replacement of a suite of MVT sulfide minerals that includes pyrite and marcasite (FeS 2 ), galena (PbS), sphalerite (ZnS), chalcopyrite (CuFeS 2 ), fluorite (CaF 2 ), and other related minerals.It is this hydrothermal system that also precipitated the arsenic, nickel, and cobalt-bearing sulfide minerals in the region, including unusually enriched sulfide cement horizons (SCH) in sandstones such as the St.Peter Sandstone of the Middle Ordovician Ancell Group (Figure 2).and replacement of a suite of MVT sulfide minerals that includes pyrite and marcasite (FeS2), galena (PbS), sphalerite (ZnS), chalcopyrite (CuFeS2), fluorite (CaF2), and other related minerals.It is this hydrothermal system that also precipitated the arsenic, nickel, and cobalt-bearing sulfide minerals in the region, including unusually enriched sulfide cement horizons (SCH) in sandstones such as the St.
Peter Sandstone of the Middle Ordovician Ancell Group (Figure 2).Pleistocene sediments and younger Silurian and Devonian strata for southeastern Wisconsin are not shown here for simplicity, but can be found elsewhere [27,28].
The sequence of mineralization in the study area is somewhat variable, but generally consists of early dolomite and quartz, followed by K-feldspar and MVT mineralization, followed by late calcite.Dolomite precipitation continues throughout the main stage of sulfide mineralization.However, unlike many MVT-related dolomites, the dolomite of eastern Wisconsin is almost exclusively planar (non-saddle).Only a few small vugs with examples of saddle dolomite have been observed in the study area, and cathodoluminescence suggests they have the same cement stratigraphy as some of the planar dolomite [23].Some late-stage euhedral planar dolomite crystals are up to 1 cm in diameter.Pleistocene sediments and younger Silurian and Devonian strata for southeastern Wisconsin are not shown here for simplicity, but can be found elsewhere [27,28].
The sequence of mineralization in the study area is somewhat variable, but generally consists of early dolomite and quartz, followed by K-feldspar and MVT mineralization, followed by late calcite.Dolomite precipitation continues throughout the main stage of sulfide mineralization.However, unlike many MVT-related dolomites, the dolomite of eastern Wisconsin is almost exclusively planar (non-saddle).Only a few small vugs with examples of saddle dolomite have been observed in the study area, and cathodoluminescence suggests they have the same cement stratigraphy as some of the planar dolomite [23].Some late-stage euhedral planar dolomite crystals are up to 1 cm in diameter.[29], B is Thornburg and Sahai [30], C is Burkel and Stoll [31], and D is Schreiber et al. [32]).Additional references for region-wide studies include Riewe et al. [33] and Johnson and Riewe [34], as well as a cancer study by Knobeloch et al. [35] that focused on Winnebago and Outagamie counties.Base map after Mudrey et al. [36] shows the following units: pЄ = Precambrian, Є = Cambrian sandstones, Opc = Ordovician Prairie du Chien Group, Oa = Ordovician Ancell Group, Os = Ordovician Sinnipee Group, Om = Ordovician Maquoketa-Neda formations, and Su = Silurian dolostones.[29], B is Thornburg and Sahai [30], C is Burkel and Stoll [31], and D is Schreiber et al. [32]).Additional references for region-wide studies include Riewe et al. [33] and Johnson and Riewe [34], as well as a cancer study by Knobeloch et al. [35] that focused on Winnebago and Outagamie counties.Base map after Mudrey et al. [36] shows the following units: p ississippi River Valley in North America.Large, well-studied deposits occur in the United Canada [1][2][3][4][5], Europe (especially Ireland and Poland) [3,[5][6][7], and elsewhere.Modern tions of the origin of MVT ore deposits generally include large-scale brine migration driven tform carbonates by large-scale tectonic events [3].pper Mississippi Valley region of southwest Wisconsin, eastern Iowa, and northwestern SA (Figure 1) contains a historic Mississippi Valley-type (MVT) ore district [8].Although ineralization in the region was discovered during the late 17th century, these deposits did major attention until the early 19th century.It was in this region that the first large "metals e United States took place in what is known as the Upper Mississippi Valley Zinc-Lead ct.Galena, and later sphalerite, were the two principal minerals of interest in small ore ated near folds and faults hosted by Paleozoic carbonate rocks [8].The likelihood of a mal origin for these deposits was recognized by the mid-20th century by Heyl et al. and [8,9].Research by Sverjensky, Garven, Bethke, Rowan and others during the 1980s and o the general acceptance that these deposits were the result of northward basin-scale brine from the Illinois basin [1,2,10,11].The economic viability of these deposits faded during 0th century, and the last mine in the district closed in 1979.ar MVT mineralization is also present throughout eastern Wisconsin, which was initially represent outlying deposits of the main ore district [14,15].Limited exploration for MVT tion took place in southeastern Wisconsin during the early 1980s through subsurface The upper Mississippi Valley region of southwest Wisconsin, eastern Iowa, and northwestern Illinois, USA (Figure 1) contains a historic Mississippi Valley-type (MVT) ore district [8].Although the lead mineralization in the region was discovered during the late 17th century, these deposits did not receive major attention until the early 19th century.It was in this region that the first large "metals rush" in the United States took place in what is known as the Upper Mississippi Valley Zinc-Lead Ore District.Galena, and later sphalerite, were the two principal minerals of interest in small ore bodies located near folds and faults hosted by Paleozoic carbonate rocks [8].The likelihood of a hydrothermal origin for these deposits was recognized by the mid-20th century by Heyl et al. and McLimans [8,9].Research by Sverjensky, Garven, Bethke, Rowan and others during the 1980s and 1990s led to the general acceptance that these deposits were the result of northward basin-scale brine migration from the Illinois basin [1,2,10,11].The economic viability of these deposits faded during the later 20th century, and the last mine in the district closed in 1979.Similar MVT mineralization is also present throughout eastern Wisconsin, which was initially thought to represent outlying deposits of the main ore district [14,15].Limited exploration for MVT mineralization took place in southeastern Wisconsin during the early 1980s through subsurface = Cambrian sandstones, Opc = Ordovician Prairie du Chien Group, Oa = Ordovician Ancell Group, Os = Ordovician Sinnipee Group, Om = Ordovician Maquoketa-Neda formations, and Su = Silurian dolostones.

Relevance to Water Quality Issues
Dissolved arsenic has been recognized as a significant water quality problem that affects millions across the world and can lead to health problems including lung, bladder, and skin cancers.Bangladesh, for example, is widely known for having high dissolved arsenic resulting from sulfide mineralization, with the situation being described as a public health emergency (e.g., [37]).Three main regions in Wisconsin exhibit elevated levels of arsenic in groundwater [13,35].Since its discovery in the late 1980s as part of a routine feasibility study for a landfill proposed in Winnebago County [33], the problem of arsenic contamination in the region's groundwater has received significant attention, with most work occurring in Winnebago and Outagamie counties (Figure 3).Eastern Wisconsin's MVT mineralization is recognized as the source of regionally extensive groundwater contamination in the form of dissolved arsenic, nickel, and other related metals [13].The focus of previous research in the region [29][30][31][32][33][34][38][39][40] was to characterize the distribution and release mechanisms of arsenic in well waters and aquifer host rocks (mainly Ancell Group), primarily in parts of the eastern Wisconsin counties of Winnebago, Outagamie, and Brown (Figure 3).Geochemical mechanisms proposed for arsenic release include oxidation of arsenic-bearing sulfide minerals (most significant), reductive dissolution of iron (hydr)oxides, and release of sorbed arsenic from mineral surfaces [29,30,39].Arsenic (As) concentrations vary from less than 1 µg/L to over 15,000 µg/L [33,41] in the study area, which contains several tens of thousands of private water wells and hundreds of municipal and other high capacity wells.In some townships (roughly 100 km 2 ), 20%-40% of the wells tested above the 10 µg/L health standard set for As in drinking water by the United States Environmental Protection Agency [13].Although wells are less commonly tested for nickel and cobalt, these metals were also observed in some wells at elevated concentrations, especially in Winnebago and Outagamie counties.An evaluation of publicly available data through the Wisconsin Department of Natural Resources [41] indicates that in four counties in the study area (Winnebago, Outagamie, Brown, and Green Lake), cobalt exceeded the Wisconsin groundwater quality enforcement standard of 40 µg/L in 5.4% to 8.8% of the wells tested.Nickel showed a similar pattern, with six counties in the study area (Winnebago, Outagamie, Brown, Marinette, Fond du Lac, and Manitowoc) exceeding the Wisconsin groundwater quality enforcement standard of 100 µg/L in 1.8% to 5.7% of the wells tested.The National Uranium Resource Evaluation (NURE) [42,43]

Relevance to Water Quality Issues
Dissolved arsenic has been recognized as a significant water quality problem that affects millions across the world and can lead to health problems including lung, bladder, and skin cancers.Bangladesh, for example, is widely known for having high dissolved arsenic resulting from sulfide mineralization, with the situation being described as a public health emergency (e.g., [37]).Three main regions in Wisconsin exhibit elevated levels of arsenic in groundwater [13,35].Since its discovery in the late 1980s as part of a routine feasibility study for a landfill proposed in Winnebago County [33], the problem of arsenic contamination in the region's groundwater has received significant attention, with most work occurring in Winnebago and Outagamie counties (Figure 3).Eastern Wisconsin's MVT mineralization is recognized as the source of regionally extensive groundwater contamination in the form of dissolved arsenic, nickel, and other related metals [13].The focus of previous research in the region [29][30][31][32][33][34][38][39][40] was to characterize the distribution and release mechanisms of arsenic in well waters and aquifer host rocks (mainly Ancell Group), primarily in parts of the eastern Wisconsin counties of Winnebago, Outagamie, and Brown (Figure 3).Geochemical mechanisms proposed for arsenic release include oxidation of arsenic-bearing sulfide minerals (most significant), reductive dissolution of iron (hydr)oxides, and release of sorbed arsenic from mineral surfaces [29,30,39].Arsenic (As) concentrations vary from less than 1 µg/L to over 15,000 µg/L [33,41] in the study area, which contains several tens of thousands of private water wells and hundreds of municipal and other high capacity wells.In some townships (roughly 100 km 2 ), 20%-40% of the wells tested above the 10 µg/L health standard set for As in drinking water by the United States Environmental Protection Agency [13].Although wells are less commonly tested for nickel and cobalt, these metals were also observed in some wells at elevated concentrations, especially in Winnebago and Outagamie counties.An evaluation of publicly available data through the Wisconsin Department of Natural Resources [41] indicates that in four counties in the study area (Winnebago, Outagamie, Brown, and Green Lake), cobalt exceeded the Wisconsin groundwater quality enforcement standard of 40 µg/L in 5.4% to 8.8% of the wells tested.Nickel showed a similar pattern, with six counties in the study area (Winnebago, Outagamie, Brown, Marinette, Fond du Lac, and Manitowoc) exceeding the Wisconsin groundwater quality enforcement standard of 100 µg/L in 1.8% to 5.7% of the wells tested.The National Uranium Resource Evaluation (NURE) [42,43], an older dataset from 1976 to 1979, also showed elevated dissolved Ni, Co, and As in the study area.Most of the wells with exceedances draw water from the Cambrian-Ordovician sandstone aquifers, especially the Ancell Group.Even municipal wells in Brown County 240-320 m deep have been impacted by elevated As, Co, and Ni during aquifer storage and recovery (ASR) tests or when valve leaks occurred that allowed oxygenated surface water into the Cambrian aquifer [44][45][46].
The most seriously impacted well studied in the region was one in northeastern Outagamie County and was drawing from the Ancell Group (St.Peter Sandstone).The well had the following chemistry (tested 1991-1992): pH = 2.05, As = 4300 µg/L, Co = 5500 µg/L, Cd = 220 µg/L, Cr = 84 µg/L, Ni = 11,000 µg/L, Al = 15,000 µg/L, and Pb = 400 µg/L [13,41].The incidence of As, Co, and Ni in drinking water wells across the study area demonstrates the need for a better understanding of the link between groundwater quality and MVT mineralization in the region.
Additional work by Knobeloch and others [35] (Figure 3) on the health implications of arsenic in the study area has indicated an increased prevalence of skin cancer in the region.In 2004, the Wisconsin Department of Natural Resources designated a "Special Well Casing Pipe Depth Area" for all of Outagamie and Winnebago counties, which included special requirements for drilling methods, disinfection methods, and grouting requirements [13].Arsenic is a major priority among regulators and local health departments in the region.Less attention has been paid to counties to the north, and they do not have special well casing requirements, despite their nearly identical geology.
The chemistry of a groundwater system is controlled in large part by the composition of its host rock [47].Sulfide minerals in eastern Wisconsin Ordovician carbonate rocks coat the walls of joints, bedding planes, faults, and macropores in these dolostones [23].In Cambrian and Ordovician sandstones, sulfides are observed as isolated nodules, stringers, and sometimes as continuous zones of cementation, especially in the Ancell Group.Oxidation of these sulfide minerals has been shown to be the major source of arsenic and other metals of concern in eastern Wisconsin [29,30,39].This occurs most often near recharge areas where the aquifers are shallow (Figure 4) or in areas of intense drawdown where the aquifer becomes partially dewatered.
While an initial stratigraphic and petrographic framework was known for bedrock mineralization in Winnebago and Outagamie counties, limited data were available regarding the whole-rock and trace element composition of the bedrock, especially in areas to the north in Marinette, Oconto, and Shawano counties (Figure 3).A better understanding of potential sources of groundwater contamination is needed, especially in these rural areas where domestic well supplies are the dominant source of drinking water.
Although the MVT minerals are unlikely to be mined within the study area, they will continue to affect groundwater quality indefinitely.Groundwater is an important source of drinking water in Wisconsin.The percentage of the population that consumes water produced by private or municipal wells ranges from 18% to 100% of the population in each county.Overall, approximately 55% of the population in the study area uses groundwater, with most surface water users limited to a few of the larger municipalities in Brown, Winnebago, and Outagamie counties [48].Many municipalities in the study area that currently use surface water have wells that penetrate the Cambrian-Ordovician sandstones in the region (Figure 4) that are kept as a backup in the event of a pipeline malfunction.
This article presents an analysis of the regional, stratigraphic, and petrographic distribution of MVT mineralization in Paleozoic rocks of northeastern Wisconsin and had two main objectives.The first was to document the origin and characteristics of the MVT-related mineralization in eastern Wisconsin.This included field and laboratory mineralogical observations along with a regional whole-rock geochemical study of bedrock from several Paleozoic stratigraphic units, with a focus on heavy metals composition.In this context, we present the first evidence in northeastern Wisconsin of discrete nickel-and cobalt-bearing mineral phases, along with a complete list of MVT-related minerals known from the region.A second objective was to provide a regional foundation for current and future groundwater studies in the region in order to better understand the causes of dissolved arsenic and other metals in the groundwater of northeastern Wisconsin.We hope this work will provide a baseline for future economic geologists and groundwater specialists in this region and in others with similar geology.

Materials and Methods
The study area includes several counties in northeastern Wisconsin, comprising approximately 13,000 km 2 (Figure 3).The sources of rock samples from the study area included quarries, outcrops, well cuttings, and limited drill cores.Design of a sampling strategy for geochemical characterization was affected by two major limitations.First, the distribution of outcrop and core samples was not random in either the stratigraphic or regional sense due to extensive Quaternary glacial deposits throughout the entire study area.Outcrops are more common near escarpments formed where Silurian and Ordovician dolostones overlie weaker shale or sandstone.Quarry locations are preferentially distributed near the base of the Silurian Mayville Formation, near the base of the Platteville Formation, and near the base of the Prairie du Chien Group carbonates (Figure 2).This outcrop distribution inevitably led to a biased stratigraphic sampling.A second limitation was that physical access to samples near quarry high walls is restricted for safety reasons.We made reasonable attempts to sample multiple stratigraphic levels wherever possible, including sets of well cuttings and drill cores.

Sample Collection
For whole-rock geochemical analysis, we collected and designated two principal categories of samples.Samples that were representative of the stratigraphic horizon or those that were taken at random are designated as "Matrix" specimens.Online Supplementary Materials Table S1 contains location numbers and sample descriptions.In contrast, sampling that was deliberately biased in an attempt to analyze the most sulfide-rich materials along fractures, bedding planes, and the SCH are designated as "Enriched".
The Paleozoic stratigraphic section (Figure 2) was divided into six principal stratigraphic units for the purpose of collection and analysis.These are the Cambrian sandstones, the Prairie du Chien Group (PDC) carbonates, the Ancell Group sandstone and shale, the Sinnipee Group carbonates, the Maquoketa Shale/Neda Ironstone, and the Silurian carbonates.Numerous samples were obtained from each of these stratigraphic units.

Scanning Electron Microscopy with Energy Dispersive Spectroscopy (SEM-EDS)
Samples of polished rock slabs, crystal-lined vugs, and well-cutting chips were carbon-coated before analysis.Samples were analyzed with a TESCAN Scanning Electron Microscope (SEM), model VEGA TS 5136SB with an excitation voltage of 30 kV.Elemental analysis was conducted using Energy Dispersive Spectroscopy (EDS) on a Genesis Spectrum EDAX detector.SEM images were obtained using backscattered electrons (BSE).This technique was used to explore the chemistry and petrographic distribution of metal sulfide minerals and associated minerals in over 30 selected samples.

Whole-Rock Geochemistry and Data Analysis
Matrix samples submitted for whole-rock analysis were typically 1-2 kg, except for drill cores and cuttings, which typically had a mass of at least 0.02 kg.Specimens indicated as "Enriched" were as large as possible, but generally were at least 0.2 kg.Digestion of 208 whole-rock samples was performed in 2007-2008 by ALS Chemex of Thunder Bay, Ontario, Canada, using Conventional 35-element Inductively Coupled Plasma-Atomic Emission Spectroscopy (ICP-AES) (ALS Method ME_ICP41) after Aqua Regia Digestion in a graphite heating block.Subsequent analysis of an additional 102 samples was performed in 2013 by ALS Minerals of Vancouver, British Columbia, Canada using two methods.Specifically, conventional 51-element ICP-AES (ALS Method ME_MS41) and ultra-low detection ICP-AES (ALS Method ME-MS41L) were performed after Aqua Regia Digestion in a graphite heating block.The results of these analyses yielded whole-rock concentrations for 35 or 51 different elements (depending on analysis), including arsenic (As), cadmium (Cd), cobalt (Co), chromium (Cr), copper (Cu), nickel (Ni), lead (Pb), zinc (Zn), molybdenum (Mo), and vanadium (V).Lower detection limits for some of the elements analyzed improved between 2008 and 2013, but this mostly affected elements that were not the focus of this study.Additional "overlimit" analyses, conducted with standard laboratory methods, were required for certain samples in which sulfur (S), iron (Fe), and phosphorous (P) exceeded instrument limits.ALS Minerals is accredited by the Standards Council of Canada and reports analytical precision of ˘10% for the methods used in this study.Quality control procedures included sample duplicates during sample crushing, and sample blanks in every analytical run.
Sample numbers, localities, analytical methods, and whole-rock geochemistry for a selected list of elements are presented in Online Supplementary Materials Table S1.Raw whole-rock geochemistry data for all samples are presented in Tables S2A and S2B.
Correlations were estimated between log 10 -transformed concentrations for selected elements.Elemental concentrations that were below detection limits were set to 0.55 times their detection limit for a given element to avoid missing data after log transformation [49].Both Pearson's r and Spearman's ρ (a non-parametric measure of correlation) were calculated as a sensitivity analysis due to some remaining non-linearity between some pairs of elements even after transformation, with both measures available in Supplementary Materials Tables S2C and S2D.R version 3.3.0for Windows was used for all analyses and the construction of plots.

Field and Core Observations
Observations made at quarries and outcrops clearly indicate that the regional distribution of MVT mineralization extends throughout the study area northward into the state of Michigan.The most abundant minerals found in vugs and fractures tended to be dolomite and calcite, with quartz present locally.Sulfide mineralization was dominated by pyrite and marcasite, but sphalerite, chalcopyrite, and galena were also observed in hand specimens.Additional macroscopic minerals included fluorite, celestine, strontianite, and barite.
The distribution of sulfide minerals (and related weathering products) in the host sedimentary rocks was heterogeneous, but stratigraphically predictable.Throughout northeastern Wisconsin, the stratigraphic unit with the most sulfide mineralization was the top of the Ordovician Ancell Group (dominated by the St.Peter Sandstone in this portion of the state).This zone of naturally enriched iron sulfide mineralization, known as the Sulfide Cement Horizon (SCH), occurs across eastern Wisconsin from the Illinois border in the south to the Michigan border in the north ( [13,23], this study).In the study area, the SCH was exposed on the quarry floors or in sump trenches at several quarries in Shawano, Oconto, Marinette, Outagamie, and Winnebago counties and areas southwest of the study area (Figure 5).No distinct ore deposits made any one part of the study area more concentrated than another.Cores and cuttings also show that the SCH is present down dip to the east of the outcrop belt, from near the surface in quarries to depths of several hundred meters (Figures 2 and 4).Nodules and intergranular cements of pyrite and marcasite were abundant where this group of rocks is rarely exposed.In one quarry near the city of Oshkosh in Winnebago County, a pink coating of amorphous material similar to erythrite (hydrated cobalt arsenate) and annabergite (hydrated nickel arsenate) was observed as a weathering product near the SCH in the St.Peter Sandstone, and was confirmed by SEM-EDS analysis.Another quarry in Green Lake County (southwest of the study area) contained primarily calcite and sphalerite cement in the SCH.A thin (~10-50 cm) layer of dark brown shale above the St.Peter Sandstone was interpreted to be the Glenwood Shale and was present at many locations.It contained abundant fine-grained sulfides, including pyritized arthropod and bryozoan fossils in some quarries.The next most mineralized stratigraphic zone occurred near the contact between the Silurian Mayville Dolomite and the underlying Ordovician Maquoketa Shale and Neda formations (Figure 6).Where present, the Upper Ordovician Neda Ironstone and related strata appeared to have sulfides associated with the upper extent of those units.In places, the mass of sulfides was striking, with large (10-40 cm thick) accumulations of pyrite easily identified, especially in central Brown County along portions of the Niagara Escarpment.This sulfide mineralization typically extended upward into the lower few meters of the Silurian Mayville Dolomite.In at least one case, a local aggregate producer The next most mineralized stratigraphic zone occurred near the contact between the Silurian Mayville Dolomite and the underlying Ordovician Maquoketa Shale and Neda formations (Figure 6).Where present, the Upper Ordovician Neda Ironstone and related strata appeared to have sulfides associated with the upper extent of those units.In places, the mass of sulfides was striking, with large (10-40 cm thick) accumulations of pyrite easily identified, especially in central Brown County along portions of the Niagara Escarpment.This sulfide mineralization typically extended upward into the lower few meters of the Silurian Mayville Dolomite.In at least one case, a local aggregate producer had problems with the "aesthetic" quality of its product due to rust stains from the weathering of pyrite fragments in asphalt and concrete products.This interval was best exposed in quarries and natural outcrops along Scray Hill in the towns of Ledgeview and Glenmore in central Brown County.Significant, but lesser amounts of pyrite were present at this contact and within the upper part of the Maquoketa Formation where the Neda Formation is not present.had problems with the "aesthetic" quality of its product due to rust stains from the weathering of pyrite fragments in asphalt and concrete products.This interval was best exposed in quarries and natural outcrops along Scray Hill in the towns of Ledgeview and Glenmore in central Brown County.
Significant, but lesser amounts of pyrite were present at this contact and within the upper part of the Maquoketa Formation where the Neda Formation is not present.The Ordovician Sinnipee Group carbonates contained significant quantities of metal sulfides, even at stratigraphic intervals farther away from the SCH in the Ancell Group.Sulfide mineralization was locally abundant along bedding planes, in vertical fractures and joints, and as vug fillings and intercrystalline cements.The mineralization was recognized throughout the Sinnipee Group, and some vertical fractures were mineralized by pyrite, K-silicates, and calcite cements throughout over 15 m of vertical exposure.
Some mineralization occurred within the Prairie du Chien Group and the Cambrian sandstones.Chalcopyrite, quartz, pyrite, and sphalerite were observed in the Prairie du Chien Group, with one notable location in western Outagamie County (Location 46) containing smoky quartz, amethyst, galena and malachite.Pyrite and marcasite nodules were observed in drill cores and cuttings from the Cambrian sandstones, with chalcopyrite, sphalerite, and galena observed locally.The top of the Jordan Sandstone (the uppermost Cambrian unit in northeastern Wisconsin) has also been observed to host zones enriched in pyrite and marcasite [34].
The carbonate units with the least amount of sulfide mineralization were the Silurian dolostones, which represent the youngest strata in the study area.Cores, abundant quarry exposures, and outcrops of the Silurian dolostones indicate that most sulfides are restricted to the lower few meters of the Mayville Formation.However, trace amounts of iron sulfides were present in various forms throughout the Silurian section.The presence of iron-stained vertical joints suggests that this unit, which is regionally karsted, has had much of its sulfide mineralization removed through near-surface weathering during the Cenozoic Era.The Ordovician Sinnipee Group carbonates contained significant quantities of metal sulfides, even at stratigraphic intervals farther away from the SCH in the Ancell Group.Sulfide mineralization was locally abundant along bedding planes, in vertical fractures and joints, and as vug fillings and intercrystalline cements.The mineralization was recognized throughout the Sinnipee Group, and some vertical fractures were mineralized by pyrite, K-silicates, and calcite cements throughout over 15 m of vertical exposure.Some mineralization occurred within the Prairie du Chien Group and the Cambrian sandstones.Chalcopyrite, quartz, pyrite, and sphalerite were observed in the Prairie du Chien Group, with one notable location in western Outagamie County (Location 46) containing smoky quartz, amethyst, galena and malachite.Pyrite and marcasite nodules were observed in drill cores and cuttings from the Cambrian sandstones, with chalcopyrite, sphalerite, and galena observed locally.The top of the Jordan Sandstone (the uppermost Cambrian unit in northeastern Wisconsin) has also been observed to host zones enriched in pyrite and marcasite [34].
The carbonate units with the least amount of sulfide mineralization were the Silurian dolostones, which represent the youngest strata in the study area.Cores, abundant quarry exposures, and outcrops of the Silurian dolostones indicate that most sulfides are restricted to the lower few meters of the Mayville Formation.However, trace amounts of iron sulfides were present in various forms throughout the Silurian section.The presence of iron-stained vertical joints suggests that this unit, which is regionally karsted, has had much of its sulfide mineralization removed through near-surface weathering during the Cenozoic Era.

Petrographic Analysis Using SEM with EDS
Samples of polished rock slabs, polished thin sections, crystal-lined vugs, and well-cuttings chips were analyzed to determine minerals present, to outline mineral paragenesis, and to identify whether or not separate mineral phases were present that could be identified as the source of certain heavy metals present in regional groundwater.
Separate mineral phases containing Fe, Zn, Pb, Cu, and other metals were observed in samples from many locations, especially where parts of the middle and lower Ordovician section are exposed (Figures 7-10).Typically, pyrite and marcasite (FeS 2 ) were the dominant sulfide phases present.Pyrite exhibited several crystal habits, including cubes, octahedrons, and pyritohedrons up to 1 cm in length.Marcasite occured as bladed-and sometimes cockscomb-habits with a greenish or bluish metallic, iridescent coating of unknown origin.Sphalerite (ZnS) was the next most commonly observed sulfide phase, with generally rounded, irregular crystals up to 5 cm.Galena (PbS) was far less common, and was typically found as small cubic crystals a few millimeters in size, but crystals up to 1 cm were also found in the southern part of the study area.Chalcopyrite crystals are 1 to 3 mm in size and display a typical disphenoid crystal form.

Petrographic Analysis Using SEM with EDS
Samples of polished rock slabs, polished thin sections, crystal-lined vugs, and well-cuttings chips were analyzed to determine minerals present, to outline mineral paragenesis, and to identify whether or not separate mineral phases were present that could be identified as the source of certain heavy metals present in regional groundwater.
Separate mineral phases containing Fe, Zn, Pb, Cu, and other metals were observed in samples from many locations, especially where parts of the middle and lower Ordovician section are exposed (Figures 7-10).Typically, pyrite and marcasite (FeS2) were the dominant sulfide phases present.Pyrite exhibited several crystal habits, including cubes, octahedrons, and pyritohedrons up to 1 cm in length.Marcasite occured as bladed-and sometimes cockscomb-habits with a greenish or bluish metallic, iridescent coating of unknown origin.Sphalerite (ZnS) was the next most commonly observed sulfide phase, with generally rounded, irregular crystals up to 5 cm.Galena (PbS) was far less common, and was typically found as small cubic crystals a few millimeters in size, but crystals up to 1 cm were also found in the southern part of the study area.Chalcopyrite crystals are 1 to 3 mm in size and display a typical disphenoid crystal form.Two copper-bearing minerals were observed in the study area.The crystal morphology and EDS spectrum of the most common Cu-bearing phase suggests that it is chalcopyrite (CuFeS2).This has been observed at several locations in Oconto, Winnebago, and Outagamie counties (Figures 3 and 9).Botryoidal and fibrous malachite has been observed in growths up to 3 mm long at one location in western Outagamie County in association with pyrite and chalcopyrite (Figure 9).Malachite has also been observed in association with chalcopyrite south of the study area at the Morris Pit in Green Lake County, Wisconsin [50].Two copper-bearing minerals were observed in the study area.The crystal morphology and EDS spectrum of the most common Cu-bearing phase suggests that it is chalcopyrite (CuFeS 2 ).This has been observed at several locations in Oconto, Winnebago, and Outagamie counties (Figures 3 and 9).Botryoidal and fibrous malachite has been observed in growths up to 3 mm long at one location in western Outagamie County in association with pyrite and chalcopyrite (Figure 9).Malachite has also been observed in association with chalcopyrite south of the study area at the Morris Pit in Green Lake County, Wisconsin [50].The vast majority of pyrite and marcasite analyzed from the region had no detectable peaks for Ni, Co, or As using SEM-EDS techniques.However, detections of Ni, Co, and As in sulfide minerals were made in some cases (Figure 10; Table 1).Detection of a discrete As peak was limited to one specimen, with values of ~2 to 4 weight percent (~0.5 to 1 mole percent), based on EDS analysis.Although this is higher than the As concentration described by Thornburg and Sahai [30] for Asbearing pyrite in which isomorphic substitution of As was interpreted, it is still likely that the As occurs in a similar form here, and not in a discrete mineral phase such as arsenopyrite.
Despite its occurrence in southern parts of Wisconsin [8,14,51,52], millerite (NiS) has not been documented in this study area of northeastern Wisconsin.Additional references to minor amounts of millerite, bravoite, violarite, honessite, cobaltite and smaltite are given for the Upper Mississippi Valley lead-zinc district in southwestern Wisconsin by Heyl et al. [8,52].
Detection of numerous microscopic crystals of sulfide minerals enriched in Ni and Co was made in samples from Brown and Oconto counties (Figure 10; Table 1).These microscopic occurrences are the first examples of Ni-and Co-bearing minerals found in northeastern Wisconsin.Crystals contain widely varying proportions of Ni-Co-Fe, based on SEM-EDS analysis.Figure 11 shows normalized mole percent values for the sulfide minerals containing Ni, Co, and Fe listed in Table 1.Nickel proportions up to 89 mole percent and cobalt proportions up to 24.6 mole percent suggest that the occurrence of these elements is not simply restricted to isomorphous substitution within pyrite, as appears to be the case with arsenic from a few locations in the study area [30].Based upon crystal The vast majority of pyrite and marcasite analyzed from the region had no detectable peaks for Ni, Co, or As using SEM-EDS techniques.However, detections of Ni, Co, and As in sulfide minerals were made in some cases (Figure 10; Table 1).Detection of a discrete As peak was limited to one specimen, with values of ~2 to 4 weight percent (~0.5 to 1 mole percent), based on EDS analysis.Although this is higher than the As concentration described by Thornburg and Sahai [30] for As-bearing pyrite in which isomorphic substitution of As was interpreted, it is still likely that the As occurs in a similar form here, and not in a discrete mineral phase such as arsenopyrite.
Despite its occurrence in southern parts of Wisconsin [8,14,51,52], millerite (NiS) has not been documented in this study area of northeastern Wisconsin.Additional references to minor amounts of millerite, bravoite, violarite, honessite, cobaltite and smaltite are given for the Upper Mississippi Valley lead-zinc district in southwestern Wisconsin by Heyl et al. [8,52].
Detection of numerous microscopic crystals of sulfide minerals enriched in Ni and Co was made in samples from Brown and Oconto counties (Figure 10; Table 1).These microscopic occurrences are the first examples of Ni-and Co-bearing minerals found in northeastern Wisconsin.Crystals contain widely varying proportions of Ni-Co-Fe, based on SEM-EDS analysis.Figure 11 shows normalized mole percent values for the sulfide minerals containing Ni, Co, and Fe listed in Table 1.Nickel proportions up to 89 mole percent and cobalt proportions up to 24.6 mole percent suggest that the occurrence of these elements is not simply restricted to isomorphous substitution within pyrite, as appears to be the case with arsenic from a few locations in the study area [30].Based upon crystal morphology and SEM-EDS composition, bravoite and vaesite are the most likely mineral species present for the Ni and Co-rich crystals.1. Red squares are for samples AP-24 (Locality 83), and green circles are for 07-OCO-MTV-2 (Locality 39).Region of bravoite, along with end members pyrite, vaesite, and cattierite are shown, as described by Kerr [53].
The precise mineralogy of these Ni and Co phases is challenging to determine due to the size of the crystals (tens of micrometers).Crystal morphology varied from indistinct to cubic to cubepyritohedron forms (Figure 10), which is possible in the solid solution series of the pyrite group.In cases where the crystal approaches the NiS2 end-member composition, the mineral vaesite (NiS2) is  1.Red squares are for samples AP-24 (Locality 83), and green circles are for 07-OCO-MTV-2 (Locality 39).Region of bravoite, along with end members pyrite, vaesite, and cattierite are shown, as described by Kerr [53].
The precise mineralogy of these Ni and Co phases is challenging to determine due to the size of the crystals (tens of micrometers).Crystal morphology varied from indistinct to cubic to cube-pyritohedron forms (Figure 10), which is possible in the solid solution series of the pyrite group.In cases where the crystal approaches the NiS 2 end-member composition, the mineral vaesite (NiS 2 ) is the most likely candidate [53].Other crystals with elevated Co and Fe are best described as bravoite ((Ni-Co-Fe)S 2 ).In the absence of X-ray data, strict identification of these minerals is not possible, but their compositions are shown in Figure 11.
In contrast to the above metals, none of the samples examined with the SEM displayed separate mineral phases for arsenic, such as arsenopyrite.While the presence of arsenopyrite cannot be ruled out in these rocks, it seems likely that As is substituting for Fe as a trace element in the iron sulfides, as has been suggested by Thornburg and Sahai [30].For a limited sample set, they concluded that As was present in isomorphous substitution with pyrite/marcasite, which can be released by oxidation and subsequently sorb onto ferric oxyhydroxides.Some fractures and vugs contained a complex microscopic array of minerals in the MVT assemblage (e.g., Figure 9B).In one 2 mm-wide fracture intercepted by a drill core, sphalerite, galena, fluorite, dolomite, marcasite, Sr-barite, and calcite were observed in close association with each other.
With only one exception near the Wisconsin-Michigan border (Locality 33), all galena observed in northeastern Wisconsin falls south an east-west regional fault that is possibly associated with a crustal boundary in the Precambrian basement rocks known as the Spirit Lake Tectonic Zone [25].The reason for this distribution is unknown, but it might be related to the metals concentrations in Precambrian bedrock through which hydrothermal deep groundwater likely passed [54,55].

Paragenetic Sequence of Mineralization
Table 2 lists the diagenetic minerals observed in Paleozoic sedimentary rocks in the region.Although no single diagram can precisely synthesize the full complexity of the mineral paragenesis in the study area, Figure 12 illustrates the general mineralization sequence in northeastern Wisconsin.
Phosphates Apatite/fluorapatite Silicates Quartz, K-feldspar, illite?, amethyst 3 , smoky quartz 3   Halides Fluorite 1 Additional rare MVT and alteration minerals not observed in this study can be found in [40] for one quarry in Winnebago County; 2 Anhydrite was identified in one sample with the assistance of David Tuschel using Raman spectroscopy at HORIBA Scientific in Edison, New Jersey, USA; 3 Amethyst and smoky quartz were observed at Locality 46 in western Outagamie County.
In general, pervasive dolomitization and silicate mineralization are the earliest events.Early quartz appears to have formed at temperatures ď ~50 ˝C, due to the presence of all-liquid fluid inclusion assemblages entrapped within quartz overgrowths.Later quartz and the majority of the dolomite were formed from a Na-Ca-Mg-Cl brine at temperatures of ~80-100 ˝C as temperatures warmed [23].Iron sulfides occurred throughout much of the sequence, but the main stage of MVT mineralization included sphalerite, galena, chalcopyrite, pyrite, marcasite, and Ni-Co-Fe-sulfide minerals (Figures 8-10).The MVT mineralization was intergrown with two episodes of dolomite, with most MVT minerals closely associated with late planar dolomite (Figures 7 and 8).Middle and late-stage dolomite, as well as MVT minerals, appear to have formed by similar brines at temperatures of 80-110 ˝C during what appears to be peak heating.Siderite was not documented in the region, although some dolomite was ferroan.The timing of K-silicates (K-feldspar and possibly illite), relative to other diagenetic phases, is difficult to determine in cases where it occurs as cements within quartz arenites.However, some of it is clearly coeval with abundant iron sulfide mineralization in vertical fractures observed in the northern part of the study area (Figure 8E).
New evidence (Figure 9) suggests that at least some of the sulfate minerals precipitated during the MVT stage of mineralization.Sr-barite, barite, and celestine were observed with fluorite, sphalerite, and pyrite in some samples, likely forming during the waning stages of mineralization.Other late-stage minerals include beef-vein and intergranular gypsum, as well as late-stage calcite.Conditions of calcite mineralization are not well known, but at least some late calcite in Silurian rocks contains all-liquid inclusions, suggesting a lower temperature of entrapment [23].Strontianite was also found in the region as a relatively late mineral, but its precise location in the sequence is not well-constrained.
The order of mineralization shown in Figure 12 should be used with caution.As a regionally derived diagram, local or stratigraphic variations may not be expressed fully.Many minerals are also found as isolated occurrences or with only one other mineral.Nevertheless, the diagram presents the only complete synopsis of the paragenetic sequence of mineralization in the study area.New evidence (Figure 9) suggests that at least some of the sulfate minerals precipitated during the MVT stage of mineralization.Sr-barite, barite, and celestine were observed with fluorite, sphalerite, and pyrite in some samples, likely forming during the waning stages of mineralization.Other late-stage minerals include beef-vein and intergranular gypsum, as well as late-stage calcite.Conditions of calcite mineralization are not well known, but at least some late calcite in Silurian rocks contains all-liquid inclusions, suggesting a lower temperature of entrapment [23].Strontianite was also found in the region as a relatively late mineral, but its precise location in the sequence is not wellconstrained.
The order of mineralization shown in Figure 12 should be used with caution.As a regionally derived diagram, local or stratigraphic variations may not be expressed fully.Many minerals are also found as isolated occurrences or with only one other mineral.Nevertheless, the diagram presents the only complete synopsis of the paragenetic sequence of mineralization in the study area.

ICP-AES Analytical Results and Discussion
Analytical results for all 310 samples are presented in Online Supplementary Materials Tables S1 and S2.The main purpose for analyzing the whole-rock chemistry of these strata was to look for regional and stratigraphic trends in the character of the mineralization, with a focus on heavy metals.Table 3 presents descriptive statistics for selected metals in all samples in the study that were designated as "Matrix" and "Enriched."Table 4 illustrates the typical (median) whole-rock concentrations for selected metals in representative "matrix" samples of bedrock from the region.Averages were not calculated because the data were heavily skewed right, as is common when measuring concentrations.Metals concentrations varied widely, from below detection limits to far above the median concentrations for the host rocks, which was consistent with observations of hand specimens and thin sections.

ICP-AES Analytical Results and Discussion
Analytical results for all 310 samples are presented in Online Supplementary Materials Tables S1  and S2.The main purpose for analyzing the whole-rock chemistry of these strata was to look for regional and stratigraphic trends in the character of the mineralization, with a focus on heavy metals.Table 3 presents descriptive statistics for selected metals in all samples in the study that were designated as "Matrix" and "Enriched."Table 4 illustrates the typical (median) whole-rock concentrations for selected metals in representative "matrix" samples of bedrock from the region.Averages were not calculated because the data were heavily skewed right, as is common when measuring concentrations.Metals concentrations varied widely, from below detection limits to far above the median concentrations for the host rocks, which was consistent with observations of hand specimens and thin sections.Based on field observations and elemental analysis, three specific stratigraphic horizons showed the most significant mineralization.These were the top of the Ancell Group, the Ordovician-Silurian contact, and, to a lesser extent, the Sinnipee Group dolostone.All three of these zones typically contained elevated metals concentrations, especially in "enriched" samples from the SCH, or along mineralized bedding planes, faults, or fractures.Figure 13 shows box-plots (on the log 10 scale) for selected metals from Table 3. Elevated metals occurred in several Cambrian and Ordovician units, but the Silurian dolostones generally exhibited lower concentrations.Additional box-plots for these and other metals are presented in Supplementary Materials Figure S3.It is important to note that all stratigraphic units show some potential for locally elevated metals concentrations, including the Prairie du Chien Group and the Cambrian sandstones, which serve as an important regional aquifer system along with the Ancell Group [56,57].
Table 5 presents Spearman's ρ correlations, estimated for a subset of log 10 -transformed metals concentrations relevant to the study.Figure 14 shows scatter plots for selected metals concentrations, again on the log 10 scale.While most pairs of elements in Table 5 show moderate correlations, the highest values occurred for Ni vs. Fe, Ni vs. Co, Pb vs. Fe, Co vs. Fe, and Cr vs. V.This result is consistent with the SEM-EDS observations of As-, Ni-, and Co-bearing sulfide minerals, suggesting that they are likely related to the same mineralization event.  2 and 4).Metals concentrations are highest in the Maquoketa-Neda and the Ancell Group, while the Silurian dolostone units tend to have the lowest concentrations for most metals.Colors correspond to those used in Figures 2 and 4. Table 5.Estimated Spearman's correlations amongst log10-transformed variables for selected elements.For all but one of the estimates (Zn vs. Fe of enriched samples), p-values were smaller than 0.02, with most < 0.0001.An expanded analysis showing Pearson's r, Spearman's , and all p-values is presented in Supplementary Materials Tables S2C and S2D.  for 309 samples (single Precambrian point excluded), separated by major stratigraphic units (Figures 2 and 4).Metals concentrations are highest in the Maquoketa-Neda and the Ancell Group, while the Silurian dolostone units tend to have the lowest concentrations for most metals.Colors correspond to those used in Figures 2 and 4.
Table 5.Estimated Spearman's ρ correlations amongst log 10 -transformed variables for selected elements.For all but one of the estimates (Zn vs. Fe of enriched samples), p-values were smaller than 0.02, with most < 0.0001.An expanded analysis showing Pearson's r, Spearman's ρ , and all p-values is presented in Supplementary Materials Tables S2C and S2D.Strong correlations for Ni vs. Fe ( ρ = 0.744), Ni vs. Co ( ρ = 0.718), and Co vs. Fe ( ρ = 0.612) suggest that these three elements are related to the same mineralization event (Figure 14).Ni and Co are well correlated in all sample subsets (Table 5), suggesting that these elements are more closely related.They are likely principally occurring as isomorphic substitutions in pyrite, or as separate Ni-Co-enriched sulfide minerals.This interpretation is consistent with the discovery of discrete, but microscopic, Ni-Co-Fe sulfide minerals vaesite and bravoite in the study area (Figures 10 and 11).Correlations between Ni, Co, and Fe are higher than correlations for Ni and As ( ρ = 0.437) or Co and As ( ρ = 0.588).This is consistent with mineralogical observations in which As is typically in separate sulfide mineral phases from Ni and Co (Figure 10).This is an important observation suggesting that water quality evaluations for wells in the study area should probably include analyses for Ni and Co more regularly, instead of only As.
Smaller correlations between Zn and Fe and between Cu and Fe are most likely due to the fact that those metals also precipitate as separate sulfide mineral phases, such as sphalerite and chalcopyrite, and are independent of the presence of pyrite or marcasite in any particular sample.The same would be expected between Pb and Fe, but there is a moderate correlation between these two elements ( ρ = 0.676 for enriched samples).Because galena is far less commonly observed in the study area than sphalerite, chalcopyrite, pyrite, or marcasite, the correlation between Pb and Fe may represent isomorphic substitution of Pb in pyrite or marcasite rather than the discrete presence of galena.
The correlation between K and B ( ρ = 0.694 for enriched samples) was unexpected, and is not well understood because boron could not be detected using SEM-EDS techniques.However, it is possible that B is substituting in the structure of authigenic feldspar or illite (e.g., Figure 8E), both of which are present in the region [23].These minerals could be the source of elevated boron levels reported for the Cambrian-Ordovician aquifer across the study area [58].
Correlations between Ni and As estimated separately by stratigraphic units (Supplementary Materials Table S2D) showed samples from the Maquoketa-Neda, the Ancell Group, and the Cambrian Sandstones were well correlated ( ρ = 0.762, 0.722, and 0.797, respectively).The plot of Ni vs.As (Figure 14C) showed that some samples contained elevated Ni concentrations relative to As.These are best represented by some Ancell Group samples from three Oconto County quarries (Duame Quarry, Chase Quarry, and Montevideo Quarry) and one Sinnipee Group sample from the Vulcan Quarry in Winnebago County.Elevated nickel relative to arsenic is consistent with the existence of Ni-bearing minerals such as vaesite and bravoite that were observed in Oconto County (Figures 10  and 11).Ni and As in the Silurian units appeared moderately correlated ( ρ = 0.501), while in the remaining units (Sinnipee and Prairie du Chien groups), they did not correlate as strongly ( ρ = 0.381 and ´0.158, respectively).
Relatively strong correlations between Ni and Co persisted even for some finer stratigraphic divisions, such as individual formations within the Sinnipee Group.Specifically, for the Galena and Platteville formations, ρ was 0.694 and 0.657, respectively (Figure 14D).As and Co were less strongly correlated for the Galena and Platteville formations, yielding ρ values of 0.572 and 0.385, respectively.Again, this is consistent with the observation of Ni and Co-bearing mineral phases as described above.
A moderate correlation ( ρ = 0.541) existed for V and Cr in all samples over all stratigraphic intervals (Figure 14E), but this correlation improved notably for the enriched sample subset ( ρ = 0.769) (Table 5).However, for a subset of samples from the mineralized zone near the Ordovician-Silurian contact in central Brown County (17 samples along a small portion of the Niagara Escarpment), a particularly high correlation between V and Cr was observed ( ρ = 0.940) (Table 5; Figure 14F).In these samples, a slight enrichment of V, Sc, Tl, and La was also observed.The precise cause is unknown, but it might relate to the presence of oolitic hematite of the Neda Formation along this part of the escarpment.Although there is a general affinity of V for iron-bearing rocks, most of the samples with high V concentrations are dolostone samples from near the contact.Another stratigraphic difference noted is that iron sulfide-rich samples from the Maquoketa-Neda contact and the Silurian dolostones just above this contact are depleted in Ni, relative to samples from the Ancell and Sinnipee groups (Figure 14A).Although depleted in Ni, they are slightly enriched in V, Sc, Tl, and La relative to the Ancell and Sinnipee Group samples (Table S1).These differences are likely the result of diagenesis by hydrothermal groundwater.The Maquoketa Shale (Figures 2 and 4) is a regional confining unit that likely isolated the Silurian dolostone units above from the Ordovician sandstones and carbonates below during mineralization.This local anomaly occurs near the only exposure of the Neda Ironstone, which suggests that mineralizing groundwater may have scavenged certain elements locally from the Neda Ironstone instead of being supplied by groundwater migrating through other strata.
Cu concentrations were highest in Oconto County (Ancell Group) and near the Ordovician-Silurian contact in central Brown County (Maquoketa-Neda units).This is consistent with the identification of chalcopyrite in field specimens from Oconto County Prairie du Chien Group samples.
Tuttle et al. [59] demonstrated that redistribution of minor and trace elements can occur in rocks with sulfide mineralization due to weathering processes and groundwater flow.While it is likely that remobilization of trace elements may have occurred in the study area, a systematic assessment of this was not performed.Extensive Late Pleistocene glaciation has removed bedrock-derived soil horizons, and most soils have developed on transported materials.In most cases, bedrock samples were obtained from relatively fresh host rocks in quarries, drill cores/cuttings, and excavated road cuts to keep sample weathering effects to a minimum.The most likely bedrock units to have been affected by weathering are the Silurian dolostones (Figure 2), which show extensive karst development [25] and relatively low trace element concentrations (Table 4).

Discussion
The origin of the SCH is not well understood, but the fact that the most prevalent sulfide mineralization in eastern Wisconsin is located at the top of the St.Peter Sandstone is well recognized.Enrichment of this zone might be the result of sulfate reduction due to organic carbon near the contact between the Ancell and Sinnipee groups, but this is absent at many locations in the study area.Alternatively, its location at the top of the St.Peter Sandstone is consistent with the possibility of a sour gas cap during mineralization.The presence of sour (H 2 S-rich) natural gas migrating with deep groundwater discharging out of the Michigan basin (e.g., [23]) could explain the observed regional and stratigraphic spatial distribution of sulfides at the top of this unit.Precipitation of sulfide minerals would be expected to occur at the site of gas exsolution from brines migrating updip or from a metals-bearing fluid mixing with a preexisting sour gas cap.Further analysis is beyond the scope of this study, but work on the sulfur-isotopic characteristics of the sulfide mineralization is ongoing and may help to resolve this.
Previous work on whole-rock chemistry in the region was limited to a small number of samples because those studies were mostly focused on detailed water quality research in a few locations (Figure 3).Results of the whole-rock analyses from this study were compared with reported metals concentrations for the Ancell Group of Winnebago and Outagamie Counties and other localities in northeastern Wisconsin [29,38,40,60,61].This procedure allowed generalized comparisons to be made between rocks in the southern part of the study area to those in the northern portion.Additional unpublished data were obtained from Dave Johnson [45] for 69 selected well cuttings from Cambrian rocks in Brown County.In those samples, nickel, arsenic, and cobalt concentrations were generally low, but values of up to 22.5 ppm for Ni, 105 ppm for As, and 93 ppm for Co were detected in different samples, with six samples containing cobalt above 25 ppm [45].This likely indicates the presence of As, Ni, or Co-bearing minerals at those locations.
Arsenic concentrations relative to iron in this study were compared to those from the studies mentioned above.In general, whole-rock arsenic concentrations varied between 0 ppm and 500 ppm for all studies.However, for three samples reported in previous studies from Winnebago and Outagamie counties, As concentrations ranged between 500 and 743 ppm.While a few samples reported by Gotkowitz et al. [29,40] had higher maximum Co concentrations than those in this study, a comparison of As vs. Co data also shows general overlap between the studies.
There were some differences noted between the studies.First, Cr concentrations in 308 of 310 samples in this study were below 40 ppm.The two exceptions were for cuttings samples in a municipal well in central Brown County with Cr concentrations of 66 and 152 ppm.In contrast, 18 of 21 samples analyzed in previous studies [29,40,61] were reported to have Cr concentrations between 100 and 764 ppm, regardless of location, with most concentrations above 200 ppm.Because the locations of some of these data overlap those in our study, these high Cr values seem anomalous and may be due to sample contamination.One value reported by Gotkowitz et al. [29,40] for an unmineralized sample of St. Peter Sandstone had a Cr concentration above 500 ppm, even though the sample had no sulfur and less than 1% Fe.Because of this anomalous behavior of Cr, we did not attempt to draw any conclusions about the difference in Cr between our studies and previous studies as they relate to rock geochemistry.It is possible that samples in the previous studies were submitted to a laboratory that crushed the samples with high Cr-steel, which is an alloy known for its toughness and which would be a suitable material for crushing rocks.In this study, contamination of the two higher Cr cuttings samples from this study by drilling equipment cannot be ruled out.
A comparison of Ni, As, and Fe shows general agreement between most samples analyzed in our study and those of Gotkowitz et al. [29,40] and Pelczar [38] (Figure 14A,B).In most samples, As is enriched relative to Ni.However, in some samples, Ni appears elevated relative to arsenic (Figure 14C), which is consistent with the presence of Ni-bearing minerals such as vaesite and bravoite.
It is unknown whether northeastern Wisconsin holds economic concentrations of MVT minerals in Paleozoic bedrock.A lack of significant local enrichment, compared to the former ore district in southwestern Wisconsin, could be due to a number of factors.It is likely that a diffuse hydrothermal fluid flow system was responsible for mineralization in the study area.This is consistent with pervasively dolomitized carbonate rocks in which 100 percent of the host-rock limestone has been replaced, along with a regionally prevalent SCH.However, occasional instances of acidic, metals-rich groundwater in water wells reported in the region suggest alteration of enriched sulfide deposits might occur in local areas.Such instances are observed mainly in areas near the Ancell Group outcrop belt (Figure 3) where oxygenated water interacts with the SCH or other sulfides due to natural conditions or aquifer drawdown [29,30,32,34,40].Farther east (i.e., down dip) from the outcrop belt (Figure 4), oxidizing conditions are less likely to occur, so limited information about the existence of mineralized zones can be inferred from groundwater chemistry alone.Extensive glacial sedimentary cover, poor knowledge of geologic structures such as folds and faults, and geochemical conditions favoring the stability of sulfide minerals, make detection of larger concentrations of sulfide minerals unfavorable in much of the region.Structural style, as well as hydrothermal fluid composition and intensity may also be factors.
It is useful to consider the potential impact of this mineralization on groundwater quality in the region outside of Winnebago and Outagamie counties where special well casing requirements are established.While the upper Ancell Group is likely to be the most significant source of metals due to the abundance of sulfides and the fact that porous media flow dominates, water quality in all units below the Silurian dolostones should be studied more carefully in the future.Because the distribution of sulfide minerals is highly heterogeneous, their locations and intensities along fractures, faults, bedding planes, and in vugs are difficult to predict locally.Fracture and bedding plane-dominated flow are likely to amplify the interaction between groundwater and sulfide minerals under the right redox conditions.In addition, certain stratigraphic zones of dolostone with high intercrystalline porosity, such as one about 6 m (20 feet) above the base of the Platteville Formation, also contain abundant Fe-sulfides locally (Figure 7A).For these reasons, the average concentrations of metals within the specific units are probably not so important.Rather, the complex distribution of sulfide minerals, redox conditions, and fluid flow pathways likely has a much larger influence on water quality in the carbonate units.raw analytical results for the whole-rock ICP-AES analyses of samples from 2007 to 2008.Table S2B: Full raw analytical results for the whole-rock ICP-AES analyses of samples from 2013.Table S2C: Expanded table showing estimated Pearson's r and Spearman's ρ correlations amongst log 10 -transformed variables for selected elements in all samples.Table S2D: Table showing estimated Pearson's r and Spearman's ρ correlations amongst log 10 -transformed Ni vs. As, listed by stratigraphic interval (Precambrian excluded).Figure S3: Box-plots for selected metals from both the full set of samples and the "matrix" subset of samples.

Figure 1 .
Figure 1.Bedrock geologic map of the western Great Lakes region of North America showing the locations of major geologic structures near the states of Wisconsin and Michigan.Geologic rock systems are as follows: pЄ = Precambrian, Є = Cambrian, O = Ordovician, S = Silurian, D = Devonian, M = Mississippian, Pn = Pennsylvanian, J = Jurassic.Red dashed lines denote the Illinois and Michigan structural basins.Black lines with double arrows indicate the positions of structural arches."UMV" indicates the location of the Upper Mississippi Valley Lead-Zinc Ore District as defined by Heyl et al. [8].Blue arrows indicate inferred flow paths of hydrothermal fluids originating in the Michigan basin.Base map modified after[12,13].

Figure 1 .
Figure 1.Bedrock geologic map of the western Great Lakes region of North America showing the locations of major geologic structures near the states of Wisconsin and Michigan.Geologic rock systems are as follows: p

Figure 1 .
Figure 1.Bedrock geologic map of the western Great Lakes region of North America showing the locations of major geologic structures near the states of Wisconsin and Michigan.Geologic rock systems are as follows: pЄ = Precambrian, Є = Cambrian, O = Ordovician, S = Silurian, D = Devonian, M = Mississippian, Pn = Pennsylvanian, J = Jurassic.Red dashed lines denote the Illinois and Michigan structural basins.Black lines with double arrows indicate the positions of structural arches."UMV" indicates the location of the Upper Mississippi Valley Lead-Zinc Ore District as defined by Heyl et al. [8].Blue arrows indicate inferred flow paths of hydrothermal fluids originating in the Michigan basin.Base map modified after[12,13].

Figure 1 .
Figure 1.Bedrock geologic map of the western Great Lakes region of North America showing the locations of major geologic structures near the states of Wisconsin and Michigan.Geologic rock systems are as follows: pЄ = Precambrian, Є = Cambrian, O = Ordovician, S = Silurian, D = Devonian, M = Mississippian, Pn = Pennsylvanian, J = Jurassic.Red dashed lines denote the Illinois and Michigan structural basins.Black lines with double arrows indicate the positions of structural arches."UMV" indicates the location of the Upper Mississippi Valley Lead-Zinc Ore District as defined by Heyl et al. [8].Blue arrows indicate inferred flow paths of hydrothermal fluids originating in the Michigan basin.Base map modified after[12,13].

Figure 2 .
Figure 2. Simplified stratigraphic column for Paleozoic rocks in northeastern Wisconsin, USA [13].Pleistocene sediments and younger Silurian and Devonian strata for southeastern Wisconsin are not shown here for simplicity, but can be found elsewhere[27,28].

Figure 2 .
Figure 2. Simplified stratigraphic column for Paleozoic rocks in northeastern Wisconsin, USA [13].Pleistocene sediments and younger Silurian and Devonian strata for southeastern Wisconsin are not shown here for simplicity, but can be found elsewhere[27,28].

Figure 3 .
Figure 3. Map with sample locations (outcrops, quarries, drill cores, well cuttings) throughout the study area in northeastern Wisconsin.Thin lines and names indicate county boundaries.Inset map shows the study area location in the state of Wisconsin, USA.Red areas with letters indicate groundwater studies performed in the region (A is reference for Gotkowitz et al.[29], B is Thornburg and Sahai[30], C is Burkel and Stoll[31], and D is Schreiber et al.[32]).Additional references for region-wide studies include Riewe et al.[33] and Johnson and Riewe[34], as well as a cancer study by Knobeloch et al.[35] that focused on Winnebago and Outagamie counties.Base map after Mudrey et al.[36] shows the following units: pЄ = Precambrian, Є = Cambrian sandstones, Opc = Ordovician Prairie du Chien Group, Oa = Ordovician Ancell Group, Os = Ordovician Sinnipee Group, Om = Ordovician Maquoketa-Neda formations, and Su = Silurian dolostones.

Figure 3 .
Figure 3. Map with sample locations (outcrops, quarries, drill cores, well cuttings) throughout the study area in northeastern Wisconsin.Thin lines and names indicate county boundaries.Inset map shows the study area location in the state of Wisconsin, USA.Red areas with letters indicate groundwater studies performed in the region (A is reference for Gotkowitz et al. [29], B is Thornburg and Sahai [30], C is Burkel and Stoll[31], and D is Schreiber et al.[32]).Additional references for region-wide studies include Riewe et al.[33] and Johnson and Riewe[34], as well as a cancer study by Knobeloch et al.[35] that focused on Winnebago and Outagamie counties.Base map after Mudrey et al.[36] shows the following units: p

1 .
Bedrock geologic map of the western Great Lakes region of North America showing the ns of major geologic structures near the states of Wisconsin and Michigan.Geologic rock s are as follows: pЄ = Precambrian, Є = Cambrian, O = Ordovician, S = Silurian, D = Devonian, ississippian, Pn = Pennsylvanian, J = Jurassic.Red dashed lines denote the Illinois and Michigan ural basins.Black lines with double arrows indicate the positions of structural arches."UMV" tes the location of the Upper Mississippi Valley Lead-Zinc Ore District as defined by Heyl et al. ue arrows indicate inferred flow paths of hydrothermal fluids originating in the Michigan basin.apmodified after[12,13].

Figure 1 .
Figure 1.Bedrock geologic map of the western Great Lakes region of North America showing the locations of major geologic structures near the states of Wisconsin and Michigan.Geologic rock systems are as follows: pЄ = Precambrian, Є = Cambrian, O = Ordovician, S = Silurian, D = Devonian, M = Mississippian, Pn = Pennsylvanian, J = Jurassic.Red dashed lines denote the Illinois and Michigan structural basins.Black lines with double arrows indicate the positions of structural arches."UMV" indicates the location of the Upper Mississippi Valley Lead-Zinc Ore District as defined by Heyl et al. [8].Blue arrows indicate inferred flow paths of hydrothermal fluids originating in the Michigan basin.Base map modified after[12,13].

Figure 4 .
Figure 4. West to east geologic cross-section through portions of Shawano, Brown, and Kewaunee counties near Green Bay, Wisconsin.Paleozoic strata typically dip southeastward between about 5 and 7.5 m/km toward the ancestral Michigan basin.Vertical lines with labels indicate well codes for selected existing and backup wells in the Green Bay region of Brown County.Vertical exaggeration is 45 times.Modified after [13].
, an older dataset from 1976 to 1979, also showed elevated dissolved Ni, Co, and As in the study area.Most of the wells with exceedances draw water from the Cambrian-Ordovician sandstone aquifers, especially the Ancell Group.Even municipal wells in Brown County 240-320 m deep have been impacted by elevated As, Co, and Ni

Figure 4 .
Figure 4. West to east geologic cross-section through portions of Shawano, Brown, and Kewaunee counties near Green Bay, Wisconsin.Paleozoic strata typically dip southeastward between about 5 and 7.5 m/km toward the ancestral Michigan basin.Vertical lines with labels indicate well codes for selected existing and backup wells in the Green Bay region of Brown County.Vertical exaggeration is 45 times.Modified after [13].

Figure 5 .
Figure 5. Exposures of the Sulfide Cement Horizon (SCH) in the upper portion of the St.Peter Sandstone (Ancell Group): (A) Discolored floor of a dolostone quarry showing the chemical alteration present at the top of the SCH at Location 34 in Oconto County.The quarry wall is the overlying Platteville Formation of the Sinnipee Group (human for scale).(B) A typical nodule of pyrite and/or marcasite-cemented quartz sandstone from the SCH at Locality 84 in Outagamie County (96-OUT-1), about 9 cm wide here.(C) Pink coating on Sulfide Cement Horizon at Locality 82 in Winnebago County (hammer for scale).

Figure 5 .
Figure 5. Exposures of the Sulfide Cement Horizon (SCH) in the upper portion of the St.Peter Sandstone (Ancell Group): (A) Discolored floor of a dolostone quarry showing the chemical alteration present at the top of the SCH at Location 34 in Oconto County.The quarry wall is the overlying Platteville Formation of the Sinnipee Group (human for scale).(B) A typical nodule of pyrite and/or marcasite-cemented quartz sandstone from the SCH at Locality 84 in Outagamie County (96-OUT-1), about 9 cm wide here.(C) Pink coating on Sulfide Cement Horizon at Locality 82 in Winnebago County (hammer for scale).

Figure 6 .
Figure 6.View of an approximately 70 cm tall section present in a bedrock sump trench at the Scott Construction Ulmen Quarry in central Brown County, Wisconsin (Location 13).Sulfides occur at the contact between the Silurian Mayville Formation (top 30 cm) and the Ordovician Neda-Maquoketa formations (lower 40 cm) and are responsible for the Fe-oxide staining.

Figure 6 .
Figure 6.View of an approximately 70 cm tall section present in a bedrock sump trench at the Scott Construction Ulmen Quarry in central Brown County, Wisconsin (Location 13).Sulfides occur at the contact between the Silurian Mayville Formation (top 30 cm) and the Ordovician Neda-Maquoketa formations (lower 40 cm) and are responsible for the Fe-oxide staining.

Figure 7 .
Figure 7. Photographs of selected MVT minerals in Ordovician carbonate strata from northeastern Wisconsin.(A) Fe-sulfide cemented dolostone breccia at Locality 48.View is approximately 12 cm across.(B) Dolostone hosted galena and sphalerite with hydrothermal dolomite from Locality 82 in Winnebago County.Vug is 2 cm across.(C) A small strike-slip fault with several centimeter-wide filling of iron sulfides) at Locality 59 in Winnebago County (note hammer for scale).

Figure 7 .
Figure 7. Photographs of selected MVT minerals in Ordovician carbonate strata from northeastern Wisconsin.(A) Fe-sulfide cemented dolostone breccia at Locality 48.View is approximately 12 cm across.(B) Dolostone hosted galena and sphalerite with hydrothermal dolomite from Locality 82 in Winnebago County.Vug is 2 cm across.(C) A small strike-slip fault with several centimeter-wide filling of iron sulfides) at Locality 59 in Winnebago County (note hammer for scale).

Figure 8 .
Figure 8. SEM-BSE images of polished sections and vug-filling cements illustrating the character of Mississippi Valley-type mineralization in northeastern Wisconsin: (A) Octahedral pyrite is intergrown with dolomite in vug from Locality 59 in Winnebago County (96-WIN-2).Left and right edges are flat, polished surfaces.(B) Bladed marcasite (bright) on dolomite (gray) from Locality 48 in Shawano County (07-SH-MCK-5).(C) Fine bladed marcasite and cubic pyrite intergrown with dolomite from Locality 48 (07-SH-Mck-0) in Shawano County.(D) Polished section showing dolomite (gray) and later pyrite (bright) filling pores and replacing parts of the dolomite at Locality 3 in Brown County (07-BN-LQ-1).(E).Fracture-filling cements from a quartz sandstone layer in the lower Platteville Formation at Locality 29 in Marinette County (Sample 07-MT-FP-2).Round gray grains are detrital quartz sand.Bright white areas in fracture are pyrite cement, whereas medium gray cements are calcite.Gray spotted areas coprecipitated with pyrite are an unknown authigenic K-Al-Silicate mineral cement, possibly illite or K-feldspar.(F) Late dolomite intergrown with chalcopyrite and sphalerite and associated with quartz from Locality 46 (07-NL-45E) in Outagamie County.

Figure 8 .
Figure 8. SEM-BSE images of polished sections and vug-filling cements illustrating the character of Mississippi Valley-type mineralization in northeastern Wisconsin: (A) Octahedral pyrite is intergrown with dolomite in vug from Locality 59 in Winnebago County (96-WIN-2).Left and right edges are flat, polished surfaces.(B) Bladed marcasite (bright) on dolomite (gray) from Locality 48 in Shawano County (07-SH-MCK-5).(C) Fine bladed marcasite and cubic pyrite intergrown with dolomite from Locality 48 (07-SH-Mck-0) in Shawano County.(D) Polished section showing dolomite (gray) and later pyrite (bright) filling pores and replacing parts of the dolomite at Locality 3 in Brown County (07-BN-LQ-1).(E) Fracture-filling cements from a quartz sandstone layer in the lower Platteville Formation at Locality 29 in Marinette County (Sample 07-MT-FP-2).Round gray grains are detrital quartz sand.Bright white areas in fracture are pyrite cement, whereas medium gray cements are calcite.Gray spotted areas coprecipitated with pyrite are an unknown authigenic K-Al-Silicate mineral cement, possibly illite or K-feldspar.(F) Late dolomite intergrown with chalcopyrite and sphalerite and associated with quartz from Locality 46 (07-NL-45E) in Outagamie County.
-EDS composition, bravoite and vaesite are the most likely mineral species present for the Ni and Co-rich crystals.

Figure 11 .
Figure 11.Ternary plot of mole percent normalized compositions indicating the relative proportions of nickel, cobalt, and iron for various members of the pyrite group.Data are for samples listed in Table1.Red squares are for samples AP-24 (Locality 83), and green circles are for 07-OCO-MTV-2 (Locality 39).Region of bravoite, along with end members pyrite, vaesite, and cattierite are shown, as described by Kerr[53].

Figure 11 .
Figure 11.Ternary plot of mole percent normalized compositions indicating the relative proportions of nickel, cobalt, and iron for various members of the pyrite group.Data are for samples listed in Table1.Red squares are for samples AP-24 (Locality 83), and green circles are for 07-OCO-MTV-2 (Locality 39).Region of bravoite, along with end members pyrite, vaesite, and cattierite are shown, as described by Kerr[53].

Figure 12 .
Figure 12.Generalized paragenetic sequence for Paleozoic rocks in northeastern Wisconsin.

Figure 12 .
Figure 12.Generalized paragenetic sequence for Paleozoic rocks in northeastern Wisconsin.

Figure 13 .
Figure13.Box-plots on the log10 scale for selected metals (As, Co, Ni, and Cu) for 309 samples (single Precambrian point excluded), separated by major stratigraphic units (Figures2 and 4).Metals concentrations are highest in the Maquoketa-Neda and the Ancell Group, while the Silurian dolostone units tend to have the lowest concentrations for most metals.Colors correspond to those used in Figures2 and 4.

Figure 13 .
Figure 13.Box-plots on the log 10 scale for selected metals ((A) As, (B) Co, (C) Ni, and (D) Cu) for 309 samples (single Precambrian point excluded), separated by major stratigraphic units (Figures2 and 4).Metals concentrations are highest in the Maquoketa-Neda and the Ancell Group, while the Silurian dolostone units tend to have the lowest concentrations for most metals.Colors correspond to those used in Figures2 and 4.

1
This row represents a subset of 17 enriched and matrix samples from Kittell Falls, Fonferek's Glen, and Ulmen Quarry along a portion of the Niagara Escarpment in central Brown County.Outcrops occur along the Ordovician-Silurian Contact.2The p-value for the correlation of Zn vs. Fe for enriched samples was 0.324.

Figure 14 .Figure 14 .
Figure 14.Scatter plots (log10 axes) of selected metals from whole-rock ICP-AES analyses on rock specimens from northeastern Wisconsin.Legends indicate the particular stratigraphic unit sampled.(A) Nickel (Ni) vs. iron (Fe) for all samples, separated by stratigraphic unit.(B) Nickel (Ni) vs. cobaltFigure 14.Scatter plots (log 10 axes) of selected metals from whole-rock ICP-AES analyses on rock specimens from northeastern Wisconsin.Legends indicate the particular stratigraphic unit sampled.(A) Nickel (Ni) vs. iron (Fe) for all samples, separated by stratigraphic unit; (B) Nickel (Ni) vs. cobalt (Co) for all samples, separated by stratigraphic unit; (C) Nickel (Ni) vs. arsenic (As) for all samples, separated by stratigraphic unit.Several samples exhibited elevated Ni, consistent with the presence of observed Ni-bearing minerals (see text); (D) Nickel (Ni) vs. cobalt (Co) for the two formations of the Sinnipee Group in northeastern Wisconsin (Platteville and Galena); (E) Vanadium (V) vs. chromium (Cr) for all samples in the study, separated by stratigraphic unit.Note the V-enriched samples (red squares and blue diamond in upper right area of plot) from central Brown County; (F) Vanadium (V) vs. chromium (Cr) for all samples from three localities along a portion of the Niagara Escarpment in central Brown County.

Table 1 .
Concentrations of elements in selected crystals of Ni-Co-Fe-As-sulfide minerals from Brown and Oconto Counties.Values are reported in weight % generated from the EDAX Software.Sample

Table 1 .
Concentrations of elements in selected crystals of Ni-Co-Fe-As-sulfide minerals from Brown and Oconto Counties.Values are reported in weight % generated from the EDAX Software.Sample AP-24 refers to the Village of Wrightstown Test Well in Brown County (Locality 83).Designations 1a, 1b, etc. refer to analysis points broadly from the core to the rim of crystal 1. Sample 07-OCO-MTV-2 refers to the Montevideo quarry in Oconto County (Locality 39).Dashes indicate no detection by the software.-24 refers to the Village of Wrightstown Test Well in Brown County (Locality 83).Designations 1a, 1b, etc. refer to analysis points broadly from the core to the rim of crystal 1. Sample 07-OCO-MTV-2 refers to the Montevideo quarry in Oconto County (Locality 39).Dashes indicate no detection by the software. AP

Table 2 .
Complete list of Mississippi Valley-type and associated diagenetic minerals observed in Paleozoic rocks of northeastern Wisconsin.

Table 3 .
Minimum, minimum detected, median, maximum, and the interquartile range of concentrations for iron and trace metals in Paleozoic whole-rock samples from northeastern Wisconsin.For Fe, the unit is percentage, while all other elements are measured in mg/kg (ppm).

Table 3 .
Minimum, minimum detected, median, maximum, and the interquartile range of concentrations for iron and trace metals in Paleozoic whole-rock samples from northeastern Wisconsin.For Fe, the unit is percentage, while all other elements are measured in mg/kg (ppm).Reported values are for all samples in the study (matrix and enriched, n = 310).ND is a value below the detection limit for that respective element.Some analyses from early in the study were below the detection limits for many metals due to the analytical methods available.

Table 4 .
Median metal concentrations (interquartile range in parentheses) for matrix samples from Paleozoic rock strata of northeastern Wisconsin (sample sizes in parentheses next to strata names).For Fe, the unit is percentage, while all other elements are measured in mg/kg (ppm).A Precambrian data point was excluded from this table due to having only one observation in that geologic unit.