Search USGS⌕ Search

Geology topics

Research about Charleston, South Carolina

Source-linked reports with geographic coverage including Charleston, South Carolina.

At least 19 recordsLinked to original sources

Reply to, “Comment on ‘The 1886 Charleston, South Carolina, earthquake: Relic railroad offset reveals rupture,’ by Roger Bilham and Susan E. Hough”

We welcome this opportunity to respond to Pratt et al. (2024) (hereinafter P24). Bilham and Hough (2023) proposed a “first-cut” elastic deformation model for the 1886 earthquake, a quantitative source model constrained by identified coseismic constraints. A key observation was the measurement of a lateral offset of a railroad line south of Summerville, leading to a model with predominately dextral slip and minor convergence, from which we concluded that active faulting had raised the Penholoway Marine Terrace >6 m since ∼770 ka. P24 questioned these constraints and proposed an alternative rupture model with predominantly reverse slip. This alternative model is neither consistent with coseismic constraints nor with other geophysical data. In a revised model presented here, we recognize that uplift of the Penholoway Terrace is confined to the eastern edge of the terrace, which we conclude results from active folding and tectonic transpression centered on the dextral fault that offset the railroad in 1886.

South Carolina↗

Static and dynamic strain in the 1886 Charleston, South Carolina, earthquake

During the 1886 Mw 7.3 Charleston, South Carolina, earthquake, three railroads emanating from the city were exposed to severe shaking. Expansion joints in segmented railroad tracks are designed to allow railroad infrastructure to withstand a few parts in 10,000 of thermoelastic strain. We show that, in 1886, transient contractions exceeding this limiting value buckled rails, and transient extensions pulled rails apart. Calculated values for dynamic strain in the meizoseismal region are in reasonable agreement with those anticipated from the relation between strain and moment magnitude proposed by Barbour et al. (2021) and exceed estimated tectonic strain released by the earthquake by an order of magnitude. Almost all of the documented disturbances of railroad lines, including evidence for shortening of the rails, can thus be ascribed to the effects of dynamic strain changes, not static strain. Little or no damage to railroads was reported outside the estimated 10 −4 dynamic strain contour. The correspondence between 10 −3 and 2×10 −4 contours of dynamic strain and Mercalli intensity 9 and 8, anticipated from the dependence of each quantity on peak ground velocity, suggests it may be possible to use railroad damage to quantitatively estimate shaking intensity. At one location, near Rantowles, ≈20 km west of Charleston, a photograph of buckled track taken one day after the earthquake has been cited as evidence for shallow dextral slip and has long focused a search for a causal fault in this region. Photogrammetric analysis reveals that the buckle was caused by transient contraction of <10 cm with no dextral offset. Our results further weaken the evidence for faulting in the swamps and forests south of the Ashley River in 1886, hitherto motivated by the photograph and limited macroseismic evidence for high‐intensity shaking.

South Carolina↗

On the provenance of field reports of the 1886 Charleston, South Carolina, earthquake: A seismo-historical whodunnit

Much of what is known about the effects of the 1886 Charleston, South Carolina, earthquake throughout the epicentral region can be attributed to meticulous field investigations by an individual with training in geology and engineering, Earle Sloan ( Clendenin, 1926 ). In a recent study, Bilham and Hough (2024) undertook a detailed analysis of the effects of the earthquake on railroads in the Charleston region, drawing heavily from Sloan’s reports. This exercise identified several inconsistencies in Sloan’s field reports, including understandable measurement imprecision, inferred data entry mistakes, and transcription errors. The study also begged the question, where was Sloan at the time of the mainshock and over the following week? And to what extent did he draw from secondhand information in compiling his reports? On this question Sloan’s reports were sometimes enigmatic, lending themselves to misinterpretation in contemporaneous as well as modern interpretations. Beyond the details that were germane for, and briefly summarized by, the studies of Bilham and Hough (2023 , 2024) , in this report we don our historical seismologist caps to chronicle Sloan’s activities following the earthquake. We summarize our inferences here for the benefit of future scholars who might attempt to retrace either Sloan’s footsteps or our own. This study also serves to highlight Sloan’s singular contributions to earthquake science, which were never published separately.

South Carolina↗

The 1886 Charleston, South Carolina, earthquake: Intensities and ground motions

The 1 September 1886 Charleston, South Carolina, earthquake was one of the largest preinstrumental earthquakes in eastern North America for which extensive contemporaneous observations were documented. The distribution of shaking was mapped shortly after the earthquake, and reconsidered by several authors in the late twentieth century, but has not been reconsidered with a modern appreciation for issues associated with macroseismic data interpretation. Detailed contemporary accounts have also never been used to map the distribution of numerical shaking intensities in the near field. In this study we reconsider macroseismic data from far‐field accounts as well as detailed accounts of damage in the near field, estimating modified Mercalli intensity values at 1297 locations including over 200 definite “not felt” reports that delineate the overall felt extent. We compare the results to the suite of ground‐motion models for eastern North America selected by the National Seismic Hazard Model, using a recently proposed mainshock rupture model and an average site condition for the locations at which intensities are estimated. The comparison supports the moment magnitude estimate, 7.3, from a recently proposed rupture model ( Bilham and Hough, 2023 ). A ShakeMap constrained by model predictions and estimated intensities further illustrates this consistency, which we show is insensitive to rupture model details. Given the uncertainty of calibration relations for magnitudes close to 7, the overall intensity distribution provides a good characterization of shaking but cannot improve the independent moment magnitude estimate. We also identify a previously unrecognized early large aftershock that occurred 9–10 min after the mainshock, for which we estimate magnitude ∼5.6.

South Carolina↗

The 1886 Charleston, South Carolina, Earthquake: Relic railroad offset reveals rupture

In the absence of documented surface rupture during the 1 September 1886 Charleston earthquake, there has been considerable speculation about the location and mechanism of the causative fault. We use an inferred coseismic offset of the South Carolina Railroad and additional numerical constraints to develop an elastic deformation model—a west‐dipping fault following strands of two previously identified faults. The constraints are consistent with a blind rupture with 6.5 ± 0.3 m of dextral slip and 2 ± 0.5 m of reverse slip below 450 m depth. We propose that repeated slip on this fault has raised the Penholoway Marine Terrace >6 m since ∼770 ka. The inferred coseismic slip on the fault in an M w 7.3 earthquake is consistent with the distribution of damage in 1886.

South Carolina↗

Shallow faulting and folding in the epicentral area of the 1886 Charleston, South Carolina, earthquake

The moment magnitude ( ⁠ M w "> M w �w ⁠ ) ∼7 earthquake that struck Charleston, South Carolina, on 31 August 1886 is the largest historical earthquake in the United States east of the Appalachian Mountains. The fault(s) that ruptured during this earthquake has never been conclusively identified, and conflicting fault models have been proposed. Here we interpret reprocessed seismic reflection profiles, reprocessed legacy aeromagnetic data, and newly collected ground penetrating radar (GPR) profiles to delineate faults deforming the Cretaceous and younger Atlantic Coastal Plain (ACP) strata in the epicentral area of the 1886 earthquake. The data show evidence for faults folding or vertically displacing ACP strata, including apparent displacements of near‐surface strata (upper ∼20 m). Aeromagnetic data show several northeast (NE)‐trending lineaments, two of which correlate with faults and folds with vertical displacements as great as 55 m on the seismic reflection and radar profiles. ACP strata show only minor thickness changes across these structures, indicating that much of the displacement postdates the shallowest well‐imaged ACP strata of Eocene age. Faults imaged on the seismic reflection profiles appear on GPR profiles to displace the erosional surface at the top of the upper Eocene to Oligocene Cooper Group, including where railroad tracks were bent during the 1886 earthquake. Some faults coincide with changes in river trends, bifurcations of river channels, and unusual river meanders that could be related to recent fault motion. In contrast to our interpreted NE fault trends, earthquake locations and some focal mechanisms in the modern seismic zone have been interpreted as defining a nearly north‐striking, west‐dipping zone of aftershocks from the 1886 earthquake. The relationship between the modern seismicity and the faults we image is therefore enigmatic. However, multiple faults in the area clearly have been active since the Eocene and deform strata in the upper 20 m, providing potential targets for field‐based geologic investigations.

South Carolina↗

Urban proximity while breeding is not a predictor of perfluoroalkyl substance contamination in the eggs of brown pelicans

Identifying sources of exposure to chemical stressors is difficult when both target organisms and stressors are highly mobile. While previous studies have demonstrated that populations of some organisms proximal to urban centers may display increased burdens of human-created chemicals compared to more distal populations, this relationship may not be universal when applied to organisms and stressors capable of transboundary movements. We examined eggs of brown pelicans ( Pelecanus occidentalis ), a nearshore seabird with daily movements ranging from local to 50 km and annual migrations ranging from year-round residency to 1500 km. Thirty-six eggs from three breeding colonies located at increasing distances to a major urban center (Charleston, South Carolina, USA) were analyzed for concentrations of per - and polyfluoroalkyl substances (PFAS). Areas of high use for each colony during the breeding season were also assessed via the tracking of adult pelicans from each colony using GPS-PTT satellite transmitters and overlapped with measures of relative urbanization via land cover data. We report potentially significant ∑PFAS concentrations in the eggs of pelicans (175.4 ± 120.1 ng/g w wt. SD), driven largely by linear perfluorooctane sulfonate (n-PFOS) (48–546 ng/g w wt.). Residues of the precursor compound perfluorooctane sulfonamide (FOSA) were also present in pelican eggs, suggesting continued exposure of local wildlife beyond implemented phaseouts of some PFAS. For most analytes, egg concentrations did not exhibit a significant spatial structure despite some differentiation in high-use areas unlike similar data for another regional apex predator, the bottlenose dolphin ( Tursiops truncatus ). We suggest that the partially migratory nature of brown pelicans during the non-breeding season, combined with daily ranges that may extend to 50 km from local point sources, may have homogenized exposure across individuals. Charleston likely remains a major source for PFAS in the overall region, however, given the high concentrations observed as well as known releases of PFAS in the nearshore environment.

South Carolina↗

Effect of phytoremediation on concentrations of benzene, toluene, naphthalene, and dissolved oxygen in groundwater at a former manufactured gas plant site, Charleston, South Carolina, USA, 1998–2014

Concentrations of benzene, toluene, naphthalene, and dissolved oxygen in groundwater at a former manufactured gas plant site near Charleston, South Carolina, USA, have been monitored since the installation of a phytoremediation system of hybrid poplar trees in 1998. Between 2000 and 2014, the concentrations of benzene, toluene, and naphthalene (BT&N) in groundwater in the planted area have decreased. For example, in the monitoring well containing the highest concentrations of BT&N, benzene concentrations decreased from 10,200 &micro;g/L to less than 4000 &micro;g/L, toluene concentrations decreased from 2420 &micro;g/L to less than 20 &micro;g/L, and naphthalene concentrations decreased from 6840 &micro;g/L to less than 3000 &micro;g/L. Concentrations of BT&N in groundwater in all wells were observed to be lower during the summer months relative to the winter months of a particular year during the first few years after installing the phytoremediation system, most likely due to increased transpiration and contaminant uptake by the hybrid poplar trees during the warm summer months; this pathway of uptake by trees was confirmed by the detection of benzene, toluene, and naphthalene in trees during sampling events in 2002, and later in the study in 2012. These data suggest that the phytoremediation system affects the groundwater contaminants on a seasonal basis and, over multiple years, has resulted in a cumulative decrease in dissolved-phase contaminant concentrations in groundwater. The removal of dissolved organic contaminants from the aquifer has resulted in a lower demand on dissolved oxygen supplied by recharge and, as a result, the redox status of the groundwater has changed from anoxic to oxic conditions. This study provides much needed information for water managers and other scientists on the viability of the long-term effectiveness of phytoremediation in decreasing groundwater contaminants and increasing dissolved oxygen at sites contaminated by benzene, toluene, and naphthalene.

South Carolina↗

Shipboard surveys track magnetic sources in marine sediments--geophysical studies of the Stono and North Edisto Inlets near Charleston, South Carolina

Magnetic field data are traditionally used to analyze igneous and metamorphic rocks, but recent efforts have shown that magnetic sources within sediments may be detectable, suggesting new applications for high-resolution magnetic field surveys. Candidates for sedimentary sources include heavy mineral sand concentrations rich in magnetite or hematite, alteration-induced glauconite, or biogenic magnetite. Magnetic field surveys can be used to map the distributions of such sources with much denser and more widespread coverage than possible by sampling. These data can then provide constraints on the composition history of local sediments. Mapping such sediments requires the sensor to be relatively close to the source, and filtering approaches may be needed to distinguish signals from both system noise and deeper basement features. Marine geophysical surveys conducted in July, 2010, over the Stono and North Edisto River inlets and their riverine inputs south of Charleston, South Carolina, showed 10- to 40-m-wide, 1- to 6-nT magnetic anomalies associated with shallow, sand-covered seabed. These anomalies are distinct from system noise but are too narrow to represent basement features. The anomalies are present mostly in shallow areas where river sediments originating from upland areas enter the inlets. Surface grab samples from the North Edisto River contain trace amounts of heavy mineral sediments including hematite, maghemite, ilmenite, and magnetite, as well as garnet, epidote, zircon, and rutile. Previous stream sediment analyses show enhanced titanium over much of the Atlantic Coastal Plain. The combined data suggest that the anomalies are generated by titanium- and iron-rich heavy mineral sands ultimately originating from the Piedmont and Blue Ridge provinces, which are then reworked and concentrated by tidal currents.

South Carolina↗

A unique approach to estimating lateral anisotropy in complex geohydrologic environments

Aquifers in fractured rock or karstic settings are likely to have anisotropic transmissivity distributions. Aquifer tests that are performed in these settings also are frequently affected by leakage from adjacent confining units. Finite-difference models such as MODFLOW are convenient tools for estimating the hydraulic characteristics of the stressed aquifer and adjacent confining units but are poor tools for the estimation of lateral anisotropy. This limitation of finite-difference methods can be overcome by application of the spin method, a technique whereby the positions of the observation wells are rotated about the production well to estimate anisotropy and orientation. Formal parameter estimation is necessary to analyze aquifer tests because of the number of parameters that are estimated. As a test, transmissivity, anisotropy, and orientation were successfully estimated for a simple hypothetical problem with known properties. The technique also was applied to estimate hydraulic properties of the Santee Limestone/Black Mingo (SL/BM) aquifer and a leaky confining unit beneath Charleston, South Carolina. A 9-day aquifer test with an average discharge of 6441/min was analyzed numerically. Drawdowns in the SL/BM aquifer and confining unit were simulated with a 12-layer MODFLOW model that was discretized into 81 rows of 81 columns. Simulated drawdowns at seven observation wells that ranged from 23 to 2700 m from the production well were matched to measured drawdowns. Transmissivity estimated along the minor axis ranged from 10 to 15 m 2 /day and along the major axis ranged from 80 to 10Om 2 /day. The major axis of transmissivity was oriented along compass heading 116° (degrees clockwise from north), which agrees with geologic interpretations. Vertical hydraulic conductivity and specific storage estimates for the overlying confining unit were 4 x 10 -5 m/day and 2 x 10 -4 1/m, respectively.

South Carolina↗

Geoelectrical evidence of bicontinuum transport in groundwater

Bicontinuum models and rate-limited mass transfer (RLMT) explain complex transport behavior (e.g., long tailing and rebound) in heterogeneous geologic media, but experimental verification is problematic because geochemical samples represent the mobile component of the pore space. Here, we present geophysical evidence of RLMT at the field scale during an aquifer-storage and recovery experiment in a fractured limestone aquifer in Charleston, South Carolina. We observe a hysteretic relation between measurements of porefluid conductivity and bulk electrical conductivity; this hysteresis contradicts advective-dispersive transport and the standard petrophysical model relating pore-fluid and bulk conductivity, but can be explained by considering bicontinuum transport models that include first-order RLMT. Using a simple numerical model, we demonstrate that geoelectrical measurements are sensitive to bicontinuum transport and RLMT parameters, which are otherwise difficult to infer from direct, hydrologic measurements.

South Carolina↗

Hydrologic and geochemical evaluation of aquifer storage recovery in the Santee Limestone/Black Mingo Aquifer, Charleston, South Carolina, 1998-2002

The hydrologic and geochemical effects of aquifer storage recovery were evaluated to determine the potential for supplying the city of Charleston, South Carolina, with large quantities of potable water during emergencies, such as earthquakes, hurricanes, or hard freezes. An aquifer storage recovery system, including a production well and three observation wells, was installed at a site located on the Charleston peninsula. The focus of this study was the 23.2-meter thick Tertiary-age carbonate and sand aquifer of the Santee Limestone and the Black Mingo Group, the northernmost equivalent of the Floridan aquifer system. Four cycles of injection, storage, and recovery were conducted between October 1999 and February 2002. Each cycle consisted of injecting between 6.90 and 7.19 million liters of water for storage periods of 1, 3, or 6 months. The volume of recovered water that did not exceed the U.S. Environmental Protection Agency secondary standard for chloride (250 milligrams per liter) varied from 1.48 to 2.46 million liters, which is equivalent to 21 and 34 percent of the total volume injected for the individual tests. Aquifer storage recovery testing occurred within two productive zones of the brackish Santee Limestone/Black Mingo aquifer. The individual productive zones were determined to be approximately 2 to 4 meters thick, based on borehole geophysical logs, electromagnetic flow-meter testing, and specific-conductance profiles collected within the observation wells. A transmissivity and storage coefficient of 37 meters squared per day and 3 x 10-5, respectively, were determined for the Santee Limestone/Black Mingo aquifer. Water-quality and sediment samples collected during this investigation documented baseline aquifer and injected water quality, aquifer matrix composition, and changes in injected/aquifer water quality during injection, storage, and recovery. A total of 193 water-quality samples were collected and analyzed for physical properties, major and minor ions, and nutrients. The aquifer and treated surface water were sodiumchloride and calcium/sodium-bicarbonate water types, respectively. Forty-five samples were collected and analyzed for total trihalomethane. Total trihalomethane data collected during aquifer storage recovery cycle 4 indicated that this constituent would not restrict the use of recovered water for drinking-water purposes. Analysis of six sediment samples collected from a cored well located near the aquifer storage recovery site showed that quartz and calcite were the dominant minerals in the Santee Limestone/Black Mingo aquifer. Estimated cation exchange capacity ranged from 12 to 36 milliequivalents per 100 grams in the lower section of the aquifer. A reactive transport model was developed that included two 2-meter thick layers to describe each of the production zones. The four layers composing the production zones were assigned porosities ranging from 0.1 to 0.3 and hydraulic conductivities ranging from 1 to 8.4 meters per day. Specific storage of the aquifer and confining units was estimated to be 1.5 x 10-5 meter-1. Longitudinal dispersivity of all layers was specified to be 0.5 meter. Leakage through the confining unit was estimated to be minimal and, therefore, not used in the reactive transport modeling. Inverse geochemical modeling indicates that mixing, cation exchange, and calcite dissolution are the dominant reactions that occur during aquifer storage recovery testing in the Santee Limestone/Black Mingo aquifer. Potable water injected into the Santee Limestone/Black Mingo aquifer evolved chemically by mixing with brackish background water and reaction with calcite and cation exchangers in the sediment. Reactive-transport model simulations indicated that the calcite and exchange reactions could be treated as equilibrium processes.

South Carolina↗

Intraplate triggered earthquakes: Observations and interpretation

We present evidence that at least two of the three 1811-1812 New Madrid, central United States, mainshocks and the 1886 Charleston, South Carolina, earthquake triggered earthquakes at regional distances. In addition to previously published evidence for triggered earthquakes in the northern Kentucky/southern Ohio region in 1812, we present evidence suggesting that triggered events might have occurred in the Wabash Valley, to the south of the New Madrid Seismic Zone, and near Charleston, South Carolina. We also discuss evidence that earthquakes might have been triggered in northern Kentucky within seconds of the passage of surface waves from the 23 January 1812 New Madrid mainshock. After the 1886 Charleston earthquake, accounts suggest that triggered events occurred near Moodus, Connecticut, and in southern Indiana. Notwithstanding the uncertainty associated with analysis of historical accounts, there is evidence that at least three out of the four known Mw 7 earthquakes in the central and eastern United States seem to have triggered earthquakes at distances beyond the typically assumed aftershock zone of 1-2 mainshock fault lengths. We explore the possibility that remotely triggered earthquakes might be common in low-strain-rate regions. We suggest that in a low-strain-rate environment, permanent, nonelastic deformation might play a more important role in stress accumulation than it does in interplate crust. Using a simple model incorporating elastic and anelastic strain release, we show that, for realistic parameter values, faults in intraplate crust remain close to their failure stress for a longer part of the earthquake cycle than do faults in high-strain-rate regions. Our results further suggest that remotely triggered earthquakes occur preferentially in regions of recent and/or future seismic activity, which suggests that faults are at a critical stress state in only some areas. Remotely triggered earthquakes may thus serve as beacons that identify regions of long-lived stress concentration.

South Carolina↗

Structural and tectonic setting of the Charleston, South Carolina, region: Evidence from the Tertiary stratigraphic record

Eleven upper Eocene through Pliocene stratigraphic units occur in the subsurface of the region surrounding Charleston, South Carolina. These units contain a wealth of information concerning the long-term tectonic and structural setting of that area. These stratigraphic units have a mosaic pattern of distribution, rather than a simple layered pattern, because deposition, erosion, and tectonic warping have interacted in a complex manner through time. By generating separate structure-contour maps for the base of each stratigraphic unit, an estimate of the original basal surface of each unit can be reconstructed over wide areas. Changes in sea level over geologic time generate patterns of deposition and erosion that are geographically unique for the time of each transgression. Such patterns fail to persist when compared sequentially over time. In some areas, however, there has been persistent, repetitive net downward of upward movement over the past 34 m.y. These repetitive patterns of persistent motion are most readily attributable to tectonism. The spatial pattern of these high and low areas is complex, but it appears to correlate well with known tectonic features of the region. This correlation suggests that the tectonic setting of the Charleston region is controlled by scissors-like compression on a crustal block located between the north-trending Adams Run fault and the northwest-trending Charleston fault. Tectonism is localized in the Charleston region because it lies within a discrete hinge zone that accommodates structural movement between the Cape Fear arch and the Southeast Georgia embayment.

South Carolina↗

Magnitude and extent of arsenic and thallium concentrations in ground water and sediments at the Charleston Naval Complex, North Charleston, South Carolina, 1994-99

Water-quality samples were collected quarterly during 1994-99 from 604 wells screened in the surficial aquifer system beneath the Charleston Naval Complex, North Charleston, South Carolina. Arsenic and thallium were selected for analysis because concentrations of these metals in some wells consistently exceeded the established (2001) drinking water maximum contaminant levels of 10 and 2 micrograms per liter, respectively. The analysis was conducted to determine the magnitude and spatial distribution of arsenic and thallium in ground water at the Charleston Naval Complex and to quantify arsenic and thallium concentrations in a dated sediment core from Shipyard Creek marsh near the southern boundary of the Naval Complex. The surficial aquifer system beneath the Charleston Naval Complex consists of an unconfined upper surficial aquifer and a confined lower surficial aquifer. Hydraulic connection between the two aquifers is limited or nonexistent throughout the system at the Naval Complex. The Charleston Naval Complex is divided into nine operational units designated as zones A through I. Arsenic and thallium concentration data were compiled and interpreted for the two surficial aquifers within each zone. Mean arsenic (n=603) and thallium (n=604) concentrations were calculated for water samples from each well screened in the upper and lower surficial aquifers. In the upper surficial aquifer, mean arsenic concentrations ranged from 0.9 to 339 micrograms per liter and exceeded 10 micrograms per liter in 29 percent of the wells. In the lower surficial aquifer, mean arsenic concentrations ranged from 1.0 to 97.4 micrograms per liter and exceeded 10 micrograms per liter in 23 percent of the wells. The greatest number of water samples with mean arsenic concentrations exceeding 10 micrograms per liter were collected from wells in the upper surficial aquifer at zone E in the northwestern part of the study area. Well clusters, defined as three or more wells in a solid-waste management unit or area of concern, where the mean arsenic concentration exceeded 10 micrograms per liter, were identified in association with 12 sites in the upper surficial aquifer-solid-waste management unit 039 (a drum-storage area) in zone A; solid-waste management units 044 (coal-storage area) and 047 (burning dump) in zone C; solid-waste management unit 065 (lead-storage area) and area of concern 556 (dry docks 3 and 4) in zone E; areas of concern 609 (building 1346 gas station) and 613 (locomotive shop) in zone F; solid-waste management units 006 (public works storage yard) and 008 (oil sludge pit), and area of concern 709 (fuel-delivery system wells 12, 13, and 14) in zone G; and solid-waste management units 009 (closed landfill) and 196 (south landfill) in zone H. One well cluster was identified in the lower surficial aquifer in association with solidwaste management unit 009 (closed landfill) in zone H. Mean thallium concentrations in water from all wells ranged from less than 1.6 to 32.6 micrograms per liter in water samples from the upper surficial aquifer, and from less than 1.6 to 67.7 micrograms per liter in water samples from the lower surficial aquifer. Mean thallium concentrations equal to or greater than 10 micrograms per liter were present in water samples from 21 of 604 wells (3.5 percent). Of the 21 wells, 14 wells were located at solid-waste management unit 009 (closed landfill) in zone H near Shipyard Creek, 8 wells in the upper aquifer, and 6 wells in the lower aquifer. One well cluster where thallium exceeded 10 micrograms per liter was identified in association with solid-waste management unit 009 (closed landfill) in the upper surficial aquifer. Mean arsenic and thallium concentrations in water were calculated for all wells screened in one aquifer and located in a single zone, and are referred to as zone mean concentration in this report. Zone mean arsenic concentrations in all nine zones ranged from 3.2 to 18 micrograms per liter in water samples from the upper surficial aquifer and from 2.7 to 22 micrograms per liter in water samples from the lower surficial aquifer. Zone mean thallium concentrations in all nine zones ranged from 3.2 to 13 micrograms per liter in water samples from the upper surficial aquifer and from 3.2 to 14 micrograms per liter in water samples from the lower surficial aquifer. Ground-water samples rarely had elevated (equal to or greater than 10 micrograms per liter) concentrations of both arsenic and thallium. Water samples had coincident elevated arsenic and thallium concentrations in 10 wells in zone H, 1 well in zone A, and 1 well in zone B. Sediment quality at Shipyard Creek marsh was investigated by collecting an 11.8-foot -long sediment core (SYC-1) adjacent to zone I. The mean arsenic concentration in sediment samples from SYC-1 (n= 160) was 3.05 milligrams per kilogram plus or minus 0.92. The mean arsenic concentration and standard deviation calculated for SYC-1 sediment samples fall within the standard error for the background mean arsenic concentration reported for South Carolina sediments (1.5 milligrams per kilogram plus or minus 2. 7). All but one sample (core depth =50 inches) was less than the threshold-effects level of 7.24 milligrams per kilogram. Acid extracts of the sediment samples were analyzed for thallium concentration, but none were detected. These data indicate no obvious change in arsenic or thallium concentrations with depth in the core.

South Carolina↗