USGS Science⌕ Search

SEARCH · USGS Science

Results for “Journal of Geophysics Research”

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 667 records · Page 37Linked to original sources
USGS2006

Distribution, 14 C chronology, and paleomagnetism of latest Pleistocene and Holocene lava flows at Haleakala volcano, Island of Maui, Hawai'i: A revision of lava flow hazard zones

New mapping and 60 new radiocarbon ages define the age and distribution of latest Pleistocene and Holocene (past 13,000 years) lava flows at Haleakalā volcano, Island of Maui. Paleomagnetic directions were determined for 118 sites, of which 89 are in lava flows younger than 13,000 years. The paleomagnetic data, in conjunction with a reference paleosecular variation (PSV) curve for the Hawaiian Islands, are combined with our knowledge of age limitations based on stratigraphic control to refine age estimates for some of the undated lava flows. The resulting volumetric rate calculations indicate that within analytical error, the extrusion rate has remained nearly constant during the past 13,000 years, in the range 0.05–0.15 km 3 /kyr, only about half the long-term rate required to produce the postshield strata emplaced in the past ∼1 Myr. Haleakalā's eruptive frequency is similar to that of Hualālai volcano on the Island of Hawai‘i, but its lava flows cover substantially less area per unit time. The reduced rates of lava coverage indicate a lower volcanic hazard than in similar zones at Hualālai.

Hawaii↗
USGS2007

Generation of fine hydromagmatic ash by growth and disintegration of glassy rinds

The deposits of mafic hydromagmatic eruptions are more fine grained and variable in vesicularity than dry magmatic deposits. Blocky, equant shapes of many hydromagmatic clasts also contrast with droplet, thread, and bubble wall morphology of dry magmatic fragments. Small (<∼180 μ m), blocky hydromagmatic pyroclasts have traditionally been interpreted to result from discrete vapor explosions, although such explosions tend to occur only under certain conditions. This paper considers a process of hydromagmatic ash formation that involves repeated growth and disintegration of glassy rinds on pyroclast surfaces as they deform within turbulent flows. This process, termed “turbulent shedding”, may occur during the expansion phase of vapor explosions or during turbulent but nonexplosive mixing of magma with water, steam, or water sprays. The occurrence of turbulent shedding and the resulting fragment sizes depend on the timescale for rind growth and the timescale between disturbances that remove or disintegrate glassy rinds. Turbulent shedding is directly observable in some small littoral jets at Kilauea. Calculations suggest that, in the presence of liquid water or water sprays, glassy rinds having a thickness of microns to millimeters should form in milliseconds to seconds. This is similar to the timescale between turbulent velocity fluctuations that can shred lava globules and remove such rinds. The fraction of a deposit consisting of fine ash should increase with the duration of this process: Large‐scale Surtseyan jets generate hundreds or thousands of shedding events; bubble bursts or tephra jets at Kilauea's coast may produce only a few.

Journal of Geophysical Research B: Solid Earth↗
USGS2006

Uplift and subsidence associated with the great Aceh-Andaman earthquake of 2004

Rupture of the Sunda megathrust on 26 December 2004 produced broad regions of uplift and subsidence. We define the pivot line separating these regions as a first step in defining the lateral extent and the downdip limit of rupture during that great Mw ??? 9.2 earthquake. In the region of the Andaman and Nicobar islands we rely exclusively on the interpretation of satellite imagery and a tidal model. At the southern limit of the great rupture we rely principally on field measurements of emerged coral microatolls. Uplift extends from the middle of Simeulue Island, Sumatra, at ??? 2.5??N, to Preparis Island, Myanmar (Burma), at ??? 14.9??N. Thus the rupture is ??? 1600 km long. The distance from the pivot line to the trench varies appreciably. The northern and western Andaman Islands rose, whereas the southern and eastern portion of the islands subsided. The Nicobar Islands and the west coast of Aceh province, Sumatra, subsided. Tilt at the southern end of the rupture is steep; the distance from 1.5 m of uplift to the pivot line is just 60 km. Our method of using satellite imagery to recognize changes in elevation relative to sea surface height and of using a tidal model to place quantitative bounds on coseismic uplift or subsidence is a novel approach that can be adapted to other forms of remote sensing and can be applied to other subduction zones in tropical regions. Copyright 2006 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS2006

Inference of postseismic deformation mechanisms of the 1923 Kanto earthquake

Coseismic slip associated with the M7.9, 1923 Kanto earthquake is fairly well understood, involving slip of up to 8 m along the Philippine Sea-Honshu interplate boundary under Sagami Bay and its onland extension. Postseismic deformation after the 1923 earthquake, however, is relatively poorly understood. We revisit the available deformation data in order to constrain possible mechanisms of postseismic deformation and to examine the consequences for associated stress changes in the surrounding crust. Data from two leveling lines and one tide gage station over the first 7-8 years postseismic period are of much greater amplitude than the corresponding expected interseismic deformation during the same period, making these data suitable for isolating the signal from postseismic deformation. We consider both viscoelastic models of asthenosphere relaxation and afterslip models. A distributed coseismic slip model presented by Pollitz et al. (2005), combined with prescribed parameters of a viscoelastic Earth model, yields predicted postseismic deformation that agrees with observed deformation on mainland Honshu from Tokyo to the Izu peninsula. Elsewhere (southern Miura peninsula; Boso peninsula), the considered viscoelastic models fail to predict observed deformation, and a model of ???1 in shallow afterslip in the offshore region south of the Boso peninsula, with equivalent moment magnitude Mw = 7.0, adequately accounts for the observed deformation. Using the distributed coseismic slip model, layered viscoelastic structure, and a model of interseismic strain accumulation, we evaluate the post-1923 stress evolution, including both the coseismic and accumulated postseismic stress changes and those stresses contributed by interseismic loading. We find that if account is made for the varying tectonic regime in the region, the occurrence of both immediate (first month) post-1923 crustal aftershocks as well as recent regional crustal seismicity is consistent with the predicted stress pattern. This suggests that the influence of the 1923 earthquake on regional seismicity is fairly predictable and has persisted for at least seven decades following the earthquake.

Journal of Geophysical Research B: Solid Earth↗
USGS2005

Thermal observations of gas pistoning at Kilauea Volcano

Data acquired by three continuously recording thermal infrared thermometers situated on the north rim of Pu'u'O' o Crater at Kilauea Volcano during 2002 revealed episodes of periodic thermal pulses originating from a degassing vent on the crater floor. These thermal pulses are interpreted as gas release (jetting events) associated with gas pistoning, a mechanism observed previously at both Mauna Ulu and Pu'u'O' o. During a 35-day-long period spanning June and July 2002, gas pistoning was frequently the dominant mode of gas release, with as many as several hundred pulses occurring in uninterrupted series. On other days, degassing alternated between periods of quasi-continuous gas jetting and intervals of gas pistoning that contained a few to a few dozen pulses. Characteristic time intervals between pistoning events ranged from 2 up to 7 min. We identify three types of pistoning. Type 1 involves emission of lava, followed by gas jetting and drain back; type 2 is the same but the elevated position of the vent does not allow postjet drain back; and type 3 involves gas jetting only with no precursory lava flow. To explain gas pistoning, we apply a model whereby a stagnant cap of degassed magma develops in the conduit below the vent. Gas bubbles rise through the magma column and collect under the cap. The collective buoyancy of these bubbles pushes the cap upward. When the cap reaches the surface, it erupts from the vent as a lava flow. Unloading of the conduit magma in this way results in an abrupt pressure drop (i.e., the overburden felt by the bubbles is reduced), causing explosive gas expansion in the form of gas jetting from the vent. This terminates the event and lava drains back into the conduit to start the cycle anew. In the case where there is no surface lava emission or drain back, the cap instead pushes into and spreads out within a subsurface cavity. Again, this unloads the conduit magma and terminates in explosive gas release. Once gas is expelled, lava in the cavity is free to drain back. We hypothesize that pistoning is a stable mode of degassing for low-viscosity basaltic magmas with appropriate conduit geometries and volatile supply rates. Copyright 2005 by the American Geophysical Union.

Hawaii↗
USGS2006

Tectonic stressing in California modeled from GPS observations

What happens in the crust as a result of geodetically observed secular motions? In this paper we find out by distorting a finite element model of California using GPS-derived displacements. A complex model was constructed using spatially varying crustal thickness, geothermal gradient, topography, and creeping faults. GPS velocity observations were interpolated and extrapolated across the model and boundary condition areas, and the model was loaded according to 5-year displacements. Results map highest differential stressing rates in a 200-km-wide band along the Pacific-North American plate boundary, coinciding with regions of greatest seismic energy release. Away from the plate boundary, GPS-derived crustal strain reduces modeled differential stress in some places, suggesting that some crustal motions are related to topographic collapse. Calculated stressing rates can be resolved onto fault planes: useful for addressing fault interactions and necessary for calculating earthquake advances or delays. As an example, I examine seismic quiescence on the Garlock fault despite a calculated minimum 0.1-0.4 MPa static stress increase from the 1857 M???7.8 Fort Tejon earthquake. Results from finite element modeling show very low to negative secular Coulomb stress growth on the Garlock fault, suggesting that the stress state may have been too low for large earthquake triggering. Thus the Garlock fault may only be stressed by San Andreas fault slip, a loading pattern that could explain its erratic rupture history.

Journal of Geophysical Research B: Solid Earth↗
USGS2006

Aftershock decay, productivity, and stress rates in Hawaii: Indicators of temperature and stress from magma sources

We examined dozens of aftershock sequences in Hawaii in terms of Gutenberg-Richter and modified Omori law parameters. We studied p, the rate of aftershock decay; A p , the aftershock productivity, defined as the observed divided by the expected number of aftershocks; and c, the time delay when aftershock rates begin to fall. We found that for earthquakes shallower than 20 km, p values >1.2 are near active magma centers. We associate this high decay rate with higher temperatures and faster stress relaxation near magma reservoirs. Deep earthquakes near Kilauea's inferred magma transport path show a range of p values, suggesting the absence of a large, deep magma reservoir. Aftershock productivity is >4.0 for flank earthquakes known to be triggered by intrusions but is normal (0.25 to 4.0) for isolated main shocks. We infer that continuing, post-main shock stress from the intrusion adds to the main shock's stress step and causes higher A p . High A p in other zones suggests less obvious intrusions and pulsing magma pressure near Kilauea's feeding conduit. We calculate stress rates and stress rate changes from pre-main shock and aftershock rates. Stress rate increased after many intrusions but decreased after large M7–8 earthquakes. Stress rates are highest in the seismically active volcano flanks and lowest in areas far from volcanic centers. We found sequences triggered by intrusions tend to have high A p , high (>0.10 day) c values, a stress rate increase, and sometimes a peak in aftershock rate hours after the main shock. We interpret these values as indicating continuing intrusive stress after the main shock.

Hawai'i↗
USGS2006

An evaluation of factors influencing pore pressure in accretionary complexes: Implications for taper angle and wedge mechanics

At many subduction zones, accretionary complexes form as sediment is off-scraped from the subducting plate. Mechanical models that treat accretionary complexes as critically tapered wedges of sediment demonstrate that pore pressure controls their taper angle by modifying basal and internal shear strength. Here, we combine a numerical model of groundwater flow with critical taper theory to quantify the effects of sediment and de??collement permeability, sediment thickness, sediment partitioning between accretion and underthrusting, and plate convergence rate on steady state pore pressure. Our results show that pore pressure in accretionary wedges can be viewed as a dynamically maintained response to factors which drive pore pressure (source terms) and those that limit flow (permeability and drainage path length). We find that sediment permeability and incoming sediment thickness are the most important factors, whereas fault permeability and the partitioning of sediment have a small effect. For our base case model scenario, as sediment permeability is increased, pore pressure decreases from near-lithostatic to hydrostatic values and allows stable taper angles to increase from ??? 2.5?? to 8??-12.5??. With increased sediment thickness in our models (from 100 to 8000 m), increased pore pressure drives a decrease in stable taper angle from 8.4??-12.5?? to <2.5-5??. In general, low-permeability and thick incoming sediment sustain high pore pressures consistent with shallowly tapered geometry, whereas high-permeability and thin incoming sediment should result in steep geometry. Our model results compare favorably with available data from active accretionary complexes. Active margins characterized by a significant proportion of fine-grained sediment within the incoming section, such as northern Antilles and eastern Nankai, exhibit thin taper angles, whereas those characterized by a higher proportion of sandy turbidites, such as Cascadia, Chile, and Mexico, have steep taper angles. Observations from active margins also indicate a strong trend of decreasing taper angle (from >15?? to <4??) with increased sediment thickness (from <1 to 7 km). One key implication is that hydrologic properties may strongly influence the strength of the crust in a wide range of geologic settings. Copyright 2006 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS2006

Dislocation pileup as a representation of strain accumulation on a strike-slip fault

The conventional model of strain accumulation on a vertical transform fault is a discrete screw dislocation in an elastic half-space with the Burgers vector of the dislocation increasing at the rate of relative plate motion. It would be more realistic to replace that discrete dislocation by a dislocation distribution, presumably a pileup in which the individual dislocations are in equilibrium. The length of the pileup depends upon the applied stress and the amount of slip that has occurred at depth. I argue here that the dislocation pileup (the transition on the fault from no slip to slip at the full plate rate) occupies a substantial portion of the lithosphere thickness. A discrete dislocation at an adjustable depth can reproduce the surface deformation profile predicted by a pileup so closely that it will be difficult to distinguish between the two models. The locking depth (dislocation depth) of that discrete dislocation approximation is substantially (???30%) larger than that (depth to top of the pileup) in the pileup model. Thus, in inverting surface deformation data using the discrete dislocation model, the locking depth in the model should not be interpreted as the true locking depth. Although dislocation pileup models should provide a good explanation of the surface deformation near the fault trace, that explanation may not be adequate at greater distances from the fault trace because approximating the expected horizontally distributed deformation at subcrustal depths by uniform slip concentrated on the fault is not justified.

Journal of Geophysical Research B: Solid Earth↗
USGS2006

A new 1649-1884 catalog of destructive earthquakes near Tokyo and implications for the long-term seismic process

In order to assess the long-term character of seismicity near Tokyo, we construct an intensity-based catalog of damaging earthquakes that struck the greater Tokyo area between 1649 and 1884. Models for 15 historical earthquakes are developed using calibrated intensity attenuation relations that quantitatively convey uncertainties in event location and magnitude, as well as their covariance. The historical catalog is most likely complete for earthquakes M ??? 6.7; the largest earthquake in the catalog is the 1703 M ??? 8.2 Genroku event. Seismicity rates from 80 years of instrumental records, which include the 1923 M = 7.9 Kanto shock, as well as interevent times estimated from the past ???7000 years of paleoseismic data, are combined with the historical catalog to define a frequency-magnitude distribution for 4.5 ??? M ??? 8.2, which is well described by a truncated Gutenberg-Richter relation with a b value of 0.96 and a maximum magnitude of 8.4. Large uncertainties associated with the intensity-based catalog are propagated by a Monte Carlo simulation to estimations of the scalar moment rate. The resulting best estimate of moment rate during 1649-2003 is 1.35 ?? 1026 dyn cm yr-1 with considerable uncertainty at the 1??, level: (-0.11, + 0.20) ?? 1026 dyn cm yr-1. Comparison with geodetic models of the interseismic deformation indicates that the geodetic moment accumulation and likely moment release rate are roughly balanced over the catalog period. This balance suggests that the extended catalog is representative of long-term seismic processes near Tokyo and so can be used to assess earthquake probabilities. The resulting Poisson (or time-averaged) 30-year probability for M ??? 7.9 earthquakes is 7-11%.

Journal of Geophysical Research B: Solid Earth↗
USGS2006

Analysis of gravity and topography in the GLIMPSE study region: Isostatic compensation and uplift of the Sojourn and Hotu Matua Ridge systems

The Gravity Lieations Intraplate Melting Petrologic and Seismic Expedition (GLIMPSE) Experiment investigated the formation of a series of non-hot spot, intraplate volcanic ridges in the South Pacific and their relationship to cross-grain gravity lineaments detected by satellite altimetry. Using shipboard gravity measurements and a simple model of surface loading of a thin elastic plate, we estimate effective elastic thicknesses ranging from ???2 km beneath the Sojourn Ridge to a maximum of 10 km beneath the Southern Cross Seamount. These elastic thicknesses are lower than predicted for the 3-9 Ma seafloor on which the volcanoes lie, perhaps due to reheating and thinning of the plate during emplacement. Anomalously low apparent densities estimated for the Matua and Southern Cross seamounts 2050 and 2250 kg m-3, respectively, probably are artifacts caused by the assumption of only surface loading, ignoring the presence of subsurface loading in the form of underplated crust and/or low-density mantle. Using satellite free-air gravity and shipboard bathymetry, we calculate the age-detrended, residual mantle Bouguer anomaly (rMBA). The rMBA corrects the free-air anomaly for the direct effects of topography, including the thickening of the crust beneath the seamounts and volcanic ridges due to surface loading of the volcanic edifices. There are broad, negative rMBA anomalies along the Sojourn and Brown ridges and the Hotu Matua seamount chain that extend nearly to the East Pacific Rise. These negative rMBA anomalies connect to negative free-air anomalies in the western part of the study area that have been recognized previously as the beginnings of the cross-grain gravity lineaments. Subtracting the topographic effects of surface loading by the ridges and seamounts from the observed topography reveals that the ridges are built on broad bands of anomalously elevated seafloor. This swell topography and the negative rMBA anomalies contradict the predictions of lithospheric cracking models for the origin of gravity lineaments and associated volcanic ridges, favoring models with a dynamic mantle component such as small-scale convection or channelized asthenospheric return flow. Copyright 2006 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS2006

Distribution of recent volcanism and the morphology of seamounts and ridges in the GLIMPSE study area: Implications for the lithospheric cracking hypothesis for the origin of intraplate, non-hot spot volcanic chains

Lithospheric cracking by remotely applied stresses or thermoelastic stresses has been suggested to be the mechanism responsible for the formation of intraplate volcanic ridges in the Pacific that clearly do not form above fixed hot spots. As part of the Gravity Lineations Intraplate Melting Petrology and Seismic Expedition (GLIMPSE) project designed to investigate the origin of these features, we have mapped two volcanic chains that are actively forming to the west of the East Pacific Rise using multibeam echo sounding and side-scan sonar. Side-scan sonar reveals the distribution of rough seafloor corresponding to recent, unsedimented lava flows. In the Hotu Matua volcanic complex, recent flows and volcanic edifices are distributed over a region 450 km long and up to 65 km wide, with an apparent, irregular age progression from older flows in the west to younger in the east. The 550-km-long Southern Cross Seamount/Sojourn Ridge/Brown Ridge chain appears to have been recently active only at its eastern end near the East Pacific Rise. A third region of recent flows is found 120 km north of Southern Cross Seamount in seafloor approximately 9 Myr old. No indication of lithospheric extension in the form of faulting or graben formation paralleling the trend of the volcanic chains is found in the vicinity of recent flows or anywhere else in the study area. Thermoelastic cracking could be a factor in the formation of a few small, very narrow volcanic ridges, but most of the volcanic activity is broadly distributed in wide swaths with no indication of formation along narrow cracks. The Sojourn and Brown chains appear to begin as distributed zones of small seamounts that later develop into segmented ridges, perhaps under the influence of membrane stresses from self-loading. We suggest that the linear volcanic chains are created by moving melting anomalies in the asthenosphere and that lithospheric cracking plays at most a secondary role. Copyright 2006 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS2006

Regional P wave velocity structure of the Northern Cascadia Subduction Zone

This paper presents the first regional three-dimensional, P wave velocity model for the Northern Cascadia Subduction. Zone (SW British Columbia and NW Washington State) constructed through tomographic inversion of first-arrival traveltime data from active source experiments together with earthquake traveltime data recorded at permanent stations. The velocity model images the structure of the subducting Juan de Fuca plate, megathrust, and the fore-arc crust and upper mantle. Beneath southern Vancouver Island the megathrust above the Juan de Fuca plate is characterized by a broad zone (25-35 km depth) having relatively low velocities of 6.4-6.6 km/s. This relative low velocity zone coincides with the location of most of the episodic tremors recently mapped beneath Vancouver Island, and its low velocity may also partially reflect the presence of trapped fluids and sheared lower crustal rocks. The rocks of the Olympic Subduction Complex are inferred to deform aseismically as evidenced by the lack of earthquakes withi the low-velocity rocks. The fore-arc upper mantle beneath the Strait of Georgia and Puget Sound is characterized by velocities of 7.2-7.6 km/s. Such low velocities represent regional serpentinization of the upper fore-arc mantle and provide evidence for slab dewatering and densification. Tertiary sedimentary basins in the Strait of Georgia and Puget Lowland imaged by the velocity model lie above the inferred region of slab dewatering and densification and may therefore partly result from a higher rate of slab sinking. In contrast, sedimentary basins in the Strait of Juan de Fuca lie in a synclinal depression in the Crescent Terrane. The correlation of in-slab earthquake hypocenters M>4 with P wave velocities greater than 7.8 km/s at the hypocenters suggests that they originate near the oceanic Moho of the subducting Juan de Fuca plate. Copyright 2006 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS2006

Damped regional-scale stress inversions: Methodology and examples for southern California and the Coalinga aftershock sequence

We present a new focal mechanism stress inversion technique to produce regional-scale models of stress orientation containing the minimum complexity necessary to fit the data. Current practice is to divide a region into small subareas and to independently fit a stress tensor to the focal mechanisms of each subarea. This procedure may lead to apparent spatial variability that is actually an artifact of overfitting noisy data or nonuniquely fitting data that does not completely constrain the stress tensor. To remove these artifacts while retaining any stress variations that are strongly required by the data, we devise a damped inversion method to simultaneously invert for stress in all subareas while minimizing the difference in stress between adjacent subareas. This method is conceptually similar to other geophysical inverse techniques that incorporate damping, such as seismic tomography. In checkerboard tests, the damped inversion removes the stress rotation artifacts exhibited by an undamped inversion, while resolving sharper true stress rotations than a simple smoothed model or a moving-window inversion. We show an example of a spatially damped stress field for southern California. The methodology can also be used to study temporal stress changes, and an example for the Coalinga, California, aftershock sequence is shown. We recommend use of the damped inversion technique for any study examining spatial or temporal variations in the stress field.

Journal of Geophysical Research B: Solid Earth↗
USGS2006

Imaging the transition from Aleutian subduction to Yakutat collision in central Alaska, with local earthquakes and active source data

In southern and central Alaska the subduction and active volcanism of the Aleutian subduction zone give way to a broad plate boundary zone with mountain building and strike-slip faulting, where the Yakutat terrane joins the subducting Pacific plate. The interplay of these tectonic elements can be best understood by considering the entire region in three dimensions. We image three-dimensional seismic velocity using abundant local earthquakes, supplemented by active source data. Crustal low-velocity correlates with basins. The Denali fault zone is a dominant feature with a change in crustal thickness across the fault. A relatively high-velocity subducted slab and a low-velocity mantle wedge are observed, and high Vp/Vs beneath the active volcanic systems, which indicates focusing of partial melt. North of Cook Inlet, the subducted Yakutat slab is characterized by a thick low-velocity, high-Vp/Vs, crust. High-velocity material above the Yakutat slab may represent a residual older slab, which inhibits vertical flow of Yakutat subduction fluids. Alternate lateral flow allows Yakutat subduction fluids to contribute to Cook Inlet volcanism and the Wrangell volcanic field. The apparent northeast edge of the subducted Yakutat slab is southwest of the Wrangell volcanics, which have adakitic composition consistent with melting of this Yakutat slab edge. In the mantle, the Yakutat slab is subducting with the Pacific plate, while at shallower depths the Yakutat slab overthrusts the shallow Pacific plate along the Transition fault. This region of crustal doubling within the shallow slab is associated with extremely strong plate coupling and the primary asperity of the Mw 9.2 great 1964 earthquake. Copyright 2006 by the American Geophysical Union.

Alaska↗
USGS2007

Microplate model for the present-day deformation of Tibet

[1] Site velocities from 349 Global Positioning System (GPS) stations are used to construct an 11-element quasi-rigid block model of the Tibetan Plateau and its surroundings. Rigid rotations of five major blocks are well determined, and average translation velocities of six smaller blocks can be constrained. Where data are well distributed the velocity field can be explained well by rigid block motion and fault slip across block boundaries. Residual misfits average 1.6 mm/yr compared to typical one standard deviation velocity uncertainties of 1.3 mm/yr. Any residual internal straining of the blocks is small and heterogeneous. However, residual substructure might well represent currently unresolved motions of smaller blocks. Although any smaller blocks must move at nearly the same rate as the larger blocks within which they lie, undetected relative motions between them could be significant, particularly where there are gaps in GPS coverage. Predicted relative motions between major blocks agree with the observed sense of slip and along-strike partitioning of motion across major faults. However, predicted slip rates across Tibet's major strike-slip faults are low, only 5–12 mm/yr, a factor of 2–3 smaller than most rates estimated from fault offset features dated by radiometric methods as ∼2000 to ∼100,000 year old. Previous work has suggested that both GPS data and low fault slip rates are incompatible with rigid block motions of Tibet. The results reported here overcome these objections.

Journal of Geophysical Research B: Solid Earth↗
USGS2006

February 2003 marine atmospheric conditions and the bora over the northern Adriatic

A winter oceanographic field experiment provided an opportunity to examine the atmospheric marine conditions over the northern Adriatic. Mean February winds are from a northeasterly direction over most of the Adriatic and a more northerly direction along the western coast. Wind speeds are fastest in jets over the NE coast during bora events and weakest in the mid-northwestern Adriatic. Diurnal air temperature cycles are smallest on the NE coast and largest in the midwestern Adriatic. The maximum sea-air difference is +10??C on the eastern coast and near zero on the midwestern Adriatic. Boras are northeasterly (from) wind events that sweep off Croatia and Slovenia, bringing slightly colder and drier air over the northern Adriatic. The main bora season is December to March. Winter 2002-2003 was normal for bora events. Synoptic-scale temporal variations are correlated over the northern Adriatic. Fastest Bora winds and highest wind stress over the northern Adriatic is concentrated in four topographically controlled jets. The strongest is the Senj Jet, while the Trieste Jet extends across the entire northern Adriatic. Between each two jets is a weak wind zone. The greatest mean net heat loss is in bora jets in the NE Adriatic, where it was -438 W m-2 and is weakest in the midwestern northern Adriatic, where it was near zero. Wind stress is concentrated over the NE half of Adriatic in four bora jets, while wind stress is weak in the NW Adriatic. There is significant variation in wind stress mean and standard deviation structure over the northern Adriatic with each bora event. Copyright 2006 by the American Geophysical Union.

Journal of Geophysical Research C: Oceans↗
USGS2006

Surface drifter derived circulation in the northern and middle Adriatic Sea: Response to wind regime and season

More than 120 satellite-tracked drifters were deployed in the northern and middle Adriatic (NMA) Sea between September 2002 and November 2003, with the purpose of studying the surface circulation at mesoscale to seasonal scale in relation to wind forcing, river runoff, and bottom topography. Pseudo-Eulerian and Lagrangian statistics were calculated from the low-pass-filtered drifter velocity data between September 2002 and December 2003. The structure of the mean circulation is determined with unprecedented high horizontal resolution by the new data. In particular, mean currents, velocity variance, and kinetic energy levels are shown to be maximal in the Western Adriatic Current (WAC). Separating data into seasons, we found that the mean kinetic energy is maximal in fall, with high values also in winter, while it is significantly weaker in summer. High-resolution Local Area Model Italy winds were used to relate the drifter velocities to the wind fields. The surface currents appear to be significantly influenced by the winds. The mean flow during the northeasterly bora regime shows an intensification of the across-basin recirculating currents. In addition, the WAC is strongly intensified both in intensity and in its offshore lateral extension. In the southeasterly sirocco regime, northward flow without recirculation dominates in the eastern half of the basin, while during northwesterly maestro the WAC is enhanced. Separating the data into low and high Po River discharge rates for low-wind conditions shows that the WAC and the velocity fluctuations in front of the Po delta are stronger for high Po River runoff. Lagrangian covariance, diffusivity, and integral time and space scales are larger in the along-basin direction and are maximal in the southern portion of the WAC. Copyright 2006 by the American Geophysical Union.

Journal of Geophysical Research C: Oceans↗