USGS Science⌕ Search

SEARCH · USGS Science

Results for “Algorithms”

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 811 records · Page 45Linked to original sources
USGS2021

Improving ChemCam LIBS long-distance elemental compositions using empirical abundance trends

The ChemCam instrument on the Curiosity rover provides chemical compositions of Martian rocks and soils using remote laser-induced breakdown spectroscopy (LIBS). The elemental calibration is stable as a function of distance for Ti, Fe, Mg, and Ca. The calibration shows small, systematically increasing abundance trends as a function of distance for Al, Na, K, and to some extent, Si. The distance effect is known to be due to a dependence with distance on the relative strengths of atomic transition lines. Emission lines representing transitions from relatively low energy levels remain intense at longer distances while emission lines representing transitions from higher energy levels decrease in intensity more rapidly as a function of distance. The multivariate algorithms used to determine elemental compositions rely on a large number of emission lines in many cases, so rather than trying to correct all emission lines, a study was made of the predicted compositions as a function of distance, in order to determine an empirical correction. Abundance trends can be well approximated by a linear trend with distance within the ranges of abundances and distances observed up to ~6 m. Data from 11 distinct geological members and data groups of the Murray formation in Gale crater, Mars, were used to form the model, selecting the members and data groups yielding the best statistics. The model was tested using data from several targets observed from two different distances, and using data from the Kimberley formation, the composition of which is significantly different from the Murray formation, showing that the model works on other compositions beyond those used to build the model. For long-distance observations up to ~6 m, corrections can be made back to an equivalent composition at the median distance of ChemCam observations (2.6 m). The model has been validated up to 6.2 m, although ChemCam is able to observe bedrock targets to >7 m, and iron meteorites to distances of >9 m.

Spectrochimica Acta Part B: Atomic Spectroscopy↗
USGS2020

The response of stream ecosystems in the Adirondack region of New York to historical and future changes in atmospheric deposition of sulfur and nitrogen

The present-day acid-base chemistry of surface waters can be directly linked to contemporary observations of acid deposition; however, pre-industrial conditions are key to predicting the potential future recovery of stream ecosystems under decreasing loads of atmospheric sulfur (S) and nitrogen (N) deposition. The integrated biogeochemical model PnET-BGC was applied to 25 forest watersheds that represent a range of acid sensitivity in the Adirondack region of New York, USA to simulate the response of streams to past and future changes in atmospheric S and N deposition, and calculate the target loads of acidity for protecting and restoring stream water quality and ecosystem health. Using measured data, the model was calibrated and applied to simulate soil and stream chemistry at all study sites. Model hindcasts indicate that historically stream water chemistry in the Adirondacks was variable, but inherently sensitive to acid deposition. The median model-simulated acid neutralizing capacity (ANC) of the streams was projected to be 55 μeq L −1 before the advent of anthropogenic acid deposition (~1850), decreasing to minimum values of 10 μeq L −1 around the year 2000. The median simulated ANC increased to 13 μeq L −1 by 2015 in response to decreases in acid deposition that have occurred over recent decades. Model projections suggest that simultaneous decreases in sulfate, nitrate and ammonium deposition are more effective in restoring stream ANC than individual decreases in sulfur or nitrogen deposition. However, the increases in stream ANC per unit equivalent decrease in S deposition is greater compared to decreases in N deposition. Using empirical algorithms, fish community density and biomass are projected to increase under several deposition-control scenarios that coincide with increases in stream ANC. Model projections suggest that even under the most aggressive deposition-reduction scenarios, stream chemistry and fisheries will not fully recover from historical acidification by 2200.

New York↗
USGS2022

Urban landcover differentially drives day and nighttime air temperature across a semi-arid city

Semi-arid urban environments are undergoing an increase in both average air temperatures and in the frequency and intensity of extreme heat events. Within cities, different composition and densities of urban landcovers (ULC) influence local air temperatures, either mitigating or increasing heat. Currently, understanding how combinations of ULC influence air temperature at the block to neighborhood scale is necessary for heat mitigation plans, and yet limited due to the complexities integrating high-resolution ULC with spatial and temporally high-resolution microclimate data. We quantify how ULC influences air temperature at 60 m resolution for day and nighttime climate normals and extreme heat conditions by integrating microclimate sensor data sensor and high-resolution (1 m 2 ) ULC for Denver, Colorado's urban core. We derive ULC drivers of air temperature using a structural equation model, then use a random forest algorithm to predict air temperatures for 30-year climate normals and an extreme heat condition. We find that, in conjunction with other ULC, urban tree canopy reduces daytime air temperatures (−0.026 °C per % cover), and the combination of impervious surfaces and buildings increases daytime air temperature (0.021 °C per % cover). Compared to daytime hours, nighttime irrigated turf temperature cooling effects are increased from being non-significant to −0.022 °C per % cover, while tree canopy effects are reduced from −0.026 °C during the day to −0.016 °C at night. Overall, ULC drives ~17% and 25% of local air temperature during the day and night, respectively. ULC influence on daytime air temperatures is altered in extreme heat events, both depending on the ULC type and time of day. Our findings inform urban planners seeking to identify potential hot and cool spots within a semi-arid city and mitigate high urban air temperatures through using ULC within larger urban climate mitigation strategies.

Colorado↗
USGS2025

Phytoplankton biomass dynamics in wet (2019) and dry (2023) years in Lake Pontchartrain estuary, Louisiana from Sentinel 2-MSI and PACE-OCI observations

This study provides a comprehensive assessment of phytoplankton biomass dynamics in Lake Pontchartrain, Louisiana, by combining monthly water quality data with multispectral and hyperspectral satellite observations using a machine learning algorithm. A machine learning model based on Variational Autoencoder (VAE), globally applicable, was used to estimate phytoplankton biomass via chlorophyll- a (Chl- a ) from Sentinel 2-MSI and NASA's new hyperspectral mission, PACE-OCI, enabling the first direct comparison between the two sensors. The model performed well in this complex estuarine system, with higher accuracy from PACE-OCI (MAE: 1.48, RMSE: 10.40, slope: 0.87) than Sentinel 2-MSI (MAE: 1.57, RMSE: 11.08, slope: 0.83). This approach enabled continuous high-resolution monitoring of phytoplankton biomass across space and time. Comparative analysis of 2019, a wet year with Bonnet Carré Spillway (BCS) openings, and 2023, a dry year with extremely low riverine inputs, revealed distinct biomass dynamics. In 2019, BCS discharge initially suppressed Chl- a within turbid waters (<5 mg Chl- a m −3 ) but later acted as a nutrient and hydrodynamic driver, transporting nutrients toward the lake outlet and Mississippi coast, promoting high biomass (25–45 mg Chl- a m −3 ) near the entrance. In contrast, dry conditions in 2023 led to more frequent-than-expected high biomass (>35 mg Chl- a m −3 ), persisting in the lake center. Similar spatial patterns were observed again in 2024, revealed for the first time by PACE-OCI. This study demonstrates the value of satellite-derived observations for capturing transient phytoplankton biomass events and highlights the potential of PACE-OCI's hyperspectral capabilities to better distinguish phytoplankton communities and improve understanding of their responses to freshwater inflows and associated processes driving pulses into estuaries.

Louisiana↗
USGS2018

AutoCNet: A Python library for sparse multi-image correspondence identification for planetary data

In this work we describe the AutoCNet library, written in Python, to support the application of Computer Vision techniques for n-image correspondence identication in remotely sensed planetary images and subsequent bundle adjustment. The library is designed to support exploratory data analysis, algorithm and processing pipeline development, and application at scale in High Performance Computing (HPC) environments for processing large data sets and generating foundational data products. We also present a brief case study illustrating high level usage for the Apollo 15 Metric camera.

SoftwareX↗
USGS2024

Spatio-temporal ecological models via physics-informed neural networks for studying chronic wasting disease

To mitigate the negative effects of emerging wildlife diseases in biodiversity and public health it is critical to accurately forecast pathogen dissemination while incorporating relevant spatio-temporal covariates. Forecasting spatio-temporal processes can often be improved by incorporating scientific knowledge about the dynamics of the process using physical models. Ecological diffusion equations are often used to model epidemiological processes of wildlife diseases where environmental factors play a role in disease spread. Physics-informed neural networks (PINN) are deep learning algorithms that constrain neural network predictions based on physical laws and therefore are powerful forecasting models useful even in cases of limited and imperfect training data. In this paper, we develop a novel ecological modeling tool using PINNs, which fits a feedforward neural network and simultaneously performs parameter identification in a partial differential equation (PDE) with varying coefficients. We demonstrate the applicability of our model by comparing it with the commonly used Bayesian stochastic partial differential equation method and traditional machine learning approaches, showing that our proposed model exhibits superior prediction and forecasting performance when modeling chronic wasting disease in deer in Wisconsin. Furthermore, our model provides the opportunity to obtain scientific insights into spatiotemporal covariates affecting spread and growth of diseases. This work contributes to future machine learning and statistical methodology development by studying spatio-temporal processes enhanced by prior physical knowledge.

Spatial Statistics↗
USGS2003

Data-based comparisons of moments estimators using historical and paleoflood data

This paper presents the first systematic comparison, using historical and paleoflood data, of moments-based flood frequency methods. Peak flow estimates were compiled from streamflow-gaging stations with historical and/or paleoflood data at 36 sites located in the United States, Argentina, United Kingdom and China, covering a diverse range of hydrologic conditions. The Expected Moments Algorithm (EMA) and the Bulletin 17B historical weighting procedure (B17H) were compared in terms of goodness of fit using 25 of the data sets. Results from this comparison indicate that EMA is a viable alternative to current B17H procedures from an operational perspective, and performed equal to or better than B17H for the data analyzed. We demonstrate satisfactory EMA performance for the remaining 11 sites with multiple thresholds and binomial censoring, which B17H cannot accommodate. It is shown that the EMA estimator readily incorporates these types of information and the LP-III distribution provided an adequate fit to the data in most cases. The results shown here are consistent with Monte Carlo simulation studies, and demonstrate that EMA is preferred overall to B17H. The Bulletin 17B document could be revised to include an option for EMA as an alternative to the existing historical weighting approach. These results are of practical relevance to hydrologists and water resources managers for applications in floodplain management, design of hydraulic structures, and risk analysis for dams. ?? 2003 Elsevier Science B.V. All rights reserved.

Journal of Hydrology↗
USGS2002

Processing, mosaicking and management of the Monterey Bay digital sidescan-sonar images

Sidescan-sonar imaging systems with digital capabilities have now been available for approximately 20 years. In this paper we present several of the various digital image processing techniques developed by the U.S. Geological Survey (USGS) and used to apply intensity/radiometric and geometric corrections, as well as enhance and digitally mosaic, sidescan-sonar images of the Monterey Bay region. New software run by a WWW server was designed and implemented to allow very large image data sets, such as the digital mosaic, to be easily viewed interactively, including the ability to roam throughout the digital mosaic at the web site in either compressed or full 1-m resolution. The processing is separated into the two different stages: preprocessing and information extraction. In the preprocessing stage, sensor-specific algorithms are applied to correct for both geometric and intensity/radiometric distortions introduced by the sensor. This is followed by digital mosaicking of the track-line strips into quadrangle format which can be used as input to either visual or digital image analysis and interpretation. An automatic seam removal procedure was used in combination with an interactive digital feathering/stenciling procedure to help minimize tone or seam matching problems between image strips from adjacent track-lines. The sidescan-sonar image processing package is part of the USGS Mini Image Processing System (MIPS) and has been designed to process data collected by any 'generic' digital sidescan-sonar imaging system. The USGS MIPS software, developed over the last 20 years as a public domain package, is available on the WWW at: http://terraweb.wr.usgs.gov/trs/software.html.

Marine Geology↗
USGS2001

Attractor reconstruction for non-linear systems: a methodological note

Attractor reconstruction is an important step in the process of making predictions for non-linear time-series and in the computation of certain invariant quantities used to characterize the dynamics of such series. The utility of computed predictions and invariant quantities is dependent on the accuracy of attractor reconstruction, which in turn is determined by the methods used in the reconstruction process. This paper suggests methods by which the delay and embedding dimension may be selected for a typical delay coordinate reconstruction. A comparison is drawn between the use of the autocorrelation function and mutual information in quantifying the delay. In addition, a false nearest neighbor (FNN) approach is used in minimizing the number of delay vectors needed. Results highlight the need for an accurate reconstruction in the computation of the Lyapunov spectrum and in prediction algorithms.

Mathematical Biosciences↗
USGS2003

Moment-tensor solutions estimated using optimal filter theory: Global seismicity, 2001

This paper is the 12th in a series published yearly containing moment-tensor solutions computed at the US Geological Survey using an algorithm based on the theory of optimal filter design ( Sipkin, 1982 and Sipkin, 1986b ). An inversion has been attempted for all earthquakes with a magnitude, m b or M S , of 5.5 or greater. Previous listings include solutions for earthquakes that occurred from 1981 to 2000 ( Sipkin, 1986b ; Sipkin and Needham, 1989 , Sipkin and Needham, 1991 , Sipkin and Needham, 1992 , Sipkin and Needham, 1993 , Sipkin and Needham, 1994a and Sipkin and Needham, 1994b ; Sipkin and Zirbes, 1996 and Sipkin and Zirbes, 1997 ; Sipkin et al., 1998 , Sipkin et al., 1999 , Sipkin et al., 2000a , Sipkin et al., 2000b and Sipkin et al., 2002 ). The entire USGS moment-tensor catalog can be obtained via anonymous FTP at ftp://ghtftp.cr.usgs.gov . After logging on, change directory to “momten”. This directory contains two compressed ASCII files that contain the finalized solutions, “mt.lis.Z” and “fmech.lis.Z”. “mt.lis.Z” contains the elements of the moment tensors along with detailed event information; “fmech.lis.Z” contains the decompositions into the principal axes and best double-couples. The fast moment-tensor solutions for more recent events that have not yet been finalized and added to the catalog, are gathered by month in the files “jan01.lis.Z”, etc. “fmech.doc.Z” describes the various fields.

Physics of the Earth and Planetary Interiors↗
USGS2001

Investigating a physical basis for spectroscopic estimates of leaf nitrogen concentration

The reflectance spectra of dried and ground plant foliage are examined for changes directly due to increasing nitrogen concentration. A broadening of the 2.1-??m absorption feature is observed as nitrogen concentration increases. The broadening is shown to arise from two absorptions at 2.054 ??m and 2.172 ??m. The wavelength positions of these absorptions coincide with the absorption characteristics of the nitrogen-containing amide bonds in proteins. The observed presence of these absorption features in the reflectance spectra of dried foliage is suggested to form a physical basis for high correlations established by stepwise multiple linear regression techniques between the reflectance of dry plant samples and their nitrogen concentration. The consistent change in the 2.1-??m absorption feature as nitrogen increases and the offset position of protein absorptions compared to those of other plant components together indicate that a generally applicable algorithm may be developed for spectroscopic estimates of nitrogen concentration from the reflectance spectra of dried plant foliage samples. ?? 2001 Published by Elsevier Science Ireland Ltd.

Remote Sensing of Environment↗
USGS2003

Mapping vegetation in Yellowstone National Park using spectral feature analysis of AVIRIS data

Knowledge of the distribution of vegetation on the landscape can be used to investigate ecosystem functioning. The sizes and movements of animal populations can be linked to resources provided by different plant species. This paper demonstrates the application of imaging spectroscopy to the study of vegetation in Yellowstone National Park (Yellowstone) using spectral feature analysis of data from the Airborne Visible/Infrared Imaging Spectrometer (AVIRIS). AVIRIS data, acquired on August 7, 1996, were calibrated to surface reflectance using a radiative transfer model and field reflectance measurements of a ground calibration site. A spectral library of canopy reflectance signatures was created by averaging pixels of the calibrated AVIRIS data over areas of known forest and nonforest vegetation cover types in Yellowstone. Using continuum removal and least squares fitting algorithms in the US Geological Survey's Tetracorder expert system, the distributions of these vegetation types were determined by comparing the absorption features of vegetation in the spectral library with the spectra from the AVIRIS data. The 0.68 μm chlorophyll absorption feature and leaf water absorption features, centered near 0.98 and 1.20 μm, were analyzed. Nonforest cover types of sagebrush, grasslands, willows, sedges, and other wetland vegetation were mapped in the Lamar Valley of Yellowstone. Conifer cover types of lodgepole pine, whitebark pine, Douglas fir, and mixed Engelmann spruce/subalpine fir forests were spectrally discriminated and their distributions mapped in the AVIRIS images. In the Mount Washburn area of Yellowstone, a comparison of the AVIRIS map of forest cover types to a map derived from air photos resulted in an overall agreement of 74.1% (kappa statistic=0.62).

Wyoming↗
USGS2000

A diurnal animation of thermal images from a day-night pair

Interpretation of thermal images is often complicated because the physical property information is contained in both the spatial and temporal variations of the data and thermal models are necessary to extract and display this information. A linearized radiative transfer solution to the surface flux has been used to derive a function that is invariant with respect to thermal inertia. This relationship makes it possible to predict the temperature variation at any time in the diurnal cycle using only two distinct measurements (e.g., noon and midnight). An animation can then be constructed from a pair of day-night images to view both the spatial and temporal temperature changes throughout the diurnal cycle. A more complete solution for the invariant function, using the method of Laplace transforms and based on the linearized solution, was introduced. These results indicate that the linear model does not provide a sufficiently accurate estimate. Using standard conditions (latitude 30??, solar declination 0??, acquisition times at noon and midnight), this new relationship was used to predict temperature throughout the diurnal cycle to an rms error of 0.2??C, which is close to the system noise of most thermal scanners. The method was further extended to include the primary effects of topographic slope with similar accuracy. The temperature was computed at 48 equally spaced times in the diurnal cycle with this algorithm using a co-registered day and night TIMS (Thermal Infrared Multispectral Scanner) data pair (330 pixels, 450 lilies) acquired of the Carlin, Nevada, area and a co-registered DEM (Digital Elevation Model). (Any reader can view the results by downloading the animation file from an identified tip site). The results illustrate the power of animation to display subtle temporal and spatial temperature changes, which can provide clues to structural controls and material property differences. This 'visual change' approach could significantly increase the use of thermal data for environmental, hazard, and resource studies. Published by Elsevier Science Inc., 2000.A linearized radiative transfer solution of determining the surface flux is proposed to predict the temperature variation at any time in the diurnal cycle using only two distinct measurements. An animation is constructed from a pair of day-night images to view the spatial and temporal temperature changes throughout the diurnal cycle. The results illustrate the effectiveness of animation to display subtle temporal and spatial temperature changes, which can provide clues to structural controls and material property differences.

Remote Sensing of Environment↗
USGS2003

The IASPEI Seismological Software Library

Since computers became widely available in the early 1960s, seismologists have been using them for data acquisition, processing, and analysis, as well as theoretical computation and modeling. For example, the book by Doornbos (1988) contains a collection of seismological algorithms with the corresponding computer programs available on tape or disk from the World Data Center A for Solid Earth Geophysics. The introduction of personal computers in the early 1980s further revolutionized the use of computers for scientific research. Instead of expensive mainframe computers that required a large staff to operate, inexpensive personal computers allowed creative applications to be implemented by individuals with a shoestring budget.

International Geophysics↗
USGS2001

1r2dinv: A finite-difference model for inverse analysis of two dimensional linear or radial groundwater flow

We have developed a program for inverse analysis of two-dimensional linear or radial groundwater flow problems. The program, 1r2dinv, uses standard finite difference techniques to solve the groundwater flow equation for a horizontal or vertical plane with heterogeneous properties. In radial mode, the program simulates flow to a well in a vertical plane, transforming the radial flow equation into an equivalent problem in Cartesian coordinates. The physical parameters in the model are horizontal or x-direction hydraulic conductivity, anisotropy ratio (vertical to horizontal conductivity in a vertical model, y-direction to x-direction in a horizontal model), and specific storage. The program allows the user to specify arbitrary and independent zonations of these three parameters and also to specify which zonal parameter values are known and which are unknown. The Levenberg-Marquardt algorithm is used to estimate parameters from observed head values. Particularly powerful features of the program are the ability to perform simultaneous analysis of heads from different tests and the inclusion of the wellbore in the radial mode. These capabilities allow the program to be used for analysis of suites of well tests, such as multilevel slug tests or pumping tests in a tomographic format. The combination of information from tests stressing different vertical levels in an aquifer provides the means for accurately estimating vertical variations in conductivity, a factor profoundly influencing contaminant transport in the subsurface. ?? 2001 Elsevier Science Ltd. All rights reserved.

Computers & Geosciences↗
USGS1997

Time-series animation techniques for visualizing urban growth

Time-series animation is a visually intuitive way to display urban growth. Animations of landuse change for the Baltimore-Washington region were generated by showing a series of images one after the other in sequential order. Before creating an animation, various issues which will affect the appearance of the animation should be considered, including the number of original data frames to use, the optimal animation display speed, the number of intermediate frames to create between the known frames, and the output media on which the animations will be displayed. To create new frames between the known years of data, the change in each theme (i.e. urban development, water bodies, transportation routes) must be characterized and an algorithm developed to create the in-between frames. Example time-series animations were created using a temporal GIS database of the Baltimore-Washington area. Creating the animations involved generating raster images of the urban development, water bodies, and principal transportation routes; overlaying the raster images on a background image; and importing the frames to a movie file. Three-dimensional perspective animations were created by draping each image over digital elevation data prior to importing the frames to a movie file.

Computers & Geosciences↗
USGS1999

UCODE, a computer code for universal inverse modeling

This article presents the US Geological Survey computer program UCODE, which was developed in collaboration with the US Army Corps of Engineers Waterways Experiment Station and the International Ground Water Modeling Center of the Colorado School of Mines. UCODE performs inverse modeling, posed as a parameter-estimation problem, using nonlinear regression. Any application model or set of models can be used; the only requirement is that they have numerical (ASCII or text only) input and output files and that the numbers in these files have sufficient significant digits. Application models can include preprocessors and postprocessors as well as models related to the processes of interest (physical, chemical and so on), making UCODE extremely powerful for model calibration. Estimated parameters can be defined flexibly with user-specified functions. Observations to be matched in the regression can be any quantity for which a simulated equivalent value can be produced, thus simulated equivalent values are calculated using values that appear in the application model output files and can be manipulated with additive and multiplicative functions, if necessary. Prior, or direct, information on estimated parameters also can be included in the regression. The nonlinear regression problem is solved by minimizing a weighted least-squares objective function with respect to the parameter values using a modified Gauss-Newton method. Sensitivities needed for the method are calculated approximately by forward or central differences and problems and solutions related to this approximation are discussed. Statistics are calculated and printed for use in (1) diagnosing inadequate data or identifying parameters that probably cannot be estimated with the available data, (2) evaluating estimated parameter values, (3) evaluating the model representation of the actual processes and (4) quantifying the uncertainty of model simulated values. UCODE is intended for use on any computer operating system: it consists of algorithms programmed in perl, a freeware language designed for text manipulation and Fortran90, which efficiently performs numerical calculations.

Computers & Geosciences↗
USGS2003

Experimental and textural constraints on mafic enclave formation in volcanic rocks

We have used experiments and textural analysis to investigate the process of enclave formation during magma mixing at Southwest Trident volcano, Alaska. Andesite enclaves are present throughout the four dacite lava flows produced by the eruption, and resemble mafic enclaves commonly found in other volcanic rocks. Our experiments replicate the pressure–temperature–time path taken by enclave-forming andesite magma as it is engulfed in dacite during magma mixing. Pressure and temperature information for the andesite and dacite are from [Coombs et al., Contrib. Mineral. Petrol. 140 (2000) 99–118]. The andesite was annealed at 1000°C, and then cooled to 890°C at rates of 110°C h −1 , 10°C h −1 , and 2°C h −1 . Once cooled to 890°C, andesite was held at this lower temperature from times ranging from 1 to 40 h. The andesite that was cooled at the slower rates of 2°C h −1 and 10°C h −1 most resembles enclave groundmass texturally and compositionally. Based on simple thermal calculations, these rates are more consistent with cooling of the andesite groundmass below an andesite–dacite interface than with cooling of enclave-sized spheres. If enclaves do crystallize as spheres, post-crystallization disaggregation must occur. Calculations using the MELTS algorithm [Ghiorso and Sack, Contrib. Mineral. Petrol. 119 (1995) 197–212] show that for incoming andesite to become less dense than the dacite ∼34 volume % of its groundmass must crystallize to undergo ∼18 volume % vesiculation; these values are similar to those determined for Southwest Trident enclaves. Thus such crystallization may lead to ‘flotation’ of enclaves and be a viable mechanism for enclave formation and dispersal. The residual melt in the cooling experiments did not evolve to rhyolitic compositions such as seen in natural enclaves due to a lack of a decompression step in the experiments. Decompression experiments on Southwest Trident dacite suggest an average ascent rate for the eruption of ∼2–3 MPa h −1 . An andesite experiment that was cooled and then decompressed at this rate contains melt that matches that of the natural enclaves. It is apparent that decompression (ascent)-induced crystallization occurs in enclaves, but not in the form of microlites as happens in the dacite host, due either to insufficient residence time at chamber temperatures or to the pre-existing microphenocrysts which act as sites for new growth.

Alaska↗