Constraining Geogenic Sources of Boron Impacting Groundwater and Wells in the Newark Basin, USA

: The Newark Basin comprises Late Triassic and Early Jurassic fluvio-lacustrine rocks (Stockton, Lockatong, Passaic, Feltville, Towaco, and Boonton Formations) and Early Jurassic diabase intrusions and basalt lava flows. Boron concentrations in private well water samples range up to 18,000 µ g/L, exceeding the U.S. Environmental Protection Agency Health Advisory of 2000 µ g/L for children and 5000 µ g/L for adults. Boron was analyzed in minerals, rocks, and water samples using FUS-ICPMS, LA-ICP-MS, and MC ICP-MS, respectively. Boron concentrations reach up to 121 ppm in sandstone of the Passaic Formation, 42 ppm in black shale of the Lockatong Formation, 31.2 ppm in sandstone of the Stockton Formation, and 36 ppm in diabase. The δ 11 B isotopic values of groundwater range from 16.7 to 32.7‰, which fall within those of the diabase intrusion (25 to 31‰). Geostatistical analysis using Principal Component Analysis (PCA) reveals that boron is associated with clay minerals in black shales and with Na-bearing minerals (possibly feldspar and evaporite minerals) in sandstones. The PCA also shows that boron is not associated with any major phases in diabase intrusion, and is likely remobilized from the surrounding rocks by the intrusion-related late hydrothermal fluids and subsequently incorporated into diabase. Calcite veins found within the Triassic rock formations exhibit relatively elevated concentrations ranging from 6.3 to 97.3 ppm and may contain micro-inclusions rich in boron. Based on the available data, it is suggested that the primary sources of boron contaminating groundwater in the area are clay minerals in black shales, Na-bearing minerals in sandstone, diabase intrusion-related hydrothermal fluids, and a contribution from calcite veins.


Introduction
Boron, a metalloid element, is widely distributed within minerals across the earth's crust, with an average concentration of 8 parts per million (ppm).Boron in groundwater can originate from both natural geogenic sources and anthropogenic activities.According to guidelines set by the World Health Organization [1], the ideal boron concentration in drinking water should not exceed 500 µg/L.Furthermore, the World Health Organization has set a guideline for boron in drinking water at 2400 µg/L and established a daily intake limit for boron for adults at 0.17 mg/kg body weight [1].A long-term threshold value of 500 µg/L has been suggested for boron concentration in irrigation water [2] and reference therein.In the Newark Basin area, boron in groundwater has concentrations exceeding the U.S. Environmental Protection Agency Health Advisory (2000 µg/L for children and 5000 µg/L for adults) [3,4] and reaching up to 24,000 µg/L in some areas [5].Despite limited testing for boron in well water, Newark Basin studies revealed that 18 out of 376 well water samples (5%) contained high concentrations of this element [5,6].
Boron can pose serious health hazards when consumed in high concentrations.It can affect the male reproductive system, reducing testosterone levels in men, and may affect fetal development [7].The presence of boron in groundwater may also pose a risk to plants when irrigated with concentrations over 1000 µg/L ( [2] and references therein).In the Newark Basin, some sensitive species of plants irrigated with well water containing boron levels exceeding 1000 µg/L sustained damage [5].
In the rural and suburban communities of the western portion of the Newark Basin, including areas in and around Hunterdon, Mercer, and Somerset Counties in New Jersey, as well as Bucks and Montgomery Counties in Pennsylvania, many communities are entirely reliant on individual private wells and are at risk for boron contamination of well water.Because well owners are responsible for testing their own water and boron is an unregulated contaminant, most well owners are unaware of the risks of boron exposure and the need to test their water.
Boron is difficult to remove from water, as testing has shown that reverse osmosis only removes about 15% of the boron in well water.The New Jersey Geological and Water Survey (NJGWS) tested and recommended a combination of reverse osmosis followed by mixed-bed deionization as an effective boron removal system for point-of-use applications, specifically for drinking and cooking water.While whole-house water treatment would be ideal, as boron is likely absorbed through the skin, a system utilizing reverse osmosis followed by mixed-bed deionization media is expensive and impractical for whole-house water treatment.
A better option is to investigate the well using geophysical methods and test the various water-bearing zones to identify the zone contributing the boron.Once identified, sealing off that zone from the well would solve the problem at its source and eliminate boron exposure throughout the home [5,6].Therefore, determining the specific rocks and minerals responsible for boron contamination in groundwater is essential for developing targeted remediation strategies and protecting public health.
Previous researchers suggested evaporites and datolite hosted in diabase, as well as Early Jurassic magmatic activity as possible sources of boron [8].We believe that sandstone and black shales of the Triassic Formations, as well as numerous Early Jurassic diabase intrusions and the associated sulfides (pyrite and chalcopyrite) and calcite, could be possible sources of boron.The numerous calcite veins hosted in the Triassic rocks are also a possible source of boron.The purpose of this project is to identify the potential source(s) of boron that contaminate groundwater and wells in the Newark Basin of New Jersey and Pennsylvania.Trace and major elements were analyzed in the investigated rocks (sandstone, black shales, diabase), and statistical analysis via Principal Component Analysis (PCA) was adopted to identify the inorganic and/or organic phases with which boron is associated.Boron and its isotopic compositions were also analyzed in groundwater.

Geology and Regional Hydrogeology of the Newark Basin
The Newark Basin, situated in the northeastern United States, stands as one of the most extensively studied rift basins within the Newark supergroup [9,10] (Figure 1).Stretching over 190 km in length and up to 50 km in width, it forms an elongated half-graben (Figure 2A) [10][11][12].To its northeast and northwest, the basin is bordered by Precambrian and Early Paleozoic rocks, a continuation of the southwestern Appalachian segment of the New England upland.Meanwhile, to its southeast lie the Paleozoic and Precambrian Appalachian highlands of the Blue Ridge and Piedmont provinces [13].These Appalachian highlands resulted from the Paleozoic collisions between the North American continent and various Gondwanan fragments, including Africa, culminating in the formation of the supercontinent Pangea [14].
The Newark Basin experienced rifting during the Pangea breakup and the initiation of the Atlantic Ocean opening [15].This continental rifting prompted Triassic extensional tectonic activity, leading to the reactivation of Paleozoic NE-SW-trending major faults and the subsequent formation of a half-graben [16,17].Sitting atop the Precambrian-Paleozoic basement are Triassic-Jurassic fluvio-lacustrine formations (6000 m to 8000 m) with occasional diabase intrusions and basaltic flows (Figure 2B) [10,11]   The Newark Basin experienced rifting during the Pangea breakup and the initiation of the Atlantic Ocean opening [15].This continental rifting prompted Triassic extensional tectonic activity, leading to the reactivation of Paleozoic NE-SW-trending major faults and the subsequent formation of a half-graben [16,17].Sitting atop the Precambrian-Paleozoic basement are Triassic-Jurassic fluvio-lacustrine formations (6000 m to 8000 m) with occasional diabase intrusions and basaltic flows (Figure 2B) [10,11], listed and briefly described below.
Lockatong Formation (1000-2000 m): Dominated by cyclic organic matter-rich black shale and grey to red mudstone, siltstone, and sandstone facies with intercalations of argillaceous carbonate.
Passaic Formation (2800 m): Predominantly consisting of red mudstone and sandstone containing evaporites with minor lacustrine gray and black shale beds.
Research conducted on the Triassic-Jurassic fluvio-lacustrine rocks by Herman [18] and investigations on cores from the Newark Basin Coring Project (NBCP) by Olsen et al. [10] have unveiled the presence of numerous calcite veins within Early Mesozoic fluviolacustrine formations of Stockton, Lockatong, and Passaic.These calcite veins were Research conducted on the Triassic-Jurassic fluvio-lacustrine rocks by Herman [18] and investigations on cores from the Newark Basin Coring Project (NBCP) by Olsen et al. [10] have unveiled the presence of numerous calcite veins within Early Mesozoic fluviolacustrine formations of Stockton, Lockatong, and Passaic.These calcite veins were formed during the Triassic-Jurassic time [18] as a result of fluid mixing between heated meteoric water and cooler formational waters [12].
Besides the geological context presented above, the hydrogeological framework is important when assessing the circulation of groundwater [19].The hydrogeological context of the Newark Basin is characterized by fractured bedrock aquifers, which include the Stockton, Lockatong, and Brunswick aquifers [20].The Brunswick aquifer encompasses any formation younger than the Lockatong Formation in the Newark Basin.Geophysical investigations, including borehole imaging, have identified and classified water-bearing features (WBFs) within these aquifers into three main categories: bedding planes and layers, fracture planes, and linear intersections of bedding and fracture planes [20].
Northern areas of the Newark Basin have been glaciated, and in these areas, the overburden may be composed of glacial till, stratified drift, or other glacial deposits.Outside the glaciated areas, the saprolite can be relatively thick.The water table is typically 3 to 6 m below land surface in overburden or shallow, highly weathered bedrock.Recharge occurs from precipitation, especially where WBFs intersect the land surface.Groundwater flow dynamics differ between shallow and deep bedrock (generally greater than 5.5 m below ground surface), with the former exhibiting uniformly distributed flow aligned with the topographic slope and the latter displaying anisotropic, confined, or semi-confined flow controlled by gently dipping strata ( [20] and references therein).
Stratigraphic bedding exerts dominant control on groundwater storage and movement throughout the basin, with vertical circulation (leakage) between water-bearing zones also influenced by joints ( [20] and references therein).Transmissivity varies significantly among Newark Basin aquifers and within discrete sections of each aquifer ( [20] and references therein).Public supply and private wells are usually in the 61 to 183 m-deep range and intersect a variety of WBFs.Geophysical logging and inflatable packer testing have shown that certain WBFs contribute significantly more boron than others, and that grouting off these zones has greatly reduced boron concentrations in several wells [6].

Materials and Methods
The samples for this study were collected from the Stockton, Lockatong, and Passaic Formations, a diabase quarry, and calcite veins.The samples were extracted from outcrops as well as the Nursery and Titusville cores of the Newark Basin Continental Drilling Project (NBCDP) encompassing a total of 6770 m of core material.
From the diabase quarry, located in Moores Station (Figure 3) (Latitude: 40 • 19 ′ 29.96 ′′ N, Longitude: 74 • 54 ′ 37.71 ′′ W), diabase samples (n = 4) were collected along with the associated chalcopyrite (n = 5), pyrite (n = 5) and calcite (n = 7).Samples of Sandstone (n = 2) of the Stockton Fm., black shale (n = 3) and Limestone (n = 1) of Lockatong Fm., and sandstone (n = 2) of the Passaic Fm. were also collected (Table 1).Calcite that occurs in numerous veins was collected from the Stockton Fm. (n = 5), the Lockatong Fm. (n = 3), and the Passaic (n = 3) (Table 1).Groundwater samples were collected for boron isotopic analysis from private we in the Passaic Formation in the Hunterdon Plateau and adjacent to diabase in the Lam bertville, New Jersey area.Two well water samples were analyzed for their B isotope com positions and boron concentrations using a Nu Plasma II multicollector inductively co pled plasma mass spectrometer (MC ICP-MS) at the Isotope Ratio Mass Spectrometry L boratory, Department of Geosciences, Stony Brook University, Stony Brook, NY.Samp aliquots were processed through column chromatography to isolate only B from the sam ple following a procedure modified from Lemarchand et al. [22].Accuracy is provided b certified SRM 951 values.The standard used to correct the sample for in-run plasma-i duced mass fractionation is a 50 ppb solution of the NIST SRM 951 boric acid.In additio background was monitored using the same nitric acid and deionized water that was use Samples of diabase (n = 4), pyrite (n = 7), chalcopyrite (n = 7), and calcite (n = 7) were carefully crushed and handpicked using a binocular microscope and microdrill.Trace element analysis in diabase were performed at Actlabs using PGNAA (Prompt Gamma Neutron Activation Analysis).Boron was analyzed at Actlabs using the Fusion-Inductively Coupled Plasma Mass Spectrometry (FUS-ICPMS) method.The reader is referred to www.Actlabs.com(accessed on 28 April 2020) for a detailed description of these methods.
To determine whether pyrite, chalcopyrite, and calcite associated with diabase contain boron-bearing inclusions, grains of these minerals underwent analysis for both trace and major elements using laser ablation inductively coupled plasma-mass spectrometry (LA-ICP-MS) at the Memorial University facility in Newfoundland, Canada.Images of specific grains were utilized to direct the LA-ICP-MS analyses.An ArF Excimer LA system operating at a wavelength of 193 nm was paired with an Element XR ICP-MS.Spot analyses were conducted to target specific locations, with single spot analyses utilizing 5 pulses/s (5 Hz) and producing a crater size of 20 µm.A helium carrier gas (flow rate of 1 L/min) was utilized to transport the ablated material, which was mixed coaxially with Ar before entering the ICP torch.Concentration and detection limit calculations followed the methodology outlined by Longerich et al. [21].Reference materials (MASS-1 and NIST SRM 610) were analyzed twice at the beginning and end of each analytical session, as well as after every 20 measurements on samples.Regarding calcite found in the veins, analyses were carried out at the Pacific Centre for Isotopic and Geochemical Research (PCIGR) in the University of British Columbia using a ArF excimer Laser (193 nm, ASI Resolution M50LR), connected to a Quadrupole ICP-MS (Agilent 7700x).Calcite samples were ablated using a 120 µm spot size and a repetition rate of 8 Hz.Energy density was set to 1.8 J/cm 2 .Helium was used as carrier gas and N 2 was admixed for higher intensity.The instrument was tuned for maximum sensitivity, low oxide rates (ThO/Th < 0.3%) optimum mass bias (238/232 ~1.05), and high stability (RSD on 139 < 3%).Calcite samples ablated better but seemed to be inhomogeneous; during ablation, several phases seemed to be mobilized.During data reduction, the phases with higher Ca indicating carbonate phases were selected.Data reduction was carried out using the software program Iolite v. 2.5.with background subtraction, data integration and calibration.For the calcite samples internal as well as external calibration was performed using Ca = 38% m/m and the synthetic glass standard SRM NIST612.The synthetic glass standard SRM NIST610 and the USGS carbonate reference material MACS-3 were used as quality control.
Groundwater samples were collected for boron isotopic analysis from private wells in the Passaic Formation in the Hunterdon Plateau and adjacent to diabase in the Lambertville, New Jersey area.Two well water samples were analyzed for their B isotope compositions and boron concentrations using a Nu Plasma II multicollector inductively coupled plasma mass spectrometer (MC ICP-MS) at the Isotope Ratio Mass Spectrometry Laboratory, Department of Geosciences, Stony Brook University, Stony Brook, NY.Sample aliquots were processed through column chromatography to isolate only B from the sample following a procedure modified from Lemarchand et al. [22].Accuracy is provided by certified SRM 951 values.The standard used to correct the sample for in-run plasma-induced mass fractionation is a 50 ppb solution of the NIST SRM 951 boric acid.In addition, background was monitored using the same nitric acid and deionized water that was used to dilute the samples and SRM 951 standards, and subtracted from all analysis.The SRM 951 boric acid standard has a certified 11 B/ 10 B value of 4.0437 ± 0.003.Results were corrected for background, drift and fractionation, and final 11 B/ 10 B ratio.Seawater was run alongside all unknowns through chemical preparation and analysis.Sample aliquots were spiked with NSIT SRM 952, thoroughly mixed and then processed through column chromatography to isolate only B. The sample was then analyzed on the MC ICP-MS and the B concentration was calculated based on the measured isotopic composition and known isotope dilution through standard calculations.The SRM 952 standard has a certified 11 B/ 10 B value of 0.0532 ± 0.0005.
Principal Component Analysis (PCA) was conducted using SPSS 23, a statistical software package developed by IBM, based in Armonk, New York.PCA is widely utilized as an effective method for studies related to water quality [23,24].The goal of PCA is to identify the specific mineral phases in different rock types (sandstone, shale, limestone, and diabase) to which boron is linked.Understanding these associations is crucial for determining the sources of boron that contaminates groundwater passing through these rock formations.In this analysis, two principal components (PC1 and PC2) were extracted, and Varimax normalized rotation was applied to derive the loading values for each parameter within these components.Each principal component encompasses a set of correlated parameters with both positive and negative loadings, facilitating the interpretation of the primary geochemical processes influencing the dataset.PCA effectively identifies clusters of related elements.The first two principal components together account for a significant portion of the variance, specifically 94.091% for shale-limestone, 80.004% for sandstone, and 93.759% for diabase.PC1 alone accounts for a substantial part of the variance in each case, highlighting its major role in explaining the dataset's variability.

Geochemical Data
The results of geochemical analysis conducted on the rocks of the Triassic formations, calcite veins, diabase, and the associated minerals are provided in Tables 2-5 and the Supplementary Materials (Tables S1 and S2).Samples of Stockton and Passaic sandstone exhibit boron concentrations ranging from 23.9 to 31.S1).Abnormally high boron contents are recorded in one sample of Stockton calcite (6,980 ppm) and one Passaic calcite (12,100 ppm).Diabase samples exhibit concentrations of boron ranging from 23.9 ppm to 35.6 ppm.
Sulfides associated with diabase display lower boron concentrations, reaching up to 7.4 ppm and 3.9 ppm in chalcopyrite and pyrite, respectively (Table S2).Arsenic, which also contaminates groundwater [9,25], exhibits high concentrations ranging from 982 to 9572 ppm in pyrite.Other trace elements exhibit very low concentrations.In calcite associated with the diabase intrusion, boron concentrations vary from 5 ppm to 14.6 ppm, while lithium content ranges from <3 ppm to 107 ppm (Table S2).High concentrations of titanium (up to 26.9%) and SiO 2 (up to 25.5%) suggest the presence of Ti-and Si-rich inclusions likely derived from the diabase.

Boron Isotope
Boron isotope compositions and boron contents are reported in Table 6.δ 11 B is expressed as The analyzed groundwater flowing through the Passaic formation displays δ 11 B isotopic values ranging from 16.7 to 32.7‰.These isotopic compositions fall within those of groundwater samples collected in the Pennsylvania portion of the Newark Basin by USGS from 40 private supply wells (δ 11 B = 1.2 to 23.8‰) and overlap those of the diabase intrusion (25 to 31‰) [8].Boron contents in the analyzed water samples vary from 7.67 ppm (7670 µg/L) to 15.3 ppm (15,300 µg/L).

Geostatistical Analysis
The result of the geostatistical (PCA) analysis conducted on the rocks of the Late Triassic formations (sandstone, shale, limestone) and the diabase are reported in Figures 4-6 and Tables 7-9.

PCA in Shale-Limestone
Principal Component Analysis (PCA) for the shale-limestone dataset reveals significant correlations among various elements (Figure 4, Table 7).PC1 captures 67.857% of the variance, while PC2 explains 26.234%, together accounting for a cumulative variance of 94.091%.In PC1, major elements such as SiO 2 (0.995), Al 2 O 3 (0.993), Na 2 O (0.990) show strong positive loadings.Boron also exhibits a notable positive loading (0.888) on PC1.Total organic carbon (TOC), iron oxides, and carbonates are clustered in different quadrants with different loading values compared to boron (Figure 4, Table 7).The PCA results for the sandstone dataset show that PC1 and PC2 explain 41.652% and 38.351% of the variance, respectively, summing to a cumulative variance of 80.004%.PC1 is characterized by high loadings of Na 2 O (0.898), with boron displaying a strong positive loading (0.916) (Figure 5, Table 8).PC2 highlights high loadings for elements such as MgO (0.919) and CaO (0.903), but boron shows a moderate negative loading (−0.268).

Discussion
Boron concentrations in the analyzed water samples range from 7800 to 15,300 µg/L, which fall within the range of <500 to 18,000 µg/L observed in most water samples determined by Rockafellow-Baldoni et al. [5] and Spayd [6] (Figure 7).
The following paragraphs highlight and discuss the possible source(s) of boron that contaminate groundwater.The PCA and geochemical data from the investigated rocks (sandstone, shale, limestone, and diabase) reveal distinct associations of boron.In sandstone, boron exhibits a significant positive loading on PC1 (loading = 0.916) alongside Na 2 O (loading = 0.898), suggesting co-occurrence with sodium-bearing minerals, which may include feldspar and evaporite minerals given that the Passaic formation contains evaporites (gypsum and anhydrite).The elevated lithium content observed in groundwater [6] further corroborate the contribution of B-rich evaporite minerals.These boron-mineral phase association along with the high concentrations of boron (23.9-121 ppm) and the major oxide phase Na 2 O recorded in sandstone samples suggest the adsorption onto and or the incorporation of boron into Na-bearing minerals (Figure 5) as suggested by several authors [8,[27][28][29][30][31]. High boron concentrations are also observed in non-marine sandstone aquifers [

Discussion
Boron concentrations in the analyzed water samples range from 7800 to 15,300 µg/L, which fall within the range of <500 to 18,000 µg/L observed in most water samples determined by Rockafellow-Baldoni et al. [5] and Spayd [6] (Figure 7).The following paragraphs highlight and discuss the possible source(s) of boron that contaminate groundwater.The PCA and geochemical data from the investigated rocks In shale-limestone, boron displays strong positive loadings (0.888) with major elements like SiO 2 (loading = 0.995) and Al 2 O 3 (loading = 0.993) on PC1.These correlations indicate potential chemical affinities of boron within aluminosilicate phases (clay minerals), reinforcing its integration into silicate structures.Conversely, boron is not associated with total organic carbon (TOC), iron oxide (Fe 2 O 3 ), or carbonate phases, as evidenced by the different loadings of these phases compared to boron (Figure 4).Varying levels of boron, ranging from about 80 to 600 ppm, are also observed in unaltered red shales of the Passaic Formation ( [8] and references therein).Relatively elevated contents of B, Cr, and Zr registered for black shales suggest fixation of these elements in the inorganic fraction (shale) during diagenesis (Figure 4) [33][34][35][36][37][38].Mather and Porteous [32] argued that the abundance of clay minerals predominantly influences the presence of boron.
Organic matter-poor limestone intercalations show lower boron concentration (15.6 ppm) compared to black shales, indicating that limestone did not adsorb boron and most other trace elements except for Sr and Sc as shown in the PCA biplot (Figure 4).Based on the aforementioned discussion, it can be inferred that clay minerals present in the Triassic rock formations might have acted as sources of boron contaminating the groundwater.
In diabase, boron with high concentrations (23.9-35.6 ppm) shows a negative loading on PC1 (loading = −0.204)and a strong positive loading on PC2 (loading = 0.884).Boron is not associated with major phases such as aluminosilicates, iron oxides, or carbonates as evidenced by different loadings values of these phases (Figure 6, Table 9).Boron's association with uranium (loading = 0.816) and zinc (loading = 0.691) on PC2 suggests secondary associations possibly linked to remobilization processes (metasomatism) rather than a direct primary association with major mineral phases like aluminosilicates or carbonates.This indicates boron's remobilization from the surrounding rocks (e.g., sandstone) by diabase-related hydrothermal fluids and its subsequent incorporation into new mineral phases in the diabase.This assertion is corroborated by the high concentrations of Zn in sandstone reaching up to 220 ppm.The contribution of diabase as a source of boron is supported by boron isotopic compositions.The δ 11 B isotopic values (16.7 to 32.7‰) of groundwater are within those of the diabase intrusion (δ 11 B = 25 to 31‰, [8]), corroborating the contribution of diabase as a potential source of boron.
The remobilization of boron from the surrounding host rocks was also proposed [8].This remobilization may have caused the accumulation of hydrothermal sulfides (chalcopyrite, pyrite) and calcite within or near the diabase.Hence, it is possible that boron-rich inclusions may be incorporated within these sulfides.However, given that boron is typically not found in sulfides, the low boron concentrations (<8 ppm) in these sulfides suggest the probable absence of boron-rich micro-inclusions.Calcite, associated with diabase, exhibits boron concentrations reaching up to 14.6 ppm and lithium content up to 107 ppm.High concentrations of titanium (up to 26.9%) and silicon (up to 25.5%) suggest the presence of Ti-and Si-rich inclusions likely originating from the diabase.
Calcite veins from the Stockton, Lockatong, and Passaic formations exhibit boron contents in the range of 6.3-97.3ppm, with very high values recorded in one Stockton calcite (6980 ppm) and one Passaic calcite (12,100 ppm).These abnormally high boron concentrations suggest either an analytical artifact or the presence of boron-rich microinclusions.Hence, it is possible that these calcite veins may contribute to boron contamination of groundwater.
Based on this discussion, a model is proposed for the fixation of boron into different potential mineral phases (Figure 8).During the Late Triassic, boron was sequestered from the Late Triassic lacustrine water and subsequently fixed in various mineral phases such as clay minerals and Na-bearing minerals (Figure 8A).In the Early Jurassic, the emplacement of the diabase intrusion led to the remobilization of boron from the surrounding rocks and its incorporation into new mineral phases within the diabase (Figure 8B).Concurrently in the Early Jurassic, calcite veins were also formed within the Late Triassic rocks.As groundwater flows through the aquifers, particularly along available open spaces such as bedding planes and fractures (WBFs), it leaches boron from the mineral phases (Figure 8C).Consequently, groundwater becomes contaminated with boron.The presence of vertical joints affecting the Triassic aquifers may facilitate the vertical flow of contaminated groundwater (leakage).
The identification of potential boron sources has significant implications for policymaking and public health measures.Recognizing these sources is crucial for determining, from a geological perspective, where to focus outreach efforts to alert private well owners that they need to test their water for boron.Currently, water treatment measures are necessary to prevent exposure and the adverse health effects of boron.The NJGWS recommends a comprehensive point-of-use boron removal system incorporating reverse osmosis, followed by mixed bed deionization (or a boron-specific resin) with a circulation pump, and concluding with a neutralizer and granular activated carbon block filter before the tap.Although this approach effectively lowers boron concentrations below the U.S. EPA health advisory level, it is costly and does not provide whole-house treatment.Alternatively, physically isolating the boron-containing water-bearing zones in the well could be the most effective and protective method for reducing boron levels in the long term [5,6].Geologists must consider the potential presence of vertical joints that could facilitate the vertical movement of contaminated groundwater to maximize boron reduction by this technique.The importance of the current study lies in identifying the sources of boron, which is crucial for determining the most appropriate zones bearing boron-rich water to be targeted and ensuring the long-term reduction in boron levels in the well water.

Concluding Remarks
The study investigates potential sources of boron contamination in groundwater and wells of the Newark Basin, focusing on Triassic rock formations of Stockton, Lockatong, and Passaic, Jurassic diabase intrusions, and associated hydrothermal minerals.Geochemical and geostatistical analyses reveal that boron is associated with clay minerals in shale and with Na-bearing minerals (e.g., feldspar, evaporites) in sandstone.These inorganic fractions (clay minerals, Na-rich minerals) are potential sources of boron.Moreover, this study shows that late hydrothermal fluids related to the emplacement of diabase intrusion remobilized boron from surrounding rocks, resulting in relatively high boron concentrations ranging from 23.9 to 35.6 ppm observed in diabase, which indicates that it is a potential source of boron.This finding is supported by the isotopic composition of boron in groundwater (δ 11 B = 16.7 to 32.7‰), which aligns with that of diabase intrusions (δ 11 B = 25 to 31‰).Calcite veins hosted in the Triassic rock formations show relatively high concentrations (6.3-97.3ppm) with the possible presence of boron-rich micro-inclusions, suggesting that these calcite veins can be considered contributors to boron contamination.The leaching of boron from these potential sources resulted in an enrichment of boron in the groundwater.The flow of boron-rich groundwater is facilitated by bedding planes and fractures in the aquifers.
The number of samples used in this project is relatively limited, and further studies with a larger sample size are required to generalize and consolidate our conclusions.Although this study provides insights into potential sources of boron such as diabase intrusions and clay and Na-bearing minerals of Triassic formations and diabase, other potential sources of boron such as mineral evaporites hosted in the Passaic Formation need to be explored, and that is the focus of our next project.With the potential sources of boron identified, the focus should then be on increasing awareness of the need to test well water for boron and implementing affordable and effective boron removal techniques, such as pointof-entry and point-of-use removal systems, as well as methods to isolate or block boroncontaining water-bearing zones.A thorough evaluation of these approaches is necessary to develop robust solutions that mitigate risks associated with elevated boron levels.

20 Figure 1 .
Figure 1.Pangea during the middle Norian showing the Early Mesozoic rifting zone (shaded) and the preserved basins of the Newark Supergroup, including the Newark basin (green) (A).Early Mesozoic rift basins of eastern North America (B) (edited by Rddad [9] after Olsen et al. [10]).For interpretation of the colors in this figure's legend, the reader is referred to the web version of this article.

Figure 1 .
Figure 1.Pangea during the middle Norian showing the Early Mesozoic rifting zone (shaded) and the preserved basins of the Newark Supergroup, including the Newark basin (green) (A).Early Mesozoic rift basins of eastern North America (B) (edited by Rddad [9] after Olsen et al. [10]).For interpretation of the colors in this figure's legend, the reader is referred to the web version of this article.

Figure 2 .
Figure 2. Location of A-B cross-section (A) and a geologic cross-section A-B in the Newark Basin showing the Triassic formations and pre-Triassic basement (B) (edited by Rddad et al. [12] after Olsen et al. [10]).

Figure 2 .
Figure 2. Location of A-B cross-section (A) and a geologic cross-section A-B in the Newark Basin showing the Triassic formations and pre-Triassic basement (B) (edited by Rddad et al. [12] after Olsen et al. [10]).

Figure 3 .
Figure 3. Selected hand-specimen samples of calcite veins hosted in the Late Triassic fluvio-lacu trine formations of Stockton (A), Lockatong (B), and Passaic (C), along with those of sulfides (Py Pyrite, Cpy = Chalcopyrite) and calcite (Ca) (D) associated with the Early Jurassic diabase.The su fides, calcite, and diabase samples were collected at the Moore station quarry (E).

Figure 3 .
Figure 3. Selected hand-specimen samples of calcite veins hosted in the Late Triassic fluvio-lacustrine formations of Stockton (A), Lockatong (B), and Passaic (C), along with those of sulfides (Py = Pyrite, Cpy = Chalcopyrite) and calcite (Ca) (D) associated with the Early Jurassic diabase.The sulfides, calcite, and diabase samples were collected at the Moore station quarry (E).

Hydrology 2024, 11 , 107 11 of 20 Figure 4 .
Figure 4. PCs biplot with loading values for 25 major and trace elements in black shale and limestone, Newark Basin.

Figure 4 .
Figure 4. PCs biplot with loading values for 25 major and trace elements in black shale and limestone, Newark Basin.

Figure 5 .
Figure 5. PCs biplot with loading values for 25 major and trace elements in sandstone, Newark Basin.

Figure 5 .
Figure 5. PCs biplot with loading values for 25 major and trace elements in sandstone, Newark Basin.

Figure 6 .
Figure 6.PCs biplot with loading values for 25 major and trace elements in diabase, Newark Basin.

Figure 6 .
Figure 6.PCs biplot with loading values for 25 major and trace elements in diabase, Newark Basin.

Figure 7 .
Figure 7. Map showing boron distribution in water in the studied area, Newark Basin (after Spayd [6]).

Figure 7 .
Figure 7. Map showing boron distribution in water in the studied area, Newark Basin (after Spayd [6]).

Figure 8 .
Figure 8. Schematic model showing the three stages of boron distribution: Late Triassic phase involves fixation of boron in Na-bearing minerals.(A) Early Jurassic phase entails the emplacement of (i) the diabase intrusion and remobilization of boron from the surrounding rocks, subsequently incorporating it into the diabase and (ii) the formation of the calcite veins.(B) Current time phase

Figure 8 .
Figure 8. Schematic model showing the three stages of boron distribution: Late Triassic phase involves fixation of boron in Na-bearing minerals.(A) Early Jurassic phase entails the emplacement of (i) the diabase intrusion and remobilization of boron from the surrounding rocks, subsequently incorporating it into the diabase and (ii) the formation of the calcite veins.(B) Current time phase includes groundwater flow through available open spaces (WBFs), leaching boron from boron-rich phases (C).

Table 1 .
The coordinates of the analyzed samples for Triassic rocks (Stockton, Lockatong, and Passaic formations) and the calcite veins, as well as the depth at which they were collected.

Table 3 .
Concentrations of trace elements, including boron, in ppm for the Triassic rocks, Newark Basin.

Table 4 .
Concentrations of major elements (%) for the Early Jurassic diabase rocks, Newark Basin.

Table 5 .
Concentrations of trace elements, including boron, in ppm, for the Early Jurassic diabase rocks, Newark Basin.

Table 6 .
[26]n concentrations and boron isotope data in groundwater flowing in the Passaic Formation.The average of the seawater B isotope values (n = 31) is also given with the published value from Foster et al.[26].

Table 7 .
Principal components (PCs) with loading values for 25 major and trace elements in shale and limestone, Newark Basin.

Table 7 .
Principal components (PCs) with loading values for 25 major and trace elements in shale and limestone, Newark Basin.

Table 8 .
Principal components (PCs) with loading values for 25 major and trace elements in sandstone, Newark Basin.

Table 8 .
Principal components (PCs) with loading values for 25 major and trace elements in sandstone, Newark Basin.

Table 9 .
Principal components (PCs) with loading values for 30 major and trace elements in diabase, Newark Basin.

Table 9 .
Principal components (PCs) with loading values for 30 major and trace elements in diabase, Newark Basin. 32].