SEARCH · USGS Science

Results for “Journal Interpretation”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 91 records · Page 5Linked to original sources

USGS2012

Spatially telescoping measurements for improved characterization of groundwater-surface water interactions

The suite of measurement methods available to characterize fluxes between groundwater and surface water is rapidly growing. However, there are few studies that examine approaches to design of field investigations that include multiple methods. We propose that performing field measurements in a spatially telescoping sequence improves measurement flexibility and accounts for nested heterogeneities while still allowing for parsimonious experimental design. We applied this spatially telescoping approach in a study of ground water-surface water (GW-SW) interaction during baseflow conditions along Lucile Creek, located near Wasilla, Alaska. Catchment-scale data, including channel geomorphic indices and hydrogeologic transects, were used to screen areas of potentially significant GW-SW exchange. Specifically, these data indicated increasing groundwater contribution from a deeper regional aquifer along the middle to lower reaches of the stream. This initial assessment was tested using reach-scale estimates of groundwater contribution during baseflow conditions, including differential discharge measurements and the use of chemical tracers analyzed in a three-component mixing model. The reach-scale measurements indicated a large increase in discharge along the middle reaches of the stream accompanied by a shift in chemical composition towards a regional groundwater end member. Finally, point measurements of vertical water fluxes -- obtained using seepage meters as well as temperature-based methods -- were used to evaluate spatial and temporal variability of GW-SW exchange within representative reaches. The spatial variability of upward fluxes, estimated using streambed temperature mapping at the sub-reach scale, was observed to vary in relation to both streambed composition and the magnitude of groundwater contribution from differential discharge measurements. The spatially telescoping approach improved the efficiency of this field investigation. Beginning our assessment with catchment-scale data allowed us to identify locations of GW-SW exchange, plan measurements at representative field sites and improve our interpretation of reach-scale and point-scale measurements.

Alaska↗

USGS2006

Influence of deflection on a fold-to-fault progression: field evidence from near Marietta, South Carolina

Four periods of deformation (D1-D4) are recognized in the Lion Park Road borrow pit near Marietta, South Carolina. Although each period is characterized by distinct structures, D3 produced two structural styles (D3a, D3b) resulting from layer-parallel shortening. D3a is characterized by detachment folding at the tip of an underlying thrust. D3b is a fold-to-fault progression that was localized by east-dipping, quartz-filled gash fractures. The fold-to-fault progression demonstrates the influence of a mechanical anisotropy on ramp development. The early stages of D3b were formed by deflection of northwest-directed, layer-parallel shortening and active, down-section propagation of folds and thrusts. Following connection with a splay of basal detachment, later D3b stages resulted from up-section movement that produced kink folding and a throughgoing thrust. This up-section movement deformed and modified the geometries of older, down-section structures. Detailed mesoscopic field observations, integrated with a combination of current thrust fault models, are used to interpret the D3b fold-to-fault progression. ?? 2006 Elsevier Ltd. All rights reserved.

Journal of Structural Geology↗

USGS2019

Lithospheric signature of late Cenozoic extension in electrical resistivity structure of the Rio Grande rift, New Mexico, USA

We present electrical resistivity models of the crust and upper mantle from two‐dimensional (2‐D) inversion of magnetotelluric (MT) data collected in the Rio Grande rift, New Mexico, USA. Previous geophysical studies of the lithosphere beneath the rift identified a low‐velocity zone several hundred kilometers wide, suggesting that the upper mantle is characterized by a very broad zone of modified lithosphere. In contrast, the surface expression of the rift (e.g., high‐angle normal faults and synrift sedimentary units) is confined to a narrow region a few tens of kilometers wide about the rift axis. MT data are uniquely suited to probing the depths of the lithosphere that fill the gap between surface geology and body wave seismic tomography, namely the middle to lower crust and uppermost mantle. We model the electrical resistivity structure of the lithosphere along two east‐west trending profiles straddling the rift axis at the latitudes of 36.2 and 32.0°N. We present results from both isotropic and anisotropic 2‐D inversions of MT data along these profiles, with a strong preference for the latter in our interpretation. A key feature of the anisotropic resistivity modeling is a broad (~200‐km wide) zone of enhanced conductivity (<20 Ωm) in the middle to lower crust imaged beneath both profiles. We attribute this lower crustal conductor to the accumulation of free saline fluids and partial melt, a direct result of magmatic activity along the rift. High‐conductivity anomalies in the midcrust and upper mantle are interpreted as fault zone alteration and partial melt, respectively.

New Mexico↗

USGS2003

Icelandic-type crust

Numerous seismic studies, in particular using receiver functions and explosion seismology, have provided a detailed picture of the structure and thickness of the crust beneath the Iceland transverse ridge. We review the results and propose a structural model that is consistent with all the observations. The upper crust is typically 7 ?? 1 km thick, heterogeneous and has high velocity gradients. The lower crust is typically 15-30 ?? 5 km thick and begins where the velocity gradient decreases radically. This generally occurs at the V p ??? 6.5 km s-1 level. A low-velocity zone ??? 10 000 km2 in area and up to ??? 15 km thick occupies the lower crust beneath central Iceland, and may represent a submerged, trapped oceanic microplate. The crust-mantle boundary is a transition zone ???5 ?? 3 km thick throughout which V p increases progressively from ???7.2 to ???8.0 km s-1. It may be gradational or a zone of alternating high- and low-velocity layers. There is no seismic evidence for melt or exceptionally high temperatures in or near this zone. Isostasy indicates that the density contrast between the lower crust and the mantle is only ???90 kg m-3 compared with ???300 kg m-3 for normal oceanic crust, indicating compositional anomalies that are as yet not understood. The seismological crust is ???30 km thick beneath the Greenland-Iceland and Iceland-Faeroe ridges, and eastern Iceland, ???20 km beneath western Iceland, and ???40 km thick beneath central Iceland. This pattern is not what is predicted for an eastward-migrating plume. Low attenuation and normal V p/V s ratios in the lower crust beneath central and southwestern Iceland, and normal uppermost mantle velocities in general, suggest that the crust and uppermost mantle are subsolidus and cooler than at equivalent depths beneath the East Pacific Rise. Seismic data from Iceland have historically been interpreted both in terms of thin-hot and thick-cold crust models, both of which have been cited as supporting the plume hypothesis. This suggests that the plume model for Iceland is an a priori assumption rather than a hypothesis subject to testing. The long-extinct Ontong-Java Plateau, northwest India and Parana??, Brazil large igneous provinces, beneath which mantle plumes are not expected are all underlain by mantle low-velocity bodies similar to that beneath Iceland. A plume interpretation for the mantle anomaly beneath Iceland is thus not required.

Geophysical Journal International↗

USGS2013

2D IR spectra of cyanide in water investigated by molecular dynamics simulations

Using classical molecular dynamics simulations, the 2D infrared (IR) spectroscopy of CN − solvated in D 2 O is investigated. Depending on the force field parametrizations, most of which are based on multipolar interactions for the CN − molecule, the frequency-frequency correlation function and observables computed from it differ. Most notably, models based on multipoles for CN − and TIP3P for water yield quantitatively correct results when compared with experiments. Furthermore, the recent finding that T 1 times are sensitive to the van der Waals ranges on the CN − is confirmed in the present study. For the linear IR spectrum, the best model reproduces the full widths at half maximum almost quantitatively (13.0 cm −1 vs. 14.9 cm −1 ) if the rotational contribution to the linewidth is included. Without the rotational contribution, the lines are too narrow by about a factor of two, which agrees with Raman and IR experiments. The computed and experimental tilt angles (or nodal slopes) α as a function of the 2D IR waiting time compare favorably with the measured ones and the frequency fluctuation correlation function is invariably found to contain three time scales: a sub-ps, 1 ps, and one on the 10-ps time scale. These time scales are discussed in terms of the structural dynamics of the surrounding solvent and it is found that the longest time scale (≈10 ps) most likely corresponds to solvent exchange between the first and second solvation shell, in agreement with interpretations from nuclear magnetic resonance measurements.

Journal of Chemical Physics↗

USGS2024

Insights on arc magmatic systems drawn from natural melt inclusions and crystallization experiments at P H2O =800 MPa under oxidizing conditions

Whole rock compositions at Buldir Volcano, western Aleutian arc, record a strong, continuous trend of iron depletion with decreasing MgO, classically interpreted as a calc-alkaline liquid line of descent. In contrast, olivine-hosted melt inclusions have higher total iron (FeO * ) than whole rocks and show little change in FeO* with decreasing MgO. To investigate this discrepancy and determine the conditions required for strong iron depletion, we conducted oxygen fugacity (ƒO 2 ) buffered, water-saturated crystallization experiments at 800 MPa and ƒO 2 = QFM + 1.6 ± 0.4 (1σ⁠) (where QFM refers to the quartz-fayalite-magnetite buffer) on a high-Al, basaltic starting material modeled after a Buldir lava. Experimental conditions were informed by olivine-hosted melt inclusions that record minimum entrapment pressures as high as 570 MPa, >6 wt % H2O, and ƒO 2 of QFM + 1.4 (±0.2), making Buldir one of the most oxidized and wettest arc volcanoes documented globally. The experiments produce melts with Si-enrichment and Fe-depletion signatures characteristic of evolved, calc-alkaline magmas at the lowest MgO, although FeO * remains roughly constant over most of the experimental temperature range. Experiments saturate CrAl-spinel and olivine at 1160°C, followed by clinopyroxene and Al-spinel at 1085°C, hornblende at 1060°C, and, finally, plagioclase and magnetite between 1040°C and 960°C. Hornblende crystallization, not magnetite, generates the largest increase in SiO2 and largest decrease in FeO * in coexisting melts. Compositions of melt inclusions are consistent with experimental melts and reflect crystallization of a basaltic parent magma at high P H2O . In contrast, the whole rock compositional trends are influenced by magma mixing and phenocryst redistribution and accumulation. The crystallization experiments and natural liquids (melt inclusions and groundmass glass) from Buldir suggest that for an oxidized, hydrous primary basalt starting composition, significant Fe depletion from the melt will not occur until intermediate to late stages of magma crystallization (< ~4.5 wt % MgO). We conclude that the Buldir whole rock trend cannot be reproduced by crystallization at arc-relevant oxygen fugacities and is not a true liquid line of descent, warranting caution when interpreting volcanic trends globally.

Journal of Petrology↗

USGS2011

A perspective on nonstationarity and water management

This essay offers some perspectives on climate-related nonstationarity and water resources. Hydrologists must not lose sight of the many sources of nonstationarity, recognizing that many of them may be of much greater magnitude than those that may arise from climate change. It is paradoxical that statistical and deterministic approaches give us better insights about changes in mean conditions than about the tails of probability distributions, and yet the tails are very important to water management. Another paradox is that it is difficult to distinguish between long-term hydrologic persistence and trend. Using very long hydrologic records is helpful in mitigating this problem, but does not guarantee success. Empirical approaches, using long-term hydrologic records, should be an important part of the portfolio of research being applied to understand the hydrologic response to climate change. An example presented here shows very mixed results for trends in the size of the annual floods, with some strong clusters of positive trends and a strong cluster of negative trends. The potential for nonstationarity highlights the importance of the continuity of hydrologic records, the need for repeated analysis of the data as the time series grow, and the need for a well-trained cadre of scientists and engineers, ready to interpret the data and use those analyses to help adjust the management of our water resources.

Journal of the American Water Resources Associatio↗

USGS2008

40Ar/39Ar data for white mica, biotite, and K-feldspar samples from low-grade metamorphic rocks in the Westminster Terrane and adjacent rocks, Maryland

This report contains reduced 40 Ar/ 39 Ar data of white mica and K-feldspar mineral separates and matrix of a whole rock phyllite, all from low-grade metamorphic rocks of the Westminster terrane and adjacent strata in central Maryland. This report presents these data in a preliminary form, but in more detail than can be accommodated in todays professional journals. Also included in this report is information on the location of the samples and a brief description of the samples. The data contained herein are not interpreted in a geological context, and care should be taken by readers unfamiliar with argon isotopic data in the use of these results; many of the individual apparent ages are not geologically meaningful. This report is primarily a detailed source document for subsequent publications that will integrate these data into a geological context.

Maryland↗

USGS1987

Regional heat flow variations in the northern Michigan and Lake Superior region determined using the silica heat flow estimator

Conventional heat flow data are sparse for northern Michigan. The groundwater silica heat flow estimator expands the database sufficiently to allow regional variations in heat flow to be examined. Heat flow shows a pattern of alternating highs and lows trending ESE across the Upper Peninsula and Lake Superior. The informal names given to these features, their characteristic heat flow and inferred causes are listed: {A table is presented} The results suggest that, for the study area, regional variations in heat flow cannot be interpreted solely in terms of regional variations of the heat generation rate of basement rocks. ?? 1987.

Journal of Volcanology and Geothermal Research↗

USGS2003

Visible/near-infrared spectra of experimentally shocked plagioclase feldspars

High shock pressures cause structural changes in plagioclase feldspars such as mechanical fracturing and disaggregation of the crystal lattice at submicron scales, the formation of diaplectic glass (maskelynite), and genuine melting. Past studies of visible/ near-infrared spectra of shocked feldspars demonstrated few spectral variations with pressure except for a decrease in the depth of the absorption feature near 1250-1300 nm and an overall decrease in reflectance. New visible/near-infrared spectra (400-2500 nm) of experimentally shocked (17-56 GPa) albite- and anorthite-rich rock powders demonstrate similar trends, including the loss of minor hydrated mineral bands near 1410, 1930, 2250, and 2350 nm. However, the most interesting new observations are increases in reflectance at intermediate pressures, followed by subsequent decreases in reflectance at higher pressures. The amount of internal scattering and overall sample reflectance is controlled by the relative proportions of micro-fractures, submicron grains, diaplectic glass, and melts formed during shock metamorphism. We interpret the observed reflectance increases at intermediate pressures to result from progressively larger proportions of submicron feldspar grains and diaplectic glass. The ensuing decreases in reflectance occur after diaplectic glass formation is complete and the proportion of genuine melt inclusions increases. The pressure regimes over which these reflectance variations occur differ between albite and anorthite, consistent with thermal infrared spectra of these samples and previous studies of shocked feldspars. These types of spectral variations associated with different peak shock pressures should be considered during interpretation and modeling of visible/near-infrared remotely sensed spectra of planetary and asteroidal surfaces.

Journal of Geophysical Research E: Planets↗

USGS2001

Resolution analysis of finite fault source inversion using one- and three-dimensional Green's functions 2. Combining seismic and geodetic data

Using numerical tests for a prescribed heterogeneous earthquake slip distribution, we examine the importance of accurate Green's functions (GF) for finite fault source inversions which rely on coseismic GPS displacements and leveling line uplift alone and in combination with near-source strong ground motions. The static displacements, while sensitive to the three-dimensional (3-D) structure, are less so than seismic waveforms and thus are an important contribution, particularly when used in conjunction with waveform inversions. For numerical tests of an earthquake source and data distribution modeled after the 1994 Northridge earthquake, a joint geodetic and seismic inversion allows for reasonable recovery of the heterogeneous slip distribution on the fault. In contrast, inaccurate 3-D GFs or multiple 1-D GFs allow only partial recovery of the slip distribution given strong motion data alone. Likewise, using just the GPS and leveling line data requires significant smoothing for inversion stability, and hence, only a blurred vision of the prescribed slip is recovered. Although the half-space approximation for computing the surface static deformation field is no longer justifiable based on the high level of accuracy for current GPS data acquisition and the computed differences between 3-D and half-space surface displacements, a layered 1-D approximation to 3-D Earth structure provides adequate representation of the surface displacement field. However, even with the half-space approximation, geodetic data can provide additional slip resolution in the joint seismic and geodetic inversion provided a priori fault location and geometry are correct. Nevertheless, the sensitivity of the static displacements to the Earth structure begs caution for interpretation of surface displacements, particularly those recorded at monuments located in or near basin environments.

Journal of Geophysical Research B: Solid Earth↗

USGS2006

Aftershock decay, productivity, and stress rates in Hawaii: Indicators of temperature and stress from magma sources

We examined dozens of aftershock sequences in Hawaii in terms of Gutenberg-Richter and modified Omori law parameters. We studied p, the rate of aftershock decay; A p , the aftershock productivity, defined as the observed divided by the expected number of aftershocks; and c, the time delay when aftershock rates begin to fall. We found that for earthquakes shallower than 20 km, p values >1.2 are near active magma centers. We associate this high decay rate with higher temperatures and faster stress relaxation near magma reservoirs. Deep earthquakes near Kilauea's inferred magma transport path show a range of p values, suggesting the absence of a large, deep magma reservoir. Aftershock productivity is >4.0 for flank earthquakes known to be triggered by intrusions but is normal (0.25 to 4.0) for isolated main shocks. We infer that continuing, post-main shock stress from the intrusion adds to the main shock's stress step and causes higher A p . High A p in other zones suggests less obvious intrusions and pulsing magma pressure near Kilauea's feeding conduit. We calculate stress rates and stress rate changes from pre-main shock and aftershock rates. Stress rate increased after many intrusions but decreased after large M7–8 earthquakes. Stress rates are highest in the seismically active volcano flanks and lowest in areas far from volcanic centers. We found sequences triggered by intrusions tend to have high A p , high (>0.10 day) c values, a stress rate increase, and sometimes a peak in aftershock rate hours after the main shock. We interpret these values as indicating continuing intrusive stress after the main shock.

Hawai'i↗

USGS2016

An analysis of water data systems to inform the Open Water Data Initiative

Improving access to data and fostering open exchange of water information is foundational to solving water resources issues. In this vein, the Department of the Interior's Assistant Secretary for Water and Science put forward the charge to undertake an Open Water Data Initiative (OWDI) that would prioritize and accelerate work toward better water data infrastructure. The goal of the OWDI is to build out the Open Water Web (OWW). We therefore considered the OWW in terms of four conceptual functions: water data cataloging, water data as a service, enriching water data, and community for water data. To describe the current state of the OWW and identify areas needing improvement, we conducted an analysis of existing systems using a standard model for describing distributed systems and their business requirements. Our analysis considered three OWDI-focused use cases—flooding, drought, and contaminant transport—and then examined the landscape of other existing applications that support the Open Water Web. The analysis, which includes a discussion of observed successful practices of cataloging, serving, enriching, and building community around water resources data, demonstrates that we have made significant progress toward the needed infrastructure, although challenges remain. The further development of the OWW can be greatly informed by the interpretation and findings of our analysis.

Journal of the American Water Resources Associatio↗

USGS2012

Determination and uncertainty of moment tensors for microearthquakes at Okmok Volcano, Alaska

Efforts to determine general moment tensors (MTs) for microearthquakes in volcanic areas are often hampered by small seismic networks, which can lead to poorly constrained hypocentres and inadequate modelling of seismic velocity heterogeneity. In addition, noisy seismic signals can make it difficult to identify phase arrivals correctly for small magnitude events. However, small volcanic earthquakes can have source mechanisms that deviate from brittle double-couple shear failure due to magmatic and/or hydrothermal processes. Thus, determining reliable MTs in such conditions is a challenging but potentially rewarding pursuit. We pursued such a goal at Okmok Volcano, Alaska, which erupted recently in 1997 and in 2008. The Alaska Volcano Observatory operates a seismic network of 12 stations at Okmok and routinely catalogues recorded seismicity. Using these data, we have determined general MTs for seven microearthquakes recorded between 2004 and 2007 by inverting peak amplitude measurements of P and S phases. We computed Green's functions using precisely relocated hypocentres and a 3-D velocity model. We thoroughly assessed the quality of the solutions by computing formal uncertainty estimates, conducting a variety of synthetic and sensitivity tests, and by comparing the MTs to solutions obtained using alternative methods. The results show that MTs are sensitive to station distribution and errors in the data, velocity model and hypocentral parameters. Although each of the seven MTs contains a significant non-shear component, we judge several of the solutions to be unreliable. However, several reliable MTs are obtained for a group of previously identified repeating events, and are interpreted as compensated linear-vector dipole events.

Alaska↗

USGS2009

Trend estimation in populations with imperfect detection

1. Trends of animal populations are of great interest in ecology but cannot be directly observed owing to imperfect detection. Binomial mixture models use replicated counts to estimate abundance, corrected for detection, in demographically closed populations. Here, we extend these models to open populations and illustrate them using sand lizard Lacerta agilis counts from the national Dutch reptile monitoring scheme. 2. Our model requires replicated counts from multiple sites in each of several periods, within which population closure is assumed. Counts are described by a hierarchical generalized linear model, where the state model deals with spatio-temporal patterns in true abundance and the observation model with imperfect counts, given that true state. We used WinBUGS to fit the model to lizard counts from 208 transects with 1–10 (mean 3) replicate surveys during each spring 1994–2005. 3. Our state model for abundance contained two independent log-linear Poisson regressions on year for coastal and inland sites, and random site effects to account for unexplained heterogeneity. The observation model for detection of an individual lizard contained effects of region, survey date, temperature, observer experience and random survey effects. 4. Lizard populations increased in both regions but more steeply on the coast. Detectability increased over the first few years of the study, was greater on the coast and for the most experienced observers, and highest around 1 June. Interestingly, the population increase inland was not detectable when the observed counts were analysed without account of detectability. The proportional increase between 1994 and 2005 in total lizard abundance across all sites was estimated at 86% (95% CRI 35–151). 5. Synthesis and applications . Open-population binomial mixture models are attractive for studying true population dynamics while explicitly accounting for the observation process, i.e. imperfect detection. We emphasize the important conceptual benefit provided by temporal replicate observations in terms of the interpretability of animal counts.

Journal of Applied Ecology↗

USGS2010

Fix success and accuracy of GPS radio collars in old-growth temperate coniferous forests

Global Positioning System (GPS) telemetry is used extensively to study animal distribution and resource selection patterns but is susceptible to biases resulting from data omission and spatial inaccuracies. These data errors may cause misinterpretation of wildlife habitat selection or spatial use patterns. We used both stationary test collars and collared free-ranging American black bears ( Ursus americanus ) to quantify systemic data loss and location error of GPS telemetry in mountainous, old-growth temperate forests of Olympic National Park, Washington, USA. We developed predictive models of environmental factors that influence the probability of obtaining GPS locations and evaluated the ability of weighting factors derived from these models to mitigate data omission biases from collared bears. We also examined the effects of microhabitat on collar fix success rate and examined collar accuracy as related to elevation changes between successive fixes. The probability of collars successfully obtaining location fixes was positively associated with elevation and unobstructed satellite view and was negatively affected by the interaction of overstory canopy and satellite view. Test collars were 33% more successful at acquiring fixes than those on bears. Fix success rates of collared bears varied seasonally and diurnally. Application of weighting factors to individual collared bear fixes recouped only 6% of lost data and failed to reduce seasonal or diurnal variation in fix success, suggesting that variables not included in our model contributed to data loss. Test collars placed to mimic bear bedding sites received 16% fewer fixes than randomly placed collars, indicating that microhabitat selection may contribute to data loss for wildlife equipped with GPS collars. Horizontal collar errors of >800 m occurred when elevation changes between successive fixes were >400 m. We conclude that significant limitations remain in accounting for data loss and error inherent in using GPS telemetry in coniferous forest ecosystems and that, at present, resource selection patterns of large mammals derived from GPS telemetry should be interpreted cautiously.

Washington↗

USGS1998

Analysis of the influence of spatial pattern in habitat selection studies

Design and analysis of wildlife habitat selection studies typically do not assess the effect of spatial pattern on the habitat selection process. Effects of landscape scale pattern on habitat selection cannot be accomplished without replicate study areas, because pattern is a single, albeit multifaceted, attribute of an area. For a single area, however, the influence of pattern-related characteristics, such as shape and edge shared with adjacent patches, can be estimated by using GLIM (McCullough and Neider 1983) procedures to model patch-specific frequency counts of animal use as a function of these parameters. This approach is evaluated and illustrated with simulated breeding-bird counts in a South Carolina study area for which a GIS land cover classification is available. A related technique for evaluating whether movement from patch to patch is selective is developed and illustrated for designs that involve collection of trajectory data from monitored individuals. These designs and analyses are feasible given current GIS and GPS technology. Statistical inferences from habitat selection studies should be interpreted within the context of a range of scales at which animals differentiate between patch attributes.

Journal of Agricultural, Biological, and Environme↗

USGS2009

High-resolution sclerochronological analysis of the bivalve mollusk Saxidomus gigantea from Alaska and British Columbia: techniques for revealing environmental archives and archaeological seasonality

The butter clam, Saxidomus gigantea, is one of the most commonly recovered bivalves from archaeological shell middens on the Pacific Coast of North America. This study presents the results of the sclerochronology of modern specimens of S. gigantea, collected monthly from Pender Island (British Columbia), and additional modern specimens from the Dundas Islands (BC) and Mink and Little Takli Islands (Alaska). The methods presented can be used as a template to interpret local environmental conditions and increase the precision of seasonality estimates in shellfish using sclerochronology and oxygen isotope analysis. This method can also identify, with a high degree of accuracy, the date of shell collection to the nearest fortnightly cycle, the time of day the shell was collected and the approximate tidal elevation (i.e., approx. water depth and distance from the shoreline) from which the shell was collected. Life-history traits of S. gigantea were analyzed to understand the timing of growth line formation, the duration of the growing season, the growth rate, and the reliability of annual increments. We also examine the influence of the tidal regime and freshwater mixing in estuarine locations and how these variables can affect both incremental structures and oxygen isotope values. The results of the sclerochronological analysis show that there is a latitudinal trend in shell growth that needs to be considered when using shells for seasonality studies. Oxygen isotope analysis reveals clear annual cycles with the most positive values corresponding to the annual winter growth lines, and the most negative values corresponding to high temperatures during the summer. Intra-annual increment widths demonstrate clear seasonal oscillations with broadest increments in summer and very narrow increments or no growth during the winter months. This study provides new insights into the biology, geochemistry and seasonal growth of S. gigantea, which are crucial for paleoclimate reconstructions and interpreting seasonality patterns of past human collection.

Alaska;British Columbia↗