USGS Science⌕ Search

Geology topics

R. A. Walters

Publications and source records attributed to R. A. Walters.

At least 19 recordsLinked to original sources

Analysis of time series of glacier speed: Columbia Glacier, Alaska

During the summers of 1984 and 1985, laser measurements were made of the distance from a reference location to markers on the surface of the lower reach of Columbia Glacier, Alaska. The data have numerous gaps, mostly because of inclement weather. The laser measurements were corrected for variations in atmospheric temperature and pressure and for the trajectory of the marker. The marker speed was calculated by using a combination of cubic splines and digital filters to arrive at a data set with 1-hour time intervals. The speed varies from 7 to 15 m/d and has three noteworthy components: (1) a low-frequency perturbation in speed with a time scale of days related to increased precipitation, (2) semidiurnal and diurnal variations related to sea tides, and (3) diurnal variations related to glacier surface melt. The low frequency and tidal period variations are separated by low-pass filtering and subtracting the filtered data from the original data. The variations from melt runoff and tides have the same frequency range and are examined using harmonic analysis and multiple regression with the various forcing functions. The high-frequency portion of the ice speed signal is dominated by tidal influences and to a lesser extent by meltwater influences. The low-frequency portion of the signal is dominated by the effects of precipitation.

Journal of Geophysical Research Solid Earth↗

Effectiveness of spinning-wing decoys varies among dabbling duck species and locations

Spinning-wing decoys are strong attractants to ducks and inc rease kill rates over traditional decoying methods. However, it is unknown whether all duck species are attracted similarly to spinning-wing decoys and whether the effectiveness of these decoys changes with latitude. We examined the effectiveness of spinning-wing decoys for 9 species of dabbling ducks during 545 experimental hunts in California (1999-2000), Minnesota (2002), Manitoba (2001-2002), Nebraska (2000-2002), Missouri (2000-2001), and Arkansas (2001-2003). During each experimental hunt, we systematically alternated between 2 paired decoy treatments every 15-30 min (depending on study site): traditional decoys only and traditional decoys with a spinning-wing decoy. Overall, 70.2% (n=1,925) of dabbling ducks were harvested (shot and retrieved) when spinning-wing decoys were turned on, ranging from 63.6% (n=187) in Missouri to 76.4% (n=356) in Minnesota. Effectiveness of spinning-wing decoys increased with latitude of study sites. Proportions of ducks shot when spinning-wing decoys were turned on differed among species, from a low of 50.0% (n=8) for cinnamon teal (Anas cyanoptera) to a high of 79.0% (n=119) for American wigeon (A. americana). The probability of being shot when spinning-wing decoys were turned on increased with annual survival rates among species; for example, spinning-wing decoys were more effective for American wigeon and mallard (A. platyrhynchos) than they were for cinnamon teal and American green-winged teal (A. crecca). Effectiveness of spinning-wing decoys did not differ consistently by age or sex of harvested ducks. Our results indicate that the effectiveness of spinning-wing decoys differs among duck species and changes with latitude; thus, consideration of these effects may be warranted when setting harvest regulations and methods of take.

Journal of Wildlife Management↗

An unstructured grid, three-dimensional model based on the shallow water equations

A semi-implicit finite difference model based on the three-dimensional shallow water equations is modified to use unstructured grids. There are obvious advantages in using unstructured grids in problems with a complicated geometry. In this development, the concept of unstructured orthogonal grids is introduced and applied to this model. The governing differential equations are discretized by means of a semi-implicit algorithm that is robust, stable and very efficient. The resulting model is relatively simple, conserves mass, can fit complicated boundaries and yet is sufficiently flexible to permit local mesh refinements in areas of interest. Moreover, the simulation of the flooding and drying is included in a natural and straightforward manner. These features are illustrated by a test case for studies of convergence rates and by examples of flooding on a river plain and flow in a shallow estuary. Copyright ?? 2000 John Wiley & Sons, Ltd.

International Journal for Numerical Methods in Flu↗

A robust, finite element model for hydrostatic surface water flows

A finite element scheme is introduced for the 2-dimensional shallow water equations using semi-implicit methods in time. A semi-Lagrangian method is used to approximate the effects of advection. A wave equation is formed at the discrete level such that the equations decouple into an equation for surface elevation and a momentum equation for the horizontal velocity. The convergence rates and relative computational efficiency are examined with the use of three test cases representing various degrees of difficulty. A test with a polar-quadrant grid investigates the response to local grid-scale forcing and the presence of spurious modes, a channel test case establishes convergence rates, and a field-scale test case examines problems with highly irregular grids.A finite element scheme is introduced for the 2-dimensional shallow water equations using semi-implicit methods in time. A semi-Lagrangian method is used to approximate the effects of advection. A wave equation is formed at the discrete level such that the equations decouple into an equation for surface elevation and a momentum equation for the horizontal velocity. The convergence rates and relative computational efficiency are examined with the use of three test cases representing various degrees of difficulty. A test with a polar-quadrant grid investigates the response to local grid-scale forcing and the presence of spurious modes, a channel test case establishes convergence rates, and a field-scale test case examines problems with highly irregular grids.

Communications in Numerical Methods in Engineering↗

Parallel iterative solution for h and p approximations of the shallow water equations

A p finite element scheme and parallel iterative solver are introduced for a modified form of the shallow water equations. The governing equations are the three-dimensional shallow water equations. After a harmonic decomposition in time and rearrangement, the resulting equations are a complex Helmholz problem for surface elevation, and a complex momentum equation for the horizontal velocity. Both equations are nonlinear and the resulting system is solved using the Picard iteration combined with a preconditioned biconjugate gradient (PBCG) method for the linearized subproblems. A subdomain-based parallel preconditioner is developed which uses incomplete LU factorization with thresholding (ILUT) methods within subdomains, overlapping ILUT factorizations for subdomain boundaries and under-relaxed iteration for the resulting block system. The method builds on techniques successfully applied to linear elements by introducing ordering and condensation techniques to handle uniform p refinement. The combined methods show good performance for a range of p (element order), h (element size), and N (number of processors). Performance and scalability results are presented for a field scale problem where up to 512 processors are used.

Advances in Water Resources↗

The sandflat habitat: Scaling from experiments to conclusions

Ecological systems are characterised by spatial and temporal variations in the density of organisms and resources, and in the intensity of processes which affect them (e.g. Watt, 1947, Dayton, 1971, Addicot et al., 1987, Kolasa and Pickett, 1991 and Giller et al., 1994). This heterogeneity represents both a difficulty for field study design and statistical testing, and a challenge to describe the spatial structuring of populations, communities and ecosystems (Legendre, 1993). Patterns in ecological and environmental variables are fundamental to developing hypotheses and in designing subsequent field studies, because they determine the spatial and temporal scales of study. Spatial and/or temporal patterns also provide useful stepping stones for addressing issues of scale. Patterns apparent at one scale can collapse to noise when viewed from other scales, indicating that patterns, processes and our perceptions vary in a scale-dependent manner. Perhaps the most important but rarely discussed issue is the question of how to scale-up from small-scale surveys and experiments to conclusions relevant at larger spatial and temporal scales.

Conference Paper↗

Characterisation of physical environmental factors on an intertidal sandflat, Manukau Harbour, New Zealand

Physical environmental factors, including sediment characteristics, inundation time, tidal currents and wind waves, likely to influence the structure of the benthic community at meso-scales (1-100 m) were characterised for a sandflat off Wiroa Island (Manukau Harbour, New Zealand). In a 500 x 250 m study site, sediment characteristics and bed topography were mostly homogenous apart from patches of low-relief ridges and runnels. Field measurements and hydrodynamic modelling portray a complex picture of sediment or particulate transport on the intertidal flat, involving interactions between the larger scale tidal processes and the smaller scale wave dynamics (1-4 s; 1-15 m). Peak tidal currents in isolation are incapable of eroding bottom sediments, but in combination with near-bed orbital currents generated by only very small wind waves, sediment transport can be initiated. Work done on the bed integrated over an entire tidal cycle by prevailing wind waves is greatest on the elevated and flatter slopes of the study site, where waves shoal over a wider surf zone and water depths remain shallow e enough for wave-orbital currents to disturb the bed. The study also provided physical descriptors quantifying static and hydrodynamic (tidal and wave) factors which were used in companion studies on ecological spatial modelling of bivalve distributions and micro-scale sediment reworking and transport.

Conference Paper↗

Application of Conjugate Gradient methods to tidal simulation

A harmonic decomposition technique is applied to the shallow water equations to yield a complex, nonsymmetric, nonlinear, Helmholtz type problem for the sea surface and an accompanying complex, nonlinear diagonal problem for the velocities. The equation for the sea surface is linearized using successive approximation and then discretized with linear, triangular finite elements. The study focuses on applying iterative methods to solve the resulting complex linear systems. The comparative evaluation includes both standard iterative methods for the real subsystems and complex versions of the well known Bi-Conjugate Gradient and Bi-Conjugate Gradient Squared methods. Several Incomplete LU type preconditioners are discussed, and the effects of node ordering, rejection strategy, domain geometry and Coriolis parameter (affecting asymmetry) are investigated. Implementation details for the complex case are discussed. Performance studies are presented and comparisons made with a frontal solver.

Advances in Water Resources↗

A three-dimensional, finite element model for coastal and estuarine circulation

This paper describes the development and application of a three-dimensional model for coastal and estuarine circulation. The model uses a harmonic expansion in time and a finite element discretization in space. All nonlinear terms are retained, including quadratic bottom stress, advection and wave transport (continuity nonlinearity). The equations are solved as a global and a local problem, where the global problem is the solution of the wave equation formulation of the shallow water equations, and the local problem is the solution of the momentum equation for the vertical velocity profile. These equations are coupled to the advection-diffusion equation for salt so that density gradient forcing is included in the momentum equations. The model is applied to a study of Delaware Bay, U.S.A., where salinity intrusion is the primary focus.

Continental Shelf Research↗

A 3D, finite element model for baroclinic circulation on the Vancouver Island continental shelf

This paper describes the development and application of a 3-dimensional model of the barotropic and baroclinic circulation on the continental shelf west of Vancouver Island, Canada. A previous study with a 2D barotropic model and field data revealed that several tidal constituents have a significant baroclinic component (the K1 in particular). Thus we embarked on another study with a 3D model to study the baroclinic effects on the residual and several selected tidal constituents. The 3D model uses a harmonic expansion in time and a finite element discretization in space. All nonlinear terms are retained, including quadratic bottom stress, advection and wave transport (continuity nonlinearity). The equations are solved as a global and a local problem, where the global problem is the solution of the wave equation formulation of the shallow water equations, and the local problem is the solution of the momentum equation for the vertical velocity profile. These equations are coupled to the advection-diffusion equation for density so that density gradient forcing is included in the momentum equations. However, the study presented here describes diagnostic calculations for the baroclinic residual circulation only. The model is sufficiently efficient that it encourages sensitivity testing with a large number of model runs. In this sense, the model is akin to an extension of analytical solutions to the domain of irregular geometry and bottom topography where this parameter space can be explored in some detail. In particular, the consequences of the sigma coordinate system used by the model are explored. Test cases using an idealized representation of the continental shelf, shelf break and shelf slope, lead to an estimation of the velocity errors caused by interpolation errors inherent in the sigma coordinate system. On the basis of these estimates, the computational grid used in the 2D model is found to have inadequate resolution. Thus a new grid is generated with increased accuracy in the region of the shelf break. However, even with increased resolution, spurious baroclinic circulation seaward of the shelf break and in the vicinity of Juan de Fuca canyon remained a significant problem when the pressure gradient terms were evaluated using the σ coordinate system and using a realistic density profile. With the new grid, diagnostic calculations of the barotropic and baroclinic residual circulation are performed using forcing from the observed σ t (density) field and from the gradient of this field.

Journal of Marine Systems↗

Columbia Bay, Alaska: An 'upside down' estuary

Circulation and water properties within Columbia Bay, Alaska, are dominated by the effects of Columbia Glacier at the head of the Bay. The basin between the glacier terminus and the terminal moraine (sill depth of about 22 m) responds as an ‘upside down’ estuary with the subglacial discharge of freshwater entering at the bottom of the basin. The intense vertical mixing caused by the bouyant plume of freshwater creates a homogeneous water mass that exchanges with the far-field water through either a two- or a three-layer flow. In general, the glacier acts as a large heat sink and creates a water mass which is cooler than that in fjords without tidewater glaciers. The predicted retreat of Columbia Glacier would create a 40 km long fjord that has characteristics in common with other fjords in Prince William Sound.

Alaska↗

Time scales of circulation and mixing processes of San Francisco Bay waters

Conceptual models for tidal period and low-frequency variations in sea level, currents, and mixing processes in the northern and southern reaches of San Francisco Bay describe the contrasting characteristics and dissimilar processes and rates in these embayments: The northern reach is a partially mixed estuary whereas the southern reach (South Bay) is a tidally oscillating lagoon with density-driven exchanges with the northern reach. The mixed semidiurnal tides are mixtures of progressive and standing waves. The relatively simple oscillations in South Bay are nearly standing waves, with energy propagating down the channels and dispersing into the broad shoal areas. The tides of the northern reach have the general properties of a progressive wave but are altered at the constriction of the embayments and gradually change in an upstream direction to a mixture of progressive and standing waves. The spring and neap variations of the tides are pronounced and cause fortnightly varying tidal currents that affect mixing and salinity stratification in the water column. Wind stress on the water surface, freshwater inflow, and tidal currents interacting with the complex bay configuration are the major local forcing mechanisms creating low-frequency variations in sea level and currents. These local forcing mechanisms drive the residual flows which, with tidal diffusion, control the water-replacement rates in the estuary. In the northern reach, the longitudinal density gradient drives an estuarine circulation in the channels, and the spatial variation in tidal amplitude creates a tidally-driven residual circulation. In contrast, South Bay exhibits a balance between wind-driven circulation and tidally-driven residual circulation for most of the year. During winter, however, there can be sufficient density variations to drive multilayer (2 to 3) flows in the channel of South Bay. Mixing models (that include both diffusive and dispersive processes) are based on time scales associated with salt variations at the boundaries and those associated with the local forcing mechanisms, while the spatial scales of variations are dependent upon the configuration of the embayments. In the northern reach, where the estuarine circulation is strong, the salt flux is carried by the mean advection of the mean salt field. Where large salinity gradients are present, the tidal correlation part of the salt flux is of the same order as the advective part. Our knowledge of mixing and exchange rates in South Bay is poor. As this embayment is nearly isohaline, the salt flux is dominated entirely by the mean advection of the mean salt field. During and after peaks in river discharge, water mixing becomes more dynamic, with a strong density-driven current creating a net exchange of both water mass and salt. These exchanges are stronger during neap tides. Residence times of the water masses vary seasonally and differ between reaches. In the northern reach, residence times are on the order of days for high winter river discharge and of months for summer periods. The residence times for South Bay are fairly long (on the order of several months) during summer, and typically shorter (less than a month) during winter when density-driven exchanges occur.

Hydrobiologia↗

Subtidal sea level and current variations in the northern reach of San Francisco Bay

Analyses of sea level and current-meter data using digital filters and a variety of statistical methods show a variety of phenomena related to non-local coastal forcing and local tidal forcing in the northern reach of San Francisco Bay, a partially mixed estuary. Low-frequency variations in sea level are dominated by non-local variations in coastal sea level and also show a smaller influence from tidally induced fortnightly sea level variations. Low-frequency currents demonstrate a gravitational circulation which is modified by changes in tidal-current speed over the spring-neap tidal cycle. Transients in gravitational circulation induce internal oscillations with periods of two to four days.

California↗

Transport and concentration controls for chloride, strontium, potassium and lead in Uvas Creek, a small cobble-bed stream in Santa Clara County, California, U.S.A.: 2. Mathematical modeling

Three models describing solute transport of conservative ion species and another describing transport of species which adsorb linearly and reversibly on bed sediments are developed and tested. The conservative models are based on three different conceptual models of the transient storage of solute in the bed. One model assumes the bed to be a well-mixed zone with flux of solute into the bed proportional to the difference between stream concentration and bed concentration. The second model assumes solute in the bed is transported by a vertical diffusion process described by Fick's law. The third model assumes that convection occurs in a selected portion of the bed while the mechanism of the first model functions everywhere. The model for adsorbing species assumes that the bed consists of particles of uniform size with the rate of uptake controlled by an intraparticle diffusion process. All models are tested using data collected before, during and after a 24-hr. pulse injection of chloride, strontium, potassium and lead ions into Uvas Creek near Morgan Hill, California, U.S.A. All three conservative models accurately predict chloride ion concentrations in the stream. The model employing the diffusion mechanism for bed transport predicts better than the others. The adsorption model predicts both strontium and potassium ion concentrations well during the injection of the pulse but somewhat overestimates the observed concentrations after the injection ceases. The overestimation may be due to the convection of solute deep into the bed where it is retained longer than the 3-week post-injection observation period. The model, when calibrated for strontium, predicts potassium equally well when the adsorption equilibrium constant for strontium is replaced by that for potassium.

California↗

Low-flow transport models for conservative and sorbed solutes — Uvas Creek, near Morgan Hill, California

Models describing low-flow transport of conservative (nonreactive) and reactive solutes, which adsorb on the streambed, are developed and tested. Temporary storage within the bed plays an important role in solute movement. Three different models of bed-storage processes are developed for conservative solutes. One model assumes the bed is a well-mixed, nondiffusing, nonreacting zone. Solute flux into the bed is then proportional to the difference between stream and bed-solute concentrations. A second model assumes that solute is transported within the bed by a vertical diffusion process. The bed-solute concentration, which matches the stream concentration at the interface, varies with depth in the bed according to Fick 's law. A third model assumes convection in the downstream direction occurs in certain parts of the bed, while the mechanism of the first model functions elsewhere. Storage of absorbing species is assumed to occur by equilibrium adsorption within streambed particles. Uptake rate is described by an intraparticle diffusion process. Model equations were solved using finite element numerical methods. Models were calibrated using data from a 24-hour injection of conservative chloride and adsorptive Sr ions at Uvas Creek near Morgan Hill, California. All models predict well except for some overestimation by the adsorption model during dieaway.

California↗