USGS Science⌕ Search

SEARCH · USGS Science

Results for “Interfaces”

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 991 records · Page 55Linked to original sources
USGS2012

Tide gage measurements of uplift along the south coast of Alaska

Annual mean sea levels along the south coast of Alaska are used to measure uplift along the Alaska‐Aleutian subduction zone. Oceanographic effects are removed from the observed annual mean sea levels by subtracting a correction that is proportional to the sea level fluctuations observed in southeast Alaska. That correction is effective in reducing fluctuations in the observed, annual mean sea level as far west as the tip of Alaska peninsula. Additional corrections to remove the eustatic rise in sea level and the apparent fall in sea level due to postglacial isostatic rebound of the land are introduced. This corrected sea level record should provide a measure of tectonic subsidence. In the area affected by the 1964 Alaska earthquake, postseismic uplift occurs where coseismic subsidence was observed, and postseismic subsidence occurs where coseismic uplift was observed. The immediate postseismic response is damped out within the first decade, and the subsequent uplift rates appear to be steady over the 1974–1989 interval. However, some of those rates seem to be too high to be sustained over the ∼1000 year earthquake recurrence interval appropriate to this area if the interseismic deformation is only to recover the coseismic displacement. Thus a long‐term ( ∼100 years) relaxation in uplift rates is postulated. The immediate (time constant ∼5 years) postseismic relaxation is attributed to postseismic slip on the plate interface directly downdip from the coseismic rupture. The long‐term (time constant ∼100 years) relaxation is attributed to flow in the asthenosphere.

Alaska↗
USGS2012

Basin and range crustal and upper mantle structure, northwest to central Nevada

We present an interpretation of the crustal and uppermost mantle structure of the Basin and Range of northwestern Nevada based on seismic refraction/wide‐angle reflection, near‐vertical reflection, and gravity data. In comparison to most previous estimates, we find that the crust is somewhat thicker (32–36 km versus 22–30 km), and the uppermost mantle velocity is somewhat higher (8.0 km/s versus 7.3–7.9 km/s). Along our transects, the crust is thinnest (32 km) in the Carson Sink‐Buena Vista Valley region and increases by 2–4 km to the west and east, respectively. There is considerable complexity throughout the crust where velocities range from of 2.5 km/s at the surface to 7.4 km/s in the lowermost crust. Variations in velocity and structure of the upper crustal layers reveal apparent basement velocity depressions (areas of lower velocities extending up to 10 km in depth) that underlie some surface ranges as well as the basins. The middle crust rises from about 20 km beneath central Nevada to within 12 km of the surface beneath the area of thinnest crust and is characterized by a modest (∼0.1 km/s) change in velocity and low‐velocity gradients. These midcrustal layers mark the onset of high crustal reflectivity and the apparent limiting depth to which Basin and Range faults can be traced in near‐vertical reflection profiles, suggesting that these midcrustal layers represent the transition between the brittle and ductile zones of the crust. The lower crust is more structurally complex, with layers thickening and thinning in a systematic manner with the upper crustal layers; generally, where there are velocity depressions in the upper crust, the lower crust is thickest and shallowest. The geometry of these lower crustal layers (derived from refraction modeling) coincides with changes in the crustal reflectivity, determined from the Consortium of Continental Reflection Profiling reflection data. The lower crustal layer is unusually high in velocity (7.4 km/s) and is likely the layer identified as mantle in some previous studies. We do not identify the 7.4 km/s layer as mantle because (1) there is an underlying layer with a velocity (8.0 km/s) that is more consistent with the worldwide average velocity for the upper mantle, and (2) the 7.4 km/s layer does not correspond to the “reflection” Moho. Gravity modeling and comparison to existing seismic models show a general consensus in many aspects with respect to crustal structure. This new model forms the basis for speculation on some of the processes associated with rifting of the Basin and Range Province. One such process, lithospheric magmatism, is inferred from the strong attenuation of transmitted seismic waves, which occurs at the same interface at which high‐amplitude, bright spot reflections originate. Unlike previous models, the overall structure and velocity of the crust and uppermost mantle of our new model are similar to other regions worldwide which have undergone high degrees of extension.

Nevada↗
USGS2010

Paleomagnetism and rotation constraints for the middle Miocene southwestern Nevada volcanic field

Middle Miocene rocks of the southwestern Nevada volcanic field (SWNVF) lie across the projection of the Walker Lane belt within the Basin and Range province and thus provide an interesting opportunity to test for late Cenozoic vertical-axis rotation. Paleomagnetic data from individual ash flow sheets document no significant relative vertical-axis rotation among localities within central SWNVF, an area of relatively low stratal tilts and widely spaced faults. A time-averaged mean paleomagnetic direction (D = 351.4°, I = 52.7°, α 95 = 4.5°) calculated from data from numerous separate rock units suggests that the central SWNVF underwent minimal counterclockwise vertical-axis rotation (R = −7.1° ± 6.6°) with respect to the North American craton. No clockwise vertical-axis rotation is found to support projection of dextral faults of the Walker Lane beneath the central SWNVF. Clockwise rotation of variable magnitude is common at numerous sites from southern and western margins of the field. These clockwise rotations probably reflect dextral shear strain developed at the interface between the little extended central SWNVF block and more strongly extended areas to the south and southwest of the field. Negligible rotation of 11.45-Ma to 13.25-Ma tuffs relative to the central SWNVF was found at the southeast margin of the field where 90° clockwise rotation at the northwest termination of the Las Vegas Valley shear zone had been postulated. Any clockwise rotation in this area must predate 13.25 Ma, and thus dextral shear within this part of the Walker Lane belt was not synchronous or connected across the southern margin of the field. Small counterclockwise vertical-axis rotation relative to the craton, as found for the central SWNVF block, might be a regional feature in the western Great Basin.

Tectonics↗
USGS2010

Local energy flux estimates for unstable conditions using variance data in semiarid rangelands

A network of meteorological stations was installed during the Monsoon '90 field campaign in the Walnut Gulch experimental watershed. The study area has a fairly complex surface. The vegetation cover is heterogeneous and sparse, and the terrain is mildly hilly, but dissected by ephemeral channels. Besides measurement of some of the standard weather data such as wind speed, air temperature, and solar radiation, these sites also contained instruments for estimating the local surface energy balance. The approach utilized measurements of net radiation ( R n ), soil heat flux ( G ) and Monin-Obukhov similarity theory applied to first- and second-order turbulent statistics of wind speed and temperature for determining the sensible heat flux ( H ). The latent heat flux ( LE ) was solved as a residual in the surface energy balance equation, namely, LE = −( R n + G + H ). This procedure (VAR-RESID) for estimating the energy fluxes satisfied monetary constraints and the requirement for low maintenance and continued operation through the harsh environmental conditions experienced in semiarid regions. Comparison of energy fluxes using this approach with more traditional eddy correlation techniques showed differences were within 20% under unstable conditions. Similar variability in flux estimates over the study area was present in the eddy correlation data. Hence, estimates of H and LE using the VAR-RESID approach under unstable conditions were considered satisfactory. Also, with second-order statistics of vertical velocity collected at several sites, the local momentum roughness length was estimated. This is an important parameter used in modeling the turbulent transfer of momentum and sensible heat fluxes across the surface-atmosphere interface.

Water Resources Research↗
USGS1996

Cyclic stressing and seismicity at strongly coupled subduction zones

We use the finite element method to analyze stress variations in and near a strongly coupled subduction zone during an earthquake cycle. Deformation is assumed to be uniform along strike (plane strain on a cross section normal to the trench axis), and periodic earthquake slip is imposed consistent with the long-term rate of plate convergence and degree of coupling. Simulations of stress and displacement rate fields represent periodic fluctuations in time superimposed on an average field. The oceanic plate, descending slab, and continental lithosphere are assumed here to respond elastically to these fluctuations, and the remaining mantle under and between plates is assumed to respond as Maxwell viscoelastic. In the first part of the analysis we find that computed stress fluctuations in space and time are generally consistent with observed earthquake mechanism variations with time since a great thrust event. In particular, trench-normal extensional earthquakes tend to occur early in the earthquake cycle toward the outer rise but occur more abundantly late in the cycle in the subducting slab downdip of the main thrust zone. Compressional earthquakes, when they occur at all, have the opposite pattern. Our results suggest also that the actual timing of extensional outer rise events is controlled by the rheology of the shallow aseismic portion of the thrust interface. The second part of the analysis shows the effects of mantle relaxation on the rate of ground surface deformation during the earthquake cycle. Models without relaxation predict a strong overall compressional strain rate in the continental plate above the main thrust zone, with the strain rate constant between mainshocks. However with significant relaxation present, a localized region of unusually low compressional, or even slightly extensional, strain rate develops along the surface of the continental plate above and somewhat inland from the downdip edge of the locked main thrust zone. The low strain rate starts in the middle or late part of the cycle, depending on position. This result suggests that the negligible or small contraction measured on the Shumagin Islands, Alaska, during 1980 to 1991, may not invalidate an interpretation of that region as being a moderately coupled subduction zone. In contrast, mantle relaxation causes only modest temporal nonuniformity of uplift rates in the overriding plate and of extensional stress rates in the subducting plate, even when the Maxwell time is an order of magnitude less than the recurrence interval.

Journal of Geophysical Research B: Solid Earth↗
USGS1996

Hydrogeology of the Hawaii Scientific Drilling Project borehole KP-1 2. Groundwater geochemistry and regional flow patterns

A series of downhole and surface water samples were taken from the 1‐km‐deep KP‐1 borehole located on the eastern flank of the island of Hawaii. Early samples from depths of more than 700 m showed salinities nearly equivalent to seawater but having anomalous cation concentrations that are attributed to ion exchange between formation fluids and residual drilling mud clays. Later deep samples found only minor variations from seawater cation chemistry that are consistent with low‐temperature weathering of basalts; δ 18 O values are equivalent to seawater values and are consistent with this interpretation. Carbon 14 activities of dissolved inorganic carbonate indicate a water age ranging from 5890 to 7170 years B.P. and fluid transport rates of 1.8 to 2.2 m/yr. Fluid samples from perforations at 310 m in the borehole demonstrate that a freshwater aquifer is present at the Mauna Kea/Mauna Loa interface; borehole resistivity logs indicate that it is ∼200 m thick. Although it has not yet been possible to obtain samples of the freshwater zone without contamination from the deep saline fluids, the chloride concentrations of the low‐salinity zone are estimated using a mixing enthalpy calculation to be less than 100 mg/L. Light stable isotope data indicate that the fresh water at 320 m is derived from recharge entering the island at an average elevation of 2000 m. Inferred 14 C activities of the dissolved bicarbonate in the freshwater zone indicate an average calibrated age of 2200 years B.P. and an average fluid velocity of at least 14 m/yr. A regional water flow model is proposed that suggests that the fresh water found at the 320‐m depth is derived from rainfall recharge from the middle elevations of Mauna Kea volcano. This rainfall is channeled beneath the Mauna Loa lavas by the thick soil layer separating the two volcanoes. A second shallow fresh‐to‐brackish water zone, derived from Mauna Loa recharge, is also inferred to exist below the carbonate formation that underlies the shallow basal lens. The results of our preliminary study of the groundwater system below the KP‐1 drill site demonstrate that intervolcano and interflow aquicludes can have a substantial impact on water circulation and discharge from young island volcanoes.

Journal of Geophysical Research B: Solid Earth↗
USGS2010

State-dependent anisotrophy: Comparison of quasi-analytical solutions with stochastic results for steady gravity drainage

Anisotropy in large-scale unsaturated hydraulic conductivity of layered soils changes with the moisture state. Here, state-dependent anisotropy is computed under conditions of large-scale gravity drainage. Soils represented by Gardner's exponential function are perfectly stratified, periodic, and inclined. Analytical integration of Darcy’s law across each layer results in a system of nonlinear equations that is solved iteratively for capillary suction at layer interfaces and for the Darcy flux normal to layering. Computed fluxes and suction profiles are used to determine both upscaled hydraulic conductivity in the principal directions and the corresponding “state-dependent” anisotropy ratio as functions of the mean suction. Three groups of layered soils are analyzed and compared with independent predictions from the stochastic results of Yeh et al. (1985b). The small-perturbation approach predicts appropriate behaviors for anisotropy under nonarid conditions. However, the stochastic results are limited to moderate values of mean suction; this limitation is linked to a Taylor series approximation in terms of a group of statistical and geometric parameters. Two alternative forms of the Taylor series provide upper and lower bounds for the state-dependent anisotropy of relatively dry soils.

Water Resources Research↗
USGS2010

Hydrologic and geochemical effects on oxygen uptake in bottom sediments of an effluent-dominated river

More than 95% of the water in the South Platte River downstream from the largest wastewater treatment plant serving the metropolitan Denver, Colorado, area consists of treated effluent during some periods of low flow. Fluctuations in effluent-discharge rates caused daily changes in river stage that promoted exchange of water between the river and bottom sediments. Groundwater discharge measurements indicated fluxes of water across the sediment-water interface as high as 18 m 3 s −1 km −1 . Laboratory experiments indicated that downward movement of surface water through bottom sediments at velocities comparable to those measured in the field (median rate ≈0.005 cm s −1 ) substantially increased dissolved oxygen uptake rates in bottom sediments (maximum rate 212 ± 10 μmol O 2 L −1 h −1 ) compared with rates obtained when no vertical advective flux was generated (maximum rate 25 ± 8.8 μmol O 2 L −1 h −1 ). Additions of dissolved ammonium to surface waters generally increased dissolved oxygen uptake rates relative to rates measured in experiments without ammonium. However, the magnitude of the advective flux through bottom sediments had a greater effect on dissolved oxygen uptake rates than did the availability of ammonium. Results from this study indicated that efforts to improve dissolved oxygen dynamics in effluent-dominated rivers might include stabilizing daily fluctuations in river stage.

Colorado↗
USGS1996

Phytoplankton bloom dynamics in coastal ecosystems: A review with some general lessons from sustained investigation of San Francisco Bay, California

Phytoplankton blooms are prominent features of biological variability in shallow coastal ecosystems such as estuaries, lagoons, bays, and tidal rivers. Long-term observation and research in San Francisco Bay illustrates some patterns of phytoplankton spatial and temporal variability and the underlying mechanisms of this variability. Blooms are events of rapid production and accumulation of phytoplankton biomass that are usually responses to changing physical forcings originating in the coastal ocean (e.g., tides), the atmosphere (wind), or on the land surface (precipitation and river runoff). These physical forcings have different timescales of variability, so algal blooms can be short-term episodic events, recurrent seasonal phenomena, or rare events associated with exceptional climatic or hydrologic conditions. The biogeochemical role of phytoplankton primary production is to transform and incorporate reactive inorganic elements into organic forms, and these transformations are rapid and lead to measurable geochemical change during blooms. Examples include the depletion of inorganic nutrients (N, P, Si), supersaturation of oxygen and removal of carbon dioxide, shifts in the isotopic composition of reactive elements (C, N), production of climatically active trace gases (methyl bromide, dimethylsulfide), changes in the chemical form and toxicity of trace metals (As, Cd, Ni, Zn), changes in the biochemical composition and reactivity of the suspended particulate matter, and synthesis of organic matter required for the reproduction and growth of heterotrophs, including bacteria, zooplankton, and benthic consumer animals. Some classes of phytoplankton play special roles in the cycling of elements or synthesis of specific organic molecules, but we have only rudimentary understanding of the forces that select for and promote blooms of these species. Mounting evidence suggests that the natural cycles of bloom variability are being altered on a global scale by human activities including the input of toxic contaminants and nutrients, manipulation of river flows, and translocation of species. This hypothesis will be a key component of our effort to understand global change at the land-sea interface. Pursuit of this hypothesis will require creative approaches for distinguishing natural and anthropogenic sources of phytoplankton population variability, as well as recognition that the modes of human disturbance of coastal bloom cycles operate interactively and cannot be studied as isolated processes.

California↗
USGS1997

Multiphase groundwater flow near cooling plutons

We investigate groundwater flow near cooling plutons with a computer program that can model multiphase flow, temperatures up to 1200°C, thermal pressurization, and temperature‐dependent rock properties. A series of experiments examines the effects of host‐rock permeability, size and depth of pluton emplacement, single versus multiple intrusions, the influence of a caprock, and the impact of topographically driven groundwater flow. We also reproduce and evaluate some of the pioneering numerical experiments on flow around plutons. Host‐rock permeability is the principal factor influencing fluid circulation and heat transfer in hydrothermal systems. The hottest and most steam‐rich systems develop where permeability is of the order of 10 −15 m 2 . Temperatures and life spans of systems decrease with increasing permeability. Conduction‐dominated systems, in which permeabilities are ≤10 −16 m 2 , persist longer but exhibit relatively modest increases in near‐surface temperatures relative to ambient conditions. Pluton size, emplacement depth, and initial thermal conditions have less influence on hydrothermal circulation patterns but affect the extent of boiling and duration of hydrothermal systems. Topographically driven groundwater flow can significantly alter hydrothermal circulation; however, a low‐permeability caprock effectively decouples the topographically and density‐driven systems and stabilizes the mixing interface between them thereby defining a likely ore‐forming environment.

Journal of Geophysical Research B: Solid Earth↗
USGS1998

Nonlinear refraction and reflection travel time tomography

We develop a rapid nonlinear travel time tomography method that simultaneously inverts refraction and reflection travel times on a regular velocity grid. For travel time and ray path calculations, we apply a wave front method employing graph theory. The first-arrival refraction travel times are calculated on the basis of cell velocities, and the later refraction and reflection travel times are computed using both cell velocities and given interfaces. We solve a regularized nonlinear inverse problem. A Laplacian operator is applied to regularize the model parameters (cell slownesses and reflector geometry) so that the inverse problem is valid for a continuum. The travel times are also regularized such that we invert travel time curves rather than travel time points. A conjugate gradient method is applied to minimize the nonlinear objective function. After obtaining a solution, we perform nonlinear Monte Carlo inversions for uncertainty analysis and compute the posterior model covariance. In numerical experiments, we demonstrate that combining the first arrival refraction travel times with later reflection travel times can better reconstruct the velocity field as well as the reflector geometry. This combination is particularly important for modeling crustal structures where large velocity variations occur in the upper crust. We apply this approach to model the crustal structure of the California Borderland using ocean bottom seismometer and land data collected during the Los Angeles Region Seismic Experiment along two marine survey lines. Details of our image include a high-velocity zone under the Catalina Ridge, but a smooth gradient zone between. Catalina Ridge and San Clemente Ridge. The Moho depth is about 22 km with lateral variations. Copyright 1998 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS1999

Coseismic slip resolution along a plate boundary megathrust: the Nankai Trough, southwest Japan

Geodetic survey measurements are used to estimate the coseismic slip distribution in the 1944 Tonankai (M w =8.1) and 1946 Nankaido (M w =8.3) earthquakes and to assess quantitatively the degree to which this slip is resolved on the plate boundary megathrust. Data used include 798 angle changes from triangulation surveys, 328 leveling section differences, and 5 coseismic tidal gage offsets. Many of the nominally coseismic triangulation data span ∼50 years, nearly half the earthquake cycle, and correction for interseismic deformation using post-1950 observations is applied. Microseismicity is used to define the configuration of the plate boundary interface and approximate it with a continuous, multisegment fault model. Because the onshore geodetic data have very limited resolving power for offshore fault segments, offshore coseismic slip was constrained by Satctke's [1993] estimation based on tsunami data. The majority of the coseismic slip occurs between 15 and 25 km depth. Although resolution declines toward the trench axis, it is sufficiently good to define two distinct high-slip regions, one off southeastern Shikoku Island (11 m maximum) and the other offshore of Kii Peninsula (3 m maximum). The slip magnitude off southeastern Shikoku, coupled with the plate convergence rate, would imply an recurrence interval of about 270 years, much-longer than the average repeat time of ∼120 years for historical great earthquakes on the Nankai Trough. However, the maximum coseismic slip is sensitive to the assumed fault geometry, and slippage on trough-parallel splay faults could significantly decrease the maximum slip to about 6 m.

Journal of Geophysical Research B: Solid Earth↗
USGS1998

Non-double-couple earthquakes. 1. Theory

Historically, most quantitative seismological analyses have been based on the assumption that earthquakes are caused by shear faulting, for which the equivalent force system in an isotropic medium is a pair of force couples with no net torque (a “double couple,” or DC). Observations of increasing quality and coverage, however, now resolve departures from the DC model for many earthquakes and find some earthquakes, especially in volcanic and geothermal areas, that have strongly non-DC mechanisms. Understanding non-DC earthquakes is important both for studying the process of faulting in detail and for identifying nonshear-faulting processes that apparently occur in some earthquakes. This paper summarizes the theory of “moment tensor” expansions of equivalent-force systems and analyzes many possible physical non-DC earthquake processes. Contrary to long-standing assumption, sources within the Earth can sometimes have net force and torque components, described by first-rank and asymmetric second-rank moment tensors, which must be included in analyses of landslides and some volcanic phenomena. Non-DC processes that lead to conventional (symmetric second-rank) moment tensors include geometrically complex shear faulting, tensile faulting, shear faulting in an anisotropic medium, shear faulting in a heterogeneous region (e.g., near an interface), and polymorphic phase transformations. Undoubtedly, many non-DC earthquake processes remain to be discovered. Progress will be facilitated by experimental studies that use wave amplitudes, amplitude ratios, and complete waveforms in addition to wave polarities and thus avoid arbitrary assumptions such as the absence of volume changes or the temporal similarity of different moment tensor components.

Reviews of Geophysics↗
USGS1998

Effect of groundwater springs on NO3− concentrations during summer in Catskill Mountain streams

Groundwater and stream water data collected at three headwater catchments in the Neversink River watershed indicate that base flow is sustained by groundwater from two sources: a shallow flow system within the till and soil and a deep flow system within bedrock fractures and bedding planes that discharges as perennial springs. Data from eight wells finished near the till/bedrock interface indicate that saturated conditions are not maintained in the shallow flow system during most summers. In contrast, the discharge of a perennial spring remained constant during two summer rainstorms, providing evidence that the deep flow system is disconnected from the shallow flow system in summer. Discharge from perennial springs was the principal source of streamflow in a headwater reach during low flow. Mean NO 3 − concentrations were 20–25 μmol L −1 in five perennial springs during the summer but only 5–10 μmol L −1 in shallow groundwater. Thus the deep flow system does not reflect typical NO 3 − concentrations in the soil during summer. A hydrologic budget at a headwater drainage reveals that March and late fall are the principal groundwater recharge periods. Residence time modeling based on analyses of 18 O and 35 S indicates that groundwater in the deep flow system is 6–22 months old. These data indicate that summer base flow largely originates from previous dormant seasons when available soil NO 3 − is greater. In these Catskill watersheds, high base flow concentrations of NO 3 − during summer do not provide sufficient evidence that the atmospheric N deposition rate exceeds the demand of terrestrial vegetation.

Water Resources Research↗
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

Influence of welded boundaries in anelastic media on energy flow, and characteristics of P, S-I, and S-II waves: Observational evidence for inhomogeneous body waves in low-loss solids

A general computer code, developed to calculate anelastic reflection-refraction coefficients, energy flow, and the physical characteristics for general P , S -I, and S -II waves, quantitatively describes physical characteristics for wave fields in anelastic media that do not exist in elastic media. Consideration of wave fields incident on boundaries between anelastic media shows that scattered wave fields experience reductions in phase and energy speeds, increases in maximum attenuation and Q −1 , and directions of maximum energy flow distinct from phase propagation. Each of these changes in physical characteristics are shown to vary with angle of incidence. Finite relaxation times for anelastic media result in energy flow due to interaction of superimposed radiation fields and contribute to energy flow across anelastic boundaries for all angles of incidence. Agreement of theoretical and numerical results with laboratory measurements argues for the validity of the theoretical and numerical formulations incorporating inhomogeneous wave fields. The agreement attests to the applicability of the model and helps confirm the existence of inhomogeneous body waves and their associated set of distinct physical characteristics in the earth. The existence of such body waves in layered, low-loss anelastic solids implies the need to reformulate some seismological models of the earth. The exact anelastic formulation for a liquid-solid interface with no low-loss approximations predicts the existence of a range of angles of incidence or an anelastic Rayleigh window, through which significant amounts of energy are transmitted across the boundary. The window accounts for the discrepancy apparent between measured reflection data presented in early textbooks and predictions based on classical elasticity theory. Characteristics of the anelastic Rayleigh window are expected to be evident in certain sets of wide-angle, ocean-bottom reflection data and to be useful in estimating Q −1 for some ocean bottom reflectors.

Journal of Geophysical Research B: Solid Earth↗
USGS2012

The crustal structure of the Wrangellia Terrane along the East Glenn Highway, eastern‐southern Alaska

Recently acquired seismic refraction data from eastern‐southern Alaska provide new information on the structure and composition of the Wrangellia and adjacent terranes. The data comprise a 160‐km‐long refraction profile along the East Glenn (Tok‐Cutoff) Highway that was collected as part of the U.S. Geological Survey's multidisciplinary Trans‐Alaska Crustal Transect program. The upper 3 km of the Wrangellia terrane and associated rocks is characterized by low compressional wave velocities (V p = 1.9, 3.3, 4.6, 5.6 km s −1 ) and high‐velocity gradients common to most onshore seismic refraction profiles. There is also clear seismic expression of the West Fork fault system as a steep, down‐to‐the‐southwest fault that separates the Peninsular terrane to the southwest and the metamorphic complex of Gulkana River to the northeast. In contrast, no seismic expression occurs for the Paxson Lake fault, which separates the Wrangellia terrane from the metamorphic complex of Gulkana River. Adjacent to the Denali fault, within the Wrangellia terrane, two high‐velocity bodies (V p = 6.6 km s −1 ) occur in the upper crust. One of these extends to ∼10‐km depth and correlates with a late Paleozoic dioritic complex, suggesting that the Wrangellia terrane is at least 10 km thick in this part of Alaska. From 5 to 23 km depth, the crust appears seismically homogeneous, with velocity increasing from V p = 6.2 to V p = 6.6 km s −1 . Beneath this level, the crust is less well resolved, although evidence exists for a low‐velocity zone between 23 and 26 km and a possible southwest dipping interface at 35 km. No identifiable mantle refraction or reflection is observed, possibly indicating a crust as thick as 55 km. The relatively low seismic velocities in the upper 23 km of the crust compare favorably with laboratory‐measured velocities on pelitic schists and intermediate‐composition plutonic rocks (granites and granodiorites), both of which are recognized in Wrangellia. We interpret the seismic velocities to indicate that silicic‐to‐intermediate‐composition rocks are important constituents of the basement of this part of Wrangellia. Geologic evidence indicates that the Alaskan part of the Wrangellia terrane is a Paleozoic and Mesozoic island arc: our seismic evidence indicates it may have been built mostly on continental crust as opposed to the fragment of Wrangellia from Vancouver Island which was probably built on oceanic crust.

Alaska↗
USGS2012

Sea ice melting in the marginal ice zone

The heat and salt flux boundary conditions together with the freezing curve relationship are a necessary component of any ice-sea water thermodynamic model. A neutral two-layer oceanic planetary boundary layer model that incorporates these boundary conditions gives the following results: The interfacial salinity is within 10% of the far-field salinity for conditions commonly encountered in the MIZ and depends only on the turbulent Lewis number and the far-field temperature and salinity. The predicted melt rates agree with the limited field observations, of the order of 1 m day −1 . The Obukov lengths, determined from the predicted interface conditions and melt rates, are generally much greater than the Ekman layer thicknesses; hence, the surface buoyancy flux has little effect on the turbulence in the planetary boundary layer.

Journal of Geophysical Research - Oceans↗