USGS ScienceSearch

SEARCH · USGS Science

Results for “Software X”

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

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

At least 19 recordsLinked to original sources

Metadata Wizard: an easy-to-use tool for creating FGDC-CSDGM metadata for geospatial datasets in ESRI ArcGIS Desktop

Creating compliant metadata for scientific data products is mandated for all federal Geographic Information Systems professionals and is a best practice for members of the geospatial data community. However, the complexity of the The Federal Geographic Data Committee’s Content Standards for Digital Geospatial Metadata, the limited availability of easy-to-use tools, and recent changes in the ESRI software environment continue to make metadata creation a challenge. Staff at the U.S. Geological Survey Fort Collins Science Center have developed a Python toolbox for ESRI ArcDesktop to facilitate a semi-automated workflow to create and update metadata records in ESRI’s 10.x software. The U.S. Geological Survey Metadata Wizard tool automatically populates several metadata elements: the spatial reference, spatial extent, geospatial presentation format, vector feature count or raster column/row count, native system/processing environment, and the metadata creation date. Once the software auto-populates these elements, users can easily add attribute definitions and other relevant information in a simple Graphical User Interface. The tool, which offers a simple design free of esoteric metadata language, has the potential to save many government and non-government organizations a significant amount of time and costs by facilitating the development of The Federal Geographic Data Committee’s Content Standards for Digital Geospatial Metadata compliant metadata for ESRI software users. A working version of the tool is now available for ESRI ArcDesktop, version 10.0, 10.1, and 10.2 (downloadable at http:/www.sciencebase.gov/metadatawizard ).

Open-File Report

A Robot Operating System (ROS) package for mapping flow fields in rivers via Particle Image Velocimetry (PIV)

Non-contact, remote sensing approaches to measuring flow velocities in river channels are widely used, but typical workflows involve acquiring images in the field and then processing data later in the office. To reduce latency between acquisition and output, with the ultimate goal of enabling real-time image velocimetry, we developed a Robot Operating System (ROS) package for Particle Image Velocimetry (PIV) that can be deployed on an embedded computer aboard an uncrewed aircraft system (UAS). The ROSPIV package consists of a series of nodes that can be run in parallel and comprise an end-to-end PIV workflow. Software development involved converting MATLAB code to C++, organizing files within a catkin workspace, and building nodes using catkin_make . The codebase is available via a repository that includes a user’s guide and demo script. This paper describes the nodes in the ROSPIV package as well as functions for preparing inputs, facilitating code generation, and visualizing PIV output. To illustrate the application of the software, we present two examples, one based on a simulated image sequence and the other based on data acquired from a UAS. For the simulated data, the velocity field derived via the ROSPIV package closely matched the known flow field used to generate the image sequence. Using real data as input demonstrated the ability of the ROSPIV package to ingest and pre-process raw images. Our initial results suggest that the ROSPIV package could become a viable approach for mapping river surface velocities in real time.

Software X

Velocity, water-quality, and bathymetric surveys of the Grays Landing and Maxwell Navigation Pools, and Selected Tributaries to the Monongahela River, Pennsylvania, 2010–11

The U.S. Geological Survey (USGS) conducted velocity, water-quality, and bathymetric surveys from spring 2010 to summer 2011 in the Grays Landing and Maxwell navigation pools of the Monongahela River, Pennsylvania, and selected tributaries in response to elevated levels of total dissolved solids (TDS) recorded in early September 2009. Velocity data were collected using an Acoustic Doppler Current Profiler. Water-quality surveys included the in-situ collection of specific-conductance, water-temperature, and turbidity data using a water-quality sonde. Additionally, discrete water samples were collected and analyzed for TDS, chloride, and sulfate. Bathymetric data were collected using an echo sounder, and the shoreline was delineated using a laser range finder and electronic compass. The data were geo-referenced using a differential global positioning system and navigational software. Horizontal (x, y) coordinates were referenced to the North American Datum of 1983. Depth (z) elevations were referenced to the North American Vertical Datum of 1988. The data are provided in electronic format (appendix 1) and may be downloaded and can be used in a geographic information system for cartographic display and data analysis.

Pennsylvania

3D Visualization of Earthquake Focal Mechanisms Using ArcScene

We created a new tool, 3D Focal Mechanisms (3DFM), for viewing earthquake focal mechanism symbols three dimensionally. This tool operates within the Environmental Systems Research Institute (ESRI®) GIS software ArcScene® 9.x. The program requires as input a GIS point dataset of earthquake locations containing strike, dip, and rake values for a nodal plane of each earthquake. Other information, such as depth and magnitude of the earthquake, may also be included in the dataset. By default for each focal point, 3DFM will create a black and white sphere or “beach ball” that is oriented based on the strike, dip, and rake values. If depth values for each earthquake are included, the focal symbol will also be placed at its appropriate location beneath the Earth's surface. In addition to the default settings, there are several other options in 3DFM that can be adjusted. The appearance of the symbols can be changed by (1) creating rings around the fault planes that are colored based on magnitude, (2) showing only the fault planes instead of a sphere, (3) drawing a flat disc that identifies the primary nodal plane, (4) or by displaying the null, pressure, and tension axes. The size of the symbols can be changed by adjusting their diameter, scaling them based on the magnitude of the earthquake, or scaling them by the estimated size of the rupture patch based on earthquake magnitude. It is also possible to filter the data using any combination of the strike, dip, rake, magnitude, depth, null axis plunge, pressure axis plunge, tension axis plunge, or fault type values of the points. For a large dataset, these filters can be used to create different subsets of symbols. Symbols created by 3DFM are stored in graphics layers that appear in the ArcScene® table of contents. Multiple graphics layers can be created and saved to preserve the output from different symbol options.

Data Series

Mineralogical characterization of strata of the Meade Peak phosphatic shale member of the Permian Phosphoria Formation: Channel and individual rock samples of measured section J and their relationship to measured sections A and B, central part of Rasmussen Ridge, Caribou County, Idaho

The Permian Phosphoria Formation of southeastern Idaho hosts one of the largest phosphate deposits in the world. Despite the economic significance of this Formation, the fine-grained nature of the phosphorite has discouraged detailed mineralogical characterization and quantification studies. Recently, selenium and other potentially toxic trace elements in mine wastes have drawn increased attention to this formation, and motivated additional study. This study uses powder X-ray diffraction (XRD), with Rietveld quantification software, to quantify and characterize the mineralogy of composite channel samples and individual samples collected from the stratigraphic sections measured by the U.S. Geological Survey in the Meade Peak Member of the Permian Phosphoria Formation at the Enoch Valley mine on Rasmussen Ridge, approximately 15 miles northeast of Soda Springs, Idaho.

Idaho

Parallel Fortran-MPI software for numerical inversion of the Laplace transform and its application to oscillatory water levels in groundwater environments

A parallel Fortran-MPI (Message Passing Interface) software for numerical inversion of the Laplace transform based on a Fourier series method is developed to meet the need of solving intensive computational problems involving oscillatory water level's response to hydraulic tests in a groundwater environment. The software is a parallel version of ACM (The Association for Computing Machinery) Transactions on Mathematical Software (TOMS) Algorithm 796. Running 38 test examples indicated that implementation of MPI techniques with distributed memory architecture speedups the processing and improves the efficiency. Applications to oscillatory water levels in a well during aquifer tests are presented to illustrate how this package can be applied to solve complicated environmental problems involved in differential and integral equations. The package is free and is easy to use for people with little or no previous experience in using MPI but who wish to get off to a quick start in parallel computing. ?? 2004 Elsevier Ltd. All rights reserved.

Environmental Modelling and Software

Incorporating temporal variation in seabird telemetry data: time variant kernel density models

A key component of the Mid-Atlantic Baseline Studies project was tracking the individual movements of focal marine bird species (Red-throated Loon [Gavia stellata], Northern Gannet [Morus bassanus], and Surf Scoter [Melanitta perspicillata]) through the use of satellite telemetry. This element of the project was a collaborative effort with the Department of Energy (DOE), Bureau of Ocean Energy Management (BOEM), the U.S. Fish and Wildlife Service (USFWS), and Sea Duck Joint Venture (SDJV), among other organizations. Satellite telemetry is an effective and informative tool for understanding individual animal movement patterns, allowing researchers to mark an individual once, and thereafter follow the movements of the animal in space and time. Aggregating telemetry data from multiple individuals can provide information about the spatial use and temporal movements of populations. Tracking data is three dimensional, with the first two dimensions, X and Y, ordered along the third dimension, time. GIS software has many capabilities to store, analyze and visualize the location information, but little or no support for visualizing the temporal data, and tools for processing temporal data are lacking. We explored several ways of analyzing the movement patterns using the spatiotemporal data provided by satellite tags. Here, we present the results of one promising method: time-variant kernel density analysis (Keating and Cherry, 2009). The goal of this chapter is to demonstrate new methods in spatial analysis to visualize and interpret tracking data for a large number of individual birds across time in the mid-Atlantic study area and beyond. In this chapter, we placed greater emphasis on analytical methods than on the behavior and ecology of the animals tracked. For more detailed examinations of the ecology and wintering habitat use of the focal species in the midAtlantic, see Chapters 20-22.

Report

Geologic map of the Chewelah 30' x 60' Quadrangle, Washington and Idaho

This data set maps and describes the geology of the Chewelah 30' X 60' quadrangle, Washington and Idaho. Created using Environmental Systems Research Institute's ARC/INFO software, the data base consists of the following items: (1) a map coverage containing geologic contacts and units, (2) a point coverage containing site-specific geologic structural data, (3) two coverages derived from 1:100,000 Digital Line Graphs (DLG); one of which represents topographic data, and the other, cultural data, (4) two line coverages that contain cross-section lines and unit-label leaders, respectively, and (5) attribute tables for geologic units (polygons), contacts (arcs), and site-specific data (points). In addition, the data set includes the following graphic and text products: (1) A PostScript graphic plot-file containing the geologic map, topography, cultural data, and two cross sections, and on a separate sheet, a Correlation of Map Units (CMU) diagram, an abbreviated Description of Map Units (DMU), modal diagrams for granitic rocks, an index map, a regional geologic and structure map, and a key for point and line symbols; (2) PDF files of the Readme text-file and expanded Description of Map Units (DMU), and (3) this metadata file. The geologic map database contains original U.S. Geological Survey data generated by detailed field observation and by interpretation of aerial photographs. The map was compiled from geologic maps of eight 1:48,000 15' quadrangle blocks, each of which was made by mosaicing and reducing the four constituent 7.5' quadrangles. These 15' quadrangle blocks were mapped chiefly at 1:24,000 scale, but the detail of the mapping was governed by the intention that it was to be compiled at 1:48,000 scale. The compilation at 1:100,000 scale entailed necessary simplification in some areas and combining of some geologic units. Overall, however, despite a greater than two times reduction in scale, most geologic detail found on the 1:48,000 maps is retained on the 1:100,000 map. Geologic contacts across boundaries of the eight constituent quadrangles required minor adjustments, but none significant at the final 1:100,000 scale. The geologic map was compiled on a base-stable cronoflex copy of the Chewelah 30' X 60' topographic base and then scribed. The scribe guide was used to make a 0.007 mil-thick blackline clear-film, which was scanned at 1200 DPI by Optronics Specialty Company, Northridge, California. This image was converted to vector and polygon GIS layers and minimally attributed by Optronics Specialty Company. Minor hand-digitized additions were made at the USGS. Lines, points, and polygons were subsequently edited at the USGS by using standard ARC/INFO commands. Digitizing and editing artifacts significant enough to display at a scale of 1:100,000 were corrected. Within the database, geologic contacts are represented as lines (arcs), geologic units as polygons, and site-specific data as points. Polygon, arc, and point attribute tables (.pat, .aat, and .pat, respectively) uniquely identify each geologic datum.

Miscellaneous Field Studies Map

ArcCN-Runoff: An ArcGIS tool for generating curve number and runoff maps

The development and the application of ArcCN-Runoff tool, an extension of ESRI@ ArcGIS software, are reported. This tool can be applied to determine curve numbers and to calculate runoff or infiltration for a rainfall event in a watershed. Implementation of GIS techniques such as dissolving, intersecting, and a curve-number reference table improve efficiency. Technical processing time may be reduced from days, if not weeks, to hours for producing spatially varied curve number and runoff maps. An application example for a watershed in Lyon County and Osage County, Kansas, USA, is presented. ?? 2004 Elsevier Ltd. All rights reserved.

Environmental Modelling and Software

Kendall-Theil Robust Line (KTRLine--version 1.0)-A Visual Basic Program for Calculating and Graphing Robust Nonparametric Estimates of Linear-Regression Coefficients Between Two Continuous Variables

The Kendall-Theil Robust Line software (KTRLine-version 1.0) is a Visual Basic program that may be used with the Microsoft Windows operating system to calculate parameters for robust, nonparametric estimates of linear-regression coefficients between two continuous variables. The KTRLine software was developed by the U.S. Geological Survey, in cooperation with the Federal Highway Administration, for use in stochastic data modeling with local, regional, and national hydrologic data sets to develop planning-level estimates of potential effects of highway runoff on the quality of receiving waters. The Kendall-Theil robust line was selected because this robust nonparametric method is resistant to the effects of outliers and nonnormality in residuals that commonly characterize hydrologic data sets. The slope of the line is calculated as the median of all possible pairwise slopes between points. The intercept is calculated so that the line will run through the median of input data. A single-line model or a multisegment model may be specified. The program was developed to provide regression equations with an error component for stochastic data generation because nonparametric multisegment regression tools are not available with the software that is commonly used to develop regression models. The Kendall-Theil robust line is a median line and, therefore, may underestimate total mass, volume, or loads unless the error component or a bias correction factor is incorporated into the estimate. Regression statistics such as the median error, the median absolute deviation, the prediction error sum of squares, the root mean square error, the confidence interval for the slope, and the bias correction factor for median estimates are calculated by use of nonparametric methods. These statistics, however, may be used to formulate estimates of mass, volume, or total loads. The program is used to read a two- or three-column tab-delimited input file with variable names in the first row and data in subsequent rows. The user may choose the columns that contain the independent (X) and dependent (Y) variable. A third column, if present, may contain metadata such as the sample-collection location and date. The program screens the input files and plots the data. The KTRLine software is a graphical tool that facilitates development of regression models by use of graphs of the regression line with data, the regression residuals (with X or Y), and percentile plots of the cumulative frequency of the X variable, Y variable, and the regression residuals. The user may individually transform the independent and dependent variables to reduce heteroscedasticity and to linearize data. The program plots the data and the regression line. The program also prints model specifications and regression statistics to the screen. The user may save and print the regression results. The program can accept data sets that contain up to about 15,000 XY data points, but because the program must sort the array of all pairwise slopes, the program may be perceptibly slow with data sets that contain more than about 1,000 points.

Techniques and Methods

Organic geochemistry and toxicology of a stream impacted by unconventional oil and gas wastewater disposal operations

Water and sediment extracts samples were analyzed for extractable hydrocarbons by gas chromatography/mass spectrometry (GC/MS) using an Agilent (Agilent Technologies, Palo Alto, CA, USA) 7890 series GC and 5975 electron ionization (EI) mass selective detector (MSD) operated in scan mode. Agilent ChemStation software was used for data acquisition and analysis (version E.02.00.493 on GC/MS computer and version F.01.03.2357 on laptop for data workup). A 30 m x 250 m x 0.25 m HP-5MS column (95% dimethyl 5% diphenyl polydimethylsiloxane) was used for GC/MS under the following conditions: 1.0 L splitless injection, constant flow of 0.7 mL/min, solvent delay of 7.5 min, injector temperature of 280C, interface at 300C, temperature program of 50-150C at 7C/min, 150-230C at 6C/min, and 230-300C at 3C/min with mass scanned from 35-500 Da.

West Virginia

Illuminating Northern California’s Active Faults

Newly acquired light detection and ranging (lidar) topographic data provide a powerful community resource for the study of landforms associated with the plate boundary faults of northern California (Figure 1). In the spring of 2007, GeoEarthScope, a component of the EarthScope Facility construction project funded by the U.S. National Science Foundation, acquired approximately 2000 square kilometers of airborne lidar topographic data along major active fault zones of northern California. These data are now freely available in point cloud (x, y, z coordinate data for every laser return), digital elevation model (DEM), and KMZ (zipped Keyhole Markup Language, for use in Google EarthTM and other similar software) formats through the GEON OpenTopography Portal (http://www.OpenTopography.org/data). Importantly, vegetation can be digitally removed from lidar data, producing high-resolution images (0.5- or 1.0-meter DEMs) of the ground surface beneath forested regions that reveal landforms typically obscured by vegetation canopy (Figure 2)

California

Development of a graphical user interface for the global land information system (GLIS)

The process of developing a Motif Graphical User Interface for the Global Land Information System (GLIS) involved incorporating user requirements, in-house visual and functional design requirements, and Open Software Foundation (OSF) Motif style guide standards. Motif user interface windows have been developed using the software to support Motif window functions war written using the C programming language. The GLIS architecture was modified to support multiple servers and remote handlers running the X Window System by forming a network of servers and handlers connected by TCP/IP communications. In April 1993, prior to release the GLIS graphical user interface and system architecture modifications were test by developers and users located at the EROS Data Center and 11 beta test sites across the country.

Pecora 12 Symposium

Bathymetric Surveys of Lake Arthur and Raccoon Lake, Pennsylvania, June 2007

In spring of 2007, bathymetric surveys of two Pennsylvania State Park lakes were performed to collect accurate data sets of lake-bed elevations and to develop methods and techniques to conduct similar surveys across the state. The lake-bed elevations and associated geographical position data can be merged with land-surface elevations acquired through Light Detection and Ranging (LIDAR) techniques. Lake Arthur in Butler County and Raccoon Lake in Beaver County were selected for this initial data-collection activity. In order to establish accurate water-surface elevations during the surveys, benchmarks referenced to NAVD 88 were established on land at each lake by use of differential global positioning system (DGPS) surveys. Bathymetric data were collected using a single beam, 210 kilohertz (kHz) echo sounder and were coupled with the DGPS position data utilizing a computer software package. Transects of depth data were acquired at predetermined intervals on each lake, and the shoreline was delineated using a laser range finder and compass module. Final X, Y, Z coordinates of the geographic positions and lake-bed elevations were referenced to NAD 83 and NAVD 88 and are available to create bathymetric maps of the lakes.

Data Series

Regional geologic map of San Andreas and related faults in Carrizo Plain, Temblor, Caliente and La Panza Ranges and vicinity, California: A digital database

This Open-File Report is a digital geologic map database. The report serves to introduce and describe the digital data. There is no paper map included in the Open-File Report. The report includes PostScript and PDF plot files that can be used to plot images of the geologic map sheet and explanation sheet. This digital map database is prepared from a previously published map by Dibblee (1973). The geologic map database delineates map units that are identified by general age, lithology, and clast size following the stratigraphic nomenclature of the U.S. Geological Survey. For descriptions of the units, their stratigraphic relations, and sources of geologic mapping, consult the explanation sheet (of99-14_4b.ps or of99-14_4d.pdf), or the original published paper map (Dibblee, 1973). The scale of the source map limits the spatial resolution (scale) of the database to 1:125,000 or smaller. For those interested in the geology of Carrizo Plain and vicinity who do not use an ARC/INFO compatible Geographic Information System (GIS), but would like to obtain a paper map and explanation, PDF and PostScript plot files containing map images of the data in the digital database, as well as PostScript and PDF plot files of the explanation sheet and explanatory text, have been included in the database package (please see the section 'Digital Plot Files', page 5). The PostScript plot files require a gzip utility to access them. For those without computer capability, we can provide users with the PostScript or PDF files on tape that can be taken to a vendor for plotting. Paper plots can also be ordered directly from the USGS (please see the section 'Obtaining Plots from USGS Open-File Services', page 5). The content and character of the database, methods of obtaining it, and processes of extracting the map database from the tar (tape archive) file are described herein. The map database itself, consisting of six ARC/INFO coverages, can be obtained over the Internet or by magnetic tape copy as described below. The database was compiled using ARC/INFO, a commercial Geographic Information System (Environmental Systems Research Institute, Redlands, California), with version 3.0 of the menu interface ALACARTE (Fitzgibbon and Wentworth, 1991, Fitzgibbon, 1991, Wentworth and Fitzgibbon, 1991). The ARC/INFO coverages are stored in uncompressed ARC export format (ARC/INFO version 7.x). All data files have been compressed, and may be uncompressed with gzip, which is available free of charge over the Internet via links from the USGS Public Domain Software page (http://edcwww.cr.usgs.gov/doc/edchome/ndcdb/public.html). ARC/INFO export files (files with the .e00 extension) can be converted into ARC/INFO coverages in ARC/INFO (see below) and can be read by some other Geographic Information Systems, such as MapInfo via ArcLink and ESRI's ArcView.

California

Three-dimensional model for multi-component reactive transport with variable density groundwater flow

PHWAT is a new model that couples a geochemical reaction model (PHREEQC-2) with a density-dependent groundwater flow and solute transport model (SEAWAT) using the split-operator approach. PHWAT was developed to simulate multi-component reactive transport in variable density groundwater flow. Fluid density in PHWAT depends not on only the concentration of a single species as in SEAWAT, but also the concentrations of other dissolved chemicals that can be subject to reactive processes. Simulation results of PHWAT and PHREEQC-2 were compared in their predictions of effluent concentration from a column experiment. Both models produced identical results, showing that PHWAT has correctly coupled the sub-packages. PHWAT was then applied to the simulation of a tank experiment in which seawater intrusion was accompanied by cation exchange. The density dependence of the intrusion and the snow-plough effect in the breakthrough curves were reflected in the model simulations, which were in good agreement with the measured breakthrough data. Comparison simulations that, in turn, excluded density effects and reactions allowed us to quantify the marked effect of ignoring these processes. Next, we explored numerical issues involved in the practical application of PHWAT using the example of a dense plume flowing into a tank containing fresh water. It was shown that PHWAT could model physically unstable flow and that numerical instabilities were suppressed. Physical instability developed in the model in accordance with the increase of the modified Rayleigh number for density-dependent flow, in agreement with previous research. ?? 2004 Elsevier Ltd. All rights reserved.

Environmental Modelling and Software

Global Earthquake Model (GEM) Risk Map

The Global Earthquake Risk Map (v2018.1) comprises four global maps. The main map presents the geographic distribution of average annual loss (USD) normalized by the average construction costs of the respective country (USD/m2 due to ground shaking in the residential, commercial and industrial building stock, considering contents, structural and non-structural components. The normalized metric allows a direct comparison of the risk between countries with widely different construction costs. It does not consider the effects of tsunamis, liquefaction, landslides, and fires following earthquakes. The loss estimates are from direct physical damage to buildings due to shaking, and thus damage to infrastructure or indirect losses due to business interruption are not included. The Global Earthquake Hazard Map depicts the geographic distribution of the Peak Ground Acceleration (PGA) with a 10% probability of being exceeded in 50 years, computed for reference rock conditions (shear wave velocity of 760-800 m/s). The Global Exposure Map depicts the geographic distribution of residential, commercial and industrial buildings. The Global Seismic Fatalities Map depicts an estimate of average annual human losses due to earthquake-induced structural collapse of buildings. The results for human losses do not consider indirect fatalities such as those from post­ earthquake epidemics. The average annual losses and number of buildings are presented on a hexagonal grid, with a spacing of 0.30 x 0.34 decimal degrees (approximately 1,000 km2 at the equator). The average annual losses were computed using the event-based calculator of the OpenQuake engine, an open-source software for seismic hazard and risk analysis developed by the GEM Foundation. The seismic hazard, exposure and vulnerability models employed in these calculations were provided by national institutions, or developed within the scope of regional programs or bilateral collaborations. These global maps and the underlying databases are based on best available and publicly accessible datasets and models. Due to possible limitations in the model, regions portrayed with low risk may experience potentially damaging earthquakes. The GEM Risk Map is intended to be a dynamic product, which will be updated when new datasets and models become available. Updated versions of the hazard, exposure, and average annual losses will be released on a regular basis. Additional metrics for each country can be explored at globalquakemodel.org/gem.

Report

A quantitative X-ray diffraction inventory of the tephra and volcanic glass inputs into the Holocene marine sediment archives off Iceland: A contribution to V.A.S.T.

This paper re-evaluates how well quantitative x-ray diffraction (qXRD) can be used as an exploratory method of the weight percentage (wt%) of volcaniclastic sediment, and to identify tephra events in marine cores. In the widely used RockJock v6 software programme, qXRD tephra and glass standards include the rhyodacite White River tephra (Alaska), a rhyolitic tephra (Hekla-4) and the basaltic Saksunarvatn tephra. Experiments of adding known wt% of tephra to felsic bedrock samples indicated that additions ≥10 wt% are accurately detected, but reliable estimates of lesser amounts are masked by amorphous material produced by milling. Volcaniclastic inputs range between 20 and 50 wt%. Primary tephra events are identified as peaks in residual qXRD glass wt% from fourth-order polynomial fits. In cores where tephras have been identified by shard counts in the > 150 µm fraction, there is a positive correlation (validation) with peaks in the wt% glass estimated by qXRD. Geochemistry of tephra shards confirms the presence of several Hekla-sourced tephras in cores B997-317PC1 and -319PC2 on the northern Iceland shelf. In core B997-338 (north-west Iceland), there are two rhyolitic tephras separated by ca. 100 cm with uncorrected radiocarbon dates on articulated shells of around 13 000 yr B.P. These tephras may be correlatives of the Borrobol and Penifiler tephras found in Scotland. The number of Holocene tephra events per 1000 yr was estimated from qXRD on 16 cores and showed a bimodal distribution with an increased number of events in both the late and early Holocene.

Polar Research