mechanics of oblique convergence the distribution and magnitude of the strike parallel component of velocity in an obliquely converging thrust wedge or accretionary prism are determined by the geometry and mechanical properties of the wedge a mechanical analysis based on the assumption of a critical or stable geometry of the wedge for which the rate of cross strike deformation is zero leads to the following conclusions for different bulk rheologies 1 in a linear viscous wedge the strike parallel motion relative to the underthrust slab decreases exponentially away from the rear and is effectively concentrated in a shear zone with a width comparable to the thickness of the wedge at the rear the wedge also deforms by corner flow producing a circulation in the cross strike plane the strike parallel and corner flow velocities depend on the thickness and viscosity of the wedge and on the shear stresses applied to its lower and rear boundaries convergence at the wedge front is normal to strike 2 a critically tapered perfect plastic wedge moves coherently without internal deformation for low and moderate obliquities of the convergence vector the wedge moves at the same velocity as the backstop upper plate for high angles of obliquity the wedge moves laterally relative to the underthrust slab at a maximum velocity dependent on its dimensions and the stress conditions on its boundaries so that it is separated from the upper plate by a strike slip fault defining a forearc sliver no geometrical configuration exists that allows the strike parallel motion to be distributed through the wedge 3 a noncohesive coulomb wedge behaves in much the same way as a plastic wedge but the geometry and velocity depend only on its mechanical properties and the shear stresses on its boundaries and they are independent of scale

geophysical geological transect and tectonic evolution of the swiss italian alps a complete alpine cross section integrates numerous seismic reflection and refraction profiles across and along strike with published and new field data the deepest parts of the profile are constrained by geophysical data only while structural features at intermediate levels are largely depicted according to the results of three dimensional models making use of seismic and field geological data the geometry of the highest structural levels is constrained by classical along strike projections of field data parallel to the pronounced easterly axial dip of all tectonic units because the transect is placed close to the western erosional margin of the austroalpine nappes of the eastern alps it contains all the major tectonic units of the alps a model for the tectonic evolution along the transect is proposed in the form of scaled and area balanced profile sketches shortening within the austroalpine nappes is testimony of a separate cretaceous age orogenic event west directed thrusting in these units is related to westward propagation of a thrust wedge resulting from continental collision along the meliata hallstatt ocean further to the east considerable amounts of oceanic and continental crustal material were subducted during tertiary orogeny which involved some 500 km of n s convergence between europe and apulia consequently only a very small percentage of this crustal material is preserved within the nappes depicted in the transect postcollisional shortening is characterized by the simultaneous activity of gently dipping north directed detachments and steeply inclined south directed detachments both detachments nucleating at the interface between lower and upper crust large scale wedging of the adriatic or apulian lower crust into a gap opening between the subduced european lower crust and the pile of thin upper crustal flakes alpine nappes indicates a relatively strong lower crust and detachment between upper and lower crust

plate coupling in the northern south island and southernmost north island new zealand as illuminated by earthquake focal mechanisms subduction of the pacific plate in the northern south island and southernmost north island of new zealand is transitional insofar as the crustal thickness of the pacific plate increases significantly along strike in the northern south island focal mechanisms of 145 events shallower than 100 km in this region have been determined using both first motion polarity data and amplitudes of seismogram envelopes the stress regime in the subducted plate appears to be dominated by slab pull t axes in both the upper and lower planes of the dipping seismic zone generally parallel the local dip of the zone and the average azimuth of these t axes is rotated some 25 clockwise out of the direction of dip of the subducted plate this can be related to the asymmetrical shape of the subducted slab in contrast the stress regime in the overlying plate appears to be dominated by subhorizontal compression low angle thrust events near the plate interface in cook strait and the southernmost north island concentrate in two areas which may mark the updip and downdip edges of a locked region identified from global positioning system gps observations an absence of low angle thrust events near the plate interface in the northern south island and the tendency of p axes of events in the subducted plate to become more horizontal suggest that plate coupling there is stronger than in the southernmost north island differential coupling at the plate interface provides a viable mechanism for producing the large tectonic rotations seen in the northern south island

lithospheric deformation during the early stages of continental collision numerical experiments and comparison with south island new zealand 1 the nature of lithospheric deformation during continental plate collision still remains unresolved while it has often been proposed that the mantle lithosphere is accommodated by distributed thickening of a viscous root an alternate hypothesis suggests that significant portions of convergent mantle lithosphere essentially undergo underthrusting or subduction to further consider this issue we model the thermochemical evolution of the lithosphere mantle system using arbitrary lagrangian eulerian finite element techniques we incorporate a mix of viscous thermally activated power law creep and plastic frictional coulomb with strain softening rheologies in the numerical experiments to treat disparate composition in the crust and mantle a range of rheological and mechanical parameters is explored to determine controls on the style of lithospheric deformation the models suggest that during the initial stages of plate collision the mantle lithosphere is characterized by plate like behavior and underthrusting subduction of the upper region in conjunction with distributed thickening and rayleigh taylor type viscous instability of the lower portion depending on the material rheology temperature regime and imposed convergence velocity the deforming mantle lithosphere demonstrates various combinations of these end member behavioral modes the modeling results are interpreted in the context of observed lithospheric deformation across south island new zealand a combined style of underthrusting and distributed thickening is consistent with the observed crustal structure in the young collisional orogen as well as seismic imaging of the geometry of the underlying lithosphere

effect of unbalanced topography and overloading on coulomb wedge kinematics insights from sandbox modeling 1 this study addresses the effect of variable unbalanced topography and overload on the kinematics of a fold and thrust belt developed within a collisional belt that underwent a subduction polarity reversal event this was done by physical modeling of doubly vergent coulomb wedges using sand as an analogue material during the experimental procedure a preexisting topography was generated by a first phase of subduction in one direction a second phase of subduction was then initiated in the opposite direction simulating a subduction flip an anomalous strong frontal growth for the wedge during the second phase was experimentally shown to be dependent upon surface slope breaks and the critical asymmetrical architecture achieved by the wedge at mature stages of deformation the latter is suggested to be the rule for doubly vergent orogens at steady state even after a subduction flip during the experiments surface processes like syntectonic erosion and sedimentation markedly altered mass transfer within the wedge in particular lowering the surface slope by syntectonic erosion favored cycling between accretion and underthrusting modes by contrast a sudden syntectonic sediment load in the prowedge region promoted prolonged phases of underthrusting retarding accretion of new imbricates at the prowedge toe however at later stages of deformation the prowedge was forced to regain its characteristic minimum critical taper as predicted by theory and did so by the sudden nucleation of long flat thrust units rapidly rebalancing the asymmetry between prowedge and retrowedge regions the experiments suggest that given a fixed critical taper related to time invariant frictional properties along the wedge bounding surfaces deformation does vary within the wedge according to time varying location of normal stress surpluses and unbalanced topography acting on potential failure surfaces these in turn alter the equilibrium between prowedge and retrowedge basal detachments that is here considered to be a major factor controlling the self regulating dynamics of collisional orogens

the western mediterranean extensional basins and the alpine orogen the western mediterranean late oligocene miocene basins alboran valencia and provencal basins are a coherent system of interrelated troughs in all basins normal faults and thermal subsidence migrated toward the east progressively moving to the miocene to pleistocene algerian and tyrrhenian basins all those troughs appear elements of the back arc opening related to the eastward roll back of the w directed apennines maghrebides subduction zone similarly to western pacific back arc settings these late oligocene early miocene basins nucleated both within the betic cordillera e g alboran sea and in its foreland valencia and provencal troughs the n40 70 direction of grabens is oblique to the coexisting n60 80 trending orogen and shows its structural independence from the orogenic roots thus as the extension cross cuts the orogen and developed also well outside the thrust belt front the westernmost basins of the mediterranean had to develop independently from the alps betics orogen therefore the alboran extension considered a classic example of a basin generated by the collapse of an orogen cannot be ascribed to the detachment or annihilation of the lithospheric root in contrast with the eastward migrating extensional basins the betic balearic thrust front was migrating westward producing interference or inversion structures

ablative subduction a two sided alternative to the conventional subduction model subduction is conventionally modelled as a convergence and overlapping of two semirigid plates such a convergence is one sided in that one lithospheric plate descends while the other remains on the surface we suggest an alternative model one based on the dynamics of fluids in this model the viscous lower lithosphere flows downward and the brittle upper lithosphere deforms in passive response this process is potentially double sided for we find that even a buoyant plate can be dragged downward by a dense descending neighbor thus an apparent overriding plate may be worn away by a process of viscous ablation with the rate of ablation a function of plate buoyancy we call this process ablative subduction ablative subduction allows us to simply interpret observations concerning slab profiles interplate seismicity back arc tectonics and complex processes such as double subduction and subduction polarity reversal in performing experiments modelling the evolution of simple fluid slabs we find that slab profile is strongly influenced by ablation in the overriding plate when ablation is weak as when a buoyant continent borders the trench deformable slabs adopt shallow andean style profiles these profiles develop over time from an initially steep shape more vigorous ablation rates as might occur in an ocean ocean convergence yield steeper marianas style profiles thus differences between marianas style and andean style slabs may result from differences in ablation rate and subduction duration the occurrence of ablation might not be easily detected at depth because material from subducting and ablating plates adhere closely as a single slab these slabs show downdip deformations that are consistent with seismic focal mechanisms plate deformations associated with ablative subduction are also consistent with observed patterns of seismicity in particular the relative aseismicity of marianas style plate boundaries is consistent with the ablative model we may thus explain how plates can converge aseismically without requiring that these plates be decoupled ablation should be weakest in andean style subduction although pulses of rapid ablation might lead to episodes of crustal shortening in oceanic regions the manifestation of ablation may depend upon tectonic regime in an extensional area like the marianas ablation may be temporally self limiting and ablative cycles may explain conflicting observations of tectonic erosion and accretion in a compressional regime vigorous ablation might lead to observations of double subduction and subduction reversal in each case the behavior of the lower lithosphere is relatively simple what varies is the rate of ablation and the response of the upper lithosphere

neogene patterns of relative plate motion for africa europe some implications for recent central mediterranean tectonics a detailed relative motion picture for the neogene africa europe plate kinematics is presented the kinematic reconstruction was carried out using the finite difference solution between the rotation parameters determined for anomalies 7 to 2 in the africa north america europe plate motion circuit the analysis shows a motion of africa with respect to europe which is nne directed during late oligocene to burdigalian times becoming nnw trending from the langhian to the early tortonian from upper tortonian times onward the motion changes to a clear north west directed convergence major late neogene tectonic features of the central mediterranean region can to a large extent be explained within the context of the reconstructed major plate motions late tortonian to recent africa europe slip vectors are compatible with a variety of geological phenomenoa such as north west directed subduction beneath calabria south east translation of calabria and extension in the tyrrhenian sea north west trending slip vectors from thrust earthquakes between gibraltar and sicily and dextral strike slip across the north african margin

dating of deformation phases using k ar and 40ar 39ar techniques results from the northern apennines abstract paleozoic to oligocene metasedimentary rocks present in the alpi apuane region of the northern apennines italy have been sequentially deformed during a tertiary progressive deformation in an attempt to date the individual deformation episodes over 50 conventional k ar and 1140 ar 39ar incremental gas release analyses have been carried out on fine grained white micas separated from samples whose structural settings were well known mineralogy x ray diffractometry and thin section analyses indicate that the constituent muscovite and phengite formed under metamorphic conditions of 3 4 kbars and 300 400 c during all deformational phases pre existing micas were variably crenulated during each subsequent deformational phase both k ar and 40ar 39ar analyses were carried out on 0 6 2 m 2 6 m and 6 20 m size separates of the phengitic white mica although the k ar apparent ages range from 11 to 27 ma and are consistent with available stratigraphic constraints the 40ar 39ar age spectra display variable internal discordancy these isotopic data indicate that 1 both the k ar and 40ar 39ar total gas ages decrease as the degree of crenulation increases 2 the k ar and 40ar 39ar total gas ages decrease as grain size decreases 3 for each sample characteristics of the 40ar 39ar age spectra depend upon grain size with fine sizes yielding discordant patterns which systematically increase in apparent age from low to high temperature and 4 phengitic micas associated with earliest structures yield generally older ages than micas associated with later structures the isotopic results are interpreted to indicate that the major deformation phase d1 occurred at approximately 27 ma with subsequent pulses ending by c 10 ma these results may be combined with finite strain data to suggest that the region was deformed at strain rates between 10 15 and 10 14 s 1 a 27 ma age indicates mid oligocene initiation of plate tectonic activity in the western mediterranean and concomitant deformation in the northern apennines

orogeny and orography the effects of erosion on the structure of mountain belts a numerical model of the coupled processes of tectonic deformation and surface erosion in convergent orogens is developed to investigate the nature of the interaction between these processes crustal deformation is calculated by a two dimensional finite element model of deformation in response to subduction and accretion of continental crust erosion operates on the uplifted surface of this model through fluvial incision which is taken to be proportional to stream power the relative importance of the tectonic and erosion processes is given by a dimensionless erosion number relating convergence velocity rock erodibility and precipitation rate this number determines the time required for a system to reach steady state and the final topographic shape and size of a mountain belt fundamental characteristics of the model orogens include asymmetric topography with shallower slopes facing the subducting plate and an asymmetric pattern of exhumation with the deepest levels of exhumation opposite to subduction these characteristics are modified when the regional climate exhibits a dominant wind direction and orographically enhanced precipitation on one side of the mountain belt the two possible cases are dominant wind in the direction of motion of the subducting plate and dominant wind direction in the opposite direction of the subducting plate velocity models of the former case predict a broad zone of exhumation with maximum exhumation in the orogen interior models of the latter case predict a focused zone of exhumation at the margin of the orogen and at high erosion number a reversal in the topographic asymmetry natural examples of these two cases are presented the southern alps of new zealand exhibits the climate and exhumation asymmetry characteristic of wind in the direction opposite to motion of the subducting plate the asymmetry of topography suggests that erosion is not efficient enough to have reversed the topographic asymmetry the contrasting example of dominant wind in the direction of subduction motion is provided by the olympic mountains of washington state in this case exhumation of deep levels of the cascadia accretionary wedge shows a broad domal pattern consistent with the observed orographic precipitation

on the strength of interplate coupling and the rate of back arc convergence in the central andes an analysis of the interseismic velocity field 1 we interpret the interseismic crustal velocity field of the central andes using a simple three plate model in which the andean mountain belt is treated as a rigid microplate located between the nazca and south american soam plates we assume that the euler vectors associated with these plates are strictly coaxial and that the surface velocity field can be decomposed into an elastic loading field driven by locking of the main plate boundary and a smaller contribution associated with back arc convergence we obtain our best fit to the geodetic velocities if the main plate boundary is fully 100 locked between depths of 10 50 km and 8 5 of nazca soam plate convergence is achieved in the back arc by underthrusting of the brazilian shield beneath the subandean zone

contemporary horizontal velocity and strain rate fields of the pacific australian plate boundary zone through new zealand we have inverted velocity solutions from nine geodetic networks distributed across new zealand to derive present day continuous horizontal velocity and strain rate fields at the earth s surface throughout the country the nine networks contain a total of 362 global positioning system gps stations that have been observed at least twice and at least a year apart between 1991 and 1998 the model velocity field is expanded as bicubic spline interpolation functions defined within a curvilinear grid that covers the country and extends into the assumed rigid australian and pacific plates to west and east the inversion jointly minimizes the magnitudes of fitted strain rates and the misfit to the observed velocity data the spline technique allows high spatial resolution of strain rate variations especially in regions with spatially dense gps data expansion of the model velocity field on the surface of a sphere allows arbitrarily large areas to be studied previously known aspects of new zealand plate boundary deformation are highlighted with improved resolution including back arc extension in the taupo volcanic zone rotation of the hikurangi forearc away from this zone in the north island and a band of high shear strain rate under the southern alps to the southeast of the alpine fault our results reveal several new features including a region of enhanced shear straining apparently associated with strike slip faults in the southern north island and a band of contractional straining subparallel and well east of the alpine fault that is similar to features found in numerical and sandbox models of continental collision

seismic image of the cantabrian mountains in the western extension of the pyrenees from integrated escin reflection and refraction data abstract integrated analysis of normal incidence and large aperture seismic reflection data collected in 1992 and 1993 within the spanish escin and complementary projects provide a first complete n s crustal transect across the northern iberian peninsula and continental margin images of the crustal structure of the cantabrian mountains and their transition to the duero basin and to the cantabrian margin are obtained from a a 65 km long vertical reflection profile escin 2 on land b a 200 km long reversed refraction profile and c wide angle recordings of the marine escin 4 profile consistent results between reflectivity pattern and velocity depth distribution reveal important lateral variations in the deep structure the reflective crust imaged in the escin 2 profile changes its attitude from sub horizontal beneath the duero basin to north dipping beyond the mountain front basement thrusts are observed in the upper crust merging into a detachment at 6 s twt and may have triggered the alpine uplift of the range the moho is identified at the bottom of the reflective lower crust and deepens from 12 to 15 s at the northern end of the profile about 35 km inland modelling of the refraction data laterally extends the seismic image and provides evidence for variscan crustal features beneath the duero basin northwards the velocity in the lower crust decreases and the moho constrained by the wide angle data from profile escin 2 deepens to about 60 km ending abruptly at the shoreline the velocity depth model is constrained along the asturian platform up to the continental slope where the crust mantle boundary is located at 24 km depth this margin moho shows a progressive deepening southwards and extends to the coast where it is found at 30 km depth the present seismic data support an important alpine reworking and thickening of the crust under the cantabrian mountains the onshore offshore transition is marked by an imbrication of two crusts of very different thicknesses this signature offers a strong parallelism with the one previously observed further east across the pyrenees in the ecors seismic profile

petrochemistry and tectonic significance of cretaceous island arc rocks cordillera oriental dominican republic abstract cretaceous island arc rocks of the caribbean island arc system have been exposed by cenozoic faulting in the cordillera oriental in eastern hispaniola high k 2 o intermediate to felsic volcanic rocks loma la vega volcanics are interbedded with marine epiclastic sedimentary rocks and tuffs las guajabas tuffs and unconformably overlie pre aptian sedimentary rocks low k 2 o volcanic rocks guamira volcanics and a granodioritic to tonalitic intrusion el valle pluton the petrology and geochemistry of these units in conjunction with regional stratigraphic data are used to speculate on the tectonics of the newly developing caribbean island arc system during early and late cretaceous time the loma la vega volcanics are characterized by the presence of large phenocrysts of sanidine and minor amounts of clinopyroxene opaque oxides and rare leucite in a devitrified matrix of chlorite and clay although the volcanic rocks have undergone some low temperature alteration metamorphism which redistributed some major elements and large ion litho phile trace elements the high field strength elements rare earth elements and radiogenic isotopes appear to have been minimally affected based on abundances of the relatively immobile elements trace element enrichment patterns and isotopic compositions the loma la vega volcanics are considered part of the high k calc alkaline ca or shoshonitic island arc volcanic series in contrast pre aptian early cretaceous volcanic and plutonic rocks of the underlying los ranchos formation have chemical characteristics similar to rocks in the island arc tholeiitic or primitive island arc pia series that form coeval and along strike sections of the early cretaceous caribbean island arc in other parts of present day hispaniola cuba puerto rico and the virgin islands an abrupt and regional change in composition from island arc tholeiites to high k calc alkaline rocks is coincident with a hypothesized reversal in subduction polarity in pre aptian time as inferred from previously published tectonic models polarity reversal may have been triggered by attempted subduction of the caribbean oceanic plateau beneath this segment of the caribbean island arc the observed magmatic and tectonic effects of the proposed cretaceous caribbean arc reversal are similar to the better documented neogene subduction reversal event in the solomon islands arc in the southwest pacific

oblique plate convergence and interseismic strain accumulation 1 a simple half space elastic model is used to investigate interseismic strain accumulation above a subduction zone associated with oblique plate convergence the horizontal surface velocity field near the leading edge of the upper plate is found to rotate as distance from the trench increases above the locked portion of the plate interface surface velocity is more oblique than the plate convergence vector farther back toward the interior of the upper plate the surface velocity is less oblique than the plate convergence vector surface velocity fields emphasize the strike slip component of plate convergence near the locked zone and emphasize the dip slip component farther back this interseismic strain partitioning is purely elastic interseismic straining will be reversed during the eventual rupture of the locked zone

sandbox model simulation of forearc evolution and noncritical wedges convergent plate boundaries exhibit a wide range of behaviors ranging from sediment accretion and underplating to large scale erosion of the forearc we have used sandbox modeling to simulate the growth of thin skinned contractional wedges over a strong backstop at an accretionary margin the backstop in our model causes the growth of two distinct wedges with quite different growth histories analogous to the accretionary wedge and the inner deformation belt retrowedge found in at least some forearcs after starting with an initially steep slope the surface flattens as the accretionary wedge grows by frontal accretion until the wedge attains and maintains a minimum possible critical taper the retrowedge driven by uplift of the outer arc high increases its surface slope as it grows until a kinematic steady state is attained it exhibits a stable taper with slope larger than the minimum critical taper but never larger than the theoretical maximum critical taper when a basal decollement is present above the backstop in models in which the retrowedge has not developed a decollement it has a taper with a slope at the angle of repose when a basal decollement is present above the backstop the main body of the retrowedge acts as a bulldozer or a backstop pushing the sediments in front of it which then form an additional wedge with a small minimum critical taper a broad outer arc high stands above the toe of the backstop where maximum principal stresses rotate to meet the stress conditions in the two contractional belts on either side of it as these belts grow in opposite directions the accretionary wedge accretes large quantities of sediment at its front but the retrowedge does not because the retrowedge can exhibit a non critical stable taper in the absence of rapid erosion or accretion it is misleading to assume that active thin skinned wedges are necessarily at maximum or minimum critical taper

new results from physical modelling of arc continent collision in taiwan evolutionary model abstract both the analysis of the geological data and the results from 2d and 3d experimental modelling of arc continent collision in taiwan reveal the following major stages in this process 1 failure and thickening of the subducting crust of the chinese continental margin under the frontal part accretionary prism of the philippine sea plate psp since ca 6 7 ma 2 uplift squeezing and exhumation of the subducted crustal slice of the chinese margin and initiation of the central range about 5 ma ago 3 flexural buckling subsidence of the luzon arc fore arc under strong horizontal compression and sedimentary infilling of this area by the material shed from the growing and exhuming central range 4 failure of the psp along the west vergent longitudinal valley fault dipping beneath the arc and closure of the fore arc 5 subduction reversal and stoppage of continental subduction which has already occurred in northern taiwan and propagates to the south the southward propagation of the collision is not stationary currently occurring initial collision in southern taiwan and offshore to the south does not follow the same described above evolution as in central and northern taiwan the evolution of the collision in time and space should be considered in the framework of a 3d model where the psp rupture is preparing in the central and northern taiwan during a long time passing through a number of evolutionary stages once it has occurred the rupture longitudinal valley fault rapidly propagates to the south into the domain which did not undergo the pre rupture evolution the rupture front is currently located near 21 30 n

oblique plate convergence slip vectors and forearc deformation slip vectors from thrust earthquakes at subduction zones where convergence is oblique to the trench often point between the directions of relative plate convergence and normal to the trench axis suggesting that oblique convergence is taken up by partial decoupling decoupling means that a component of arc parallel motion of the leading edge of the upper plate results in less oblique thrusting at the trench partial decoupling is modeled by partitioning of oblique convergence into slip on thrust and strike slip faults that are parallel to the trench and to each other and starting with a force equilibrium condition a relationship between the obliquity and the earthquake slip vector orientation is derived assuming that either fault slips when shear stress on it reaches a yield stress oblique slip parallel to the plate vector should occur on the thrust fault when obliquity is smaller than a critical angle for obliquity at or greater than this angle the stress on the strike slip fault is large enough to start it slipping and when both faults are active the arc parallel motion of the forearc deflects the slip vector back toward the trench normal if we assume that continued slip on either fault occurs at constant stress but the two faults can be at different stresses the slip vector will maintain a constant angle relative to the trench normal even when obliquity is larger than the critical angle this limiting angle of the slip vector called max measured relative to the trench normal is simply the arcsine of the ratio of the shear forces resisting slip on the strike slip and thrust faults a consequence is that when the obliquity exceeds max the slip vectors on the thrust fault are sensitive only to the thrust fault orientation and contain no information about the convergence direction between the plates slip vectors at the java trench southwest of sumatra show the relationship clearly with max 20 5 while slip vectors at the aleutian trench show the relationship less clearly with max 25 to 45 the greater angle at the aleutian trench suggests that the upper plate is stronger in the aleutian arc relative to the thrust fault than in the sumatran arc consistent with the sumatran arc being continental and having a well developed strike slip fault while the aleutian arc is oceanic and without a clear transcurrent fault slip vectors at the philippine trench which like sumatra has a large strike slip fault inboard of it tend to stay within 25 of the trench normal when obliquity is as large as 50 if obliquity exceeds max and continues to increase along a subduction zone the rate of motion of the forearc relative to the upper plate will vary with obliquity in which case the forearc sliver should extend or contract parallel to the arc from the geometry of modern island arcs arc parallel extension should be the more common and has been hypothesized for both sumatra and the aleutians on the basis of earthquake slip vectors and for these and other arcs from geological observations from estimates of max and the arc parallel gradients in obliquity arc parallel strain rates are estimated to be 1 to 3 10 8 yr for the sumatran forearc 2 to 6 10 8 yr for the aleutian forearc and 0 3 to 3 10 8 yr for the philippine forearc oblique convergence and subsequent arc parallel extension if accompanied by crustal thinning may provide an important yet little appreciated mechanism for bringing high grade metamorphic rocks to the surface of subduction complexes

the determination of glacier speed by time lapse photography under unfavorable conditions two practical problems in the use of time lapse photography for the measurement of speed were encountered during the recent surge of west fork glacier in the central alaska range alaska u s a the first is severe rotational camera instability we show how natural unsurveyed features of the valley wall can be used to make the necessary corrections the second problem is the computation of absolute speed when many different unsurveyed glacier surface features are used as targets

an elastoplastic model of frost deformation for the porous rock under freeze thaw abstract freeze thaw action has caused many engineering geology disasters in cold regions such as frost cracking of tunnel and pavement rockfall of rock slope falling off of building materials the freeze thaw damage of porous rock is mainly induced by the frost heaving pressure due to pore water freezing although some elastic crystallization theories have been proposed to model the freezing process of pore water the long term frost heave characteristics under cyclic freeze thaw are rarely studied the ice crystallization and frost heave theory under freeze thaw have been developed in the elastoplastic regime in this study a novel elastoplastic model has been proposed to estimate the long term frost heave and frost damage under freeze thaw this model has considered the influence of the pore size yield stress of matrix and degree of water saturation on the development of the frost heaving pressure in pores during freezing a plastic region will produce inside the matrix around the pore this plastic region gradually expands outwards with the increase of this pressure during thawing the elastic deformation will be recovered but a considerable residual plastic strain is accumulated however the growth rate of the residual plastic strain will decrease under freeze thaw due to the production of new empty voids therefore the freezing of pore water in porous rock is an elastoplastic problem besides the water saturation and mechanical properties of the porous rock are very important to quantify the frost heave this study provides a better understanding of the frost damage of the porous rock and a reference for the stability of rock engineering under freeze thaw in cold regions

fault development and interaction in distributed strike slip shear zones an experimental approach analogue model experiments using both brittle and viscous materials were performed to investigate the development and interaction of strike slip faults in zones of distributed shear deformation at low strain bulk dextral shear deformation of an initial rectangular model is dominantly accommodated by left stepping en echelon strike slip faults riedel shears r that form in response to the regional bulk stress field push up zones form in the area of interaction between adjacent left stepping riedel shears in cross sections faults bounding push up zones have an arcuate shape or merge at depth adjacent left stepping r shears merge by sideways propagation or link by short synthetic shears that strike subparallel to the bulk shear direction coalescence of en echelon r shears results in major through going faults zones master faults several parallel master faults develop due to the distributed nature of deformation spacing between master faults is related to the thickness of the brittle layers overlying the basal viscous layer master faults control to a large extent the subsequent fault pattern with increasing strain relatively short antithetic and synthetic faults develop mostly between old but still active master faults the orientation and evolution of the new faults indicate local modifications of the stress field in experiments lacking lateral borders closely spaced parallel antithetic faults cross faults define blocks that undergo clockwise rotation about a vertical axis with continuing deformation fault development and fault interaction at different stages of shear strain in our models show similarities with natural examples that have undergone distributed shear

active detachment of taiwan illuminated by small earthquakes and its control of first order topography we use 110 000 small earthquakes to locate and map active faults in three dimensions within the taiwan arc continent collision the structure is dominated by a nearly horizontal band of small earthquakes at 10 km depth that is interpreted to seismically illuminate the main detachment zone of the mountain belt with other illuminated fault zones abutting the detachment zone the zone steepens below eastern taiwan to 30 90 and reaches depths of 30 60 km the three dimensional shape of the detachment zone in relation to topography allows a new test of critical taper wedge mechanics and suggests that the reversal of topographic slope across taiwan is controlled by the shape of the detachment

deformation partitioning in flat subduction setting case of the andean foreland of western argentina 28 s 33 s 1 the andean foreland of western argentina 28 s 33 s corresponds to back arc deformations associated with the ongoing flat subduction of the nazca plate beneath the south american lithosphere in this paper short and long term tectonic regimes for this andean region are investigated and compared to nuvel 1 and the new gps derived velocity field for the andes mountains from this comparison it appears that 1 the slight convergence obliquity at the plate boundary is accommodated by the subduction zone itself precluding any deformation partitioning at lithospheric scale 2 clockwise rotation of 1 and p axes within the andean foreland of western argentina suggests that plio quaternary deformation partitioning may occur in the upper plate at crustal scale and 3 comparison of seismic rates of shallow deformations determined from moderate to large earthquakes localized between 22 s and 36 s shows that the amount of shortening is about the same 2 4 cm yr 1 between the forearcs of the 30 dipping slab segments and the back arc of the flat segment all together those observations are in close agreement with the idea that interplate coupling in flat subduction settings may play a key role in the present day transfer of plate boundary forces consequently the andean back arc of western argentina might be regarded as an obliquely converging foreland where plio quaternary deformations are partitioned between strike slip and thrust motions that are localized in the east verging thin skinned argentine precordillera and west verging thick skinned pampean ranges respectively

mechanical model for the tectonics of doubly vergent compressional orogens a mechanical model of crustal shortening and deformation driven by the relative convergence of rigid underlying mantle plates explains many features of convergent orogens results based on numerical models and supported by sandbox models show that a coulomb crustal layer subject to basal velocity boundary conditions corresponding to asymmetric detachment and subduction of the underlying mantle passes through three stages of orogenic growth 1 block uplift bounded by step up shear zones 2 development of a low taper wedge over the underthrusting mantle plate and 3 development of a low taperwedge overlying the overthrusting mantle plate and verging in the opposite direction when modified by isostasy basal viscous flow surface erosion and denudation and sedimentation the resultant model orogens exhibit a variety of styles with characteristics in common with small rapidly denuded orogens large orogens with plateaus and extensional characteristics and active subduction margins with doubly vergent accretionary wedges and deformed fore arc basins

three dimensional critical wedges tectonics and topography in oblique collisional orogens in orogens deforming in oblique compression both the long wavelength mean topography of averaged slopes and the higher frequency topography first harmonic of valleys and ridges are related to the mechanics of deformation in three dimensions and cannot be adequately modeled in two dimensions a method of analysis of the principal stresses failure planes and topography is presented for three dimensional critical wedges in a critical coulomb orogen which is at failure throughout the components of the stress tensor are interrelated through the failure condition defined in terms of the principal stresses variation of any of the cartesian components must be accommodated by adjustment of one or more of the other components so that the orogen remains at failure this coupling of cartesian components to maintain the wedge at failure permits exploration of the behavior of a three dimensional critical orogen in which shear stresses arise outside of the plane of convergence internal adjustment of shear stresses to externally imposed stress conditions results in the formation of characteristic orogen shapes for an oblique collision zone such as the southern alps of new zealand the resultant mountain belt is one of an obliquely deforming outboard fold thrust belt with subdued relief rising abruptly into a steep backslope to the main divide which falls steeply to the indentor forming the inboard slope the outboard topographic slope is reduced in relation to that of the two dimensional orogen the rigid indentor has a fixed strike and therefore imposes a failure orientation adjacent to the indentor that differs from that of the ideal three dimensional critical wedge solution this produces several different slip orientations and minor partitioning of strain near the indentor first harmonic topography of valleys and ridges arises due to nonlinear heterogeneous behavior of both erosional and mechanical processes and is characterized by uplift controlled ridges in the outboard and incision controlled valleys in the inboard

role of preexisting topography and overburden on strain partitioning of oblique doubly vergent convergent wedges 1 scaled sandbox modeling was used to analyze the three dimensional aspects of strain partitioning along obliquely convergent margins that underwent a subduction polarity reversal event in the experiments the first phase of shortening produced a preexisting topography that affected wedge kinematics during the second phase increasing degrees of obliquity during the first shortening phase produced different amounts of inherited topography and overburden giving way to different effects on wedge response during the second phase of shortening during the second phase the models with an inherited heavy overburden show an anomalous orientation for riedel shears predating the development of a margin parallel strike slip fault these shears are oriented at an angle that is higher than kinematically predicted with respect to the subduction slot suggesting a transient rotation of the stress applied by the obliquity of subduction in contrast experiments with a lighter overburden do not show this inconsistency in addition all the models displayed an unexpectedly high accretion rate during the second phase of shortening which we interpret to be dependent on topographical slope breaks inherited after the first phase of shortening and the critical asymmetrical architecture of doubly vergent wedges

implications of shortening in the himalayan fold thrust belt for uplift of the tibetan plateau 1 recent research in the himalayan fold thrust belt provides two new sets of observations that are crucial to understanding the evolution of the himalayan tibetan orogenic system first u pb zircon ages and sm nd isotopic studies demonstrate that 1 greater himalayan medium to high grade metasedimentary rocks are much younger than true indian cratonic basement and 2 these rocks were tectonically mobilized and consolidated with the northern margin of gondwana during early paleozoic orogenic activity these observations require that greater himalayan rocks be treated as supracrustal material in restorations of the himalayan fold thrust belt rather than as indian cratonic basement in turn this implies the existence of greater himalayan lower crust that is not exposed anywhere in the fold thrust belt second a regional compilation of shortening estimates along the himalayan arc from pakistan to sikkim reveals that 1 total minimum shortening in the fold thrust belt is up to 670 km 2 total shortening is greatest in western nepal and northern india near the apex of the himalayan salient and 3 the amount of himalayan shortening is equal to the present width of the tibetan plateau measured in an arc normal direction north of the indus yalu suture zone this information suggests that a slab of greater indian lower crust composed of both indian cratonic lower crust and greater himalayan lower crust with a north south length of 700 km may have been inserted beneath the tibetan crust during the cenozoic orogeny we present a modified version of the crustal underthrusting model for himalayan tibetan orogenesis that integrates surface geological data recent results of mantle tomographic studies and broadband seismic studies of the crust and upper mantle beneath the tibetan plateau previous studies have shown that incremental mesozoic and early cenozoic shortening had probably thickened tibetan crust to 45 55 km before the onset of the main cenozoic orogenic event thus the insertion of a slab of greater indian lower crust with maximum thickness of 20 km tapering northward could explain the cenozoic uplift of the tibetan plateau the need for tibetan crust to stretch laterally as the greater indian lower crust was inserted may explain the widespread east west extension in the southern half of the plateau our reconstruction of the himalayan fold thrust belt suggests that indian cratonic lower crust of presumed mafic composition and high strength should extend approximately halfway across the tibetan plateau to the banggong suture from there northward we predict that the tibetan plateau is underlain by more felsic and therefore weaker lower crust of greater himalayan affinity two slab break off events are predicted by the model the first involved neotethyan oceanic lithosphere that foundered 45 35 ma and the second consisted of greater indian lithosphere most likely composed of greater himalayan material that delaminated and foundered 20 10 ma asthenospheric upwelling associated with the break off events may explain patterns of cenozoic volcanism on the tibetan plateau although the model predicts a northward migrating topographic front due solely to insertion of greater indian lower crust the actual uplift history of the plateau was complicated by early middle tertiary shortening of tibetan crust

experimental insights into the partitioning of motion within zones of oblique subduction abstract in zones of oblique convergence the motion is often partitioned into subduction and an associated strike slip component fitch 1972 in this paper we first describe the geodynamic setting of the philippine region as a case study of deformation by strike slip faulting behind parallel subduction zones then we present six small scale laboratory experiments with analogue materials to investigate how partitioning may be controlled by internal boundaries between oceanic lithosphere and continental or island arc lithosphere during the early stages of a typical experiment deformation is taken up on oblique thrusts subduction zones which root into internal boundaries paired strike slip faults develop at later stages thereafter subduction zones and paired strike slip faults operate simultaneously the partitioning occurs without any need for complex constraints on external boundaries

3d evolution of fold and thrust belts formed by oblique convergence abstract this paper presents the results of a series of scaled analogue models designed to simulate thrust belt development formed by orthogonal and oblique convergence the models were constructed to produce doubly vergent orogenic wedges with a pro wedge and a retro wedge separated by an uplifted orogenic core models were constructed from homogeneous layered sandpacks in a 3 m 1 3 m deformation rig the progressive evolution of the models was recorded using digital photography and the photographs were animated for analysis the experiments investigated convergence obliquities from orthogonal 90 to as low as 15 obliquity orthogonal models produced long linear critically tapered pro wedge fold and thrust belts parallel to the convergence margin together with a narrow uplifted core and a steep retro wedge thrust system the taper angle of the pro wedge was typically 11 12 whereas the retro wedge maintained a steeper taper of 38 42 models where the convergence vector was oblique at angles of 60 and 45 to the margin produced doubly vergent thrust wedges with thrust faults trending parallel to the margin there was little evidence of linked penetrative strike slip faulting but en echelon riedel like shear systems developed in the 45 oblique model subordinate oblique slip motion was observed on some thrust faults in contrast for 30 and 15 oblique convergence models strong deformation partitioning developed with steeply sloping pro wedges that did not have critical tapers through going linked strike slip faulting was developed that cut the uplifted axial zones of these models analysis of the models by animation of the digital photographs shows how the thrust systems initiated and propagated in particular it is clear that at any one time several thrusts were moving simultaneously serial cross sections through the moderately oblique convergence models 60 45 are almost indistinguishable from the orthogonal model cross sections the results of orthogonal convergence experiments compare well to the map patterns of faults and folds in the salt range pakistan 60 and 45 oblique convergence models compare well to the zagros where some strain partitioning on strike slip faults occurs the highly oblique models show very strong strain partitioning with penetrative margin parallel strike slip faulting similar to that found in ne venezuela

crustal extension in the northern apennines the transition from compression to extension in the alpi apuane core complex an interpretation of northern apennine geology is presented which relates the temporal and spatial occurrence of both compressional and extensional deformation features in terms of the changing dynamic evolution within an accretionary wedge after a model proposed by platt 1986 followed by the initiation and development of continental rifting during cretaceous to eocene time an accretionary wedge formed as remnant tethyan oceanic crust subducted beneath the rotating corsica sardinia microplate microplate collision during the oligocene was characterized by the rapid imbrication of buoyant continental crust of the italian continental margin the record of which is preserved within the duplex structure geometry of the alpi apuane region the overthickened wedge geometry returned to a more stable configuration by developing extensional features during the miocene both listric normal faults at upper crustal levels and shear zones indicating evidence of distributed ductile extensional strain at mid crustal levels are recorded it is proposed that large scale regional extension with associated volcanism beginning in the messinian was represented by the intrusion of asthenospheric material from the subducted plate into the already attenuated accretionary complex further rifting perhaps aided by subduction and back arc processes in the southern apennines led to the formation of the tyrrhenian sea as an oceanic basin both the apennines and north american core complexes record evidence of crustal thickening followed by crustal thinning and finally of continental rifting this suggests that the similar histories of these regions with vastly different plate tectonic settings may both be explained by processes linked to the changing internal dynamics of accretionary wedges

computerized x ray tomography analysis of three dimensional fault geometries in basement induced wrench faulting abstract computerized x ray tomography applied to sandbox experiments makes it possible to analyze the kinematic evolution as well as the three dimensional geometry of faults in basement controlled wrench faulting with increasing basement displacement the evolution of shear bands developed in the model sand bed is observed as follows 1 the right stepping shears that have a cirque or shell or ship body shape develop on both sides of the basement fault 2 the shears on the one side of the basement fault join those on the other side and therefore each shear surface forms a helicoidal shape shears reach the surface of the sand near or above the basement fault thereby en echelon riedel shears are observed at the surface of the sand 3 lower angle shears and p shears lying within the zone defined by the first riedels have steeper dips and a less pronounced helicoidal aspect than the riedels such three dimensional shape and kinematic evolution of shears in the fault experiments agree well with the similar faults which propagate up through the unconsolidated sediments from a wrench fault in bedrock as observed in the excavated trenches at some active faults in the fault experiments the strike of the riedel shear to the basement fault in dense sand is larger than that in loose sand in reasonable agreement with the directions predicted by the coulomb criterion in terms of the angle of internal friction likewise the geometries of the riedel shears are considered to depend on the overburden materials in natural cases of wrench faulting

first deep seismic reflection images of the eastern alps reveal giant crustal wedges and transcrustal ramps 1 the alps are considered as a classical example for an orogen created by continental plate collision in this study we present new images obtained from deep seismic reflection profiling in the eastern alps between munich and venice which give rise to examine and revise existing concepts the seismic sections exhibit a prominent bi verging reflection pattern at crustal scale a major ramp like structure outcropping at the inn valley fault can be traced southward over 80 km into the mountain root where relics of the subducted penninic ocean are expected new models of the evolution of the eastern alps show an upper lower crustal decoupling along transcrustal thrust faults with opposite thrust directions of both the european and the adriatic african continents

arc continent collision and subduction reversal in the cenozoic evolution of the northwest pacific an example from kamchatka ne russia abstract the principal new feature in the cenozoic evolution of the northwest pacific is based on the fact that two arc terranes of upper cretaceous paleogene age were incorporated in the kamchatka orogen during the early tertiary and at the end of the miocene respectively these two arc terranes so called achaivayam valaginskaya campanian lower paleocene and kronotskaya coniacian eocene arcs are separated by a suture zone that contains deformed slices of young paleocene lower eocene oceanic crust paleomagnetic data on the original paleolatitudes of the arc positions as well as geochemical data characterized the arcs as intraoceanic originated in the northwest pacific the arcs were separated by the vetlovka oceanic plate suggested to be a trapped fragment of the kula plate in which a new spreading zone developed in the early tertiary the principal event in the reconstructed evolution is the eocene collision of the extended achaivayam valaginskaya arc terrane with the continental margin of asia the margin is interpreted as the southeastern facing the pacific cretaceous accretionary boundary of the okhotsk microplate docked against the asian continent at 55 65 ma the arc continent collision developed progressively northward from the end of the paleocene early eocene in southern kamchatka to the middle eocene in the olutorka region the different stages of the collision occurred simultaneously along strike of the achaivayam valaginskaya arc according to the geodynamic model based on geological data for southern kamchatka collision included 1 the entrance of the continental margin in the subduction zone deformation of the overriding plate and subduction of fore arc block 2 arc obduction on the continental margin blocking of continental margin subduction and 3 subduction reversal with onset of oceanic subduction at the eastern side of the deformed arc this subduction reversal produced deformation and reverse faulting along the eastern side of the arc thus allowing accretion of tectonic slices of young oceanic crust of the vetlovka plate which began to subduct during the subduction reversal the subducted plate dipping to the southeast was broken off and detached at depth under the arc area the subduction polarity reversal the most characteristic feature of the collision model is strongly supported by structural and geologic data of southern kamchatka and explains the main tectonic metamorphic and magmatic events related to and postdating the collision

critical wedges in three dimensions analytical expressions from mohr coulomb constrained perturbation analysis we present simple analytical expressions for evaluating three dimensional stress fields in critical wedges coulomb failure throughout that are deforming within an oblique convergent zone we assume that the load and geometry variations in the lateral direction and the topographic slope are relatively small the stress state of a two dimensional critical wedge is perturbed by admitting a small lateral shear which is accommodated to maintain criticality by a reduction in the dominant normal compressive stress cohesionless cases are further simplified with taylor series expansions in the topographic slope the analytical expressions distinguish the behavior of strong base wedges associated with steep slopes from weak base wedges with more modest slopes and are reliable for lateral shear values as large as half the lithostatic load in the weak base case where the vertical shear is small the orientations of the principal stresses are very sensitive to the perturbation as the lateral shear increases through relatively small values the roles of the lateral and vertical coordinate axes are interchanged with a resultant transition from predominantly thrust failure to strike slip failure transition to strike slip in the strong base case is more gradual the particular value of horizontal shear for which transition occurs also depends on the difference of the minor normal stresses we produce an analytical expression predicting displacement patterns compatible with those observed in natural and analog orogens displacement partitioning is predicted for weak base wedges due to minor changes in lateral boundary conditions associated with variable orogen geometry

thrust emplacement of the hispaniola peridotite belt orogenic expression of the mid cretaceous caribbean arc polarity reversal articulo sobre el emplazamiento fallero en la zona de peridotitas de la hispaniola republica dominicana estudiando si es la expresion orogenica de la polaridad reversa en la zona del caribe en el cretacico medio

three dimensional numerical experiments of strain partitioning at oblique plate boundaries implications for contrasting tectonic styles in the southern coast ranges california and central south island new zealand a three dimensional numerical model is used to investigate the nature and distribution of crustal scale deformation and strain partitioning across obliquely convergent transpressive plate boundaries a uniform elastic plastic crustal layer obeying murrell s extension to three dimensions of the griffith failure criterion was deformed under basal kinematic boundary conditions that correspond to the oblique convergence between two mantle lithospheres one of which detaches and subducts beneath the other the results show deformation concentrated in relatively narrow shear zones the character planar versus curved and dip of these shear zones or crustal scale faults depend on the relative amount of convergence versus transcurrent motion imposed across the plate boundary in models dominated by transcurrent boundary conditions two separate sets of coexisting structures develop 1 a pair of planar relatively shallow dipping thrust zones arranged in a v structure that accommodate convergent motion and 2 a set of curved near vertical strike slip zones arranged in a flower structure that accommodates the transcurrent motion in models dominated by convergence both contractional and transcurrent deformations are accommodated by oblique slip across a single pair of structures that have the same geometry as the thrust shear zones there is a gradational behavior between these end member partitioned and nonpartitioned modes for intermediate boundary conditions in all cases the dip of the shear zones predicted by the numerical model is in good agreement with the dip of ideal faults predicted by mohr s incipient faulting hypothesis the two basic modes are interpreted to underlie the contrasting tectonic styles of the southern coast ranges california and central south island new zealand factors other than obliquity of convergence which play a role in strain partitioning are discussed

subduction of asian lithospheric mantle beneath tibet inferred from models of continental collision the relative northward motion of the indian subcontinent that followed the onset of continental collision with asia has produced extensive deformation of the earth s crust giving rise to the world s highest mountains in the himalayan chain and the world s largest high elevation region the tibetan plateau the formation of the broad mountain belt implies that contrary to the original tenets of plate tectonics the lithospheric plates have experienced widespread deformation far from the plate boundary1 several models have been proposed2 6 to explain the manner in which this post collisional deformation is distributed within the continental lithosphere of the indian and asian plates here we propose an alternative model in which subduction of the asian lithospheric mantle develops following the collision of india our model is supported by numerical calculations of crustal deformation and thickening and is consistent with available geological and geophysical data7 9 this picture suggests that lithospheric mantle is not deformed along with the crust and would imply that continental collision zones are more analogous to oceanic subduction zones than was previously believed