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 91 records · Page 5Linked to original sources

Causal analysis of fire regime drivers in California

Background Understanding the relative contribution of climate and human factors to wildfires is critical for managing risk across California’s diverse ecosystems, in the United States (US). Aims We propose a model that distinguishes between proximate and ultimate drivers of fire regimes and apply it to a century of fire and climate data to assess regional variation in causal mechanisms. Methods We analyzed fire statistics (1910–2021) alongside climate and weather data, stratifying the state by 10 ecoregions. Key results Northern forests had the strongest correlation with the proximate factor fuel aridity, ultimately due to climate. Fire rotation intervals exceeded 100 years, implicating woody fuel accumulation as an additional factor. Lightning ignitions occurred in decadal bursts, with dense strike events potentially overwhelming fire-fighting resources. Lower elevation/latitude foothill ecoregions experienced highest fire activity following wet winters and springs, implicating control by herbaceous fuel loads and a negative effect of global warming on future fires. Human ignitions dominate in these ecoregions, and population growth contributes to expansion of powerlines, a major ignition source. Conclusions While climate change may increase fire activity in forested ecoregions, its role is less pronounced in non-forested ecoregions, where human ignition sources are the dominant factor. Implications Different areas within ecoregions may require different management actions that reflect the specific proximate and ultimate factors at play.

California↗

An eddy-resolving numerical model to study turbulent flow, sediment and bed evolution using detached eddy simulation in a lateral separation zone at the field-scale

Turbulence-resolving simulations elucidate key elements of fluid dynamics and sediment transport in fluvial environments. This research presents a feasible strategy for applying state-of-the-art computational fluid mechanics to the study of sediment transport and morphodynamic processes in lateral separation zones, which are common features in canyon rivers where massive lateral flow separation causes large-scale turbulence that controls sediment erosion and deposition. An eddy-resolving model was developed and tested at the field-scale, coupling a viscous flow and sediment transport solver using Detached Eddy Simulation techniques. A morphodynamic model was applied to the viscous flow/sediment solver to calculate erosion and deposition. A simulation of turbulence was performed at the grid resolution for a straight channel to determine the relative contributions of modeled and resolved diffusivity. The time-dependent, energetically important, correlative, non-stationary signals of the simulated quantities were captured at the lateral separation zone. Strong periodic signals featured by high amplitude were found at the separation zone, while low frequency pulsations were observed at the reattachment zone of the lateral separation zone. Interactions between the eddies and the loose bed boundaries resulted in erosion of sediment at the main channel followed by deposition at the primary eddy and eddy bars. tions elucidate key elements of fluid dynamics and sediment transport in fluvial environments. This research presents a feasible strategy for applying state-of-the-art computational fluid mechanics to the study of sediment transport and morphodynamic processes in lateral separation zones, which are common features in canyon rivers where massive lateral flow separation causes large-scale turbulence that controls sediment erosion and deposition. An eddy-resolving model was developed and tested at the field scale, coupling viscous flow and sediment transport solver using Detached Eddy Simulation (DES) techniques. A morphodynamic model was applied to the viscous flow/sediment solver to calculate erosion and deposition. A simulation of turbulence was performed at the grid resolution for a straight channel to determine the relative contributions of modeled and resolved diffusivity. The time-dependent, energetically important, correlative non-stationary signals of the simulated quantities were captured at the lateral separation zone. Strong periodic signals featured by high amplitude were found at the separation zone, while low frequency pulsations were observed at the reattachment zone of the lateral separation zone. Interactions between the eddies and the loose bed boundaries resulted in massive erosion of sediment at the main channel followed by deposition at the primary eddy and eddy bars.

Arizona↗

W phase source inversion for moderate to large earthquakes (1990-2010)

Rapid characterization of the earthquake source and of its effects is a growing field of interest. Until recently, it still took several hours to determine the first-order attributes of a great earthquake (e.g. M w ≥ 7.5), even in a well-instrumented region. The main limiting factors were data saturation, the interference of different phases and the time duration and spatial extent of the source rupture. To accelerate centroid moment tensor (CMT) determinations, we have developed a source inversion algorithm based on modelling of the W phase, a very long period phase (100–1000 s) arriving at the same time as the P wave. The purpose of this work is to finely tune and validate the algorithm for large-to-moderate-sized earthquakes using three components of W phase ground motion at teleseismic distances. To that end, the point source parameters of all M w ≥ 6.5 earthquakes that occurred between 1990 and 2010 (815 events) are determined using Federation of Digital Seismograph Networks, Global Seismographic Network broad-band stations and STS1 global virtual networks of the Incorporated Research Institutions for Seismology Data Management Center. For each event, a preliminary magnitude obtained from W phase amplitudes is used to estimate the initial moment rate function half duration and to define the corner frequencies of the passband filter that will be applied to the waveforms. Starting from these initial parameters, the seismic moment tensor is calculated using a preliminary location as a first approximation of the centroid. A full CMT inversion is then conducted for centroid timing and location determination. Comparisons with Harvard and Global CMT solutions highlight the robustness of W phase CMT solutions at teleseismic distances. The differences in M w rarely exceed 0.2 and the source mechanisms are very similar to one another. Difficulties arise when a target earthquake is shortly (e.g. within 10 hr) preceded by another large earthquake, which disturbs the waveforms of the target event. To deal with such difficult situations, we remove the perturbation caused by earlier disturbing events by subtracting the corresponding synthetics from the data. The CMT parameters for the disturbed event can then be retrieved using the residual seismograms. We also explore the feasibility of obtaining source parameters of smaller earthquakes in the range 6.0 ≤M w < 6.5. Results suggest that the W phase inversion can be implemented reliably for the majority of earthquakes of M w = 6 or larger.

Geophysical Journal International↗

Bluff evolution along coastal drumlins: Boston Harbor Islands, Massachusetts

A series of partially drowned drumlins forms the backbone of the inner islands within Boston Harbor. The shoreline of these rounded glacial deposits is composed of actively retreating bluffs formed by continual wave attack. Comparisons of bluffs reveal variability in their height and lateral extent, as well as in the dominant mechanism causing their retreat. Two processes are responsible for bluff erosion and yield distinct bluff morphologies: (1) wave attack undercuts the bluff and causes episodic slumping, yielding planar bluff slopes, and (2) subaerial processes such as rainfall create irregular slopes characterized by rills and gullies. We propose a model of drumlin bluff evolution that is based on processes of erosion and physical characteristics such as bluff height, slope morphology, and the orientation of the bluff with respect to the long axis of the drumlin and its topographic crest. The four phases of drumlin bluff evolution consist of (1) initial formation of bluff, with retreat dominated by wave notching and slumping processes; (2) rill and gully development as bluff heights exceed 10 m and slumped sediment at bluff base inhibits wave attack; (3) return of wave notching and slumping as bluff heights decrease; and (4) final development of boulder retreat lag as last remnants of drumlin are eroded by wave action. These phases capture the important physical processes of drumlin evolution in Boston Harbor and could apply to other eroding coastal drumlin deposits.

Massachusetts↗

Source parameters of microearthquakes at Mount St Helens (USA)

We estimate the source parameters for a selection of microearthquakes that occurred at Mount St Helens in the period 1995–1998. Excluding the activity of 2004 September, this time period includes the most intense episode of earthquake activity since the last dome-building eruption in 1986 October. 200 seismograms were processed to obtain seismic moments, source radii, stress drops and average fault slip. The source parameters were determined from the spectral analysis of P waves, after correction for attenuation and site effects. In particular, P -wave quality ( Q p ) and site ( S ) factors have been previously calculated in the frequency ranges 2–7 Hz and 18–30 Hz. Because it was impossible to perform corrections for Q p and S over the whole spectrum we applied a new approach, based on the notion of ‘ holed spectrum ’, to estimate spectral parameters. The term ‘ holed spectrum ’ indicates a spectrum lacking corrected spectral amplitude values at certain frequencies. We carried out a statistical study to verify that dealing with the ‘ holed spectrum ’ does not lead to significant differences in the estimates of spectral parameters. We also investigated the dependence of spectral parameters (low-frequency level, corner frequency and high-frequency decay) on the bandwidth of spectral hole, and defined the threshold values for three different spectral models. Displacement ‘ holed spectra ’, corrected by attenuation and site response, are then used to determine spectral parameters in order to calculate seismic source parameters. Seismic moments range from 10 17 to 10 19 dyne-cm, source dimensions from 100 to 350 m, and average fault slip from 0.003 to 0.1 cm. Self-similarity seems to break down in that stress drops are very low (0.1–1 bars). We postulate that seismicity is associated with a brittle shear failure mechanism occurring in a highly heterogeneous material under a relatively low stress regime.

Washington↗

Identification of seismo-volcanic regimes at Whakaari/White Island (New Zealand) via systematic tuning of an unsupervised classifier

We present an algorithm based on Self-Organizing Maps (SOM) and k-means clustering to recognize patterns in a continuous 12.5-year tremor time series recorded at Whakaari/White Island volcano, New Zealand (hereafter referred to as Whakaari). The approach is extendable to a variety of volcanic settings through systematic tuning of the classifier. Hyperparameters are evaluated by statistical means, yielding a combination of “ideal” SOM parameters for the given data set. Extending from this, we applied a Kernel Density Estimation approach to automatically detect changes within the observed seismicity. We categorize the Whakaari seismic time series into regimes representing distinct volcano-seismic states during recent unrest episodes at Whakaari (2012/2013, 2016, and 2019). There is a clear separation in classification results between background regimes and those representing elevated levels of unrest. Onset of unrest is detected by the classifier 6 weeks before the August 2012 eruption, and ca. 3.5 months before the December 2019 eruption, respectively. Regime changes are corroborated by changes in commonly monitored tremor proxies as well as with reported volcanic activity. The regimes are hypothesized to represent diverse mechanisms including: system pressurization and depressurization, degassing, and elevated surface activity. Labeling these regimes improves visualization of the 2012/2013 and 2019 unrest and eruptive episodes. The pre-eruptive 2016 unrest showed a contrasting shape and nature of seismic regimes, suggesting differing onset and driving processes. The 2016 episode is proposed to result from rapid destabilization of the shallow hydrothermal system, while rising magmatic gases from new injections of magma better explain the 2012/2013 and 2019 episodes.

Whakaari, White Island↗

Rapid intrusion of magma into wet rock: Groundwater flow due to pore pressure increases

Analytical and numerical solutions are developed to simulate the pressurization, expansion, and flow of groundwater contained within saturated, intact host rocks subject to sudden heating from the planar surface of an igneous intrusion. For most rocks, water diffuses more rapidly than heat, assuring that groundwater is not heated along a constant-volume pressure path and that thermal expansion and pressurization adjacent to the intrusion drives a flow that extends well beyond the heated region. The forcing parameter for pressurization and flow is α Δ T , where α is a thermal expansion coefficient reflecting the overall expansion of water heated through the temperature difference Δ T between the initial ambient and intrusive values. Pore pressure increases due to heating are greatest when the intrusion is emplaced rapidly and where the intrusive contact is impervious to groundwater contained in stiff, impermeable rocks with high thermal diffusivities and porosities. The maximum velocity of water flowing in pores decays with the inverse square root of time and is insensitive to hydraulic properties of the host rocks. Pressures are lessened and flow directions are reversed with the onset of hydrothermal convection. This occurs at times ranging from hours to weeks after onset of intrusion. As magma rises into near-surface rocks, steam can be generated. Solutions indicate that pressure increases and velocities are sensitive to the overall amount of expansion rather than the behavior of the water-steam transition. Both the overall thermal expansion coefficient α and the temperature difference Δ T are greater in shallow (<1 km) environments than in deep (∼5 km) ones. Thus, for rocks with similar transport properties, pressure increases due to heating are greatest in shallow environments. Although solutions can be applied to rocks with a wide variety of properties, pressure increases are calculated for compliant quartz-rich sedimentary rocks with a porosities between 1 and 20% and permeabilities between 1 darcy and 1 μdarcy, subject to temperature increases of 500 and 1000 K at depths ranging from 0.1 to 5 km in a region of hydrostatic pressures and normal geothermal gradient. Such rocks, with porosities greater than 5%, permeabilities less than a 0.1 mdarcy, and drained hydrostatic compressibilities of 10 −4 /MPa, undergo pressure increases greater than 10 MPa (100 bars)for conditions typical of water table depths of 2.5 km and heating to 500 K above ambient. Similar rocks, but with permeabilities less than 1 mdarcy, undergo pressure increases of 10 MPa for conditions typical of 1 km water table depth. Rocks commonly considered to be good aquifers undergo pressure increases of less than 1 MPa, primarily because of their high permeability. Although these estimates neglect the effects of fracturing and brecciation that may accompany such pressure increases, calculations indicate that pressure increases due to heating of cool groundwater can lead to failure of host rocks by a phreatic mechanism.

Journal of Geophysical Research Solid Earth↗

Morphological plasticity following species-specific recognition and competition in two perennial grasses

Morphological characteristics and biomass allocation of two perennial grasses, Pseudoroegneria spicata (Pursh) A. Löve ssp. spicata (bluebunch wheatgrass) and Agropyron desertorum (Fisch. ex Link) Schult. (crested wheatgrass), were compared under different competition and nutrient treatments. The competitive responses of two plants grown in containers under field conditions were assessed in monocultures and mixtures in two experiments using different scales of nutrient application. In the Small-Scale Experiment, a localized fertilization was applied in the rooting zone between two plants; in the Large-Scale Experiment the entire container was supplied with nutrients. Agropyron responded more vigorously to fertilization than did Pseudoroegneria , but based on the relative performance of Agropyron in monoculture and mixture, it was not superior to Pseudoroegneria in resource competition. Pseudoroegneria was apparently able to recognize neighboring plants as either conspecifics or individuals of the other species. The responses included changes in shoot architecture, root morphology, and allocation between roots and shoots. Agropyron generally did not exhibit such morphological flexibility. In field plot plantings of 4-yr-old tussocks similar shoot differences were seen in Pseudoroegneria. There was, however, no indication of superior resource competition for Agropyron. Thus, any early advantage of Agropyron in vigorous growth of young plants in response to nutrients was apparently lost by the time the plants had reached this stage of development. Morphological and allocation flexibility of Pseudoroegneria may have compensated for slower, less vigorous growth. If species-specific recognition and morphological plasticity are common in nature, this complicates our attempts to understand mechanisms of competition.

American Journal of Botany↗

Generation and verification of theoretical rating curves in the Whitewater River basin, Kansas

[1] A new method for generating stage-discharge relations (rating curves) for geomorphically stable channels is presented and applied to two streams in the Whitewater River basin, Kansas. The approach converts measurements of stage into discharge using a fluid mechanically based model. The model does not use empirical roughness coefficients, such as Manning coefficients, but rather determines channel roughness from field measurements of the (1) channel geometry, (2) the physical roughness of the bed, banks, and floodplain, and (3) the vegetation density on the banks and floodplain. These measurements are used to calculate explicitly the drag on the small-scale topographic features on the boundary, the drag on the stems and branches of woody vegetation, and the friction on the bed, banks, and floodplain. The theoretical rating curves produced by the model for two study reaches, which are near U.S. Geological Survey (USGS) streamflow-gauging stations, are in good agreement with direct measurements of discharge made by the USGS. Our method has the potential of providing accurate estimates of stream flows less expensively than conventional gauging methods. In addition, the method can be used to obtain more accurate discharge estimates than conventional indirect methods for determining discharge, which are based on estimates of Manning's coefficient.

Journal of Geophysical Research F: Earth Surface↗

Seismicity around Parkfield correlates with static shear stress changes following the 2003 M w 6.5 San Simeon earthquake

Earthquakes trigger other earthquakes, but the physical mechanism of the triggering is currently debated. Most studies of earthquake triggering rely on earthquakes listed in catalogs, which are known to be incomplete around the origin times of large earthquakes and therefore missing potentially triggered events. Here we apply a waveform matched-filter technique to systematically detect earthquakes along the Parkfield section of the San Andreas Fault from 46 days before to 31 days after the nearby 2003 M w 6.5 San Simeon earthquake. After removing all possible false detections, we identify ~8 times more earthquakes than in the Northern California Seismic Network catalog. The newly identified events along the creeping section of the San Andreas Fault show a statistically significant decrease following the San Simeon main shock, which correlates well with the negative static stress changes (i.e., stress shadow) cast by the main shock. In comparison, the seismicity rate around Parkfield increased moderately where the static stress changes are positive. The seismicity rate changes correlate well with the static shear stress changes induced by the San Simeon main shock, suggesting a low friction in the seismogenic zone along the Parkfield section of the San Andreas Fault.

California↗

Postwildfire soil‐hydraulic recovery and the persistence of debris flow hazards

Deadly and destructive debris flows often follow wildfire, but understanding of changes in the hazard potential with time since fire is poor. We develop a simulation‐based framework to quantify changes in the hydrologic triggering conditions for debris flows as postwildfire infiltration properties evolve through time. Our approach produces time‐varying rainfall intensity‐duration thresholds for runoff‐ and infiltration‐generated debris flows with physics‐based hydrologic simulations that are parameterized with widely available hydroclimatic, vegetation reflectance, and soil texture data. When we apply our thresholding protocol to a test case in the San Gabriel Mountains (California, USA), the results are consistent with existing regional empirical thresholds and rainstorms that caused runoff‐ and infiltration‐generated debris flows soon after and three years following a wildfire, respectively. We find that the hydrologic triggering mechanisms for the two observed debris flow types are coupled with the effects of fire on the soil saturated hydraulic conductivity. Specifically, the rainfall intensity needed to generate debris flows via runoff increases with time following wildfire while the rainfall duration needed to produce debris flows via subsurface pore‐water pressures decreases. We also find that variations in soil moisture, rainfall climatology, median grain size, and root reinforcement could impact the median annual probability of postwildfire debris flows. We conclude that a simulation‐based method for calculating rainfall thresholds is a tractable approach to improve situational awareness of debris flow hazard in the years following wildfire. Further development of our framework will be important to quantify postwildfire hazard levels in variable climates, vegetation types, and fire regimes.

Journal of Geophysical Research: Earth Surface↗

Semiautomated estimates of directivity and related source properties of small to moderate southern California earthquakes using second seismic moments

We develop a semiautomated method for estimating with second seismic moments the directivity, rupture area, duration, and centroid velocity of earthquakes. The method is applied to 41 southern California earthquakes with magnitude in the range 3.5–5.2 and provides stable results for 28 events. Apparent source time functions (ASTFs) of P and S phases are derived using deconvolution with three stacked empirical Green's functions (seGf). The use of seGf suppresses nongeneric source effects, improves the focal mechanism correspondence to the analyzed earthquakes, and typically allows inclusion of 5 to 15 more ASTFs compared with analysis using a single eGf. Most analyzed earthquakes in the Trifurcation area of the San Jacinto Fault have directivities toward the northwest, while events around Cajon Pass and San Gabriel Mountain tend to propagate toward the southeast. These results are generally consistent with predictions for dynamic rupture on bimaterial interfaces associated with the imaged velocity contrasts in the area. The second moment inversions also provide constraints on the upper and lower bounds of rupture areas in our data set. Stress drops and uncertainties are estimated for elliptical ruptures using the derived characteristic rupture length and width. The semiautomated second moment method with seGfs can be used for routine application to moderate earthquakes in locations with good station coverage.

California↗

Implications of regional gravity for state of stress in the earth's crust and upper mantle

Topography is maintained by stress differences within the earth. Depending on the distribution of the stress we classify the support as either local or regional compensation. In general, the stresses implied in a regional compensation scheme are an order of magnitude larger than those corresponding to local isostasy. Gravity anomalies, a measure of the earth's departure from hydrostatic equilibrium, can be used to distinguish between the two compensation mechanisms and thus to estimate the magnitude of deviatoric stress in the crust and upper mantle. Topography created at an oceanic ridge crest or in a major continental orogenic zone appears to be locally compensated. Such features were formed on weak crust incapable of maintaining stress differences much greater than the stress from the applied load. Oceanic volcanoes formed on an already cooled, thickened lithosphere are regionally supported with elastic stresses. Similarly, the broad topographic rise seaward of subduction zones is elastically supported as the lithosphere is bent near the plate margin. Although the implied stress is to some degree dependent on the rheological model assumed, the gravity anomalies and surface deformation produced by these features demonstrate that the upper 30–40 km of the oceanic lithosphere is capable of regionally supporting stress differences in the 100-MPa range. Given certain conditions of load emplacement, continental crust can also support loads regionally over 100-m.y. time scales, but the effects of erosion only allow an estimate of a lower bound on stress. Data from space probes indicate that the upper layers of other terrestrial planets also support topographic-induced stress differences in excess of 100 MPa.

Journal of Geophysical Research Solid Earth↗

Numerical model of frazil ice and suspended sediment concentrations and formation of sediment laden ice in the Kara Sea

A one-dimensional (vertical) numerical model of currents, mixing, frazil ice concentration, and suspended sediment concentration has been developed and applied in the shallow southeastern Kara Sea. The objective of the calculations is to determine whether conditions suitable for turbid ice formation can occur during times of rapid cooling and wind- and wave-induced sediment resuspension. Although the model uses a simplistic approach to ice particles and neglects ice-sediment interactions, the results for low-stratification, shallow (∼20-m) freeze-up conditions indicate that the coconcentrations of frazil ice and suspended sediment in the water column are similar to observed concentrations of sediment in turbid ice. This suggests that wave-induced sediment resuspension is a viable mechanism for turbid ice formation, and enrichment mechanisms proposed to explain the high concentrations of sediment in turbid ice relative to sediment concentrations in underlying water may not be necessary in energetic conditions. However, salinity stratification found near the Ob' and Yenisey Rivers damps mixing between ice-laden surface water and sediment-laden bottom water and probably limits incorporation of resuspended sediment into turbid ice until prolonged or repeated wind events mix away the stratification. Sensitivity analyses indicate that shallow (≤20 m), unstratified waters with fine bottom sediment (settling speeds of ∼1 mm s −1 or less) and long open water fetches (>25 km) are ideal conditions for resuspension.

Kara Sea↗

Analysis of atmospheric paniculate matter; application of optical and selected geochemical techniques

An increase in particulate matter in the atmosphere has been shown to be linked to increased mortality but this relationship is poorly understood. Light microscopy, electron microscopy , electron microprobe , and micro-FTIR techniques have been applied to study atmospheric particulates in Brisbane, Australia as a part of a study on asthma. The particulate matter samples were collected daily from April to August 1992, and the sampling covered the autumn period which is typically a time of high asthma incidence in Brisbane. Volumetrically, most atmospheric particulate matter is less than 2 μ m in size. The microscopic analysis reveals that this material is composed mainly of combusted and incompletely burned hydrocarbons from motor vehicle exhaust emissions, quiescent spores of Mucorales, soil bacteria, and inorganic matter in the form of quartz and other silicates. Elemental and functional group analyses confirm microscope identification, documenting carbon-rich, aromatic exhaust material, more aliphatic pollen and spore material and inorganic matter. Fungal spores dominate bioaerosol and are very abundant from the end of April through May to mid-June. The cytoplasmic content of pollens or fungal spores is commonly regarded as allergenic. Particulates from the exhaust emissions and crustal material in a sub-micrometer size range may act as carriers or dispersive mechanisms for cytoplasmic material from fungal spores and pollens, perhaps causing periods of the highest exhaust emission to be the most allergenic.

International Journal of Coal Geology↗

An accurate method for measuring triploidy of larval fish spawns

A standard flow cytometric protocol was developed for estimating triploid induction in batches of larval fish. Polyploid induction treatments are not guaranteed to be 100% efficient, thus the ability to quantify the proportion of triploid larvae generated by a particular treatment helps managers to stock high-percentage spawns and researchers to select treatments for efficient triploid induction. At 3 d posthatch, individual Grass Carp Ctenopharyngodon idella were mechanically dissociated into single-cell suspensions; nuclear DNA was stained with propidium iodide then analyzed by flow cytometry. Following ploidy identification of individuals, aliquots of diploid and triploid cell suspensions were mixed to generate 15 levels (0–100%) of known triploidy ( n = 10). Using either 20 or 50 larvae per level, the observed triploid percentages were lower than the known, actual values. Using nonlinear regression analyses, quadratic equations solved for triploid proportions in mixed samples and corresponding estimation reference plots allowed for predicting triploidy. Thus, an accurate prediction of the proportion of triploids in a spawn can be made by following a standard larval processing and analysis protocol with either 20 or 50 larvae from a single spawn, coupled with applying the quadratic equations or reference plots to observed flow cytometry results. Due to the universality of triploid DNA content being 1.5 times the diploid level and because triploid fish consist of fewer cells than diploids, this method should be applicable to other produced triploid fish species, and it may be adapted for use with bivalves or other species where batch analysis is appropriate.

North American Journal of Aquaculture↗

Estimating recharge at Yucca Mountain, Nevada, USA: Comparison of methods

Obtaining values of net infiltration, groundwater travel time, and recharge is necessary at the Yucca Mountain site, Nevada, USA, in order to evaluate the expected performance of a potential repository as a containment system for high-level radioactive waste. However, the geologic complexities of this site, its low precipitation and net infiltration, with numerous mechanisms operating simultaneously to move water through the system, provide many challenges for the estimation of the spatial distribution of recharge. A variety of methods appropriate for arid environments has been applied, including water-balance techniques, calculations using Darcy's law in the unsaturated zone, a soil-physics method applied to neutron-hole water-content data, inverse modeling of thermal profiles in boreholes extending through the thick unsaturated zone, chloride mass balance, atmospheric radionuclides, and empirical approaches. These methods indicate that near-surface infiltration rates at Yucca Mountain are highly variable in time and space, with local (point) values ranging from zero to several hundred millimeters per year. Spatially distributed net-infiltration values average 5 mm/year, with the highest values approaching 20 mm/year near Yucca Crest. Site-scale recharge estimates range from less than 1 to about 12 mm/year. These results have been incorporated into a site-scale model that has been calibrated using these data sets that reflect infiltration processes acting on highly variable temporal and spatial scales. The modeling study predicts highly non-uniform recharge at the water table, distributed significantly differently from the non-uniform infiltration pattern at the surface.

Hydrogeology Journal↗

Mechanisms of wave‐driven water level variability on reef‐fringed coastlines

Wave‐driven water level variability (and runup at the shoreline) is a significant cause of coastal flooding induced by storms. Wave runup is challenging to predict, particularly along tropical coral reef‐fringed coastlines due to the steep bathymetric profiles and large bottom roughness generated by reef organisms, which can violate assumptions in conventional models applied to open sandy coastlines. To investigate the mechanisms of wave‐driven water level variability on a reef‐fringed coastline, we performed a set of laboratory flume experiments on an along‐shore uniform bathymetric profile with and without bottom roughness. Wave setup and waves at frequencies lower than the incident sea‐swell forcing (infragravity waves) were found to be the dominant components of runup. These infragravity waves were positively correlated with offshore wave groups, signifying they were generated in the surf zone by the oscillation of the breakpoint. On the reef flat and at the shoreline, the low‐frequency waves formed a standing wave pattern with energy concentrated at the natural frequencies of the reef flat, indicating resonant amplification. Roughness elements used in the flume to mimic large reef bottom roughness reduced low frequency motions on the reef flat and reduced wave run up by 30% on average, compared to the runs over a smooth bed. These results provide insight into sea‐swell and infragravity wave transformation and wave setup dynamics on steep‐sloped coastlines, and the effect that future losses of reef bottom roughness may have on coastal flooding along reef‐fringed coasts.

Journal of Geophysical Research C: Oceans↗