Search USGS⌕ Search

SEARCH · Search USGS

Results for “Solutions Journal”

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 613 records · Page 34Linked to original sources

The enigma of the Arthur's Pass, New Zealand, earthquake: 1. Reconciling a variety of data for an unusual earthquake sequence

The 1994 Arthur's Pass earthquake ( M W 6.7) is the largest in a recent sequence of earthquakes in the central South Island, New Zealand. No surface rupture was observed, the aftershock distribution was complex, and routine methods of obtaining the faulting orientation of this earthquake proved contradictory. We use a range of data and techniques to obtain our preferred solution, which has a centroid depth of 5 km, M 0 =1.3 × 10 19 N m, and a strike, dip, and rake of 221°, 47°, 112°, respectively. Discrepancies between this solution and the Harvard centroid moment tensor, together with the Global Positioning System (GPS) observations and unusual aftershock distribution, suggest that the rupture may not have occurred on a planar fault. A second, strike slip, subevent on a more northerly striking plane is suggested by these data but neither the body wave modeling nor regional broadband recordings show any complexity or late subevents. We relocate the aftershocks using both one-dimensional and three-dimensional velocity inversions. The depth range of the aftershocks (1–10 km) agrees well with the preferred mainshock centroid depth. The aftershocks near the hypocenter suggest a structure dipping toward the NW, which we interpret to be the mainshock fault plane. This structure and the Harper fault, ∼15 km to the south, appear to have acted as boundaries to the extensive aftershock zone trending NNW-SSE. Most of the M L ≥ 5 aftershocks, including the two largest ( M L 6.1 and M L 5.7), clustered near the Harper fault and have strike slip mechanisms consistent with motion on this fault and its conjugates. Forward modeling of the GPS data suggests that a reverse slip mainshock, combined with strike slip aftershock faulting in the south, is able to match the observed displacements. The occurrence of this earthquake sequence implies that the level of seismic hazard in the central South Island is greater than previous estimates.

Arthur's Pass↗

Movement of sediment through a burned landscape: Sediment volume observations and model comparisons in the San Gabriel Mountains, California, USA

Post-wildfire changes to hydrologic and geomorphic systems can lead to widespread sediment redistribution. Understanding how sediment moves through a watershed is crucial for assessing hazards, developing debris flow inundation models, engineering sediment retention solutions, and quantifying the role that disturbances play in landscape evolution. In this study, we used terrestrial and airborne lidar to measure sediment redistribution in the 2016 Fish Fire, in the San Gabriel Mountains in southern California, USA. The lidar areas are in two adjacent watersheds, at spatial scales of 900 m 2 to 4 km 2 , respectively. Terrestrial lidar data were acquired prior to rainfall, and two subsequent surveys show erosional change after rainstorms. Two airborne lidar flights occurred (1) 7 months before, and (2) 14 months after the fire ignition, capturing the erosional effects after rainfall. We found hillslope erosion dominated the overall sediment budget in the first rainy season after wildfire. Only 7% of the total erosion came from the active channel bed and channel banks, and the remaining 93% of eroded sediment was derived from hillslopes. Within the channelized portion of the watershed erosion/deposition could be generally described with topographic metrics used in a stream power equation. Observed sediment volumes were compared with four empirical models and one process-based model. We found that the best predictions of sediment volume were obtained from an empirical model developed in the same physiographic region. Moreover, this study showed that post-wildfire erosion rates in the San Gabriel Mountains attain the same magnitude as millennial time scale bedrock erosion rates.

California↗

An improved method for quantifying soil macroporosity

Quantitative information on macroporosity is needed to predict water flow and solute transport in field soils. A method was developed for determining the number, shape, and size distribution of soil macropores. Horizontal serial sections sawed from paraffin-impregnated soil cores were photographed under ultraviolet (UV) light. Anthracene, mixed with the paraffin, fluoresces a bright bluish white under UV light and provides a sharp contrast between the soil matrix and the paraffin-filled pore space. Section photographs were converted to 256 level, grey-scale digital images using a flat-bed scanner. Image processing was used to classify each pixel in a digital image as pore space or soil matrix, to group the pore space pixels into pores, and to measure the area and perimeter of each pore. The method was able to measure pores with an equivalent radius ≥85 µm. Macroporosity in soil cores sampled form a tillage path and from an adjacent, undisturbed (notill) region was quantified. Tillage sections contained, on the average, 9.4 macropores/cm 2 with an equivalent macroporosity of 8%, while no-till sections contained 0.8 macropores/cm 2 with an equivalent macroporosity of 0.3%. Computed intrinsic permeabilities for tillage sections that included macropore information were significantly larger than values computed using micropore information alone, suggesting that macropores must be included in permeability calculations when the number of macropores is large. The developed method should be useful for quantifying macroporosity in nonskeletal soils.

Soil Science Society of America Journal↗

Direct measurement of asperity contact growth in quartz at hydrothermal conditions

Earthquake recurrence requires interseismic fault restrengthening which results from solid state deformation in room-temperature friction and indentation experiments. In contrast exhumed fault zones show solution-transport processes such as pressure solution and contact overgrowths influence fault zone properties . In the absence of fluid flow, overgrowths are driven by gradients in surface curvature where material is dissolved, diffuses, and precipitates at the contact without convergence normal to the contact. To determine the rate of overgrowth for quartz, we conducted single contact experiments in an externally heated pressure vessel. Convergence was continuously monitored using reflected-light interferometry through a long-working-distance microscope. Contact normal force was constant with an initial effective normal stress of 1.7 MPa, temperature was between 350 and 530{degree sign}C, and water pressure was constant at 150 MPa. Two control experiments were conducted: one dry at 425{degree sign}C and one bi-material (sapphire) at 425{degree sign}C and 150 MPa water pressure. No contact growth or convergence was observed in the controls. For wet single-phase contacts, growth was initially rapid and then decreased with time. No convergence was observed. Fluid inclusions indicate that the contact is not uniformly wetted. The contact is bounded by small regions of high aperture, reflecting local free-face dissolution as the source for the overgrowth. The apparent activation energy is ~125 kJ/mol. Extrapolation predicts rates of contact area increase orders of magnitude faster than in dry, room-temperature and hydrothermal friction experiments, suggesting that natural strength recovery near the base of the seismogenic zone could be dominated by contact overgrowth.

Journal of Geophysical Research B: Solid Earth↗

Estimation of suspended-sediment rating curves and mean suspended-sediment loads

Suspended-sediment loads are often estimated from an empirical relation between suspended-sediment load ( L ) and streamflow ( S ). This relation is usually defined as a power function, L = aS h , and is referred to as a suspended-sediment rating curve. This function can be formulated as either a linear or non-linear model to find the solution of the rating-curve parameters ( a and b ). Formulation of the power function as a linear model requires a logarithmic transformation to linearize the function and a subsequent correction for transformation bias. Rating-curve parameter estimates for both the bias-corrected, transformed-linear or non-linear models can be obtained by the method of least squares. Each model has distinct advantages and disadvantages. A unique solution of the parameters of the transformed-linear model may be obtained algebraically. These parameter estimates have some optimal properties when certain attainable conditions are met. However, the parameter estimates must be corrected for transformation bias when obtained this way. Parameter estimates obtained for the non-linear model do not require a correction for transformation bias. However, these estimates must be obtained by iterative methods which do not always converge to a solution. In addition, the residual errors of the non-linear model typically are not identically distributed throughout the range of streamflow values. This problem adversely affects the precision of the parameter estimates. Weighted non-linear least squares can be used to improve the parameter estimates for the non-linear model, but the weights must be approximated and their appropriate form may be difficult to determine. A simulation study was done to evaluate: (1) the accuracy and precision of parameter estimates for the bias-corrected, transformed-linear and non-linear models obtained by the method of least squares; (2) the accuracy of mean suspended-sediment loads calculated by the flow-duration, rating-curve method using model parameters obtained by the alternative methods. Parameter estimates obtained by least squares for the bias-corrected, transformed-linear model were considerably more precise than those obtained for the non-linear or weighted non-linear model. The accuracy of parameter estimates obtained for the biascorrected, transformed-linear and weighted non-linear model was similar and was much greater than the accuracy obtained by non-linear least squares. The improved parameter estimates obtained by the biascorrected, transformed-linear or weighted non-linear model yield estimates of mean suspended-sediment load calculated by the flow-duration, rating-curve method that are more accurate and precise than those obtained for the non-linear model.

Journal of Hydrology↗

Solving the dynamic rupture problem with different numerical approaches and constitutive laws

We study the dynamic initiation, propagation and arrest of a 2-D in-plane shear rupture by solving the elastodynamic equation by using both a boundary integral equation method and a finite difference approach. For both methods we adopt different constitutive laws: a slip-weakening (SW) law, with constant weakening rate, and rate- and state-dependent friction laws (Dieterich-Ruina). Our numerical procedures allow the use of heterogeneous distributions of constitutive parameters along the fault for both formulations. We first compare the two solution methods with an SW law, emphasizing the required stability conditions to achieve a good resolution of the cohesive zone and to avoid artificial complexity in the solutions. Our modelling results show that the two methods provide very similar time histories of dynamic source parameters. We point out that, if a careful control of resolution and stability is performed, the two methods yield identical solutions. We have also compared the rupture evolution resulting from an SW and a rate- and state-dependent friction law. This comparison shows that despite the different constitutive formulations, a similar behaviour is simulated during the rupture propagation and arrest. We also observe a crack tip bifurcation and a jump in rupture velocity (approaching the P-wave speed) with the Dieterich-Ruina (DR) law. The rupture arrest at a barrier (high strength zone) and the barrier-healing mechanism are also reproduced by this law. However, this constitutive formulation allows the simulation of a more general and complex variety of rupture behaviours. By assuming different heterogeneous distributions of the initial constitutive parameters, we are able to model a barrier-healing as well as a self-healing process. This result suggests that if the heterogeneity of the constitutive parameters is taken into account, the different healing mechanisms can be simulated. We also study the nucleation phase duration Tn, defined as the time necessary for the crack to reach the half-length Ic. We compare the Tn values resulting from distinct simulations calculated using different constitutive laws and different sets of constitutive parameters. Our results confirm that the DR law provides a different description of the nucleation process than the SW law adopted in this study. We emphasize that the DR law yields a complete description of the rupture process, which includes the most prominent features of SW.

Geophysical Journal International↗

Infrared spectroscopic examination of the interaction of urea with the naturally occurring zeolite clinoptilolite

Infrared spectroscopy has shown for the first time that the naturally occurring zeolite clinoptilolite can absorb urea, (NH 2 ) 2 CO, under ambient conditions from either aqueous or ethanolic solutions. The two strongest NH stretching bands at 3441 and 3344 cm −1 in pure, solid urea shift to higher frequency (about 3504 and 3401 cm −1 ) after absorption. Two of the four urea bands in the 1800-1300 cm −1 range (at 1683 and 1467 cm −1 ) undergo marked downward shifts to about 1670 and 1445 cm −1 . The other two bands show little change in frequency. The strong band at 1602 cm −1 , however, diminishes in intensity to little more than an ill-defined shoulder on the 1626-cm −1 peak. When clinoptilolite is heated to 450 °C and then treated with molten urea (ca. 140 °C) for several minutes, and finally washed twice with ethanol to remove excess unreacted urea, further changes become apparent in the spectrum of the urea-treated clinoptilolite. The two NH stretching bands broaden without significant change in frequency. Two new bands appear in the midfrequency range at 1777 (weak) and 1719 (medium strong) cm −1 . Of the four original midfrequency peaks, the one at 1602 cm −1 is now absent. Two others (1627 and 1440 cm −1 ) exhibit little change, while the fourth has broadened and shifted down to 1663 cm −1 , where it appears as a shoulder on the band at 1627 cm −1 . Both treatments clearly induce interaction between urea and the zeolite which seems to result in significant modifications in the nature of the hydrogen bonding of the substrate.

Microchemical Journal↗

Mechanism of smoke-induced seed germination in a post-fire chaparral annual

1 Smoke-stimulated germination in the post-fire flora of California chaparral does not appear to be triggered by nitrate. Application of freshly prepared unbuffered KNO 3 solutions (pH c. 6.2) failed to enhance germination of five populations of Emmenanthe penduliflora or one Phacelia grandiflora population, regardless of light or stratification conditions. 2 KNO 3 buffered at acidic pH (or unbuffered solutions equilibrated with atmospheric CO 2 ) did induce germination, but KNO 3 solutions at pH 7 failed to induce germination. Induction of germination is therefore not due to the nitrate ion per se , but rather to high [H + ], although buffered controls gave weak germination at low pH, suggesting a role for H + plus nitrate. However, other anions such as sulphate were equally as effective as nitrate at breaking dormancy. 3 The germination response to KNO 3 was affected by the type of filter paper used and this may be linked to differences in pH. 4 NO 2 , at concentrations present in biomass smoke, was highly effective at inducing germination, and other oxidizing agents also induced germination. 5 Several growth regulators, including nitrite and gibberellin, were stimulatory only at acidic pH, but KCN was stimulatory across a broad pH range. 6 Germination decreased at smoke exposures longer than a few minutes. Also, smoked water samples effective at breaking dormancy were acidic and were less effective when buffered to pH >7. 7 Physical scarification of the seed coat induced germination but the effect was not due to penetration of a water barrier, or to enhanced oxygen uptake or to wound responses such as CO 2 or ethylene production. 8 Different effects of the gibberellin inhibitor CCC (chlorocholine chloride) suggested that the mechanisms of scarification-induced and smoke-induced germination may differ. 9 We conclude that either oxidizing gases in smoke and/or acids generated on burnt sites play a role in germination of post-fire annuals in chaparral.

Journal of Ecology↗

Lower temperature terminations of the three-phase region plagioclase-alkali feldspar-liquid

Geological and experimental evidence indicate that the three-phase field, plagioclase-alkalifeldspar-liquid, may terminate in several different ways. The possible terminations have been developed from Schreinemakers' rules governing the disappearance of three-phase fields. In igncous rocks, these different terminations may arise from variations in the relative amounts of additional components in magmas, or from changes of total pressure, or from structural changes that effect the extent of solid solution in feldspar.The three-phase region originates from the intersection of the solidus and the feldspar solvus. The available evidence regarding this intersection is reviewed, as is the evidence for the existence and form of the boundary curve on the feldspar liquidus. The data are used to project a series of isobaric polythermal and isobaric isothermal diagrams for each possible termination . Subsequent discussions relate the theoretical arguments to the natural evidence, suggest the more probable geological environments of some of the terminations , and indicate characteristic features of each termination .It may be possible to identify the type of termination involved in the crystallization of some rocks. The necessary data are the amounts and compositions of both kinds of coexisting feldspars and of the feldspar components of the coexisting liquid at one or more stages of the crystallization process. Volcanic rocks are most suitable for such studies.The terminations explain many of the compositional relationships possible between feldspar crystals and silicate melt under geological and experimental conditions. These compositional relationships are fundamental to understanding the crystallization of igneous rocks and the formation of melts by partial fusion. Other geological applications include a simple explanation for some resorbed feldspars, the separation, correlation, and comparison of porphyritic rock units, interpretation of compositional changes of successive zones of zoned feldspars, and mantling of one feldspar by another. With all types of termination , extensive fractionation yields a liquid rich in Or and Ab that may crystallize to alkali feldspar or feldspars.

Journal of Petrology↗

Effects of karst features on circulation of water in carbonate rocks in coastal areas

The normal balance between fresh water in coastal aquifers and sea water applies also to carbonate-rock aquifers that have been karstified, but there are local modifications in the balance that need to be considered. Uneven distribution of permeability, expressed by a network of solution channels bounded by relatively impermeabler rock, causes an uneven distribution of head of the environmental water along the seacoast. Where sinkholes and (or) vertical solution shafts below sea level penetrate the aquifer, the fresh ground water may discharge through these karst features if the fresh-water head is greater than that of the salt water. However, under some conditions the salt-water head may exceed that of the fresh water, and the direction of movement is reversed as sea water flows into the aquifer. This sea-water flow into the aquifer occurs (1) where sinkholes, acting as “cased wells,” penetrate less permeable rock before reaching a lateral solution channel and (2) where (or when) the fresh-water head is less than that required to balance the salt water. On Andros Island, Bahamas, the range in tide (as much as 5 feet) from low tide to high tide is sufficient to cause such a reversal locally. During low tide the salt-water head becomes sufficiently low that the ground-water head exceeds that of the sea water, and the ground water flows through the sinkholes to the ocean floor; during high tide sea water flows in the sinkholes. In the Adriatic Sea along the coast of Yugoslavia, apparently the fresh-water head is sufficient to produce perennial springs in some localities, but in other areas, as in the Bay of Kastela near Split, the fresh-water head becomes low enough during some seasons that the flow is reversed and salt water enters the aquifer through the karst features. The development of sinkholes and other karst features near the present coast and extending below sea level occurred generally during a low stand of the Pleistocene sea when the top of the saturated zone stood lower than the bottom of the deepest sinkholes or natural wells. Integrated evaluations of (1) the distribution of permeability in coastal karst regions and (2) the principles relating to the dynamic balance between fresh aquifer water and sea water are leading to better knowledge of methods that may salvage much karst water which is lost to the sea.

Split↗

Factors controlling the geochemical evolution of fumarolic encrustations, Valley of Ten Thousand Smokes, Alaska

Factor and canonical correlation analysis of geochemical data from eight fossil fumaroles suggest that six major factors controlled the formation and evolution of fumarolic encrustations on the 1912 ash-flow sheet in the Valley of Ten Thousand Smokes (VTTS). The six-factor solution model explains a large proportion (low of 74% for Ni to high of 99% for Si) of the individual element data variance. Although the primary fumarolic deposits have been degraded by secondary alteration reactions and up to 75 years of weathering, the relict encrustations still preserve a signature of vapor-phase element transport. This vapor-phase transport probably occurred as halide or oxyhalide species and was significant for As, Sb and Br. At least three, and possibly four, varied temperature leaching events affected the fumarolic deposits. High-temperature gases/liquids heavily altered the ejecta glass and mineral phases adjacent to the fumarolic conduit. As the fumaroles cooled. Fe-rich acidic condensate leached the ejecta and primary fumarolic deposits and resulted in the subsequent precipitation of Fe-hydroxides and/or Fe-oxides. Low- to ambient-temperature leaching and hydration reactions generated abundant hydrated amorphous phases. Up to 87% of the individual element data variance is apparently controlled by the chemistry of the ejecta on which the relict encrustations are found. This matrix chemistry factor illustrates that the primary fumarolic minerals surrounding the active VTTS vents observed by earlier workers have been effectively removed by the dissolution reactions. Element enrichment factors calculated for the VTTS relict encrustations support the statistical factor interpretations. On the average, the relict encrustations are enriched, relative to visibly unaltered matrix protolith, in As, Br, Cr, Sb, Cu, Ni, Pb, Fe, and LOI (an indirect measure of sample H2O content). ?? 1993.

Journal of Volcanology and Geothermal Research↗

A distribution-based parameterization for improved tomographic imaging of solute plumes

Difference geophysical tomography (e.g. radar, resistivity and seismic) is used increasingly for imaging fluid flow and mass transport associated with natural and engineered hydrologic phenomena, including tracer experiments, in situ remediation and aquifer storage and recovery. Tomographic data are collected over time, inverted and differenced against a background image to produce ‘snapshots’ revealing changes to the system; these snapshots readily provide qualitative information on the location and morphology of plumes of injected tracer, remedial amendment or stored water. In principle, geometric moments (i.e. total mass, centres of mass, spread, etc.) calculated from difference tomograms can provide further quantitative insight into the rates of advection, dispersion and mass transfer; however, recent work has shown that moments calculated from tomograms are commonly biased, as they are strongly affected by the subjective choice of regularization criteria. Conventional approaches to regularization (Tikhonov) and parametrization (image pixels) result in tomograms which are subject to artefacts such as smearing or pixel estimates taking on the sign opposite to that expected for the plume under study. Here, we demonstrate a novel parametrization for imaging plumes associated with hydrologic phenomena. Capitalizing on the mathematical analogy between moment-based descriptors of plumes and the moment-based parameters of probability distributions, we design an inverse problem that (1) is overdetermined and computationally efficient because the image is described by only a few parameters, (2) produces tomograms consistent with expected plume behaviour (e.g. changes of one sign relative to the background image), (3) yields parameter estimates that are readily interpreted for plume morphology and offer direct insight into hydrologic processes and (4) requires comparatively few data to achieve reasonable model estimates. We demonstrate the approach in a series of numerical examples based on straight-ray difference-attenuation radar monitoring of the transport of an ionic tracer, and show that the methodology outlined here is particularly effective when limited data are available.

Geophysical Journal International↗

Thermodynamic properties and crystal structure refinement of ferricopiapite, coquimbite, rhomboclase, and Fe2(SO4)3(H2O)5

Enthalpies of formation of ferricopiapite [nominally Fe4.67(SO4)6(OH)2 (H2O)20]. coquimbite [Fe2(SO4)3(H2O)9], rhomboclase [(H3O)Fe(SO4)2 (H2O)3], and Fe2(SO4)3(H2O)5 were measured by acid (5 N HCl) solution calorimetry. The samples were characterized by wet chemical analyses and synchrotron powder X-ray diffraction (XRD). The refinement of XRD patterns gave lattice parameters, atomic positions, thermal factors, and occupancies of the sites. The calculated formulae differ slightly from the nominal compositions: Fe4.78(SO4)6 (OH)2.34(H2O)20.71 (ferricopiapite), (Fe1.47Al0.53)(SO4)3 (H2O)9.65 (coquimbite), (H3O)1.34Fe(SO4)2.17 (H2O)3.06 (rhomboclase), and Fe2(SO4)3 (H2O)5.03. All thermodynamic data are given per mole of these formulae. The measured standard enthalpies (in kJ/mol) of formation from the elements (crystalline Fe, Al, S, and ideal gases O2 and H2) at T = 298.15 K are -4115.8??4.1 [Fe2(SO4)3 (H2O)5.03], -12045.1??9.2 (ferricopiapite), -5738.4??3.3 (coquimbite), and -3201.1??2.6 (rhomboclase). Standard entropy (S??) was estimated as a sum of entropies of oxide, hydroxide, and sulfate components. The estimated S?? (in J/mol.K) values for the iron sulfates are 488.2 [Fe2(SO4)3 (H2O)5.03], 1449.2 (ferricopiapite), 638.3 (coquimbite), and 380.1 (rhomboclase). The calculated Gibbs free energies of formation (in kJ/mol) are -3499.7??4.2 [Fe2(SO4)3 (H2O)5.03], -10089.8??9.3 (ferricopiapite), -4845.6??3.3 (coquimbite), and -2688.0??2.7 (rhomboclase). These results combined with other available thermodynamic data allow construction of mineral stability diagrams in the FeIII2(SO4)3-FeII SO4-H2O system. One such diagram is provided, indicating that the order of stability of ferric sulfate minerals with decreasing pH in the range of 1.5 to -0.5 is: hydronium jarosite, ferricopiapite, and rhomboclase. ?? 2006 E. Schweizerbart'sche Verlagsbuchhandlung.

European Journal of Mineralogy↗

Stress field variations in the Swiss Alps and the northern Alpine foreland derived from inversion of fault plane solutions

This study is devoted to a systematic analysis of the state of stress of the central European Alps and northern Alpine foreland in Switzerland based on focal mechanisms of 138 earthquakes with magnitudes between 1 and 5. The most robust feature of the results is that the azimuth of the minimum compressive stress, S 3 , is generally well constrained for all data subsets and always lies in the NE quadrant. However, within this quadrant, the orientation of S 3 changes systematically both along the structural strike of the Alpine chain and across it. The variation in stress along the mountain belt from NE to SW involves a progressive, counterclockwise rotation of S 3 and is most clear in the foreland, where it amounts to 45°–50°. This pattern of rotation is compatible with the disturbance to the stress field expected from the indentation of the Adriatic Block into the central European Plate, possibly together with buoyancy forces arising from the strongly arcuate structure of the Moho to the immediate west of our study area. Across the Alps, the variation in azimuth of S 3 is defined by a progressive, counterclockwise rotation of about 45° from the foreland in the north across the Helvetic domain to the Penninic nappes in the south and is accompanied by a change from a slight predominance of strike-slip mechanisms in the foreland to a strong predominance of normal faulting in the high parts of the Alps. The observed rotation can be explained by the perturbation of the large-scale regional stress by a local uniaxial deviatoric tension with a magnitude similar to that of the regional differential stress and with an orientation perpendicular to the strike of the Alpine belt. The tensile nature and orientation of this stress is consistent with the “spreading” stress expected from lateral density changes due to a crustal root beneath the Alps.

European Alps↗

Cost-benefit analysis for evacuation decision-support: Challenges and possible solutions for applications in areas of distributed volcanism

During a volcanic crisis, evacuation is the most effective mitigation measure to preserve life. However, the decision to call an evacuation is typically complex and challenging, in part due to uncertainties related to the behaviour of the volcano. Cost-benefit analysis (CBA) can support decision-makers: this approach compares the cost of evacuating versus the expected loss from not evacuating, expressed as a ‘break-even’ probability of fatality. Here we combine CBA with a Bayesian Event Tree for Short-term Volcanic Hazard (BET_VHst) to create an evacuation decision-support tool to identify locations that are cost-beneficial to evacuate in the event of volcanic unrest within a distributed volcanic field. We test this approach with the monogenetic Auckland Volcanic Field (AVF), situated beneath the city of Auckland, New Zealand. We develop a BET_VHst for the AVF, extending a recently revised Bayesian Event Tree for Eruption Forecasting (BET_EF) to consider the eruptive style, phenomena produced, and the impact exceedance probability as a function of distance. The output of the BET_VHst is a probability of volcanic hazard impact at a given location. Furthermore, we propose amending the weight of the monitoring component within the BET_VHst framework to a transitional parameter, addressing limitations identified in a previous study. We examine how three possible transitional monitoring component weights affect the spatial vent likelihood and subsequent BET_VHst outputs, compared to the current default weight. For the CBA, we investigate four thresholds, based on two evacuation durations and two different estimates for the value of life that determine the cost of not evacuating. The combinations of CBA and BET_VHst are tested using a synthetic unrest dataset to define an evacuation area for each day. While suitable evacuation areas were identified, there are further considerations required before such an approach can be applied operationally to support crisis management.

Auckland Volcanic Field vent↗

Cryopreservation of Sperm from the Endangered Colorado Pikeminnow

We developed methods for the cryopreservation of sperm of the endangered Colorado pikeminnow Ptychocheilus lucius. Sperm were collected from a captive broodstock population of Colorado pikeminnow reared and maintained at the Dexter National Fish Hatchery and Technology Center. Our objectives were to (1) evaluate the effects on sperm motility of 24-h storage in Hanks' balanced salt solution (HBSS); (2) characterize sperm motility and duration; (3) examine the relationship between sperm motility and osmotic pressure; (4) examine the effect of four cryoprotectants (dimethyl sulfoxide [DMSO], dimethyl acetamide [DMA], glycerol, and methanol [MeOH] at two concentrations [5% and 10%]) on postthaw motility; and (5) compare the effect of two cooling rates (40??C/ min and 4??C/min) on postthaw motility. The sperm samples diluted with HBSS retained higher motility (mean ??SD, 77 ?? 22%; n = 9) than did undiluted samples (12 ?? 30%; n = 9) after 24 h of storage. When exposed to HBSS at 274 mosmols/kg or more, few sperm became motile (???1%). Exposure to HBSS at 265 mosmols/kg elicited threshold activation (defined as 10% motility), and maximum motility (>95%) was observed at 93 mosmols/ kg. The maximum motility of sperm was observed within 10 s after activation with deionized water, and sperm remained motile for 57 s. The sperm that were cooled at a rate of 40??C/min and cryopreserved with 5% MeOH retained higher postthaw motility (56 ?? 13%) than did sperm cryopreserved with DMSO, DMA, or glycerol (at 5% and 10%). When the sperm samples were cooled at a rate of 4??C/min, sperm cryopreserved with MeOH (5% or 10%) or DMSO (5% or 10%) retained the highest postthaw motilities (???14%). The use of cryopreserved sperm can assist hatchery managers in the production of fish, provide for the long-term conservation of genetic resources, and assist in the recovery of endangered species such as the Colorado pikeminnow.

North American Journal of Aquaculture↗

Water-level changes induced by local and distant earthquakes at Long Valley caldera, California

Distant as well as local earthquakes have induced groundwater-level changes persisting for days to weeks at Long Valley caldera, California. Four wells open to formations as deep as 300 m have responded to 16 earthquakes, and responses to two earthquakes in the 3-km-deep Long Valley Exploratory Well (LVEW) show that these changes are not limited to weathered or unconsolidated near-surface rocks. All five wells exhibit water-level variations in response to earth tides, indicating they can be used as low-resolution strainmeters. Earthquakes induce gradual water-level changes that increase in amplitude for as long as 30 days, then return more slowly to pre-earthquake levels. The gradual water-level changes are always drops at wells LKT, LVEW, and CH-10B, and always rises at well CW-3. At a dilatometer just outside the caldera, earthquake-induced strain responses consist of either a step followed by a contractional strain-rate increase, or a transient contractional signal that reaches a maximum in about seven days and then returns toward the pre-earthquake value. The sizes of the gradual water-level changes generally increase with earthquake magnitude and decrease with hypocentral distance. Local earthquakes in Long Valley produce coseismic water-level steps; otherwise the responses to local earthquakes and distant earthquakes are indistinguishable. In particular, water-level and strain changes in Long Valley following the 1992 M7.3 Landers earthquake, 450 km distant, closely resemble those initiated by a M4.9 local earthquake on November 22, 1997, during a seismic swarm with features indicative of fluid involvement. At the LKT well, many of the response time histories are identical for 20 days after each earthquake, and can be matched by a theoretical solution giving the pore pressure as a function of time due to diffusion of a nearby, instantaneous, pressure drop. Such pressure drops could be produced by accelerated inflation of the resurgent dome by amounts too small to be detected by the two-color electronic distance-measuring network. Opening-mode displacement in the south moat, inferred to have followed a M4.9 earthquake on November 22, 1997, could also create extensional strain on the dome and lead to water-level changes similar to those following dome inflation. Contractional strain that could account for earthquake-induced water-level rises at the CW-3 well is inconsistent with geodetic observations. We instead attribute these water-level rises to diffusion of elevated fluid pressure localized in the south moat thermal aquifer. For hydraulic diffusivities appropriate to the upper few hundred meters at Long Valley, an influx of material at temperatures of 300°C can thermally generate pressure of 6 m of water or more, an order of magnitude larger than needed to account for the CW-3 water-level rises. If magma or hot aqueous fluid rises to within 1 km of the surface in the eastern part of the south moat, then hydraulic diffusivities are high enough to allow fluid pressure to propagate to CW-3 on the time scale observed. The data indicate that seismic waves from large distant earthquakes can stimulate upward movement of fluid in the hydrothermal system at Long Valley.

California↗

Viability and fertilizing capacity of cryopreserved sperm from three North American acipenseriform species: A retrospective study

Populations of sturgeon across the globe are threatened due to unregulated harvest and habitat loss, and the status varies among species across North America. Ready access to viable and functional sperm would contribute to recovery programmes for these species. In this study, we examined the motility, viability (cell membrane integrity) of cryopreserved sperm from three North American acipenseriform species and fertilizing capacity. Milt samples were collected from captive shortnose sturgeon (Acipenser brevirostrum), wild paddlefish (Polyodon spathula) and pallid sturgeon (Scaphirhynchus albus) and cryopreserved using combinations of Modified Tsvetkova's (MT) extender, Original Tsvetkova's extender, and modified Hanks' balanced salt solution, along with the cryoprotectants methanol (MeOH) or dimethyl sulfoxide (DMSO). A dual-staining technique using the fluorescent stains SYBR-14 and propidium iodide was employed with flow cytometry to determine the percentages of spermatozoa that were viable by virtue of having intact membranes. The percentage of viable spermatozoa ranged from 5% to 12% in shortnose sturgeon, 30-59% in paddlefish, and 44-58% in pallid sturgeon. In the first experiment with shortnose sturgeon sperm, methanol allowed for higher values for dependent variables than did DMSO, and sperm viability generally correlated with post-thaw motility. However, fertilization rate, neurulation, or hatching rates were independent from these factors. In the second experiment with shortnose sturgeon, 5% MeOH combined with MT yielded higher values for all parameters tested than the other combinations: viability was correlated with motility, fertilization rate, and hatching rate. Overall, viability and post-thaw motility was not affected by the use of hyperosmotic extenders (OT) or cryoprotectants (DMSO), but their use decreased fertilization percentages. For paddlefish sperm (experiment 3), MT combined with 10% MeOH was clearly a good choice for cryopreservation; viability and motility results were correlated, but independent of fertilization. For pallid sturgeon sperm (experiment 4), MT with 5-10% MeOH showed significantly higher sperm quality and fertilization parameters. Membrane integrity can be used as a predictor of fertilization by cryopreserved sperm, however additional sperm quality parameters, supplementary to motility and membrane integrity, would be useful in the refining and optimizing cryopreservation protocols with acipenseriform sperm. ?? 2008 Blackwell Verlag, Berlin.

Journal of Applied Ichthyology↗