Geology ReportsSearch

SEARCH · Geology Reports

Results for “Nature Methods”

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 397 records · Page 22Linked to original sources

Analysis of spurious oscillation modes for the shallow water and Navier-Stokes equations

The origin and nature of spurious oscillation modes that appear in mixed finite element methods are examined. In particular, the shallow water equations are considered and a modal analysis for the one-dimensional problem is developed. From the resulting dispersion relations we find that the spurious modes in elevation are associated with zero frequency and large wave number (wavelengths of the order of the nodal spacing) and consequently are zero-velocity modes. The spurious modal behavior is the result of the finite spatial discretization. By means of an artificial compressibility and limiting argument we are able to resolve the similar problem for the Navier-Stokes equations. The relationship of this simpler analysis to alternative consistency arguments is explained. This modal approach provides an explanation of the phenomenon in question and permits us to deduce the cause of the very complex behavior of spurious modes observed in numerical experiments with the shallow water equations and Navier-Stokes equations. Furthermore, this analysis is not limited to finite element formulations, but is also applicable to finite difference formulations.

Computers and Fluids

Resolving structural uncertainty in natural resources management using POMDP approaches

In recent years there has been a growing focus on the uncertainties of natural resources management, and the importance of accounting for uncertainty in assessing management effectiveness. This paper focuses on uncertainty in resource management in terms of discrete-state Markov decision processes (MDP) under structural uncertainty and partial observability. It describes the treatment of structural uncertainty with approaches developed for partially observable resource systems. In particular, I show how value iteration for partially observable MDPs (POMDP) can be extended to structurally uncertain MDPs. A key difference between these process classes is that structurally uncertain MDPs require the tracking of system state as well as a probability structure for the structure uncertainty, whereas with POMDPs require only a probability structure for the observation uncertainty. The added complexity of the optimization problem under structural uncertainty is compensated by reduced dimensionality in the search for optimal strategy. A solution algorithm for structurally uncertain processes is outlined for a simple example in conservation biology. By building on the conceptual framework developed for POMDPs, natural resource analysts and decision makers who confront structural uncertainties in natural resources can take advantage of the rapid growth in POMDP methods and approaches, and thereby produce better conservation strategies over a larger class of resource problems.

Ecological Modelling

Simultaneous determination of Cr(iii) and Cr(vi) using reversed-phased ion-pairing liquid chromatography with dynamic reaction cell inductively coupled plasma mass spectrometry

A method for the simultaneous determination of Cr(iii) and Cr(vi) species in waters, soil leachates and synthetic bio-fluids is described. The method uses reversed-phase ion-pairing liquid chromatography to separate the chromium species and a dynamic reaction cell (DRC??) equipped ICP-MS for detection of chromium. Separation of the chromium species is carried out in less than 2 min. Cr(iii) is complexed with ethylenediaminetetraacetic acid (EDTA) prior to separation by mixing samples with the mobile phase containing 2.0 mM tetrabutylammonium hydroxide (TBAOH), 0.5 mM EDTA (dipotassium salt), and 5% (vol/vol) methanol, adjusted to pH 7.6. The interfering 40Ar 12C+ background peak at mass 52 was reduced by over four orders of magnitude to less than 200 cps by using 0.65 mL min-1 ammonia as a reaction gas and an RPq setting on the DRC of 0.75. Method detection limits (MDLs) of 0.09 ??g L-1 for Cr(iii) and 0.06 ??g L-1 for Cr(vi) were obtained based on peak areas at mass 52 for 50 ??L injections of low level spikes. Reproducibility at 2 ??g L-1 was 3% RSD for 5 replicate injections. The tolerance of the method to various levels of common cations and anions found in natural waters and to matrix constituents found in soil leachates and simulated gastric and lung fluids was tested by performing spike recovery calculations for a variety of samples. ?? The Royal Society of Chemistry.

Journal of Analytical Atomic Spectrometry

Occam's shadow: levels of analysis in evolutionary ecology - where to next?

Evolutionary ecology is the study of evolutionary processes, and the ecological conditions that influence them. A fundamental paradigm underlying the study of evolution is natural selection. Although there are a variety of operational definitions for natural selection in the literature, perhaps the most general one is that which characterizes selection as the process whereby heritable variation in fitness associated with variation in one or more phenotypic traits leads to intergenerational change in the frequency distribution of those traits. The past 20 years have witnessed a marked increase in the precision and reliability of our ability to estimate one or more components of fitness and characterize natural selection in wild populations, owing particularly to significant advances in methods for analysis of data from marked individuals. In this paper, we focus on several issues that we believe are important considerations for the application and development of these methods in the context of addressing questions in evolutionary ecology. First, our traditional approach to estimation often rests upon analysis of aggregates of individuals, which in the wild may reflect increasingly non-random (selected) samples with respect to the trait(s) of interest. In some cases, analysis at the aggregate level, rather than the individual level, may obscure important patterns. While there are a growing number of analytical tools available to estimate parameters at the individual level, and which can cope (to varying degrees) with progressive selection of the sample, the advent of new methods does not reduce the need to consider carefully the appropriate level of analysis in the first place. Estimation should be motivated a priori by strong theoretical analysis. Doing so provides clear guidance, in terms of both (i) assisting in the identification of realistic and meaningful models to include in the candidate model set, and (ii) providing the appropriate context under which the results are interpreted. Second, while it is true that selection (as defined) operates at the level of the individual, the selection gradient is often (if not generally) conditional on the abundance of the population. As such, it may be important to consider estimating transition rates conditional on both the parameter values of the other individuals in the population (or at least their distribution), and population abundance. This will undoubtedly pose a considerable challenge, for both single- and multi-strata applications. It will also require renewed consideration of the estimation of abundance, especially for open populations. Thirdly, selection typically operates on dynamic, individually varying traits. Such estimation may require characterizing fitness in terms of individual plasticity in one or more state variables, constituting analysis of the norms of reaction of individuals to variable environments. This can be quite complex, especially for traits that are under facultative control. Recent work has indicated that the pattern of selection on such traits is conditional on the relative rates of movement among and frequency of spatially heterogeneous habitats, suggesting analyses of evolution of life histories in open populations can be misleading in some cases.

Journal of Applied Statistics

Detection and attribution of temperature changes in the mountainous Western United States

Large changes in the hydrology of the western United States have been observed since the mid-twentieth century. These include a reduction in the amount of precipitation arriving as snow, a decline in snowpack at low and midelevations, and a shift toward earlier arrival of both snowmelt and the centroid (center of mass) of streamflows. To project future water supply reliability, it is crucial to obtain a better understanding of the underlying cause or causes for these changes. A regional warming is often posited as the cause of these changes without formal testing of different competitive explanations for the warming. In this study, a rigorous detection and attribution analysis is performed to determine the causes of the late winter/early spring changes in hydrologically relevant temperature variables over mountain ranges of the western United States. Natural internal climate variability, as estimated from two long control climate model simulations, is insufficient to explain the rapid increase in daily minimum and maximum temperatures, the sharp decline in frost days, and the rise in degree-days above 0??C (a simple proxy for temperature driven snowmelt). These observed changes are also inconsistent with the model-predicted responses to variability in solar irradiance and volcanic activity. The observations are consistent with climite simulations that include the combined effects of anthropogenic greenhouse gases and aerosols. It is found that, for each temperature variable considered, an anthropogenic signal is identifiable in observational fields. The results are robust to uncertainties in model-estimated fingerprints and natural variability noise, to the choice of statistical down-scaling method, and to various processing options in the detection and attribution method. ?? 2008 American Meteorological Society.

Journal of Climate

Multivariate model of female black bear habitat use for a Geographic Information System

Simple univariate statistical techniques may not adequately assess the multidimensional nature of habitats used by wildlife. Thus, we developed a multivariate method to model habitat-use potential using a set of female black bear ( Ursus americanus ) radio locations and habitat data consisting of forest cover type, elevation, slope, aspect, distance to roads, distance to streams, and forest cover type diversity score in the Ozark Mountains of Arkansas. The model is based on the Mahalanobis distance statistic coupled with Geographic Information System (GIS) technology. That statistic is a measure of dissimilarity and represents a standardized squared distance between a set of sample variates and an ideal based on the mean of variates associated with animal observations. Calculations were made with the GIS to produce a map containing Mahalanobis distance values within each cell on a 60- × 60-m grid. The model identified areas of high habitat use potential that could not otherwise be identified by independent perusal of any single map layer. This technique avoids many pitfalls that commonly affect typical multivariate analyses of habitat use and is a useful tool for habitat manipulation or mitigation to favor terrestrial vertebrates that use habitats on a landscape scale.

Arkansas

Thorium deposits of the United States — Energy resources for the future?

Many nations are exploring new ways to meet their growing energy supply needs, with a particular focus upon methods that produce lower carbon dioxide emissions compared to traditional oil, natural gas, and coal power plants. As a result, thorium-based nuclear power has experienced renewed attention as a potential energy source. Thus, it benefits the United States and other countries to identify and evaluate their indigenous thorium resources. This report describes the geology and resources of the principal thorium districts of the United States.

Circular

Spatial distribution and trends in trace elements, polycyclic aromatic hydrocarbons, organochlorine pesticides, and polychlorinated biphenyls in Lake Worth sediment, Fort Worth, Texas

In spring 2000, the Texas Department of Health issued a fish consumption advisory for Lake Worth in Fort Worth, Texas, because of elevated concentrations of polychlorinated biphenyls (PCBs) in fish. In response to the advisory and in cooperation with the U.S. Air Force, the U.S. Geological Survey collected 21 surficial sediment samples and three gravity core sediment samples to assess the spatial distribution and historical trends of selected hydrophobic contaminants, including PCBs, and to determine, to the extent possible, sources of hydrophobic contaminants to Lake Worth. Compared to reference (background) concentrations in the upper lake, elevated PCB concentrations were detected in the surficial sediment samples collected in Woods Inlet, which receives surface runoff from Air Force facilities and urban areas. Gravity cores from Woods Inlet and from the main part of the lake near the dam indicate that the concentrations of PCBs were three to five times higher in the 1960s than in 2000. A regression method was used to normalize sediment concentrations of trace elements for natural variations and to distinguish natural and anthropogenic contributions to sediments. Concentrations of several trace elements—cadmium, chromium, copper, lead, and zinc—were elevated in sediments in Woods Inlet, along the shoreline of Air Force facilities, and in the main lake near the dam. Concentrations of these five trace elements have decreased since 1970. Polycyclic aromatic hydrocarbons also were elevated in the same areas of the lake. Concentrations of total polycyclic aromatic hydrocarbons, normalized with organic carbon, were mostly stable in the upper lake but steadily increased near the dam, except for small decreases since 1980. The Woods Inlet gravity core showed the largest increase of the three core sites beginning about 1940; total polycyclic aromatic hydrocarbon concentrations in post-1940 sediments from the core showed three apparent peaks about 1960, 1984, and 2000. The concentrations of organochlorine pesticides were low relative to consensus-based sediment-quality guidelines and either decreased or remained constant since 1970. The two likely sources of hydrophobic contaminants to the lake are urban areas around the lake and the drainage area of Meandering Road Creek that contributes runoff to Woods Inlet and includes Air Force facilities.

Texas

A demonstration of the instream flow incremental methodology, Shenandoah River, Virginia

Current and projected demands on the water resources of the Shenandoah River have increased concerns for the potential effect of these demands on the natural integrity of the Shenandoah River system. The Instream Flow Incremental Method (IFIM) process attempts to integrate concepts of water-supply planning, analytical hydraulic engineering models, and empirically derived habitat versus flow functions to address water-use and instream-flow issues and questions concerning life-stage specific effects on selected species and the general well being of aquatic biological populations. The demonstration project also sets the stage for the identification and compilation of the major instream-flow issues in the Shenandoah River Basin, development of the required multidisciplinary technical team to conduct more detailed studies, and development of basin specific habitat and flow requirements for fish species, species assemblages, and various water uses in the Shenandoah River Basin. This report presents the results of an IFIM demonstration project, conducted on the main stem Shenandoah River in Virginia, during 1996 and 1997, using the Physical Habitat Simulation System (PHABSIM) model. Output from PHABSIM is used to address the general flow requirements for water supply and recreation and habitat for selected life stages of several fish species. The model output is only a small part of the information necessary for effective decision making and management of river resources. The information by itself is usually insufficient for formulation of recommendations regarding instream-flow requirements. Additional information, for example, can be obtained by analysis of habitat time-series data, habitat duration data, and habitat bottlenecks. Alternative-flow analysis and habitat-duration curves are presented.

Virginia

Estimating riparian and agricultural evapotranspiration by reference crop evapotranspiration and MODIS Enhanced Vegetation Index

Dryland river basins frequently support both irrigated agriculture and riparian vegetation and remote sensing methods are needed to monitor water use by both crops and natural vegetation in irrigation districts. We developed an algorithm for estimating actual evapotranspiration (ET a ) based on the Enhanced Vegetation Index (EVI) from the Moderate Resolution Imaging Spectrometer (MODIS) sensor on the EOS-1 Terra satellite and locally-derived measurements of reference crop ET (ET o ). The algorithm was calibrated with five years of ETa data from three eddy covariance flux towers set in riparian plant associations on the upper San Pedro River, Arizona, supplemented with ETa data for alfalfa and cotton from the literature. The algorithm was based on an equation of the form ET a = ET o [a(1 − e −bEVI ) − c], where the term (1 − e −bEVI ) is derived from the Beer-Lambert Law to express light absorption by a canopy, with EVI replacing leaf area index as an estimate of the density of light-absorbing units. The resulting algorithm capably predicted ET a across riparian plants and crops (r 2 = 0.73). It was then tested against water balance data for five irrigation districts and flux tower data for two riparian zones for which season-long or multi-year ET a data were available. Predictions were within 10% of measured results in each case, with a non-significant (P = 0.89) difference between mean measured and modeled ET a of 5.4% over all validation sites. Validation and calibration data sets were combined to present a final predictive equation for application across crops and riparian plant associations for monitoring individual irrigation districts or for conducting global water use assessments of mixed agricultural and riparian biomes.

Remote Sensing

Advancing best practices for the analysis of the vulnerability of military installations in the Pacific Basin to coastal flooding under a changing climate – RC-2644

Coastal flooding takes many forms, ranging from major flooding associated with storms to minor flooding associated with exceptionally high tides and other oceanic and atmospheric phenomena on storm-free days. A major societal challenge is to understand and predict how flood magnitude and frequency will manifest at particular places and times, now and in the future. Of particular interest here is how coastal flooding will impact Department of Defense (DoD) installations. In response to this need, this work aims to advance the practical application of statistical and other analytical techniques that can be used to assess the exposure, and ultimately the vulnerability, of built and natural environments to the impacts of coastal flooding. A variety of methods are described and applied to assess exposure. This includes tide gauge station-based diagnosis and prognosis of patterns and trends of Still Water Level, techniques to characterize the expression of ‘lesser extremes’ (e.g., sub-annual to subdecadal event probabilities), and region-wide analysis that improves upon results obtained from conventional single-tide gauge analyses. A novel hybrid statistical and dynamical modeling approach is applied to the analysis of Total Water Levels, necessary for exposure assessment along shorelines exposed to wave action. The hybrid exposure assessment modeling approach is incorporated into a broader mission-based protocol for the assessment of resilience to coastal flooding at the installation level. Demonstrated via an exemplar assessment, which takes into account functional (lost day) as well as financial impacts (lost dollars), the protocol meets the demand for an actionable characterization of how DoD installations will be affected by coastal flooding and improves DoD’s ability to make informed decisions about how to adapt to its effects. The methods described, evaluated, and applied here, including innovative approaches and proof-of-concept products developed through this work, are incorporated into and considered within an analytical framework that serves as guidance as to their relative merits with respect to coastal flood exposure assessment in various circumstances and settings, and illustrates best practices. This will provide engineers, scientists and other practitioners with an enhanced capability to generate information that can be used to support area-wide assessment related to climate adaptation planning and disaster risk reduction as well as site-specific analysis related to design and maintenance of facilities and infrastructure. While the focus is on a select set of DoD sites in the Pacific Basin, the results have broad applicability nationally as well as globally.

Final Report

The presence of silicate melt may enhance rates of cation diffusion in olivine

Olivine is commonly used as a ‘crystal clock’ to extract timescales relevant to pre-eruptive perturbations within mafic magmatic systems. Diffusion chronometry applications require accurate calibrations for the rates at which Fe-Mg or other commonly measured elements like Ni, Mn, and Ca diffuse through the crystal lattice. In the past, these rates have been mainly characterized using solid-solid diffusion couple experiments involving olivine single crystals, thin films, or powder sources. Despite the presence of melt surrounding olivine in natural magmatic systems, very few experiments involving magma have been performed, largely because controlling interface reactions is difficult. For this study, we carried out olivine-melt diffusion experiments as a test of the diffusion chronometry method, and to determine whether the presence of melt influences the calculated timescales. To approximate a natural system, we incorporated small natural Kīlauea and San Carlos olivine seeds within a natural Kīlauea basalt and tracked diffusive re-equilibration through time. To better control interface reactions, after some equilibration period at an initial superliquidus temperature of 1290 °C, the runs were rapidly cooled to form a rim and left to dwell at various final temperatures (1200, 1220, 1240, 1255 °C) for 6–84 h. Concentration gradients for Fe-Mg, Mn, Ni, Ca were measured, and the step-wise nature of the core-rim transition was ascertained using slow diffusing elements like P or Al. When these gradients are modeled using published diffusivities, the timescales retrieved are typically 10 times longer than the actual experiment durations. Thus, measured diffusivities are an order of magnitude faster than those previously obtained in olivine-solid source experiments, but they are in excellent agreement with the only two other melt-olivine datasets. We explore reasons for why melt-bearing olivine diffusion experiments tend to yield faster rates. The possible effects of (1) growth during diffusion, (2) diffusion during any initial dissolution step, and (3) extended tube or planar defects at the interface on calculated diffusivities are all considered but found to be inconsequential. Instead, we argue that additional point defects (vacancies) are likely created at the interface by higher concentrations in elements like Al or H in the basalt melt compared to other solid couple diffusant sources. Future applications of diffusion chronometry in olivine may require a complete re-evaluation of published diffusivities using melt-bearing experimental configurations.

Earth and Planetary Science Letters

Direct and indirect estimates of natural mortality for Chesapeake Bay blue crab

Analyses of the population dynamics of blue crab Callinectes sapidus have been complicated by a lack of estimates of the instantaneous natural mortality rate ( M ). We developed the first direct estimates of M for this species by solving Baranov's catch equation for M given estimates of annual survival rate and exploitation rate. Annual survival rates were estimated from a tagging study on adult female blue crabs in Chesapeake Bay, and female-specific exploitation rates for the same stock were estimated by comparing commercial catches with abundances estimated from a dredge survey. We also used eight published methods based on life history parameters to calculate indirect estimates of M for blue crab. Direct estimates of M for adult females in Chesapeake Bay for the years 2002–2004 ranged from 0.42 to 0.87 per year and averaged 0.71 per year. Indirect estimates of M varied considerably depending on life history parameter inputs and the method used. All eight methods yielded values for M between 0.99 and 1.08 per year, and six of the eight methods yielded values between 0.82 and 1.35 per year. Our results indicate that natural mortality of blue crab is higher than previously believed, and we consider M values between 0.7 and 1.1 per year to be reasonable for the exploitable stock in Chesapeake Bay. Remaining uncertainty about M makes it necessary to evaluate a range of estimates in assessment models.

Transactions of the American Fisheries Society

Airborne geophysics for geologic mapping of critical mineral systems in the United States southern midcontinent

The increased demand for clean energy technology and a significant reliance on foreign supply chains have given impetus to understanding critical mineral systems and locating potential resources within the United States. At least thirteen critical mineral-bearing systems have been identified throughout the U.S. southern Midcontinent (Hofstra and Kreiner, 2020) but much of the region’s geologic framework is concealed by vegetation and sedimentary cover that hinder traditional geologic mapping efforts. Airborne geophysical data provide an effective way to overcome these obstacles and to provide additional insight into the deeper structures that underlie shallow mineralization. However, legacy airborne magnetic and radiometric data were collected using now-outdated instruments and methods, inconsistent survey parameters, and large flight-line spacings resulting in low-resolution data that present challenges to regional-scale study and interpretation. Over the last decade, the U.S. Geological Survey Earth Mapping Resources Initiative (EMRI) and National Cooperative Geologic Mapping Program have conducted a series of high-resolution airborne magnetic and radiometric surveys across the southern Midcontinent (Fig. 1) as part of an effort to improve understanding of the geophysical framework and natural resource potential in the region. These surveys are designed using modern survey methods and instruments with consistent parameters for flight-line spacing and flight height relative to magnetic sources. The EMRI airborne surveys are planned in collaboration with State geological surveys based on focus areas (Dicken et al., 2022) according to the presence of or potential for critical mineral deposits. High-resolution airborne magnetic and radiometric data cover focus areas such as the southeast Missouri iron metallogenic province and South-Central iron-oxide-apatite (IOA) – iron-oxide-copper-gold (IOCG) province, the Magnet Cove alkaline-carbonatite complex, the Midwest Permian ultramafic dike district, the Illinois-Kentucky fluorspar district, and several Mississippi Valley-type lead-zinc deposits and districts (Fig. 1). These focus areas represent known deposits or prospective host systems of critical minerals including rare earth elements (REEs), platinum-group elements (PGEs), cobalt, lithium, fluorspar, niobium, titanium, vanadium, lead, zinc, gallium, germanium, and many more. Other significant geologic and geophysical features covered include the Reelfoot rift, the New Madrid seismic zone, the Illinois basin, the Arkoma basin, the South-Central magnetic lineament, and the Kentucky-Tennessee magnetic anomaly (Fig. 1). This presentation focuses on new airborne magnetic and radiometric data with continuous coverage across parts of six states, preliminary interpretations, examples of geologic mapping applications, and discussion of newly discovered magnetic anomalies and follow-up investigations.

Conference Paper

Stability and uncertainty of finite-fault slip inversions: Application to the 2004 Parkfield, California, earthquake

The 2004 Parkfield, California, earthquake is used to investigate stability and uncertainty aspects of the finite-fault slip inversion problem with different a priori model assumptions. We utilize records from 54 strong ground motion stations and 13 continuous, 1-Hz sampled, geodetic instruments. Two inversion procedures are compared: a linear least-squares subfault-based methodology and a nonlinear global search algorithm. These two methods encompass a wide range of the different approaches that have been used to solve the finite-fault slip inversion problem. For the Parkfield earthquake and the inversion of velocity or displacement waveforms, near-surface related site response (top 100 m, frequencies above 1 Hz) is shown to not significantly affect the solution. Results are also insensitive to selection of slip rate functions with similar duration and to subfault size if proper stabilizing constraints are used. The linear and nonlinear formulations yield consistent results when the same limitations in model parameters are in place and the same inversion norm is used. However, the solution is sensitive to the choice of inversion norm, the bounds on model parameters, such as rake and rupture velocity, and the size of the model fault plane. The geodetic data set for Parkfield gives a slip distribution different from that of the strong-motion data, which may be due to the spatial limitation of the geodetic stations and the bandlimited nature of the strong-motion data. Cross validation and the bootstrap method are used to set limits on the upper bound for rupture velocity and to derive mean slip models and standard deviations in model parameters. This analysis shows that slip on the northwestern half of the Parkfield rupture plane from the inversion of strong-motion data is model dependent and has a greater uncertainty than slip near the hypocenter.

California

Methods for peak-flow frequency analysis for streamgages in or near Montana, North Dakota, South Dakota, and Wyoming

The U.S. Geological Survey, in cooperation with the Montana Department of Natural Resources and Conservation, North Dakota Department of Water Resources, South Dakota Department of Transportation, and the Wyoming Water Development Office, has developed standard methods of peak-flow frequency analysis for studies in Montana, North Dakota, South Dakota, and Wyoming. These methods describe the implementation of national flood frequency guidelines described in Bulletin 17C ( https://doi.org/10.3133/tm4B5 ) for the four States and deviations from Bulletin 17C standard procedures to accommodate unusual hydrologic conditions. A U.S. Geological Survey data release accompanying this report ( https://doi.org/10.5066/P1WHRK8H ) provides example peak-flow frequency analyses for selected streamgages in the study area. The methods described in this report can be used to publish similar data releases for other streamgages in the study area.

Montana, North Dakota, South Dakota, Wyoming

Methods for determining manning's coefficients for Illinois streams

Determination of Manning's coefficient, n, for natural streams remains a challenge in practices. One source for determining the n-values that has received practitioners' attention is presenting the n-values determined from field data (measured discharge and water-surface slope) in combination of photographs and site descriptions (ancillary information). Further improvements in the visual approach can be made in presenting site characteristics and describing site ancillary information. In this manner, users can use the presented information for sites of interest with similar features. This approach in a current project on the subject for Illinois streams is discussed.

Conference Paper

Simple, efficient allocation of modelling runs on heterogeneous clusters with MPI

In scientific modelling and computation, the choice of an appropriate method for allocating tasks for parallel processing depends on the computational setting and on the nature of the computation. The allocation of independent but similar computational tasks, such as modelling runs or Monte Carlo trials, among the nodes of a heterogeneous computational cluster is a special case that has not been specifically evaluated previously. A simulation study shows that a method of on-demand (that is, worker-initiated) pulling from a bag of tasks in this case leads to reliably short makespans for computational jobs despite heterogeneity both within and between cluster nodes. A simple reference implementation in the C programming language with the Message Passing Interface (MPI) is provided.

Environmental Modelling and Software