A finite difference code for the calculation of Fe-Mg interdiffusion timescales in orthopyroxene with initial condition optimization (version 1.0)
September 10, 2025
These Python scripts calculate Fe-Mg diffusion timescales in orthopyroxene crystals for use in studies of pre-eruption processes. The model is based on a centered space, forward time implicit finite difference solution to the diffusion equation with with a compositionally depended diffusivity. The model includes an initial condition optimization step and a Monte Carlo approach to determine the confidence interval of the calculated timescale.
Citation Information
| Publication Year | 2025 |
|---|---|
| Title | A finite difference code for the calculation of Fe-Mg interdiffusion timescales in orthopyroxene with initial condition optimization (version 1.0) |
| DOI | 10.5066/P1WWP3QM |
| Authors | Nathan L Andersen |
| Product Type | Software Release |
| Record Source | USGS Asset Identifier Service (AIS) |
| USGS Organization | USGS Volcano Science Center |
| Rights | This work is marked with CC0 1.0 Universal |
Related
The transition from melt accumulation to eruption initiation recorded by orthopyroxene Fe-Mg diffusion timescales in late Holocene rhyolites, South Sister volcano, Oregon Cascade Range The transition from melt accumulation to eruption initiation recorded by orthopyroxene Fe-Mg diffusion timescales in late Holocene rhyolites, South Sister volcano, Oregon Cascade Range
South Sister volcano, Oregon Cascade Range, USA, has repeatedly erupted rhyolite since ca. 40 ka. The youngest such eruptions are the ca. 2 ka Rock Mesa and Devils Chain rhyolites, erupted several hundred years apart from two multi-vent complexes separated by 3–6 km. Fe-Mg interdiffusion models of orthopyroxene rims from both rhyolites produce timescales up to several-thousand years, but
Authors
Nathan Lee Andersen, Annika E. Dechert, Dawn Catherine Sweeney Ruth, May (Mai) Sas, Julie Chouinard, Josef Dufek
Related
The transition from melt accumulation to eruption initiation recorded by orthopyroxene Fe-Mg diffusion timescales in late Holocene rhyolites, South Sister volcano, Oregon Cascade Range The transition from melt accumulation to eruption initiation recorded by orthopyroxene Fe-Mg diffusion timescales in late Holocene rhyolites, South Sister volcano, Oregon Cascade Range
South Sister volcano, Oregon Cascade Range, USA, has repeatedly erupted rhyolite since ca. 40 ka. The youngest such eruptions are the ca. 2 ka Rock Mesa and Devils Chain rhyolites, erupted several hundred years apart from two multi-vent complexes separated by 3–6 km. Fe-Mg interdiffusion models of orthopyroxene rims from both rhyolites produce timescales up to several-thousand years, but
Authors
Nathan Lee Andersen, Annika E. Dechert, Dawn Catherine Sweeney Ruth, May (Mai) Sas, Julie Chouinard, Josef Dufek