Dolomitization of Paleozoic Successions, Huron Domain of Southern Ontario, Canada: Fluid Flow and Dolomite Evolution

: Integrated petrographic, isotopic, ﬂuid inclusion microthermometry, and geochemical analyses of Paleozoic carbonate successions from multiple boreholes within the Huron Domain, southern Ontario were conducted to characterize the diagenetic history and ﬂuid composition, on a regional scale, and evaluate the nature and origin of dolomitized beds. Multiple related diagenesis during the regional tectonic evolution of the Michigan Basin.


Introduction
As part of an ongoing investigation into the occurrence of strata-bound dolomite layers within sedimentary, Paleozoic-age formations of southern Ontario an initial study using core acquired from the site of a formerly proposed deep geological repository (Bruce nuclear site) was conducted [1]. This initial study provided a site-specific analog to examine fluid migration and rock-formation barrier integrity at geological timescales. Results showed that hydrothermal alterations were minimal, burial-related, and lacked isotopic signatures associated with fault-controlled diagenesis, as observed in other regions nature and origin of dolomitized beds within Paleozoic carbonates of the Huron Domain with the intention on deciphering the nature of dolomitization processes, evolution of diagenetic fluids, the nature of paleohydrogeological system(s) and whether they are connected on not, the source of heat for hydrothermal fluids (if any) and the driving force for diagenetic fluid(s) migration, in contrast to previous studies in southwestern Ontario.

Regional Geologic, Tectonics, and Stratigraphic Framework
The study area is located on the northeastern margin of the Michigan Basin (Figures 1 and 2). It represents a portion of the northwestern flank of the Algonquin Arch that consists of Paleozoic sediments overlying the subsurface basement [34]. The thickness of the Paleozoic sequence is estimated to range from a maximum of 4800 m at the center of the Michigan Basin to 850 m on the flank of the Algonquin Arch.  [35,36]); (B) Tectonic history of the Phanerozoic successions in southern Ontario (modified from [37]). Southern Ontario separates two major Paleozoic sedimentary basins, the Appalachian Basin to the southeast and Michigan Basin to the west ( Figure 1A). The Algonquin Arch, which separates the two basins, extends from northeast of the Canadian shield to the southwest near the city of Chatham. Near the Windsor area, another structural high named the Findlay Arch begins and extends to the southwest into Michigan and Ohio states [39]. The Michigan basin is a nearly circular intracratonic basin dominated by carbonate and evaporite sediments, whereas the Appalachian Basin is an elongate foreland basin developed because of the collisional tectonic event along the eastern margin of the North American continent during the Paleozoic and is mainly characterized by deposition of siliciclastic sediments [40].
During Precambrian to Cambrian, a rifting event that separated North America and Africa initiated subsidence (approximately 1100 Ma) and deposition within the Michigan Basin ( Figure 1B; [37]). Thermal subsidence of the Precambrian basement followed approximately 580 Ma to 500 Ma after the cooled lithosphere increased in density [41]. The subsequent and continuous deposition of sediment into the Michigan Basin caused the bending of the basement [42][43][44]. However, the development of the basin is not the result of a continuous subsidence but represents the outcome of a series of tectonic events that characterized the Paleozoic [45,46].
The Taconic and Acadian Orogenies ( Figure 1B) had a dominant control on the Paleozoic strata of the eastern flank of Michigan [47]. In fact, during the Taconic Orogeny (Ordovician), a large-scale eastward tilting of the Laurentian margin has been produced with the consequent disappearance of the concentric depositional geometry of the Michigan Basin and with the subsequent deposition of Ordovician limestone and overlying shale successions [46,48,49].
The Middle Devonian Acadian Orogeny ( Figure 1B) determined the reactivation of the Algonquin-Findlay Arch system with consequent return of the concentric geometry of the Michigan Basin [45]. During the Upper Silurian where restricted marine conditions dominated, the deposition was mainly characterized by the evaporites of the Salina Group [50,51]. In addition, the presence of an unconformity reveals emergent conditions at the end of the Silurian [47].
The Caledonian (Devonian) and Alleghenian (Carboniferous) Orogenies ( Figure 1B) played a major role in terms of diagenetic fluid migration [5]. The regional maximum principal stress at that time was northwesterly oriented in accordance with the main direction of the thrust motion along the Appalachian tectonic front [47,52].
The Paleozoic sequence is separated from underlying Precambrian rocks of the ca. 1.1 Ga old Grenville orogeny by an unconformity overlain by a relatively thin Cambrian sandstone, which represents an aquifer [53,54]. The depositional environment of the Ordovician and Cambrian formations encompasses various facies assemblages, including supratidal, lagoonal to subtidal environments (e.g., [55,56]). The diverse facies changes observed in these formations were controlled by relative sea-level changes and sources of detritus, related to the regional tectonic setting during the Paleozoic [37,46]. In the study area, many of the Cambrian units along the Algonquin Arch ( Figure 1) were absent due to arch rejuvenation and uplift during early Paleozoic times (e.g., [56]). Ordovician formations consist of thin, decimeter-scale, alternating shale, and carbonate beds; some of which have been dolomitized [47]. The Michigan Basin has been characterized by a typical tropical climate during mid to late Silurian because it was located at 25 • south of the equator [57] leading to the development of shallow waters reef complexes. Moreover, the Michigan Basin has been isolated and is becoming an evaporative basin in which alternation of evaporites and carbonates is deposited during the rest of the Silurian [58]. Conversely, in the Middle Devonian, the Michigan basin has been characterized by an extremely arid climate with consequent deposition of shallow intertidal and supratidal lithofacies such as Amherstburg and Lucas formations [59]. The Devonian sequences are represented by limestones, shales, dolostones, and sandstones [39] with arid climate and deposition of shallow intertidal and supratidal lithofacies (e.g., Amherstburg and Lucas formations).

Sampling and Analytical Methods
This section may be divided into subheadings. It should provide a concise and precise description of the experimental results, their interpretation, as well as the experimental conclusions that can be drawn.
Cores (n = 91) from specific sections from deep boreholes (n = 10) from Bruce Nuclear Site (Tiverton, Ontario) and across the Huron Domain of southern Ontario were sampled ( Figure 2). These samples encompass formations ranging from Cambrian to Devonian age successions and representing a wide range of rocks including dolomitized limestones, dolostones, sandy dolostones, and sandstones ( Figure 2). These core samples were examined and photographed prior to their selection for making thin sections and fluid inclusion wafers. Petrographic examination of 145 thin sections was performed under a standard petrographic microscope to compliment core descriptions and by cathodoluminescence microscopy (CL) using a Technosyn cold cathodoluminescence stage with a 12-15 kV beam and a current intensity of 0.42-0.43 mA.
Samples for oxygen and carbon isotopic analyses were extracted from calcite and dolomite samples (n = 259) from polished slabs using a microscope-mounted drill assembly. The samples were reacted in a vacuum with 100% pure phosphoric acid for four hours at 25 • and 50 • C for calcite and dolomite, respectively. CO 2 gas from samples that contained a mixture of calcite and dolomite was extracted using the chemical separation method of Al-Aasm et al. [60]. The evolved CO 2 gas was analyzed for isotopic ratios on a Delta Plus mass spectrometer. Delta (δ) values for oxygen and carbon are reported per mil (‰) relative to the PeeDee Belemnite (VPDB) standard. Precision was better than 0.05‰ for both δ 18 O and δ 13 C.
Twenty-six core samples were selected for Sr isotopic analyses from various calcitic and dolomitic components representing host rock, skeletal components, and diagenetic phases. Prior to 87 Sr/ 86 Sr analysis, powdered samples (n = 78) were extracted employing a microscope-mounted drill assembly to extract 1-10 mg of powdered calcite or dolomite from polished slabs. Care was taken to avoid contamination from other components such as later cement fill. These samples were initially dissolved in 2.5 N HCl suprapure for 25 h in sealed PFA vessels at room temperature. They were then passed through cation exchange quartz columns. To remove the eventual pore salts that may contain Sr, these samples were, after crushing and pulverization were washed twice with MilliQ water. Strontium isotopic ratios were measured on a Finnigan MAT 261 mass spectrometer. All analyses were performed in the static multi-collector mode using Re filaments. NBS and ocean water were used as standard references and 87 Sr/ 86 Sr ratios were normalized to 87 Sr/ 86 Sr = 8.375209. The mean standard error was 0.00003 for NBS-987.
Thirty-one samples were analyzed for fluid inclusion microthermometric studies. Petrographic work was done with a conventional petrographic microscope. Careful and detailed observations of size, occurrence habit, and paragenetic contexts of the fluid inclusions were conducted in the double polished wafers. Distinctions were made between primary, secondary, and pseudo-secondary fluid inclusions. This was followed by a determination of the number of phases (i.e., the liquid to vapor ratio) and their ratios entrapped in each of the fluid inclusions. Then determination of the composition (aqueous or hydrocarbon fluids) of the inclusion using fluorescence microscopy; most oil-filled inclusions fluoresce, while the aqueous inclusions do not. Finally, microthermometric measurements of the homogenization (Th) and final ice melting temperatures (Tmice) of fluid inclusions were measured. These measurements were conducted using a Linkam THMGS600 heating-freezing stage, which was calibrated using synthetic fluid inclusions of known compositions. Homogenization temperatures (Th), and ice-melting temperatures (Tm-ice) were measured with a precision (reproducibility) of ±1 • C and ±0.1 • C, respectively. The salinities for inclusions in which ice was the last melting phase were calculated using the program of Chi and Ni [61], and those for inclusions with hydrohalite as the last melting phase was done with the program of Steele-MacInnis et al. [62].
Deformation of inclusion walls, necking-down of single inclusions into smaller inclusions, or leakage of fluid can all lead to changes in homogenization temperature. To limit the possibility of measuring deformed aqueous inclusions, only primary inclusions from the same field of view were measured during a single heating or freezing run. By restricting measurements to inclusions within the same field of view, any sudden changes in liquid/vapor ratios due to inclusion deformation could be observed and removed from consideration. This does not exclude the possibility of subtle deformation in the z-axis, which is difficult to detect. Heating runs were conducted before freezing runs to reduce the possibility of inclusion stretching by freezing [63].
Major (Ca, Mg), minor, trace elements (Sr, Na, Mn, and Fe) and rare-earth element (REE) analyses (n = 95) were performed using ICP-MS. The procedure used for the sample digestion has been adapted to digest carbonates only. Around 50-100 mg of powdered carbonate samples from different calcite and dolomite phases were digested utilizing 5% v/v acetic acid. Measurements were carried employing an Agilent 8800 inductively coupled mass spectrometer at the University of Waterloo. All results we reported in part per million (ppm).
For rare-earth (REE) element analysis, samples were analyzed using external calibration, internal standard, and standard addition methods, following the procedure described in Jenner et al. [64].  [65]. The proportions of LREE(La-Nd) over HREE (Ho-Lu) were calculated via (La/Yb) SN ratios as proposed by Kucera et al. [66]. The proportions of MREE (Sm-Dy) was calculated using (Sm/La) SN and (Sm/Yb) SN ratios [67] over LREE and HREE, respectively.
An environmental scanning electron microscope (ESEM) coupled with an EDAX detector was utilized to analyze selected carbon-coated thin section samples to investigate the nature of dolomitization and subsequent recrystallization.

Petrographic Analysis
Cambrian strata are dominated by mixed carbonate-siliciclastic lithofacies ( Figure 3). The carbonate lithofacies is composed of dolomitized and silicified crinoidal and bryozoan components. The siliciclastic facies consist of silicified ooids and other coated grains as well as detrital carbonates with large detrital quartz, cemented by carbonates.
The Middle Ordovician carbonates in the subsurface are divided into the Black River Group and the overlying Trenton Group ( [47]; Figure 3). The original depositional facies of the investigated Ordovician formations (Shadow Lake, Gull River, Coboconk, Kirkfield, Sherman Fall, Cobourg, and Collingwood) consists of nodular, bioturbated, micritic, and bioclastic facies-now partly dolomitized. In the Coboconk Formation, a thin fossiliferous layer (grainstone facies) occurs with a sharp contact with the nodular facies. This layer is also pervasively dolomitized. Skeletal components include bryozoans and echinoderms. In the Gull River Formation, the matrix is highly dolomitized but fossils such as brachiopods, crinoids, and trilobites are dolomitized but their textural integrity is preserved.
The diagenetic history of the Cambrian, Ordovician, Silurian, and Devonian successions in the Huron Domain encompasses multiple episodes of calcite cementation, dolomitization, anhydrite cementation, mechanical and chemical compaction, dissolution, and silicification. It is understood (e.g., [2,3]) that these diagenetic processes occurred during deposition and continued upon burial. The paragenetic sequence shown in Figure 4 is based on petrographic relationships, spatial distribution, and geochemical evidence. Due to the significance of dolomitization in these carbonates, the discussion in this paper will be concentrated on the most significant and relevant diagenetic events.  [2]) and paragenesis of the Silurian-Devonian formations (modified from Tortola et al. [3]).

Dolomite Types
In this study, multiple generations of dolomite have been observed ( Figure 4). The dolomite types occur as both replacements (D) and cement (SD) (cf. [2,3]). According to Sibley and Gregg's classification [68] and based on petrographic characteristics (e.g., crystal size, shape, extinction, fabric, and paragenetic relationships), several types of dolomite development are distinguished. Table 1 summarizes the main types of dolomite and their petrographic and CL characteristics. In the mixed carbonate-siliciclastic lithofacies of the Cambrian strata, the presence of dolomite is restricted to the dolomitized bioclastic lithofacies and scattered dolomite crystals within the silicified lithofacies ( Figure 5). Two replacive dolomite types occur in the Cambrian rocks (microcrystalline subhedral-D1 and medium to coarse crystalline subhedral to anhedral, zoned-D2; Figure 5) and minor amounts of saddle dolomite (SD) cement present in vugs (Table 1; Figure 5). This dolomite cement is succeeded by large blocky calcite cement ( Figure 5I) in fractures and vugs. D2 mimically replaces allochems, such as fossils and ooids. Ordovician: Two types of replacive dolomite are recognized in the Gull River, Coboconk, and Cobourg formations (Table 1; Figure 6): (1) microcrystalline dolomite replacing skeletal components and matrix (D1), and (2) medium crystalline replacive matrix dolomite (D2), which is fabric destructive, partly zoned and replaces the intranodular micritic matrix ( Figure 6). D2 dolomite is volumetrically more abundant than D1 and it is also observed along dissolution seams in partially dolomitized limestones, usually characterized by cloudy-dark cores and clear rims ( Figure 6C). A minor presence of saddle dolomite (SD) occludes pores and fractures and postdates D1 and D2 ( Figure 6D,E). Anhydrite cement postdates saddle dolomite ( Figure 6E,F). In other Ordovician formations, such as the Shadow Lake, Kirkfield and Sherman Fall, Silurian, which are mostly dominated by limestones of varied lithologies, no significant dolomitization is observed. However, the scattered occurrence of D1 dolomite is observed in Sherman Fall Formation replacing the calcitic muddy matrix.  Figure 7). These are: D1 which is a pervasive replacive micro to fine crystalline matrix dolomite and D2, which consists of medium crystalline, subhedral, zoned dolomite ( Figure 7B,C). It also occurs as a less abundant selective replacive medium (>50-100 µm) euhedral crystalline dolomite, only observed along dissolution seams in partially dolomitized limestones and usually characterized by cloudy-dark cores and clear rims ( Figure 7C). Ferroan saddle dolomite (SD; Figure 7D-F) occludes fractures and vugs predating late blocky calcite cement ( Figure 7E,F). Devonian: Two types of dolomites occur in the Devonian formations ( Figure 3). These are: a pervasive replacive euhedral to subhedral micro to fine crystalline matrix dolomite (D1) and a pervasive replacive euhedral to subhedral medium crystalline, zoned matrix dolomite (D2) that is also associated sometimes associated with dissolution seams ( Figure 8E). They are primarily observed in the Lucas and Amherstburg formations (Table 1; Figure 8).

Calcite Cementation
Early and late calcite type of cement are present in the investigated Paleozoic formations. These types of cement display various crystal habits and sizes and occlude both primary as well as secondary porosity, pre-and post-dolomitization (Figures 4 and 9). The types of calcite cement identified include isopachous, syntaxial overgrowth, equant, dog-tooth, and blocky cement (Table 1).
Within Cambrian strata, fractures are occluded by a coarse, equant/ blocky calcite cement with strong, bright CL ( Figure 9). In contrast, calcite cementation in the Ordovician formations is represented by an early equant calcite cement surrounding fossil allochems, isopachous rims around fossils and occluding hairline fractures (Figure 9), and later blocky calcite in fractures and vugs. The blocky calcite (BKC) in fractures postdates saddle dolomite cement ( Figure 9).

Other Diagenetic Phases
Together with the calcite and dolomite diagenetic phases reported above, silica-bearing minerals and sulfate minerals are also present in the studied successions ( Figure 10). Quartz is common in the Cambrian rocks typically replacing skeletal and nonskeletal components such as corals in coral-bearing facies and coated grains, such as ooids ( Figure 5G). Silicification pre-dates dolomitization but postdates early isopachous calcite cementations ( Figure 4). In the Ordovician strata, silicification replaced part of the muddy calcitic matrix, with microcrystalline dolomite (D1) replacing the silicified micrite and fossils. Anhydrite and gypsum are present in Cambrian, Ordovician, and Silurian rocks as fracture-fill, which postdates replacive dolomite ( Figure 10). The rare occurrence of authigenic K-feldspar is observed by SEM filling some pores in Ordovician and Silurian rock but predated D2. Also, very minor pyrite and fluorite are observed in Silurian dolostones. Table 2 summarizes the stable and radiogenic isotope results of dolomite and calcite cement types. Figures 11-15 show the relationships among these isotopes with respect to their petrographic characteristics.      The fine crystalline replacive matrix dolomite (D1) in the Ordovician carbonates has δ 18 O values vary from −12.4 to −5.3‰, and less negative δ 13 C values compared to the Cambrian D1 dolomite (ranging from −1.2 to 1.1‰; Figure 11). Its 87 Sr/ 86 Sr ratios range from 0.70808 to 0.72346, respectively. The highly enriched ratios are samples that occur along dissolution seams ( Figure 13). D2 overlaps with D1 with δ 18 O values ranging from −9.6 to −4.1‰, δ 13 C values from −0.4 to 2.2‰, and 87 Sr/ 86 Sr ratios from 0.70814 to 0.71255, respectively. The enriched ratio also represents a sample occurring along dissolution seams ( Figure 13)

Major, Minor, Trace and REE Elements Content
Analytical results of major, trace, and rare-earth (REE) element concentrations for the dolomites and calcites from various age groups are summarized in Table 3. All analyzed replacement and cement dolomite types are non-stoichiometric, consistently calcium-rich (average CaCO 3 ranging from 54.3 to 61.8 mol% for D1 and D2, and from 57.7 to 63.3 mol% for SD, respectively; Figure 16). There are no significant differences in Mg content in zoned D2 and SD between the cores and rims, as evidenced from the semiquantitative EDX analysis. In general, Mn and Fe contents in the Cambrian and Ordovician dolomites are much higher than the Silurian and Devonian dolomites (Table 3). There is also a covariant trend between Mn and Fe for these dolomites ( Figure 17). However, the average of Sr content in all dolomites does not show large variations among the investigated dolomite types (vary between 83 and 253 ppm, Table 3), but many samples from the Ordovician dolomites show high Sr content and covariant trends with Fe and Mn (Figure 17). Such a relationship is not apparent in the dolomites from other age groups.     Samples of dolomite and fracture-and vug-filling calcite cement from different age groups were analyzed for REE content ( Table 4). The ∑REE of Cambrian dolomite types and calcite cement is significantly higher than those in Ordovician, Silurian, and Devonian strata (Table 4) compared with the average of warm water brachiopods shown by Azmy et al. [70]. ∑REE of Devonian samples record the lowest values among the rest.
Replacive dolomite D1 and D2 sampled from Cambrian formations show different trends ( Figure 18) compared to the pattern of warm water brachiopods shown by Azmy et al. [70], with a dominant-negative La anomaly and positive Ce anomalies ( Figure 19).    Samples of dolomite and fracture-and vug-filling calcite cement from different age groups were analyzed for REE content ( Table 4). The ∑REE of Cambrian dolomite types and calcite cement is significantly higher than those in Ordovician, Silurian, and Devonian strata (Table 4) compared with the average of warm water brachiopods shown by Azmy et al., [70]. ∑REE of Devonian samples record the lowest values among the rest.
Replacive dolomite D1 and D2 sampled from Cambrian formations show different trends ( Figure 18) compared to the pattern of warm water brachiopods shown by Azmy et al. [70], with a dominant-negative La anomaly and positive Ce anomalies ( Figure 19).
Both dolomite D1 and D2 ( Figure 18) sampled from Ordovician formations show different trends compared to the pattern of warm water brachiopods shown by Azmy et al. [70], with a negative La anomaly and both cases of positive and negative Ce anomalies ( Figure 19). Both D1 and D2 exhibit higher average ΣREE (11.63 ± 21.56 ppm, and 8.54 ± 7.36 ppm, respectively) compared with the average of warm water brachiopods shown by Azmy et al., [70]. Saddle dolomite cement (SD) patterns show lower average ΣREE (9.55 ± 4.36 ppm) than D1 but slightly higher than D2. It also shows, similarly to D1 and D2, Both dolomite D1 and D2 ( Figure 18) sampled from Ordovician formations show different trends compared to the pattern of warm water brachiopods shown by Azmy et al. [70], with a negative La anomaly and both cases of positive and negative Ce anomalies ( Figure 19). Both D1 and D2 exhibit higher average ΣREE (11.63 ± 21.56 ppm, and 8.54 ± 7.36 ppm, respectively) compared with the average of warm water brachiopods shown by Azmy et al. [70]. Saddle dolomite cement (SD) patterns show lower average ΣREE (9.55 ± 4.36 ppm) than D1 but slightly higher than D2. It also shows, similarly to D1 and D2, slightly different patterns compared to warm water brachiopod's pattern with both negative La and Ce anomalies.
Replacive dolomite D1 and D2 sampled from Silurian formations show different trends compared to the pattern of warm water brachiopods shown by Azmy et al. [70] ( Figure 18), with a minor negative La anomaly and both cases of positive and negative Ce anomalies ( Figure 19). Brachiopods, D1, and D2 patterns with both negative La and Ce anomalies as well as a slight negative Eu anomaly ( Figure 19).
Both fine and medium dolomite matrix samples (D1 and D2, respectively) from Devonian formations have shale-normalized patterns sub-parallel to those of modern warm water brachiopods from Azmy et al. [70]. They both show the typical characteristics of seawater patterns such as depletion of LREE over HREE and slightly positive La and Gd anomalies ( Figure 19).
There is a subtle difference in Cambrian and Ordovician dolomite vs. Silurian and Devonian samples in terms of Sm/Yb vs. Eu/Sm ( Figure 20A

Fluid Inclusions Microthermometry
Microthermometric measurements of primary fluid inclusions from D2, SD, and BKC from each age group provide melting and homogenization temperatures (Th) and estimates of salinity (cf. [61,71] for each phase (Table 5 and Figure 21). Primary fluid inclusions within D1 crystals were too small to measure (cf. [71]). However, monophase fluid inclusions are present in this dolomite. For calcite, measurements were taken for coarse crystalline cements (BKC) in vugs and fractures. The measured Th provides an estimate of minimum entrapment temperatures. The low first melting temperatures and ice-melting temperatures suggest that the fluid within the fluid inclusions may approximate a H2O-NaCl-CaCl2 system. The examined inclusions are parallel to crystal facets and the vapor bubble size for the two-phase (liquid-vapor) fluid inclusions. There

Fluid Inclusions Microthermometry
Microthermometric measurements of primary fluid inclusions from D2, SD, and BKC from each age group provide melting and homogenization temperatures (T h ) and estimates of salinity (cf. [61,71] for each phase (Table 5 and Figure 21). Primary fluid inclusions within D1 crystals were too small to measure (cf. [71]). However, monophase fluid inclusions are present in this dolomite. For calcite, measurements were taken for coarse crystalline cements (BKC) in vugs and fractures. The measured T h provides an estimate of minimum entrapment temperatures. The low first melting temperatures and ice-melting temperatures suggest that the fluid within the fluid inclusions may approximate a H 2 O-NaCl-CaCl 2 system. The examined inclusions are parallel to crystal facets and the vapor bubble size for the two-phase (liquid-vapor) fluid inclusions. There is a notable overlap in T h and salinity values among all dolomite types in Cambrian, Ordovician and Silurian samples (Figure 21

Data Analysis and Discussion
Petrographic investigations, stable-isotope, elemental analyses, and fluid-inclusion microthermometry provide important clues regarding dolomitization and other diagenetic processes, fluid chemistry, and the history of fluids migration in the sedimentary basin. The Paleozoic successions in the Huron Domain of Michigan Basin were affected by early (marine to shallow eogenetic zone; [72]) and deeper (mesogenetic zone) diagenetic processes [72]. These include early calcite cementation, silicification, partial dissolution, early dolomitization and later dolomitization, recrystallization, fracturing, stylolitization, late calcite, and sulphate cementation. There are similarities in the diagenetic processes in Cambrian and Ordovician carbonates and Silurian and Devonian carbonates which allow us to consider and discuss these carbonates together (Figure 4).

Constraints from Petrography-Paragenetic Sequence
The paragenetic sequence of the main diagenetic events in the studied Paleozoic successions (Figure 4), which was constructed originally based on textural relationship and geochemical and fluid inclusion evidence, provides a useful platform for further discussion of the data obtained in this study.
As reviewed in Section 2 above, the depositional environments of the Cambrian and Ordovician carbonates encompass various lithofacies assemblages, including supratidal, lagoonal to subtidal environments (e.g., [47,55]). The diverse lithofacies changes observed in these formations are directly controlled by relative sea-level changes and sources of detritus. These were evident especially in the Cambrian strata where mixed siliciclasticcarbonate facies are common. Ultimately, these depositional processes also affected the diagenetic changes encountered in these sediments during shallow and deeper burial stages. In carbonate-rich units, the bioclastic grainstone facies has been influenced by early stabilization of calcitic/aragonitic components and dolomite replacement at shallow burial (D1; Figure 5B) to be succeeded by medium to coarse crystalline dolomite (D2; Figure 5B-D,I), possibly at an intermediate burial. Minor presence of saddle dolomite (SD) occluding vugs postdates D2 but predates BKC, which is occluding fractures; both of which probably formed later than D2 in the intermediate burial stage (Figure 4; Figure 5E). In the ooid grainstone lithofacies, which is extensively silicified (Figure 5G), the early isopachous calcite cement was silicified. However, the silicified ooids are well preserved and show evidence of very early diagenetic, the pre-compaction origin of these allochems. The siliciclastic-dominated facies were cemented earlier by calcite cement later to be compacted with the creation of some minor fractures occluded by equant calcite cement, which was partly dolomitized. The basal Cambrian sandstone is assumed to have been a major migratory path for diagenetic fluids [73,74]. Strata in the eastern region of the basin pinches out towards the margins. These thin fingers of clastic material which have large amounts of bioclastic material increase the likelihood for extensive dolomitization of a thin layer.
In the Ordovician formations (e.g., Coboconk Formation) diagenesis was initiated by early cementation that produced the nodular facies on the seafloor. Burial and compaction affected the formation of these nodules and may have created hairline fractures, later occluded by dolomite. Dolomitization (D1; Figure 6) occurred by selective early replacement of skeletal components such as bryozoans and echinoderms and the micritic matrix, as well as, pervasively replacing the intra-nodule matrix during shallow burial. These dolomites were crosscut by hairline fractures and stylolites during a burial at a later stage, suggesting the formation of D1 at a shallow burial. Silicification, which occurs in some facies of the Coboconk Formation occurred prior to dolomitization as evidenced by silica replacing the matrix and fossils that later were selectively dolomitized. In the Gull River Formation, a similar paragenetic sequence to the Coboconk Formation is observed (Figure 4). Initially, dolomitization proceeded with fine crystalline (D1) and medium to coarse crystalline matrix replacement dolomite (D2). The minor saddle dolomite (SD) cement postdates earlier dolomite formation and occludes vugs and fractures and is also postdated by late, blocky calcite cement (BKC), which is occluding vugs and fractures ( Figure 6). Other compaction features such as stylolites postdate or sometimes appear to be contemporaneous with the medium to coarse crystalline dolomite (D2; Figure 6C). In both the Coboconk and Gull River formations, the intra and internodule dolomites typically present as petrographically distinct. These distinctions reflect differences in the underlying substrate that was originally dolomitized and relate to permeability, not fluid chemistry, or timing (notwithstanding some evolution of brine chemistry over time). An early, shallow burial origin is suggested for the fine crystalline (D1) dolomite as evidenced by the selective replacement of this dolomite, its shape and size, and it predates stylolite formation, consistent with a shallow burial diagenetic environment. In contrast, the medium and coarse crystalline (D2) dolomite may have formed at medium burial stages as evidenced by its increasing crystal size, subhedral and anhedral shape (e.g., [26]), and close association sometimes with compaction features, such as stylolites. The stylolites observed in the samples all occur at~90 • to the fractures. This relationship may be explained by compressive stress, imposed by the Taconic Orogeny, which was occurring concurrently with deposition. These stylolites may have acted as conduits for Mg-rich fluids (e.g., [75]).
In contrast, the depositional environments of the Silurian and Devonian successions encompass various lithofacies assemblages, including shallow marine ramp and deeper water subtidal lithofacies dominated by limestone and dolostone lithofacies with intercalations of shales and evaporite deposits in both age groups (e.g., [47,76,77]). In the Silurian and Devonian carbonates early calcite cementation (e.g., ISC and SXC; Table 1) and micritization initiated on the seafloor to be followed during burial by occlusion or porosity by later calcite cementing vugs and fractures (e.g., DC, and BKC; [78]). BKC postdates saddle dolomite (SD; Figure 4). Evaporite minerals, such as gypsum and anhydrite represent the final stage of cementation in the Silurian rocks. Dolomitization commenced early (D1) formed replacing calcitic and aragonitic allochems and muddy calcitic matrix [3]. This dolomite may have formed at a shallow burial, in a reflux evaporative setting (e.g., [33,79,80]). This is supported by its fine crystalline crystal size [68] and the cross-cutting relationship with compactional features, such as dissolution seams and stylolites which crosscut this dolomite. This dolomitization phase is followed by the formation of medium to coarse crystalline replacive dolomite (D2) at the medium burial stage. The increase in crystal size of D2 may reflect dolomitization of coarser lithofacies (such as grainstones) and/or recrystallization of precursor dolomite matrix (D1; [23]). In the Silurian carbonates, only saddle dolomite (SD) is observed occluding fractures, postdating D1 and D2 but predating blocky calcite (BKC; Figure 7). This dolomite formed possibly in an intermediate burial setting.
Evaporite minerals (gypsum and anhydrite; Figure 10) represent the last diagenetic phases and form in a deeper burial realm. These minerals occlude fractures and vugs in Cambrian, Ordovician, and Silurian, carbonates.

Stable and Radiogenic Isotopes
Cambrian: The stable isotopic composition of D1, D2, and SD in Cambrian rocks ( Figure 11) are notably lighter than the postulated dolomite precipitated in equilibrium with marine waters (e.g., [69]). These values are also distinct from other age groups ( Figure 11). There is little isotopic fractionation of 13 C/ 12 C with temperature, and hence the δ 13 C values of these dolomite reflect isotopic fractionation of C related to fluids that are enriched in lighter CO 2 originated from oxidized carbon, which can be related to seawater-sourced, reduced, bacterial sulfate reduction process for D1 and hydrocarbon sources for D2 and SD [2].
However, the δ 18 O values of the dolomites reflect both the temperature and the composition of the parent fluids, and thus the departure from equilibrium values reflects precipitation of these dolomite under a different diagenetic regime [81]. The negative shift in δ 18 O of D1 and D2 (Figure 11) signifies the effect of temperature increase during burial, change in fluid composition, or both [81]. Recrystallization of D1 with increasing temperature/change in fluid chemistry can also result in a negative shift in D1 (e.g., [23]). This recrystallization is also observed petrographically, as shown by the increasing crystal size of D1 and the presence of overgrowth zonation in D2 ( Figure 5). The co-variant trend between carbon and oxygen isotopes in these dolomites reflect varied water-rock interaction with increasing burial (cf. [82]). Water-rock interaction is also evidenced by the presence of more radiogenic 87 Sr/ 86 Sr ratios in D1 and D2 (Figures 13 and 15) in comparison to the postulated range for Cambrian seawater. This can be related to the presence of more siliciclastics, rich in clays and feldspars, in the Cambrian rocks as well as can be sourced from the Precambrian shield [74,83]. The same trend and comparable negative isotopic values are evident in late calcite cement (BKC; Figure 12) and thus reflecting similar diagenetic conditions during later fluid flow events for these carbonates unique to Cambrian paleohydrologic flow system. Similarly, BKC also shows more radiogenic Sr isotope ratios than the Cambrian seawater values (Figure 14).
Ordovician: Most of the oxygen and carbon isotopic values of the Ordovician replacive dolomites (D1 and D2) show overlapping signatures with a departure by at least 3‰ from their coeval postulated seawater equilibrium values with respect to their age (Figure 11), especially for oxygen isotopes. This can be related to the recrystallization of D1 to D2 at elevated temperatures during the intermediate burial stage (Figure 4). C isotope values are closer to their seawater equilibrium values reflecting the buffering of this isotope with original carbonate values (e.g., [26,29]). However, δ 13 C values of saddle dolomite show slight negative shifts compared to the postulated equilibrium values. This may reflect an oxidized organic C source during burial. The 87 Sr/ 86 Sr ratios of many samples of D1 and D2 show comparable and/or slight radiogenic signature to seawater values ( Figure 13; [69]) with the exception of few dolomites that show highly radiogenic rations (Figures 13 and 14). These dolomite samples occur along dissolution seams and thus reflecting a radiogenic fluid source of Sr, perhaps from Cambrian siliciclastic fluid. In contrast, the isotopic composition of early calcite cement shows preservation of the original seawater signature (Figure 12), while late calcite cement shows some negative trend in both oxygen and carbon isotopes relative to the postulated values in Ordovician calcite. This reflects the diagenetic progression and precipitation of this cement from fluids modified at the burial. This is also shown by 87 Sr/ 86 Sr ratios of BKC which show seawater signature (Figure 14).
Silurian and Devonian: The isotopic data of the overlying Silurian and Devonian carbonates show overlaps between δ 13 C and δ 18 O values ( Figure 11). Dolomite samples from both age groups show similar δ 13 C values compared to those of carbonates deposited in equilibrium with seawater of respective age except for few D1 samples from the Silurian carbonates that show a slight negative shift, which can be related to a depleted carbon source related to bacterial sulfate reduction at the shallow burial [84]. This is supported by the presence of minute framboidal pyrite crystals embedded within RD1 (this study and Tortola et al. [3]). Also, few BKC samples from the Silurian carbonates show a negative shift from the equilibrium values ( Figure 12). This can be related to an oxidized organic carbon source, possibly related to interaction with hydrocarbons (e.g., [2]). However, δ 18 O values of D1 show evidence of dolomite recrystallization as these values depart from their respective age equilibrium values. This negative shift in oxygen isotopes also reflects increasing temperature and/or changes in the diagenetic fluid composition with increasing burial (cf., [3,26]). 87 Sr/ 86 Sr ratios of many samples, both dolomite and calcite, from these age groups cluster close to their coeval seawater composition ( Figure 13). However, a few dolomite samples show slight enrichment of radiogenic Sr ratios. This can be related to the presence of modified seawater fluids upon burial.

Major, Minor, and REE Elements
The average values of CaCO 3 in all analyzed replacement and cement dolomite types show that all dolomite types are non-stoichiometric, consistently calcium-rich (Table 4, average CaCO 3 range from 54.3 to 61.8 mol% for D1 and D2, and from 57.7 to 63.3 mol% for SD, respectively; Figure 16). This suggests that the precursor mineralogy of these dolomites is of high-Mg calcites and aragonite and the stabilization of these minerals to dolomite occur in semi-closed diagenetic systems with less fluid mobility and limited supply of Mg ions (cf., [18,85]). Hence, dolomitizing fluids reflect a lowering Mg/Ca ratio in these fluids which encourages Ca-rich, non-stoichiometric dolomite to form. In general, Mn and Fe contents in the Cambrian and Ordovician dolomites, which are indicators of redox conditions during dolomitization, are much higher than the Silurian and Devonian dolomites (Table 3) reflecting a decrease in oxidizing conditions in the older successions. There is also a covariant trend between Mn and Fe for these dolomites (Figure 17) reflecting the geochemical affinity of these two elements [86]. However, the average of Sr content in all dolomites does not show large variations among the investigated dolomite types (vary between 83 and 253 ppm, Table 3) and are comparable to those of other ancient dolomites but lower than the estimated values for dolomite in equilibrium with seawater (~500-800 ppm, (e.g., [87]). However, many samples from the Ordovician dolomites show high Sr content and covariant trends with Fe and Mn (Figure 17), perhaps reflecting the aragonitic precursor mineralogy of these dolomites. Such a relationship is not apparent in dolomites from other age groups.
Rare earth elements (REEs) have been commonly utilized for the reconstruction of a paleoenvironmental history of ancient seawater [88][89][90], diagenetic alteration, and the temperature of fluids precipitating carbonate minerals [70]. REE can also provide information about the relative water depth and the proximity to hydrothermal sources [91]. Some workers indicate that carbonates might preserve their REE compositions and/or patterns during diagenetic processes [89,92]. However, many studies have revealed that consequent diagenetic alteration modifies REE concentrations and/or patterns [70,93]. Therefore, the redistribution of REE in carbonates is an indicator that can be utilized to ascertain the potential contributing factor to observed REE concentrations in diagenetic minerals, especially in high water/rock ratio diagenetic systems [70]. The rare earth elements generally exist in the stable trivalent (3+) oxidation state, whereas Eu (Eu 2+ , Eu 3+ ) and Ce (Ce 3+ , Ce 4+ ) exhibit distinct behavior due to the redox state of the natural system (e.g., [93]). Ce (Ce 3+ ) can be oxidized to Ce (Ce 4+ ) in marine systems under surface conditions, whereas high temperature and/or high reducing conditions cause Eu (III) to be reduced to Eu (II) [65,94]. These characteristics provide the utilization of REEs as tracers of chemical processes during the diagenetic history of carbonates [89,[95][96][97].
The ∑REE of Cambrian dolomite types and calcite cement (Table 3) is significantly higher than those in Ordovician, Silurian, and Devonian strata when compared with the average of warm water brachiopods shown by Azmy et al. [70]. ∑REE of Devonian samples record the lowest values among those examined in this study. The higher concentration of ∑REE of Cambrian dolomites signify that the dolomitizing fluids contained higher concentrations of REE, likely related to the interaction of these fluids with basement rocks.
Replacive dolomite D1 and D2 sampled from Cambrian formations show different trends ( Figure 18) compared to the pattern of warm water brachiopods shown by Azmy et al. [70], with a dominant-negative La anomaly and positive Ce anomalies ( Figure 19). These dolomites may have formed under suboxic to anoxic conditions [70]) from early crustal fluids highly enriched in REE particularly MREE and HREE relative to dolomites from other age groups. Saddle dolomite from the Cambrian rocks show evidence of different fluids that produced D1 and D2 are likely related to hydrothermal fluids that originated from deeper sources and transported via fractures and faults [2]. The flat Ce* anomaly, slightly positive Eu*, flat LREE over HREE, and sight MREE enrichment with low La/Sm in saddle dolomite reflect higher temperature fluids of seawater origin under low oxic conditions modified by siliciclastic contamination.
Both D1 and D2 ( Figure 18) sampled from Ordovician formations show different trends compared to the pattern of warm water brachiopods shown by Azmy et al. [70], with a flat to slight Eu anomaly, negative La anomaly, and both cases of positive and negative Ce anomalies ( Figure 19). Both, D1 and D2 exhibit higher average ΣREE (11.59 ± 22.44 ppm, and 8.54 ± 7.36 ppm, respectively) compared with the average of warm water brachiopods shown by Azmy et al. [70]. Saddle dolomite cement (SD) patterns show lower average ΣREE (9.55 ± 4.36 ppm) than D1 but slightly higher than D2. It also shows, similarly to D1 and D2, slightly different patterns compared to Warm water Brachiopod's pattern with both negative La and Ce anomalies. High Eu*, high La/YbN* and Y/Ho can indicate interaction with hydrothermal fluids [90].
Silurian replacive dolomites (D1 and D2) show different trends compared to the pattern of warm water brachiopods shown by Azmy et al. ([70]; Figure 18), with a minor negative La anomaly and both cases of positive and negative Ce anomalies ( Figure 19). Brachiopods, D1 and D2 show patterns with both negative La and Ce anomalies as well as a slight negative Eu anomaly ( Figure 18). These data suggest a seawater-dominated source of diagenetic fluids evolved during burial and increasing water/rock interaction [3]. The higher ΣREE and a more pronounced enrichment in LREE in SD and BK reflect a different source of REEs and/or involvement of saline fluids.
Both, fine and medium dolomite matrix samples (D1 and D2, respectively) from Devonian formations have shale-normalized patterns sub-parallel to those of modern warm water brachiopods from Azmy et al. [70]. They both show the typical characteristics of seawater patterns such as slight depletion of LREE over HREE and slightly positive La and Gd anomalies ( Figure 19) potentially modified from the interaction of meteorically derived fluids and fluids of marine parentage. This can occur at near-surface, slightly reducing conditions. Although there is no evidence of karstification in the samples collected and analyzed in this study, a previous study [98] suggested that dissolution processes were mostly active in the shallow subsurface, usually <200 m depth in south-western Ontario. Thus, the possibility of meteoric water influence cannot be discounted completely.
There is a subtle difference in Cambrian and Ordovician dolomite vs. Silurian and Devonian samples in terms of Sm/Yb vs. Eu/Sm ( Figure 20A), La/Yb vs. La/Sm ( Figure 20B), Y/Ho vs. Sm/Yb ( Figure 20C) and Y/Ho vs. Eu/Sm ( Figure 20D) ratios. These differences may reflect changes in fluid composition, temperature, and/or redox conditions [99,100]. For example, Y/Ho vs. Sm/Yb ratios showing two possible fluids (marine and hydrothermal), and their interactions more evident in Cambrian and Ordovician dolomites. In contrast, the effect of seawater on the formation of Silurian and Devonian dolomites is more important. Other ratios document the presence of modified seawater as the primary fluid type in the formation of these dolomites [100].

Fluid Inclusions
Previous studies on fluid inclusions in Cambrian to Devonian dolomite and calcite in both central and marginal parts of the Michigan basin demonstrated fairly a wide range of T h and salinity values (e.g., [2,3,32,33,[101][102][103]). These investigations reported high T h and salinity values in the central part of the basin compared to the margin of the basin, which some of them attributed to hydrothermal fluids which originated from the central part of the basin. Haeri-Ardakani et al. [33] further suggest a thermal anomaly related to the mid-continent rift under the Michigan Basin during late Devonian-Mississippian is the major source of heat for hydrothermal fluids.
The presence of single-phase fluid inclusions in D1 suggests very low temperatures, up to 50 • C [71]. However, the salinity-homogenization temperature plots ( Figure 21A,B) for both dolomite (D2 and SD) and calcite (BKC) demonstrated a clear and subtle difference between age groups. For instance, higher T h and salinity values are recorded in SD and some D2 in Ordovician and Silurian samples compared to Devonian dolomite. In Cambrian, D2 overlaps with Ordovician samples. In contrast, BKC in Devonian samples show the highest T h and lowest salinity values compared to the rest ( Figure 21B). High salinity and homogenization temperatures values are shown by Silurian dolomite and calcite. This reflects the effect of fluids sourced from saline brines related to evaporite dissolution in the Silurian rocks (e.g., [3,32]). The differences in homogenization temperature and salinity values in BKC ( Figure 21) between Silurian and Devonian may reflect divergent fluid systems, such as a fluid that is characterized by high temperature and high salinity in the Silurian rocks and much higher temperature and lower salinity fluids in the Devonian rocks. These phenomena can also reflect the effect of Alleghenian Orogeny in Carboniferous time ( Figure 1) as a factor that affected the fluid composition and sources.
In all analyzed D2 and SD dolomite and late BKC in the investigated age groups the high T h values are significantly higher than the maximum burial temperatures suggested by other studies (e.g., [2,32,33,101,102]). This may suggest an involvement of hydrothermal fluids during the formation/recrystallization of dolomite and precipitation of late calcite cement. Increasing temperature of formation/recrystallization, as well as water/rock interaction, can explain the negative departures from equilibrium δ 18 O values of dolomite.

Dolomite Recrystallization
Dolomite crystal shape, its boundaries, and its characteristics are directly related to the crystallization environment in which they formed [104]. Recrystallization of earlier formed dolomite (e.g., [23,[105][106][107]) can involve change in texture, such as coarsening of crystal size or increase in nonplanar boundaries, change in structure (progressive ordering and strain), change in composition, which include isotopes, trace elements, stoichiometry, fluid inclusions, and zonation. It can also include a change in paleomagnetic properties [11].
Recrystallization of D1and D2 dolomites in the studied rocks has been invoked because of petrographic and geochemical data. Alteration of early formed dolomite is a common diagenetic process during progressive burial and/or increasing water/rock interactions (e.g., [8,108,109]). SEM and petrographic evidence for recrystallization of D1 include: (1) increase of crystal size ( Figure 22); (2) lack of zonation (e.g., overgrowth) under CL, especially in D2; (3) lack of stoichiometry in D1 and D2 can be related to recrystallization of D1 to D2 by Ca-rich, warmer basinal fluids, perhaps under a semi-isolated diagenetic system (e.g., [3,13,110]). Geochemical and isotopic evidence for recrystallization with increasing burial include: (1)

Evolution and Origin of the Diagenetic Fluids
Fluid flow in sedimentary basins can be linked to the tectonic evolution of the basin which can involve tectonic thrusting, sedimentary loading, uplift, and compression [4]. These fluids can transfer dolomitizing fluids as well as metals and hydrocarbons. Recent studies on diagenetic fluid flow in southern Ontario provided models on the fluid com-position and its source(s). Coniglio and Williams-Jones [111], Middleton et al. [55] and Coniglio et al. [112] advocated that lateral and vertical movement of Mg-bearing fluids in Michigan Basin resulting from the transformation of clay minerals, originating from younger, older, or correlative argillaceous sediments from deeper in the Michigan and Appalachian basins towards the basin margin was considered as a mechanism for migration of dolomitizing fluids. Yoo et al. [113] also considered compaction driven flow to explain the mass and heat from the Michigan Basin. Haeri-Ardakani et al. [32,33] suggested that hot brines in the central part of the basin migrated through the basal Cambrian sandstone and ascended through a fracture network precipitating dolomite and late-stage calcite while carrying hydrocarbons.
Hobbs et al. [114] suggested two geochemical systems present on a regional scale; a shallow system (<200 m) below ground surface containing fresh through brackish waters and an intermediate to deep (>200 m) system containing brines associated with hydrocarbon reservoirs. Petts et al. [115] using microthermometric and isotopic data of secondary minerals in Cambrian and Ordovician strata suggested the presence of hydrothermal brines mixed with connate water. Al-Aasm and Crowe [2] in their investigations of dolomitization in the Cambrian and Ordovician carbonates in the Huron Domain identified two somewhat isolated digenetic fluid systems; an earlier Cambrian system and a later Ordovician system, both displaying unique geochemical attributes. Bouchard et al. [116] suggested that the fluids in the Ordovician low permeability sequences in the eastern part of Michigan Basin have originated as Ordovician seawater that was chemically overprinted by evaporated seawater that infiltered during the late Silurian and mixed with pre-existing brines from a deeper part of the hydrostratigraphic section. More recently, Tortola et al. [3] suggested the presence of a diagenetic fluid system that affected Silurian carbonates and was altered by salt dissolution post-Silurian and a Devonian fluid system that affected by the Alleghanian orogeny and resulted in a different fluid source and composition.
Petrographic relationships and isotope data can be used to infer changes in the composition of pore fluids in the Paleozoic successions during progressive diagenesis. Using the mean Th values and the range of δ 18 O for D2 and SD from the investigated time frames, the δ 18 O composition of the precipitating fluids was determined. The relationships among mineral δ 18 O values, pore water δ 18 O values, and temperature for diagenetic phases from the investigated carbonates are shown in Figure 23A,B. The suggested pathways for pore water evolution reflect changes in temperature, changes in fluid chemistry, or both. For dolomite types/generations ( Figure 23A), there is a subtle difference in temperature and hence fluid composition between D2 and SD in Cambrian strata. The dolomitizing fluids for D2 reflect marine parentage, while for SD it is more enriched in δ 18 O values suggesting the modified brines formed under higher temperatures. In contrast, there is a significant overlap in the δ 18 O composition of dolomitizing fluids for Ordovician D2, SD, and Devonian D2 and a comparable range of formation/recrystallization temperatures. However, there is a measurable difference between D2 and SD in Silurian samples, as D2 formed under lower temperatures and comparable seawater δ 18 O values but higher temperatures and more enriched values for SD formation. This can be related to an increase in salinity with burial, which can be linked to the Salina salt (dissolution) in the Upper Silurian, and the incursion of hydrothermal fluids during the formation of SD (e.g., [117]).  The formation of blocky calcite (BKC) later in the diagenetic sequence reflects variations in formation temperatures and fluid composition in the investigated age groups ( Figure 23B). All calcites formed at higher temperatures and deviate from the seawater composition of corresponding ages. However, similarities in δ 18 O fluid composition between BKC and D2 in the Cambrian samples suggest the fluids could have originated from the Precambrian rocks below and migrated to the overlying porous Cambrian strata (cf. [74]). In contrast, the composition and temperatures of fluids that formed the late calcite in other age groups changed, becoming more enriched and thus reflecting evolved basinal brines [2]. The pathways of these fluids are shown by increasing both δ 18 O values and temperature with time (from older to younger systems). Variations in pore fluid isotopic composition in the Paleozoic age groups suggest compartmentalization of these diagenetic fluids resulting in the formation and/or recrystallization of dolomite and calcite cementation. This phenomenon is also supported by radiogenic Sr isotopic ratios of these diagenetic phases as well as REE distribution of these carbonates (Figures 14 and 23).
Fluid compartmentalization in carbonate strata in the Huron Domain is also illustrated in Figure 24, which shows some of the important geochemical parameters used in this study vs. their stratigraphic position. There is a subtle difference in stable isotopic composition between different age groups, especially in the Cambrian samples. In all dolomite types there is a significant departure in δ 18 O values from the postulated range of dolomite that was precipitated in equilibrium with seawater of respective ages. This can be related to increasing temperature during burial as well as modification of pore fluid isotopic composition, in SD this negative shift can be related to the effects of hydrothermal fluid interaction with host rocks; perhaps limited in extent and distribution. Salinity and Th measurements in the Silurian also show a distinctive pattern with higher and salinity temperatures of formation comparing to other rock units, which suggest the effect of evaporite dissolution and migration of fluids in these rocks. Figure 24. Stratigraphy vs. isotopic and microthermometric data for the investigate age groups.
The conceptual model for fluid flow mechanism in the Huron Domain shown in Figure 25 illustrates the concept of fluid compartmentalization in different stratigraphic rock units suggested in this study. The current burial depth is shown also in the model based on various studies [80,111,120]. As mentioned earlier this region of the Michigan Basin has experienced multiple stages of tectonic evolution, both passive and active tectonic events which affected fluid composition and migration pathways (e.g., [2,5,32,37]). Migration of warm, basinal brines from the center of Michigan Basin towards the margin [32] transported dolomitizing fluids and other solutes. The source of heat of these fluids could be related to the reactivation of the buried mid-continental rifts during late Devonian to Mississippian time [121]. A limited influx of hydrothermal fluids moved upwards from the Precambrian basement via fractures and faults [2,47] results in minor precipitation of saddle dolomite in fractures and vugs and altered earlier formation of matrix dolomite. This likely occurred during Taconian and Acadian Orogenies. Interaction of connate fluids of marine parentage with warm, basinal brines controlled the formation of dolomite and calcite cement at a shallow to intermediate burial depth. We believe that diagenesis is horizontally (stratabound) uniform across the Huron domain, displaying few signs of significant vertical connectivity beyond formation tops, except in areas with local heterogeneity (faults; [122]). Downward migration of saline fluids that originated from the dissolution of Silurian evaporites and possible interaction with basinal brines moving up-dip for the center of the basin was responsible for the modification of diagenetic phases [33]. In the Devonian hydrogeochemical system, the interaction of basinal brines with possible infiltration of younger fluids, perhaps of meteoric water origin, resulted in the formation of late dolomite and calcite cement (e.g., [55,111,112]). This probably occurred during the Alleghenian Orogeny (e.g., [3]). Sutcliff et al. [124] in their study of late calcite cement recovered from sub-horizontal fractures near the base of the Silurian sequence dated show U-Pb ages of 318 ± 10 Ma determined by LA-ICP-MS and 313 ± 1 Ma by ID-TIMS. Hence, these authors also suggested that the fluid flow of hydrothermal brines was influenced by the Alleghanian Orogeny which caused the migration of hydrothermal fluids westwards from high-level Alleghanian plutons through the basal Cambrian rocks.

Conclusions
Based on petrographic, geochemical, and fluid inclusion data presented, the following conclusions have been drawn with regards to the diagenetic evolution of the Huron Domain of southern Ontario: • Dolomitization of Cambrian to Devonian carbonates occurred through multiple events, with the formation of the earliest replacive dolomite phase (D1) occurring at shallow burial depths. This was followed by the precipitation of a later replacive dolomite phase (D2) and the formation of dolomite cement (SD) occluding vugs and fractures, which were formed at an intermediate burial depth. Evaporite and late blocky calcite cement followed the formation of these dolomite phases. • These types of dolomite types display unique geochemical and microthermometric signatures suggesting the flow of diagenetic fluids that were driven by tectonic processes over the history of the Michigan Basin. The diagenetic fluid composition was dominated initially by seawater and later modified upon burial and interaction with basinal brines. For Cambrian and Ordovician carbonates two possibly isolated diagenetic fluid systems affected them. These are: (i) an earlier fluid system that is characterized by a pronounced negative shift in oxygen and carbon isotopic composition, more radiogenic ratios, warm and saline signature with higher average ∑REE, negative La anomaly and positive Ce anomaly; and (ii) a later Ordovician system, characterized by less negative shifts in oxygen and carbon isotopes, comparable homogenization temperature, hypersaline, a less radiogenic fluid system, less negative La anomaly and positive and negative Ce anomaly, but also higher average ∑REE compared to warm water marine brachiopods. • There is a significant overlap of δ 13 C and δ 18 O values and 87 Sr/ 86 Sr ratios in Silurian and Devonian carbonates reflecting both a modified marine source of fluids as well recrystallization of earlier formed dolomite. High homogenization temperatures and salinity in saddle dolomite (SD) in the Silurian sequences and its REE shale-normalized patterns suggest that in both age groups the diagenetic fluids were originally of coeval seawater composition and were subsequently modified via water-rock interaction. The different evolution of the diagenetic fluids is more prominent in REE SN patterns from Silurian samples and it is possibly related to brines, which were modified by the dissolution of Silurian evaporites from the Salina series. Data Availability Statement: The data will be available with the authors.