Search USGS⌕ Search

SEARCH · Search USGS

Results for “Journal of Applied Mechanics”

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 109 records · Page 6Linked to original sources
USGS2002

Waveform inversion of oscillatory signatures in long-period events beneath volcanoes

The source mechanism of long-period (LP) events is examined using synthetic waveforms generated by the acoustic resonance of a fluid-filled crack. We perform a series of numerical tests in which the oscillatory signatures of synthetic LP waveforms are used to determine the source time functions of the six moment tensor components from waveform inversions assuming a point source. The results indicate that the moment tensor representation is valid for the odd modes of crack resonance with wavelengths 2 L / n , 2 W / n , n = 3, 5, 7, …, where L and W are the crack length and width, respectively. For the even modes with wavelengths 2 L / n , 2 W / n , n = 2, 4, 6, …, a generalized source representation using higher-order tensors is required, although the efficiency of seismic waves radiated by the even modes is expected to be small. We apply the moment tensor inversion to the oscillatory signatures of an LP event observed at Kusatsu-Shirane Volcano, central Japan. Our results point to the resonance of a subhorizontal crack located a few hundred meters beneath the summit crater lakes. The present approach may be useful to quantify the source location, geometry, and force system of LP events, and opens the way for moment tensor inversions of tremor.

Journal of Geophysical Research B: Solid Earth↗
USGS2016

Three-dimensional surface deformation derived from airborne interferometric UAVSAR: Application to the Slumgullion Landslide

In order to provide surface geodetic measurements with “landslide-wide” spatial coverage, we develop and validate a method for the characterization of 3-D surface deformation using the unique capabilities of the Uninhabited Aerial Vehicle Synthetic Aperture Radar (UAVSAR) airborne repeat-pass radar interferometry system. We apply our method at the well-studied Slumgullion Landslide, which is 3.9 km long and moves persistently at rates up to ∼2 cm/day. A comparison with concurrent GPS measurements validates this method and shows that it provides reliable and accurate 3-D surface deformation measurements. The UAVSAR-derived vector velocity field measurements accurately capture the sharp boundaries defining previously identified kinematic units and geomorphic domains within the landslide. We acquired data across the landslide during spring and summer and identify that the landslide moves more slowly during summer except at its head, presumably in response to spatiotemporal variations in snowmelt infiltration. In order to constrain the mechanics controlling landslide motion from surface velocity measurements, we present an inversion framework for the extraction of slide thickness and basal geometry from dense 3-D surface velocity fields. We find that the average depth of the Slumgullion Landslide is 7.5 m, several meters less than previous depth estimates. We show that by considering a viscoplastic rheology, we can derive tighter theoretical bounds on the rheological parameter relating mean horizontal flow rate to surface velocity. Using inclinometer data for slow-moving, clay-rich landslides across the globe, we find a consistent value for the rheological parameter of 0.85 ± 0.08.

Colorado↗
USGS2019

Effects of infiltration characteristics on the spatial-temporal evolution of stability of an interstate highway embankment

Infiltration-induced landslides are among the most common natural disasters threatening modern civilization, but conventional methods for studying the triggering mechanisms and predicting the occurrence of these slides are limited by incomplete consideration of underlying physical processes and the lack of precision inherent in limit-equilibrium analyses. To address this problem the spatial-temporal evolution of failure is investigated in a seasonally unstable section of interstate highway embankment, known as the Straight Creek landslide, Colorado. The study includes multi-year site investigation, monitoring, and numerical simulation using a rigorous hydromechanical framework along with a field of local factor of safety method. The sensitivity of episodic landslide reactivation to infiltration characteristics is evaluated. Results indicate that annual cumulative snowmelt infiltration, which typically accounts for approximately 75% of total annual cumulative infiltration and occurs over a short period in the spring, has the most substantial impact on slide activation. The rate of snowmelt infiltration varies independently of annual cumulative snowmelt infiltration and cumulative infiltration in the previous year, but still affects antecedent soil moisture conditions at the onset of snowmelt infiltration and therefore also the level of slide activation. These findings are used to establish specific thresholds for exacerbated slide movement using annual snowpack accumulation, forecasted snowmelt rate, and the previous year’s snowmelt, an approach which may be applied for predicting movement at this and other recurring or potential slide sites.

Colorado↗
USGS2014

Mechanisms driving recruitment variability in fish: comparisons between the Laurentian Great Lakes and marine systems

In his seminal work, Hjort (in Fluctuations in the great fisheries of Northern Europe. Conseil Parmanent International Pour L'Exploration De La Mar. Rapports et Proces-Verbaux, 20: 1–228, 1914) observed that fish population levels fluctuated widely, year-class strength was set early in life, and egg production by adults could not alone explain variability in year-class strength. These observations laid the foundation for hypotheses on mechanisms driving recruitment variability in marine systems. More recently, researchers have sought to explain year-class strength of important fish in the Laurentian Great Lakes and some of the hypotheses developed for marine fisheries have been transferred to Great Lakes fish. We conducted a literature review to determine the applicability of marine recruitment hypotheses to Great Lakes fish. We found that temperature, interspecific interactions, and spawner effects (abundance, age, and condition of adults) were the most important factors in explaining recruitment variability in Great Lakes fish, whereas relatively fewer studies identified bottom-up trophodynamic factors or hydrodynamic factors as important. Next, we compared recruitment between Great Lakes and Baltic Sea fish populations and found no statistical difference in factors driving recruitment between the two systems, indicating that recruitment hypotheses may often be transferable between Great Lakes and marine systems. Many recruitment hypotheses developed for marine fish have yet to be applied to Great Lakes fish. We suggest that future research on recruitment in the Great Lakes should focus on forecasting the effects of climate change and invasive species. Further, because the Great Lakes are smaller and more enclosed than marine systems, and have abundant fishery-independent data, they are excellent candidates for future hypothesis testing on recruitment in fish.

Laurentian Great Lakes↗
USGS2003

Characterization of the time-dependent strain field at seismogenic depths using first-motion focal mechanisms: Observations of large-scale decadal variations in stress along the San Andrea fault system

We present a method for summing moment tensors derived from first-motion focal mechanisms to study temporal dependence in features of the subsurface regional strain field. Time-dependent processes are inferred by comparing mechanisms summed over differing time periods. We apply this methodology to seismogenic zones in central and southern California using focal mechanisms produced by the Northern and Southern California Seismograph Networks for events during 1980-1999. We find a consistent pattern in both the style of deformation (strike-slip versus compressional) and seismicity rate across the entire region. If these temporal variations are causally related, it suggests a temporal change in the regional-scale stress field. One change consistent with the observations is a rotation in the regional maximum horizontal compressive stress direction, followed by a reversal to the original direction. Depending upon the dominant style of deformation locally, this change in orientation of the regional stress will tend to either enhance or hinder deformation. The mode of enhanced deformation can range from increased microseismicity and creep to major earthquakes. We hypothesize that these temporal changes in the regional stress field are the result of subtle changes in apparent relative plate motion between the Pacific and North American plates, perhaps due to long-range postseismic stress diffusion. Others have hypothesized that small changes in plate motion over thousands of years, and/or over decades, are responsible for changes in the style of deformation in southern California. We propose that such changes, over the course of just a few years, also affect the style of deformation.

Journal of Geophysical Research B: Solid Earth↗
USGS2002

Topographic stress perturbations in southern Davis Mountains, west Texas 1. Polarity reversal of principal stresses

We have applied a previously developed analytical stress model to interpret subsurface stress conditions inferred from acoustic televiewer logs obtained in two municipal water wells located in a valley in the southern Davis Mountains near Alpine, Texas. The appearance of stress-induced breakouts with orientations that shift by 90° at two different depths in one of the wells is explained by results from exact solutions for the effects of valleys on gravity and tectonically induced subsurface stresses. The theoretical results demonstrate that above a reference depth termed the hinge point, a location that is dependent on Poisson's ratio, valley shape, and magnitude of the maximum horizontal tectonic stress normal to the long axis of the valley, horizontal stresses parallel to the valley axis are greater than those normal to it. At depths below this hinge point the situation reverses and horizontal stresses normal to the valley axis are greater than those parallel to it. Application of the theoretical model at Alpine is accommodated by the fact that nearby earthquake focal mechanisms establish an extensional stress regime with the regional maximum horizontal principal stress aligned perpendicular to the valley axis. We conclude that the localized stress field associated with a valley setting can be highly variable and that breakouts need to be examined in this context when estimating the orientations and magnitudes of regional principal stresses.

Texas↗
USGS2015

A new method for discovering behavior patterns among animal movements

Advanced satellite tracking technologies enable biologists to track animal movements at fine spatial and temporal scales. The resultant data present opportunities and challenges for understanding animal behavioral mechanisms. In this paper, we develop a new method to elucidate animal movement patterns from tracking data. Here, we propose the notion of continuous behavior patterns as a concise representation of popular migration routes and underlying sequential behaviors during migration. Each stage in the pattern is characterized in terms of space (i.e., the places traversed during movements) and time (i.e. the time spent in those places); that is, the behavioral state corresponding to a stage is inferred according to the spatiotemporal and sequential context. Hence, the pattern may be interpreted predictably. We develop a candidate generation and refinement framework to derive all continuous behavior patterns from raw trajectories. In the framework, we first define the representative spots to denote the underlying potential behavioral states that are extracted from individual trajectories according to the similarity of relaxed continuous locations in certain distinct time intervals. We determine the common behaviors of multiple individuals according to the spatiotemporal proximity of representative spots and apply a projection-based extension approach to generate candidate sequential behavior sequences as candidate patterns. Finally, the candidate generation procedure is combined with a refinement procedure to derive continuous behavior patterns. We apply an ordered processing strategy to accelerate candidate refinement. The proposed patterns and discovery framework are evaluated through conceptual experiments on both real GPS-tracking and large synthetic datasets.

International Journal of Geographical Information ↗
USGS2000

Stress resistance strategy in an arid land shrub: interactions between developmental instability and fractal dimention

This paper investigates allocation of energy to mechanisms that generate and preserve architectural forms (i.e. developmental stability, complexity of branching patterns) and productivity (growth and reproduction) in response to environmental disturbances (i.e. grazing and resource availability). The statistical error in translational symmetry was used to detect random intra-individual variability during development. This can be thought of as a measure of developmental instability caused by stress. Additionally, we use changes in fractal complexity and shoot distribution of branch structures as an alternate indicator of stress. These methods were applied to Anthyllis cytisoides L., a semi-arid environment shrub, to ascertain the effect of grazing and slope exposure on developmental traits in a 2×2 factorial design. The results show that A. cytisoides maintains developmental stability at the expense of productivity. Anthyllis cytisoides was developmentally more stable when grazed and when on south-facing, as opposed to north-facing slopes. On the contrary, shoot length, leaf area, fractal dimension and reproductive-to-vegetative allocation ratio were larger in north- than in south-facing slopes. As a consequence, under extreme xeric conditions, shrub mortality increased in north-facing slopes, especially when not grazed. The removal of transpiring area and the reduction of plant competition favoured developmental stability and survival in grazed plants. Differences between grazed and ungrazed plants were most evident in more mesic (north-facing) areas.

Journal of Fish Diseases↗
USGS2005

Applying petrophysical models to radar travel time and electrical resistivity tomograms: Resolution-dependent limitations

[1] Geophysical imaging has traditionally provided qualitative information about geologic structure; however, there is increasing interest in using petrophysical models to convert tomograms to quantitative estimates of hydrogeologic, mechanical, or geochemical parameters of interest (e.g., permeability, porosity, water content, and salinity). Unfortunately, petrophysical estimation based on tomograms is complicated by limited and variable image resolution, which depends on (1) measurement physics (e.g., electrical conduction or electromagnetic wave propagation), (2) parameterization and regularization, (3) measurement error, and (4) spatial variability. We present a framework to predict how core‐scale relations between geophysical properties and hydrologic parameters are altered by the inversion, which produces smoothly varying pixel‐scale estimates. We refer to this loss of information as “correlation loss.” Our approach upscales the core‐scale relation to the pixel scale using the model resolution matrix from the inversion, random field averaging, and spatial statistics of the geophysical property. Synthetic examples evaluate the utility of radar travel time tomography (RTT) and electrical‐resistivity tomography (ERT) for estimating water content. This work provides (1) a framework to assess tomograms for geologic parameter estimation and (2) insights into the different patterns of correlation loss for ERT and RTT. Whereas ERT generally performs better near boreholes, RTT performs better in the interwell region. Application of petrophysical models to the tomograms in our examples would yield misleading estimates of water content. Although the examples presented illustrate the problem of correlation loss in the context of near‐surface geophysical imaging, our results have clear implications for quantitative analysis of tomograms for diverse geoscience applications.

Journal of Geophysical Research B: Solid Earth↗
USGS2020

Conceptualizing alternate regimes in a large floodplain-river ecosystem

Regime shifts –persistent changes in the structure and function of an ecosystem - are well-documented in many ecosystems but remain poorly understood in floodplain-river ecosystems. We apply a resilience perspective to large floodplain-river ecosystems by presenting three examples of plausible sets of alternate regimes that are relevant to natural resource management interests within the Upper Mississippi River and Illinois River. These alternate regimes include: 1) a clear water and abundant vegetation regime vs. a turbid water and sparse vegetation regime in lentic, off-channel areas, 2) a diverse native fish community regime vs. an invasive-dominated fish community regime, and 3) a regime characterized by a diverse and dynamic mosaic of floodplain vegetation types vs. one characterized as a persistent invasive wet meadow monoculture. For each set of potential alternate regimes, we synthesize known or hypothesized feedback mechanisms that reinforce regimes, controlling variables that drive regime transitions, and restoration pathways. The conceptual models presented here provide a framework for synthesizing our understanding of the dynamics of this ecosystem and are relevant to other large floodplain-river ecosystems that face similar human pressures across the world. The models are currently being used to prioritize future research, test hypotheses, and inform restoration and management on the Upper Mississippi River and Illinois River. Through sharing our approach, we provide a case study in which we document an important step in operationalizing resilience concepts for the management of natural resources.

Minnesota, Wisconsin, Iowa, Illinois, Missouri↗
USGS1985

Three-dimensional trend mapping from wire-line logs

Mapping of lithofacies and porosities of stratigraphic units is complicated because these properties vary in three dimensions. The method of moments was proposed by Krumbein and Libby (1957) as a technique to aid in resolving this problem. Moments are easily computed from wireline logs and are simple statistics which summarize vertical variation in a log trace. Combinations of moment maps have proved useful in understanding vertical and lateral changes in lithology of sedimentary rock units. Although moments have meaning both as statistical descriptors and as mechanical properties, they also define polynomial curves which approximate lithologic changes as a function of depth. These polynomials can be fitted by least-squares methods, partitioning major trends in rock properties from finescale fluctuations. Analysis of variance yields the degree of fit of any polynomial and measures the proportion of vertical variability expressed by any moment or combination of moments. In addition, polynomial curves can be differentiated to determine depths at which pronounced expressions of facies occur and to determine the locations of boundaries between major lithologic subdivisions. Moments can be estimated at any location in an area by interpolating from log moments at control wells. A matrix algebra operation then converts moment estimates to coefficients of a polynomial function which describes a continuous curve of lithologic variation with depth. If this procedure is applied to a grid of geographic locations, the result is a model of variability in three dimensions. Resolution of the model is determined largely by number of moments used in its generation. The method is illustrated with an analysis of lithofacies in the Simpson Group of south-central Kansas; the three-dimensional model is shown as cross sections and slice maps. In this study, the gamma-ray log is used as a measure of shaliness of the unit. However, the method is general and can be applied, for example, to suites of neutron, density, or sonic logs to produce three-dimensional models of porosity in reservoir rocks. ?? 1985 Plenum Publishing Corporation.

Journal of the International Association for Mathe↗
USGS2025

Concentration-discharge relations and transient metal loads reveal spatiotemporal variability in solute-generation mechanisms in a mine-affected watershed

Concentration-discharge (CQ) relations are commonly used to understand geochemical and hydrologic controls on the generation of solutes in watersheds. Despite the widespread application of CQ relations, this technique has been infrequently applied to acid mine drainage (AMD) sites, but the CQ framework may allow mechanistic understanding of remedial outcomes such as impoundment of water within underground mines. Results of CQ analyses and changes in metal loads in an AMD affected watershed in Colorado, USA indicate that dissolved loads increased at many individual locations following water impoundment within mine workings. Although increased loads were observed at most individual locations, these increases were offset by a large decrease in loading from the largest mine. A loading analysis that included data from an instream monitoring location showed a statistically significant decrease in Fe and Zn after bulkhead emplacement, indicating a net positive effect of bulkheads. Streams generally displayed dilution CQ patterns whereas mines and springs showed either flushing or chemostatic patterns prior to bulkheading, which transitioned to chemostatic patterns following bulkheading, indicating a transition from dynamic to equilibrium geochemical processes. Saturation indices for sulfide and secondary minerals indicated that mines and springs were near equilibrium for phases including schwertmannite, fluorite, and gypsum. Saturation indices vary through time for mines suggesting progressive leaching of sulfide minerals as the mass of available minerals in the mine workings decreases. Together, these diverse analyses provide an integrated understanding of the variability in solute generating processes in this watershed and may inform remediation plans for similarly affected sites by indicating the nature of mineralogic controls on water quality.

Colorado↗
USGS2003

Application of two hydrologic models with different runoff mechanisms to a hillslope dominated watershed in the northeastern US: A comparison of HSPF and SMR

Differences in the simulation of hydrologic processes by watershed models directly affect the accuracy of results. Surface runoff generation can be simulated as either: (1) infiltration-excess (or Hortonian) overland flow, or (2) saturation-excess overland flow. This study compared the Hydrological Simulation Program - FORTRAN (HSPF) and the Soil Moisture Routing (SMR) models, each representing one of these mechanisms. These two models were applied to a 102 km2 watershed in the upper part of the Irondequoit Creek basin in central New York State over a seven-year simulation period. The models differed in both the complexity of simulating snowmelt and baseflow processes as well as the detail in which the geographic information was preserved by each model. Despite their differences in structure and representation of hydrologic processes, the two models simulated streamflow with almost equal accuracy. Since streamflow is an integral response and depends mainly on the watershed water balance, this was not unexpected. Model efficiency values for the seven-year simulation period were 0.67 and 0.65 for SMR and HSPF, respectively. HSPF simulated winter streamflow slightly better than SMR as a result of its complex snowmelt routine, whereas SMR simulated summer flows better than HSPF as a result of its runoff and baseflow processes. An important difference between model results was the ability to predict the spatial distribution of soil moisture content. HSPF aggregates soil moisture content, which is generally related to a specific pervious land unit across the entire watershed, whereas SMR predictions of moisture content distribution are geographically specific and matched field observations reasonably well. Important is that the saturated area was predicted well by SMR and confirmed the validity of using saturation-excess mechanisms for this hillslope dominated watershed. ?? 2003 Elsevier B.V. All rights reserved.

Journal of Hydrology↗
USGS2003

Automated calibration of a stream solute transport model: Implications for interpretation of biogeochemical parameters

The hydrologic processes of advection, dispersion, and transient storage are the primary physical mechanisms affecting solute transport in streams. The estimation of parameters for a conservative solute transport model is an essential step to characterize transient storage and other physical features that cannot be directly measured, and often is a preliminary step in the study of reactive solutes. Our study used inverse modeling to estimate parameters of the transient storage model OTIS (One dimensional Transport with Inflow and Storage). Observations from a tracer injection experiment performed on Uvas Creek, California, USA, are used to illustrate the application of automated solute transport model calibration to conservative and nonconservative stream solute transport. A computer code for universal inverse modeling (UCODE) is used for the calibrations. Results of this procedure are compared with a previous study that used a trial-and-error parameter estimation approach. The results demonstrated 1) importance of the proper estimation of discharge and lateral inflow within the stream system; 2) that although the fit of the observations is not much better when transient storage is invoked, a more randomly distributed set of residuals resulted (suggesting non-systematic error), indicating that transient storage is occurring; 3) that inclusion of transient storage for a reactive solute (Sr2+) provided a better fit to the observations, highlighting the importance of robust model parameterization; and 4) that applying an automated calibration inverse modeling estimation approach resulted in a comprehensive understanding of the model results and the limitation of input data.

Journal of the North American Benthological Societ↗
USGS2003

Source mechanism of long-period events at Kusatsu-Shirane Volcano, Japan, inferred from waveform inversion of the effective excitation functions

We investigate the source mechanism of long-period (LP) events observed at Kusatsu-Shirane Volcano, Japan, based on waveform inversions of their effective excitation functions. The effective excitation function, which represents the apparent excitation observed at individual receivers, is estimated by applying an autoregressive filter to the LP waveform. Assuming a point source, we apply this method to seven LP events the waveforms of which are characterized by simple decaying and nearly monochromatic oscillations with frequency in the range 1-3 Hz. The results of the waveform inversions show dominant volumetric change components accompanied by single force components, common to all the events analyzed, and suggesting a repeated activation of a sub-horizontal crack located 300 m beneath the summit crater lakes. Based on these results, we propose a model of the source process of LP seismicity, in which a gradual buildup of steam pressure in a hydrothermal crack in response to magmatic heat causes repeated discharges of steam from the crack. The rapid discharge of fluid causes the collapse of the fluid-filled crack and excites acoustic oscillations of the crack, which produce the characteristic waveforms observed in the LP events. The presence of a single force synchronous with the collapse of the crack is interpreted as the release of gravitational energy that occurs as the slug of steam ejected from the crack ascends toward the surface and is replaced by cooler water flowing downward in a fluid-filled conduit linking the crack and the base of the crater lake. ?? 2003 Elsevier Science B.V. All rights reserved.

Journal of Volcanology and Geothermal Research↗
USGS2026

Groundwater nutrient contributions to surface water in the upper Maumee River basin

Harmful algal blooms triggered by nutrient loading have negatively affected Lake Erie water quality for decades, but the role groundwater plays in transporting nutrients to tributary streams is unclear. Here, we apply a multi-scale approach to investigate groundwater nitrogen (N) and phosphorus (P) dynamics at the field scale, along a groundwater flow path, and at the Maumee River basin scale. Unsaturated zone monitoring indicates that macropore flow is the primary mechanism for N and P leaching with nitrate and orthophosphate (OP) concentrations from macropore samplers more than double those from soil-matrix samplers. Shallow (< 4.5 m) groundwater nitrate concentrations indicate the timing of saturated-zone nitrate fluxes is crop dependent with more nitrate mobilized during the maize growing season whereas concentrations below soybeans can remain elevated (> 3 mg/L) through winter months. A three-end-member hydrograph analysis for the upper Maumee River basin indicates that tile drains transport most subsurface N and P with greater than 15,000 kg per day of nitrate transported between January and May. Orthophosphate concentrations in tile drainage at an edge-of-field site correlate well with concentrations in a nearby tributary of the Maumee River. While runoff is the dominant source of OP loading, tile drainage and groundwater discharge contributed an estimated 36% of OP load during a 2018–2022 analysis period with most subsurface loading occurring between January and March. These findings highlight how best management practices (such as cover crops and tile drainage control structures) can be used to help limit leaching and tile effluent during the nongrowing season.

Indiana↗
USGS1992

Mechanisms and function of school formation in subyearling American shad (Alosa sapidissima)

Summary The behavioural repertoire of subyearling American shad ( Alosa sapidissima ) was studied in the laboratory to examine pathways leading to Aggregation and School formation at various group sizes (3–100 fish per 500-liter tank). Three hypotheses were tested: (1) schooling is the dominant behavioural activity of subyearling American shad, (2) development and maintenance of cohesive social behaviour is affected b group size, and (3) survival is a function of group size. Seventeen discrete behaviour patterns were observed. In relation to group size, the amount of time spent Schooling varied directly whereas Following and Aggregating varied inversely. The number of bouts of Parallel orient, Parallel swim, and Pivot also varied inversely with group size. A conceptual model was developed and tested to show how normal schooling behaviour is achieved in small or fragmented groups through intermediary modal action patterns. Sequence analysis showed behavioural diversity to diminish with group size, thereby increasing group cohesion and uniformity of response to stimuli. Probability of survival, determined from weekly mortality in experimental tanks, increased significantly with group size (0.87–0.98 for 3–100 fish). It is hypothesized that survival improves with an increase in the proportion of time spent Schooling, which may increase feeding efficiency, reduce stress, or reduce metabolic expenditures.

Journal of Applied Ichthyology↗
USGS2012

Patterns and controlling factors of species diversity in the Arctic Ocean

Aim The Arctic Ocean is one of the last near-pristine regions on Earth, and, although human activities are expected to impact on Arctic ecosystems, we know very little about baseline patterns of Arctic Ocean biodiversity. This paper aims to describe Arctic Ocean-wide patterns of benthic biodiversity and to explore factors related to the large-scale species diversity patterns. Location Arctic Ocean. Methods We used large ostracode and foraminiferal datasets to describe the biodiversity patterns and applied comprehensive ecological modelling to test the degree to which these patterns are potentially governed by environmental factors, such as temperature, productivity, seasonality, ice cover and others. To test environmental control of the observed diversity patterns, subsets of samples for which all environmental parameters were available were analysed with multiple regression and model averaging. Results Well-known negative latitudinal species diversity gradients (LSDGs) were found in metazoan Ostracoda, but the LSDGs were unimodal with an intermediate maximum with respect to latitude in protozoan foraminifera. Depth species diversity gradients were unimodal, with peaks in diversity shallower than those in other oceans. Our modelling results showed that several factors are significant predictors of diversity, but the significant predictors were different among shallow marine ostracodes, deep-sea ostracodes and deep-sea foraminifera. Main conclusions On the basis of these Arctic Ocean-wide comprehensive datasets, we document large-scale diversity patterns with respect to latitude and depth. Our modelling results suggest that the underlying mechanisms causing these species diversity patterns are unexpectedly complex. The environmental parameters of temperature, surface productivity, seasonality of productivity, salinity and ice cover can all play a role in shaping large-scale diversity patterns, but their relative importance may depend on the ecological preferences of taxa and the oceanographic context of regions. These results suggest that a multiplicity of variables appear to be related to community structure in this system.

Journal of Biogeography↗