Software

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...

The code in this repository generates and analyzes directed acyclic graphs (DAGs) representing hypotheses about variables affecting the occurrence and surface availability of Plethodon salamanders in northern California; formats data from field surveys and GIS for analysis; and analyzes multi-scale occupancy models, including model selection, to evaluate the model and variables with the...

Welcome to the README for the Bull Trout Population Dynamics Simulation Tool. This tool can be used directly in an R console (including IDEs such as RStudio) or it can be used via the web and an R Shiny interface (https://rconnect.usgs.gov/BullTroutPopulationDynamicsSimulationTool/).

shakingstats is a Python package for calculating earthquake statistics from continuous seismic amplitudes. It calculates the time-varying amplitude completeness (A_C) and amplitude-based b‑value (b_A​) by fitting the joint distribution of seismic noise and earthquake shaking. It also provides an implementation of the b-positive (β+) method for comparison to catalog-based statistics.

This release adds functionality for building debris-flow likelihood lookup tables, and uploading the results to the National Weather Service (NWS) FTP endpoint. The release: * Adds the `authenticate nws` command, which collects authentication credentials for uploading to the NWS FTP endpoint * Updates the `run` command to compute debris-flow likelihood lookup tables, and save the results...

Effective future fire planning and decision-making in Alaska requires accessible, interpretable, and locally relevant climate data to anticipate fire danger. To meet this need, we utilized an existing dataset of projected fire danger metrics from the Canadian Forest Fire Danger Rating System (Young et al. 2025). This dataset contains daily projections from 1980 through 2099 for key...

R code and Stan models to simulate sample data set, perform parameter estimation, project future scenarios, and produce plots for accompanying publication examining the integrated dynamics of waterfowl and hunter populations through harvest.

This software was developed using the statistical program R (R Core Team 2024) to compile multiple wildlife tracking datasets that contain locations of animals recorded through time. The software is organized through a series of scripts that run hierarchically through script sourcing. To demonstrate how the software works, we simulated four wildlife tracking datasets, each representing...

This script is for analyzing raw eDNA qPCR data. It is designed to efficiently concatenate qPCR data from multiple plates exported as *.csv files. It will assist with completing several QA/QC checks and output documentation of these checks as well as helpful files for further review.

The Operational Aftershock Forecasting server runs continuously in the cloud, monitoring the USGS ComCat earthquake catalog. It automatically generates aftershock forecasts and publishes them on the earthquake event pages. The analyst utilities are GUI programs that are used by analysts to evaluate the operation of the OAF servers, and make adjustments to the statistical model used by...

"Max spots" or "maxspots" refers to the locations of the local maxima of the horizontal gradient of potential fields such as gravity or magnetic anomalies. Because maxspots mark locations of abrupt change in the anomalies, they can be used to identify the edges of subsurface features with significant contrasts in density and/or magnetization, such as faults, intrusions, or lithologic...

This software contains 3 separate R Markdown files and 1 R script file that comprise an analysis to 1. collate freshwater fish survey data from within the northeastern United States, 2. develop separate models to predict the probability of occurrence and the proportional abundance of 53 freshwater fish species, and 3. combine model results by fish species traits and habitat...