USGS Science⌕ Search

SEARCH · USGS Science

Results for “Journal of Geophysical Research - Oceans”

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 1,531 records · Page 85Linked to original sources
USGS2012

Magmatic inclusions in silicic and intermediate volcanic rocks

Fine‐grained ellipsoidal inclusions from a few millimeters to over l m in size are present in many intermediate to silicic lava flows and domes. Only recently has it become widely accepted that such inclusions are chilled blobs of magma. Their magmatic origin is manifested by vesicularity and high groundmass porosity, by ellipsoidal shapes, by mingling at contacts with the host, and by textural evidence for groundmass crystallization in an undercooled state. Groundmass textures reflect the degree of undercooling, which is a function of thermal contrast (Δ T ) between inclusions and host magmas before mingling. Compositions of inclusions are invariably more mafic than their hosts and indicate that the magmatic systems were compositionally zoned. Many inclusions were formed from hybrid magmas. Four examples are given in order to document effects of decreasing Δ C (compositional contrast) and Δ T , as hosts for mafic andesitic to andesitic inclusions vary from high‐silica rhyolite (HSR) to silicic andesite. In the HSR example from the Coso volcanic field, California, inclusions track the evolution of zonation in a long‐lived magmatic system. At the other extreme, andesitic inclusions in silicic andesite from Crater Lake, Oregon, demonstrate zonation of a system of low Δ T and Δ C and show effects of in situ differentiation by gas filter‐pressing of residual liquid. Inclusions form when a relatively small proportion of mafic magma comes into contact with a larger proportion of silicic magma. Several mechanisms by which the magmas mix have been proposed: buoyant rise of vesiculated mafic magma or convective stirring at the interface between mafic and overlying silicic magma, forcible injection of mafic magma, diking of fractured silicic magma, and mixing during ascent in a conduit. No single mechanism explains all field observations. Many inclusions in plutons are analogous to magmatic inclusions in volcanic rocks but are likely to have suffered greater chemical and textural change during comparatively slow postentrapment cooling. Magmatic inclusions are evidence of open system behavior of the differentiated parts of magma bodies, and the inclusions themselves may behave as open systems where Δ T is low.

Journal of Geophysical Research B: Solid Earth↗
USGS2012

Paleomagnetism of Middle Tertiary volcanic rocks from the Western Cascade Series, northern California

The Western Cascade Series (WCS) is a 3.5‐km‐thick, crudely homoclinal (east dipping) calcalkaline volcanic sequence of mid‐Oligocene to early Miocene age that crops out near the southern tip of the Cascade Range in northern California. The mean direction of remanent magnetization in the WCS is D, 4.9°; I, 57.6° (N, 53; k, 14.4; α 95 , 5.3°). When compared to a reference direction for the North American craton, the WCS direction indicates that the southern Cascade Range has rotated 14.0° ± 9.0° since the WCS accumulated. A difference in mean direction between the lower and upper halves of the WCS suggests that much of this rotation occurred during the late Oligocene. Six other paleomagnetic studies of rock units of an age roughly comparable to the WCS also are available from western Oregon, northwestern California, and southwestern Washington. All show statistically significant clockwise rotation (inclinations of all but one are concordant). Comparison of directions indicates that the Pacific Northwest did not rotate as a rigid body, but neither did it behave as a collection of small, independently rotating domains. The amount of rotation found throughout the area increases to the west or northwest, suggesting a driving force for rotation operating at the continental margin.

California↗
USGS2012

Paleomagnetic evidence for the timing of collapse and resurgence of the Lake City Caldera, San Juan Mountains, Colorado

Rocks of the 23.1‐m.y.‐old Lake City caldera consist of the compositionally zoned Sunshine Peak Tuff, postcollapse intracaldera lava flows, and resurgent quartz syenite intrusions. Declinations of reversely magnetized ( I = −45° to −75°) Sunshine Peak Tuff change from easterly ( D = 93°–130°) throughout most of the tuff to southerly ( D = 195°–207°) within the late eruptive phases. The postcollapse lava flows are also reversely magnetized, and their emplacement appears to have spanned resurgence of the caldera. In contrast, the uppermost levels of the resurgent intrusion record two distinct components of magnetization: a reversed component isolated by thermal demagnetization typically at temperatures greater than 400°C and a normal component defined by thermal demagnetization at lower temperatures. A baked‐contact test implies that emplacement and initial cooling of the intrusion occurred during a time of reversed polarity and that continued cooling at lower blocking temperatures took place during normal polarity. Assuming that the reversed directions of the extrusive and intrusive rocks record a single period of reversed polarity as suggested by the paleomagnetic and radiometric age data, the time for caldera development (from eruption of the ash flow tuffs to emplacement of the resurgent intrusion) was less than about 300,000 years on the basis of best estimates of durations of reversed polarities in the early Miocene. The estimated 300,000‐year time span for the development of the Lake City caldera is very close in duration to the radiometrically determined time spans of much younger (Pleistocene) calderas having similar dimensions and similar volumes of erupted material.

Colorado↗
USGS2012

Short-period strain (0.1–105 s): Near-source strain field for an earthquake (ML 3.2) near San Juan Bautista, California

Measurements of dilational earth strain in the frequency band 25–10 −5 Hz have been made on a deep borehole strainmeter installed near the San Andreas fault. These data are used to determine seismic radiation fields during nuclear explosions, teleseisms, local earthquakes, and ground noise during seismically quiet times. Strains of less than 10 −10 on these instruments can be clearly resolved at short periods (< 10 s) and are recorded with wide dynamic range digital recorders. This permits measurement of the static and dynamic strain variations in the near field of local earthquakes. Noise spectra for earth strain referenced to 1 (strain) 2 /Hz show that strain resolution decreases at about 10 dB per decade of frequency from −150 dB at 10 −4 Hz to −223 dB at 10 Hz. Exact expressions are derived to relate the volumetric strain and displacement field for a homogeneous P wave in a general viscoelastic solid as observed on colocated dilatometers and seismometers. A rare near-field recording of strain and seismic velocity was obtained on May 26, 1984, from an earthquake ( M L 3.2) at a hypocentral distance of 3.2 km near the San Andreas fault at San Juan Bautista, California. While the data indicate no precursory strain release at the 5 × 10 −11 strain level, a coseismic strain release of 1.86 nanostrain was observed. This change in strain is consistent with that calculated from a simple dislocation model of the event. Ground displacement spectra, determined from the downhole strain data and instrument-corrected surface seismic data, suggest that source parameters estimated from surface recordings may be contaminated by amplification effects in near-surface low-velocity materials.

California↗
USGS2012

Paleomagnetism of the Tertiary Clarno Formation of central Oregon and its significance for the tectonic history of the Pacific Northwest

The Clarno Formation, a mostly Eocene and partly early Oligocene sequence of andesitic lavas and volcaniclastic rocks, is the oldest Tertiary formation exposed in north central Oregon. Remanent magnetization directions at 46 sites in the lavas provide a paleomagnetic pole at 84°N, 278°E with a 95% confidence cone of 7°. Comparison of this pole with the North American reference pole for Eocene time indicates that the Clarno Formation has rotated 16° clockwise with an uncertainty of 10° but has undergone no significant latitudinal displacement. Further comparison with paleomagnetic data from the Miocene Columbia River Basalt Group and Steens Basalt shows that large‐scale block rotations in the eastern Columbia Embayment ceased sometime between 38 m.y. and 15 m.y. A paleogeographic reconstruction for 38 m.y. is offered which is consistent with the observed rotation of the Clarno and with the other paleomagnetic data from the Pacific Northwest, which are briefly reviewed. This reconstruction shows that it is possible to account for virtually all of the paleomagnetically indicated rotations in pre‐Miocene Tertiary rocks of the Pacific Northwest by an extensional tectonic model. As a consequence, only part of the rotations in the pre‐Tertiary rocks of the eastern Columbia Embayment need to be the result of accretionary tectonics. The initiation of magmatism in the eastern Columbia Embayment that is represented by the Clarno Formation itself may have been contemporaneous with the beginning of crustal extension in the Pacific Northwest.

Journal of Geophysical Research B: Solid Earth↗
USGS2012

Usefulness of weak bands in midinfrared remote sensing of particulate planetary surfaces

Midinfrared (2.5–25 μm) reflectance spectra of minerals are often used to predict emittance qualitatively. These spectra display weak overtone and combination tone bands, which may be as diagnostic of composition as the strong fundamental molecular vibration bands usually considered for remote sensing applications, but which have been widely ignored. However, unlike the strong bands, the contrast of weak bands relative to the continuum in the midinfrared usually does not decrease with decreasing particle size, but typically increases. To illustrate this behavior, transmittance and reflectance spectra of calcite and quartz are presented for the wavelength range from 4000 to 400 cm −1 (2.5–25 μm). It is the purpose of this paper to point out that these weak bands are potentially useful for compositional remote sensing of particulate planetary surfaces; this will require completion of supporting laboratory studies to document the occurrence of weak bands.

Journal of Geophysical Research Solid Earth↗
USGS2012

Age of tilted reefs, Hawaii

Submerged carbonate reefs are preserved as a series of submarine terraces between Molokai and Hawaii along a 200-km span of the southeastern Hawaiian Ridge. Limestones from two of the terraces have been sampled from submersibles and dated radiometrically at 13 and 120 ka. Recognition that the terraces are tilted permits assignment of about a dozen terraces from 150 to 1300 m depth to eight general reef platforms. These reefs were drowned by the combined effects of island subsidence and sea level rise at the end of successive glacial stages from 13 to 647 ka. The platforms are tilted 5 m/km southeast toward the locus of volcanic loading centered on the island of Hawaii. The end of the active period of volcanic shield building is represented on each nonactive volcano by a submerged change in slope. Dating of this slope change by reference to the dated reefs provides independent evidence that shield building ended at about the same time that the volcano changed from eruption of tholeiitic to alkalic basalt. The variation in depth of the dated tilted reefs along the ridge indicates that the end of subsidence follows the end of shield building by about one half million years. The present northwest margin of major subsidence probably crosses southeast Maui. During the last half million years the rate of southeast migration of the beginning of both subsidence and of tholeiitic volcanism does not appear to have increased. However, the rate of southeast migration of the end of both subsidence and of tholeiitic volcanism has apparently increased. This shortening of the life history of the younger volcanoes implies an increase in eruption rates.

Journal of Geophysical Research Solid Earth↗
USGS2012

Petrology and geochemistry of alkali gabbronorites from Lunar Breccia 67975

Clasts of an unusual type of lunar highlands igneous rock, alkali gabbronorite, have been found in Apollo 16 breccia 67975. The alkali gabbronorites form two distinct subgroups, magnesian and ferroan. Modes and bulk compositions are highly varied. The magnesian alkali gabbronorites are composed of bytownitic plagioclase (OrAn), hypersthene (WoEn), augite (WoEn), a silica mineral, and trace Ba‐rich K‐feldspar. The ferroan alkali gabbronorites are composed of ternary plagioclase (OrAn), pigeonite (WoEn), augite (WoEn), Ba‐rich K‐feldspar, and a silica mineral. Trace minerals in both subgroups are apatite, REE‐rich whitlockite, and zircon. The magnesian and ferroan alkali gabbronorites appear to have formed by progressive differentiation of the same, or closely related, parent magmas; the compositional data indicate that these magmas were REE‐rich. The ternary plagioclase is probably a high‐temperature metastable phase formed during crystallization. In composition and mineralogy, the 67975 alkali gabbronorites show many similarities to Apollo 12 and 14 alkali norites, alkali gabbronorites, and alkali anorthosites, and all these rocks together constitute a distinctive alkali suite. In addition, the alkali gabbronorites show some similarities to KREEP basalts, Mg‐norites, and some felsites. These data suggest genetic links between some or all of these types of pristine rocks. Two types of relationships are possible. The first is that alkali‐suite rocks crystallized in plutons of KREEP basalt magma, and KREEP basalts are their extrusive equivalents. The second is that the alkali‐suite rocks and some felsites all crystallized in plutons of Mg‐norite parent magmas, and KREEP basalt magmas formed by remelting of these plutons. Additional studies are needed to resolve which of these hypotheses is correct.

Journal of Geophysical Research B: Solid Earth↗
USGS2012

Origin of lead from green glass of Apollo 15426: A search for primitive lunar lead.

A major obstacle in lunar U-Pb chronology is the elusiveness of the primary Pb isotopic composition and U/Pb ratio and therefore the Pb evolution for the early history of the moon. In an attempt to seek the primitive lunar Pb isotopic composition, green glass from lunar clod 15426,49 was studied for U-Th-Pb systematics because it is extremely Mg-rich and known to be the most primitive among sampled lunar volcanic rocks. Green glass and yellow impact glass were handpicked, washed with water, then washed in an acid mixture (0.3 N HBr + 0.2 N HF), and the washes and residues were analyzed. In contrast to the Pb from Apollo 17 orange glass, Pb from the surface of the green glass spherules is totally different from that of the spherule interior, but is similar to lunar surface lead. Data for the interior Pb yield a 207 Pb/ 206 Pb age of 3.41 Ga. Because of the low Pb concentration and high U/Pb ratio observed for the interior, the initial Pb was poorly defined. Nevertheless, the data indicate that lunar Pb evolved in an environment with 238 U/ 204 Pb = 19–55, which is considerably lower than those for mare basalts (around 300) but higher than values for the earth (6–8).

Journal of Geophysical Research Solid Earth↗
USGS2012

Remagnetization of the Coast Range Ophiolite and Lower Part of the Great Valley Sequence in Northern California and Southwest Oregon

Overprinted magnetizations have been found at four localities in the Middle Jurassic Coast Range ophiolite and the overlying Upper Jurassic and Lower Cretaceous Great Valley sequence in northern California and at one locality in the partially correlative Lower Cretaceous Days Creek Formation in southwest Oregon. At Del Puerto Canyon, on the east side of the Diablo Range, a pilot study of the upper Jurassic Lotta Creek Formation gives in situ paleomagnetic directions grouped around the present geomagnetic field, suggesting magnetic overprinting long after deposition and folding. We suggest that the loss of original magnetization could be the result of long burial at about 7 km depth, followed by Late Tertiary uplift; the possibility of chemical remagnetization, however, can not be excluded. Paleomagnetic data from sandstones from the Lower Cretaceous Great Valley sequence in the Wilbur Springs area suggest Cenozoic remagnetization that could be related to Pliocene and Pleistocene volcanic and hydrothermal activity in the area. The results from calcareous concretions in the Wilbur Springs area, from 21 sites from basalts and overlying sedimentary rocks of the Great Valley sequence at Stonyford, and from several sites in various rock types in the Coast Range ophiolite and Great Valley sequence near Paskenta suggest unblocking of magnetization during the long period of burial at depth of 7 km or more and remagnetization during tectonic uplift in the Late Cretaceous or Tertiary. Chemical remagnetization, however, can not be excluded. The paleomagnetic data from six sites in Lower Cretaceous sedimentary rocks of the Days Creek Formation in southwest Oregon, which overlies a dismembered ophiolite, fail the fold test. In situ paleomagnetic directions group near the expected Tertiary field directions. These sedimentary rocks may have been remagnetized during a major Eocene tectonic event known to have occurred in this region. The paleomagnetic results confirm that the Coast Range ophiolite and the seemingly little-deformed Great Valley sequence in northern California, as well as correlative rocks in southwest Oregon have had a complex tectonic history. Once the nature and timing of these events are better understood, the timing of remagnetization may be further constrained, which, in turn, could give further insight into the nature of the tectonic events

California, Oregon↗
USGS2012

Evidence for a shallow pluton beneath the Goat Rocks Wilderness, Washington, from gravity and magnetic data

A gravity and aeromagnetic study was conducted in and adjacent to the Goat Rocks Wilderness area, Washington. This work was done in conjunction with geologic and geochemical studies which were critically important to our interpretation. The Bouguer gravity anomaly map displays a large central positive gravity anomaly and three superimposed smaller positive anomalies. The main anomaly, about 26 by 11 km with an amplitude of about 16 mGal seems to delineate a dense large pluton or sheeted dike complex intruding less dense, older volcanic and sedimentary rocks lying directly below the volcano. The pluton or sheeted dike complex may have fed the Goat Rocks volcano throughout its 1–2 m.y. history. Three small feeders, mapped on the surface as andesite intrusions, may be the cause of three smaller gravity anomalies superimposed on the main anomaly. The magnetic anomaly maps also exhibit a large central anomaly whose source is the Goat Rocks pluton. The magnetic data suggest that the three shallow intrusive cupolas have undergone a significant amount of hydrothermal alteration. In addition, to some extent the magnetic data reflect the terrain of the area. Topographic highs produce magnetic highs, and topographic lows produce magnetic lows. The deviations from this pattern delineate subsurface geology. Magnetic lows may be of particular interest in locating zones of mineralization where magnetic minerals have been altered and have become significantly less magnetic.

Washington↗
USGS2012

Character and distribution of borehole breakouts and their relationship to in situ stresses in deep Columbia River basalts

The character and distribution of borehole breakouts in deeply buried basalts at the Hanford Site in south central Washington State are examined in light of stress indicator data and hydraulic-fracturing stress data by means of acoustic televiewer and acoustic waveform logging systems. A series of boreholes penetrating the Grande Ronde Basalt of the Columbia River Basalt Group were logged to examine the extent of breakouts at depths near 1000 m. Breakouts occur discontinuously throughout the interiors of most flows. In some boreholes the distribution of borehole wall breakouts closely correlates with the incidence of core disking. Differences in the distribution of breakouts and disking are attributed to differences in failure mechanisms. A thin interval of breakout-free basalt occurs near the upper and lower limits of flow interiors, with many intervals of breakouts terminating at the intersection of oblique fractures with the borehole. Hydraulic-fracturing stress measurement results obtained from four deep boreholes indicate anisotropic horizontal principal stresses, with maximum principal stress along an approximate north trending axis, consistent with the east trending orientation of breakouts. Acoustic waveform logs indicate that there is no measurable difference between the seismic properties of breakout-free flow tops and flow interiors. The highly coherent waveforms obtained in almost all flow interiors indicate that damage to the borehole wall associated with breakout formation remains confined to the thin annulus of stress concentration.

Journal of Geophysical Research Solid Earth↗
USGS2012

Detection of a locked zone at depth on the Parkfield, California, segment of the San Andreas fault

The Parkfield, California, segment of the San Andreas fault is transitional in character between the creeping segment of the fault to the northwest and the locked Carrizo Plain segment to the southeast. The rate of shallow fault slip decreases from 25–30 mm/yr northwest of the epicenter of the 1966 Parkfield earthquake to zero at the southeastern end of the 1966 rupture zone. Data from a network of trilateration lines spanning the San Andreas fault near Parkfield and extending to the Pacific coast near San Luis Obispo shed light on the rate of fault slip at depth since the 1966 earthquake. In this study, average rates of line length change and shallow fault slip were inverted to determine the slip rate at depth on the Parkfield fault segment. The fault is taken to be a vertical surface with unknown distribution of strike-slip displacement in an elastic half-space. A striking result of the inversions is that all solutions providing acceptable fits to the data exhibit a locked zone essentially coincident with the rupture surface of the 1966 Parkfield earthquake. The data require that the locked zone extend nearly as far north as the 1966 epicenter; however, the vertical extent of the locked zone is not well resolved. Over much of the Parkfield segment the fault is slipping faster at the earth's surface than it is at seismogenic depths. In order to fit the trilateration measurements it is necessary to include a component of contraction normal to the trend of the San Andreas. The inversion results suggest a spatially uniform normal strain of −0.06 μstrain/yr. The orientation of the contraction is compatible with geologic and seismic evidence of active folding and reverse faulting in the region. The magnitude of the contraction is consistent with convergence rates inferred from global plate motion models.

Journal of Geophysical Research Solid Earth↗
USGS2012

Stress rotation during the Coalinga aftershock sequence

This study considers spatial and temporal changes of the stress regime during the 1983 Coalinga aftershock sequence. In both cases the observed changes manifest themselves as rotations of the most compressive stress axis. Inversion of 165 M ≥3 aftershocks shows that in the southern aftershock zone the azimuth of the most compressional stress axis is oriented 24°E of north, while in the northern and central zones it is oriented 51°E of north. Statistical analysis by use of non-parametric bootstrap resampling shows that this rotation is significant at above the 95% confidence level. Inverting the focal mechanisms of 122 M ≥3 aftershocks in the central and northern aftershock zones reveals a rotation of the principal stress axes during the first 500 days after the main shock. The compressional axis rotates from an azimuth of 62°E of north to 47°E of north. Statistical analysis by use of nonparametric bootstrap resampling shows that this rotation is just significantly different from zero at the 95% confidence level, subject to the assumption that 80% of the fault planes can be correctly picked from the two nodal planes. As the number of fault planes picked correctly varies from 50% (random choices) to 100% the significance of the rotation varies from 87% to 98%. Further investigation fails to find the signature of the stress changes predicted by dislocation theory in the focal mechanisms of the aftershocks. Active folds in the area suggest that the pre-main shock stress regime had a compressional axis with a NE-SW orientation. Thus the rotation during the aftershock sequence is thought to be a rebound from a stress anomaly induced by the main shock. This stress anomaly may be causally related to the occurrence of the aftershocks.

Journal of Geophysical Research Solid Earth↗
USGS2012

Fracture characterization by means of attenuation and generation of tube waves in fractured crystalline rock at Mirror Lake, New Hampshire

Results are presented from experiments carried out in conjunction with the U. S. Geological Survey at the Hubbard Brook Experimental Forest near Mirror Lake, New Hampshire. The study focuses on our ability to obtain orientation and transmissivity estimates of naturally occurring fractures. The collected data set includes a four-offset hydrophone vertical seismic profile, full waveform acoustic logs at 5, 15, and 34 kHz, borehole televiewer, temperature, resistivity, and self-potential logs, and borehole-to-borehole pump test data. Borehole televiewer and other geophysical logs indicate that permeable fractures intersect the Mirror Lake boreholes at numerous depths, but less than half of these fractures appear to have significant permeability beyond the annulus of drilling disturbance on the basis of acoustic waveform log analysis. The vertical seismic profiling (VSP) data indicate a single major permeable fracture near a depth of 44 m, corresponding to one of the most permeable fractures identified in the acoustic waveform log analysis. VSP data also indicate a somewhat less permeable fracture at 220 m and possible fractures at depths of 103 and 135 m; all correspond to major permeable fractures in the acoustic waveform data set. Pump test data confirm the presence of a hydraulic connection between the Mirror Lake boreholes through a shallow dipping zone of permeability at 44 m in depth. Effective fracture apertures calculated from modeled transmissivities correspond to those estimated for the largest fractures indicated on acoustic waveform logs but are over an order of magnitude larger than effective apertures calculated from tube waves in the VSP data set. This discrepancy is attributed to the effect of fracture stiffness. A new model is presented to account for the mechanical strength of asperities in resisting fracture closure during the passage of seismic waves during the generation of VSPs.

New Hampshire↗
USGS2012

The Steens Mountain (Oregon) geomagnetic polarity transition: 3. Its regional significance

Study of the variations of direction and intensity of the geomagnetic field as recorded by the Miocene lava flows on Steens Mountain, southeastern Oregon, has resulted in a detailed description of total field behavior during a reversal in polarity. In addition to information about the polarity reversal itself, the detailed paleomagnetic record includes several thousand years of geomagnetic history preceding and following the polarity transition at 15.5 Ma. In order to test the feasibility of using this record as a means of correlation in this part of the western United States, comparisons are made of reconnaissance and previously published paleomagnetic records obtained from what has been thought to be the Steens Basalt or rocks of equivalent age. Despite the fact that many of these earlier studies were not done in detail and were not intended for correlation purposes, convincing similarities among some of the records are evident. The Steens Basalt paleomagnetic record does, indeed, have potential as a correlation tool during this time of widespread basaltic volcanism. Additionally, paleomagnetic data from flows that were sampled in detail yield a middle Miocene paleomagnetic pole at 88.3°N, 209.0° (α 95 = 6.3°) for the High Lava Plains of Oregon. This pole position is statistically indistinguishable from the earth's rotational axis and implies that no tectonic rotation of this region has occurred since these lava flows were erupted. Data from selected sites within the coeval part of the Columbia River Basalt Group yield a paleomagnetic pole at 88.7°N, 171.6°E (α 95 = 4.0°). The Columbia River Basalt Group pole is statistically indistinguishable from either the rotational axis or from the High Lava Plains pole. These findings indicate no post‐20 Ma differential rotation between south‐eastern Washington and south‐central Oregon, in contrast to previous interpretations.

Journal of Geophysical Research B: Solid Earth↗
USGS2012

Columbia Glacier, Alaska: Changes in velocity 1977-1986

The Columbia Glacier, a grounded, iceberg-calving tidewater glacier near Valdez, Alaska, began to retreat about 1977. Drastic retreat occurred in 1984, and by early 1986, retreat amounted to 2 km. The glacier has thinned more than 100 m since 1974 at a point 4 km behind the 1974 terminus position. Between 1977 and 1985 the lower glacier ice velocity increased from 3–8 m/d to 10–15 m/d. Ice velocity in the region 0.5 km above the terminus was highest near the time the glacier was most receded (late fall), and lowest near the time of maximum length (early summer), for years 1977–1982. Velocity in the region 52–57 from the head of the glacier was highest in mid-spring, and lowest in early fall from 1977 to 1985. Through the years 1983–1985, the dates of maximum and minimum velocities within 0.5 km of the receding terminus tended toward the dates of the 52–57 km maximum and minimums. This occurred because as the terminus receded, it was no longer strongly influenced by the reverse slope of the terminal moraine shoal. Velocities near the terminus fluctuated by 2–3 m/d during summer and fall, when liquid water input was variable, and were relatively constant during winter. Hourly variations in ice velocities are controlled by liquid water input to the glacier hydraulic system and tide stage. Velocity increases near periods of high surface water input and decreases during periods of high tide as a result of hydrostatic back pressure.

Journal of Geophysical Research Solid Earth↗
USGS2012

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↗