wave run up on slender piles in design conditions model tests and design rules for offshore wind wave run up on foundations is a very important factor in the design of entrance platforms for offshore wind turbines when the horns reef 1 wind turbine park in denmark was designed the vertical wave run up phenomenon was not well known in the industry hence not sufficiently considered in the design of horns reef 1 as a consequence damage was observed on the platforms this has been the situation for several sites and design tools for platform loads are lacking as a consequence a physical model test study was initiated at aalborg university to clarify wave run up on cylindrical piles for different values of diameter to water depth ratios d h and different wave heights to water depth ratios h h for both regular and irregular waves a calculation model is calibrated based on stream function theory for crest kinematics and velocity head stagnation theory due to increased velocities close to the pile an empirical factor is included on the velocity head the evaluation of the calculation model shows that an accurate design rule can be established even in breaking wave conditions however calibration of a load model showed that it was necessary to increase the run up factor on the velocity head by 40 to take into account the underestimation of run up for breaking or nearly breaking waves given that they produce thin run up wedges and air entrainment two factors not coped with by the measurement system

new global marine gravity model from cryosat 2 and jason 1 reveals buried tectonic structure gravity models are powerful tools for mapping tectonic structures especially in the deep ocean basins where the topography remains unmapped by ships or is buried by thick sediment we combined new radar altimeter measurements from satellites cryosat 2 and jason 1 with existing data to construct a global marine gravity model that is two times more accurate than previous models we found an extinct spreading ridge in the gulf of mexico a major propagating rift in the south atlantic ocean abyssal hill fabric on slow spreading ridges and thousands of previously uncharted seamounts these discoveries allow us to understand regional tectonic processes and highlight the importance of satellite derived gravity models as one of the primary tools for the investigation of remote ocean basins

geodetic constraints on postseismic deformation following the 1990 ms 7 8 luzon earthquake and implications for luzon tectonics and philippine sea plate motion 1 we have obtained three epochs of high quality dual frequency gps data at 15 stations in the luzon philippines region since the july 16 1990 ms 7 8 earthquake from which we generate horizontal velocity fields at the earth s surface for the 1993 1996 and 1996 1998 periods we use these velocities plus 1996 1998 gps velocity results reported elsewhere from 13 additional stations to show that present day deformation in luzon is dominated by strike slip motion along the philippine fault system the measured strike slip rate across central luzon is faster than the expected long term slip rate on the fault which we attribute to postseismic deformation following the 1990 earthquake the 1993 1996 and 1996 1998 deformation patterns cannot be statistically distinguished within the 10 uncertainties of the data if the observed deformation is interpreted using an elastic half space model with uniform slip below a subsurface locking depth we find a fault slip rate of 40 mm yr and a locking depth of 15 km this is too fast to be a steady interseismic rate as it disagrees with what is known from paleoseismic and other evidence about earthquake recurrence and long term slip rate on the philippine fault system in central and northern luzon also a 15 km locking depth for steady interseismic slip is probably too shallow compared to the depth of rupture in the 1990 earthquake if the deformation is instead interpreted with two types of viscoelastic model we find that the observed velocities can be fit well for a range of values of lower lithosphere viscosity and long term slip rate using a two dimensional 2 d viscoelastic coupling model the minimum allowable lower lithosphere viscosity is 0 5 1019 pa s with an associated long term slip rate of 15 22 mm yr faster long term slip rates also fit well with correspondingly higher values of viscosity with a possibly more realistic model that includes 3 d viscoelastic effects we find preferred long term velocities in the range 20 35 mm yr and viscosities in the range 2 6 1019 pa s our preferred viscosities correspond to maxwell times on the order of 7 20 of the typical earthquake recurrence interval and are several times higher than values found for a number of californian strike slip earthquakes other implications of our 1993 1998 surface velocities are that the philippine sea eurasia ph eu euler vector of seno et al 1993 is a better description of the ph eu convergence rate than a euler vector recently estimated from gps velocities in the northern ph and that the great majority of the normal component of ph eu plate convergence is taking place west of luzon presumably largely at the manila trench we also find that luzon is rotating counterclockwise relative to ph at 1 2 myr only 25 of the rate recently suggested by other workers on the basis of limited gps data

a new branch of the philippine fault system as observed from aeromagnetic and seismic data abstract the philippine fault is one of the major strike slip faults of the world extending over 1200 km from the lingayen la union region of northwestern luzon to the davao gulf south of mindanao nevertheless its total slip history of displacement and its exact location have been uncertain and even controversial significant new insight has been provided by aeromagnetic and seismic surveys made as part of a basin evaluation project of the philippines bureau of energy development funded by the world bank this paper outlines evidence from these surveys for a previously unrecognized major branch of the philippine fault system and tentatively discusses its paleogeographic implications we have interpreted aeromagnetic data covering most of the length of the fault and seismic profiles crossing it at various locations the philippine fault is well displayed in the aeromagnetic data surprisingly the data also show a major previously unrecognized branch of the philippine fault here named the sibuyan sea branch this branch is almost entirely offshore or is covered by quaternary or recent volcanic rocks it leaves the known philippine fault northeast of masbate island passes along the northeast edge of the sibuyan sea southwest of the bondoc peninsula passes northeast of marinduque island apparently passes through the taal volcano and then offshore west of the batangas peninsula the sibuyan sea branch is 350 km long and is presently active it is more obvious on seismic and aeromagnetic data than the northern part of the known philippine fault and thus may have accommodated more motion we have made a preliminary evaluation of the total slip on the sibuyan sea branch by attempting to restore long wavelength geophysical and geologic features including forearc basins ophiolite belts volcanic belts and regional bouguer free air and magnetic anomalies which suggests the possibility of 200 300 km of slip since about the early or middle miocene in this reconstruction the central luzon basin and luzon volcanic arc are continuous with the forearc basin and volcanic belt on panay as a continuous arc and forearc extending from northern ilocos to southern panay after collision with the palawan mindoro continental block continued motion along the manila trench north of mindoro has ruptured the central philippines along the sibuyan sea branch of the philippine fault thus the sibuyan sea branch could account for most of the post early to middle miocene relative motion between luzon and the visayas

damped regional scale stress inversions methodology and examples for southern california and the coalinga aftershock sequence 1 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

earthquake induced crustal gravitational potential energy change in the philippine area abstract the crustal gravitational potential energy change gpe caused by earthquakes in the philippine area from january 1976 to november 2011 was estimated in this study the active convergence between the philippine sea plate and the sundaland eurasian margin is reflected by the greatest gains in gpe along the philippine negros and cotabato trenches whereas the manila trench is covered by a gpe loss pattern although the philippine mobile belt pmb itself is actually affected by the ongoing collision and subduction processes almost the entire philippine fault zone is dominated by gpe loss revealing a slightly extensional environment along the fault the time evolution of the cumulated gpe for different segments along the philippine archipelago shows distinct patterns due to the numerous large underthrusting events that have occurred along the philippine trench the cumulated gpe is regularly increasing in its most southern segment however in the middle segments where the palawan block enters into collision with the pmb the increase in cumulated gpe is relatively small in the most northern segment where the north luzon is located a decrease of cumulated gpe demonstrates that the seismic characteristic of the manila trench is dissimilar from other subduction systems in the world we suggest that the collision of both the palawan block and the benham rise with the pmb promotes the rotation of the pmb and facilitates the northward escape of the northeastern luzon resulting in a decrease of cumulated gpe in the northern philippines

plate convergence measured by gps across the sundaland philippine sea plate deformed boundary the philippines and eastern indonesia summary the western boundary of the philippine sea ph plate in the philippines and eastern indonesia corresponds to a wide deformation zone that includes the stretched continental margin of sundaland the philippine mobile belt pmb extending from luzon to the molucca sea and a mosaic of continental blocks around the ph australia sunda triple junction the gps geodyssea data are used to decipher the present kinematics of this complex area in the philippines the overall scheme is quite simple two opposing rotations on either side of the left lateral philippine fault clockwise to the southwest and counterclockwise to the northeast transfer 55 per cent of the ph sundaland convergence from the manila trench to the northwest to the philippine trench to the southeast further south 80 per cent of the ph sunda convergence is absorbed in the double subduction system of the molucca sea and less than 20 per cent along both continental margins of northern borneo finally within the triple junction area between the sundaland ph and australia plates from sulawesi to irian jaya preferential subduction of the celebes sea induces clockwise rotation of the sulu block which is escaping toward the diminishing celebes sea oceanic space from the eastwardadvancing ph plate to the south we identify an undeformed banda block that rotates counterclockwise with respect to australia and clockwise with respect to sundaland the kinematics of this block can be defined and enable us to compute the rates of southward subduction of the banda block within the flores trench and of eastward convergence of the makassar straits with the banda block the analysis made in this paper confirms that this deformation is compatible with the eastward motion of sundaland with respect to eurasia determined by the geodyssea programme but is not compatible with the assumption that sundaland belongs to eurasia as was often assumed prior to this study

present day crustal deformation along the philippine fault in luzon philippines abstract the philippine fault results from the oblique convergence between the philippine sea plate and the sunda block eurasian plate the fault exhibits left lateral slip and transects the philippine archipelago from the northwest corner of luzon to the southeast end of mindanao for about 1200 km to better understand fault slip behavior along the philippine fault eight gps surveys were conducted from 1996 to 2008 in the luzon region we combine the 12 yr survey mode gps data in the luzon region and continuous gps data in taiwan along with additional 15 international gnss service sites in the asia pacific region and use the gamit globk software to calculate site coordinates we then estimate the site velocity from position time series by linear regression our results show that the horizontal velocities with respect to the sunda block gradually decrease from north to south along the western luzon at rates of 85 49 mm yr in the west northwest direction this feature also implies a southward decrease of convergence rate along the manila trench significant internal deformation is observed near the philippine fault using a two dimensional elastic dislocation model and gps velocities we invert for fault geometries and back slip rates of the philippine fault the results indicate that the back slip rates on the philippine fault increase from north to south with the rates of 22 37 and 40 mm yr respectively on the northern central and southern segments the inferred long term fault slip rates of 24 40 mm yr are very close to back slip rates on locked fault segments suggesting the philippine fault is fully locked the stress tensor inversions from earthquake focal mechanisms indicate a transpressional regime in the luzon area directions of 1 axes and maximum horizontal compressive axes are between 90 and 110 consistent with major tectonic features in the philippines the high angle between 1 axes and the philippine fault in central luzon suggests a weak fault zone possibly associated with fluid pressure

spatial variations of the crustal stress field in the philippine region from inversion of earthquake focal mechanisms and their tectonic implications abstract this paper presents the spatial variations in the crustal stress field for the philippine region based on the stress configuration we divided the stress regime in the area of the philippines into four parts 1 trench perpendicular compressive stress axes 1 with an intermediate plunge were observed along the eastern subduction systems of the philippine mobile belt pmb which suggests that the subduction of the philippine sea plate psp is the primary factor controlling the stress distribution 2 characterized by systematically rotated 1 with very shallow plunge the central western portion of the philippines appears to be largely affected by the collision between the palawan block and the pmb 3 the 1 distribution in the northern part of luzon island north of 14 n is compatible with the plate motion of the psp whereas the extensional axes 3 has a fan shaped stress signature with n s direction at the south and a more spreading direction to the north which implies its north escapement 4 bounded by the scarborough seamount the northern and southern portions of the manila trench reveal a very distinct stress pattern with a trench parallel 3 to the north and a trench perpendicular 3 to the south this observation may infer a potential effect caused by the presence of the oceanic bathymetric highs on the behavior of the subduction process in addition to the identification of the four provinces with different stress regimes we suggest that the philippine fault zone may play an important role to alter the stress state in the region of the philippines

large scale model test investigation on wave run up in irregular waves at slender piles an experimental large scale study on wave run up generated loads on entrance platforms for offshore wind turbines was performed the experiments were performed at grosser wellenkanal gwk forschungszentrum kuste fzk in hannover germany the present paper deals with the run up heights determined from high speed video recordings based on the measured run up heights different types of prediction formulae for run up in irregular waves were evaluated in conclusion scale effects on run up levels seem small except for differences in spray however run up of individual waves is difficult to predict due to high importance of wave form

onland signatures of the palawan microcontinental block and philippine mobile belt collision and crustal growth process a review abstract the collision of the palawan microcontinental block with the philippine mobile belt had significantly influenced the geological evolution of the philippines multiple collisions involving several fragments through space and time resulted into the collage of terranes of varying origin exposed in this part of central philippines cusping of the overriding plate volcanic arc gap ophiolite emplacement incipient back arc rifting island rotation and tilting raised coastal terraces metamorphism intrusion of igneous rocks and steepened subducted slab as seen in focal mechanism solutions are some of the manifestations of this collision a late early miocene to early middle miocene age 20 16 ma is proposed for the major collision between the palawan indenter and the philippine mobile belt the collision boundary is located from the northern part of mindoro through the central mountain range swinging east of sibuyan island in the romblon island group and finally threading along the buruanga peninsula and eastern side of the antique ophiolite complex before exiting and connecting with the negros trench the collision through accretion and crustal thickening has contributed to the crustal growth of the philippine archipelago

spatial variation of the crustal stress field along the ryukyu taiwan luzon convergent boundary 1 we applied an improved stress inversion method to a comprehensive data set of earthquake focal mechanisms to depict the pattern of crustal stress along the western convergent boundary of the philippine sea plate our results indicate that the crustal stress along the ryukyu fore arc is segmented with boundaries at or near the places of seamount subduction including the tokara channel an extensional stress regime is observed along the entire ryukyu back arc implying that back arc rifting may have extended northward to kyushu a triangular area near the southernmost terminus of the ryukyu arc is characterized by a unique stress signature the eastern boundary of this ryukyu taiwan stress transition coincides with the 123 e meridian where the gagua ridge intercepts the ryukyu trench whereas its western boundary agrees remarkably well with the border between the postcollision and waning collision domains in northern taiwan the taiwan collision zone is dominated by compression that rotates locally according to the structural configuration of the lukang magnetization high lmh suggesting that the lmh may be critical in controlling the local stress distribution the stress signature of the luzon arc taiwan collision reaches as far south as 19 5 n the tectonic stress along the manila trench luzon fore arc is dominated by a complex regime of extension that cannot be explained by simple plate bending or in slab membrane stress since this extensional regime is observed only south of 22 n it probably marks the northern limit of the contemporary boundary between the subduction along the manila trench and the collision in taiwan

plate boundary forces are not enough second and third order stress patterns highlighted in the world stress map database 1 the world stress map project compiles a global database of contemporary tectonic stress information of the earth s crust early releases of the world stress map project demonstrated the existence of first order plate scale stress fields controlled by plate boundary forces and second order regional stress fields controlled by major intraplate stress sources such as mountain belts and zones of widespread glacial rebound the 2005 release of the world stress map project database provides for some areas high data density that enables us to investigate third order local stress field variations and the forces controlling them such as active faults local inclusions detachment horizons and density contrasts these forces act as major controls on the stress field orientations when the magnitudes of the horizontal stresses are close to isotropic we present and discuss examples for venezuela australia romania brunei western europe and southern italy where a substantial increase of data records demonstrates some of the additional factors controlling regional and local stress patterns

plate convergence at the westernmost philippine sea plate abstract to understand the convergent characteristics of the westernmost plate boundary between the philippine sea plate psp and eurasian plate ep we have calculated the stress states of plate motion by focal mechanisms cataloged by the harvard centroid moment tensor solutions harvard cmt and the broadband array in taiwan bats moment tensor 251 focal mechanisms are used to determine the azimuths of the principal stress axes we first used all the data to derive the mean stress tensor of the study area the inversion result shows that the stress regime has a maximum compression along the direction of azimuth n299 this result is consistent with the general direction of the rigid plate motion between the psp and ep in the study area in order to understand the spatial variation of the regional stress pattern we divided the study area into six sub areas blocks a to f based on the feature of the free air gravity anomaly we compare the compressive directions obtained from the stress inversion with the plate motions calculated by the euler pole and the global positioning system gps analysis as a result the azimuth of the maximum stress axis 1 generally agrees with the directions of the theoretical plate motion and gps velocity vectors except block c lanhsu region and block f ilan plain region the discrepancy of convergent direction near the ilan plain region is probably caused by the rifting of the okinawa trough the deviation of the 1 azimuth in the lanhsu region could be attributed to a southwestward extrusion of the luzon arc la block between 21 n and 22 n whose northern boundary may be associated with the right lateral ne sw trending fault i e huatung fault hf along the taitung canyon comparing the 1 stress patterns between block c and block d great strain energy along hf may not be completely released yet alternatively the upper crust of block c may significantly have decoupled from its lower crust or uppermost mantle

swell transformation across the continental shelf part i attenuation and directional broadening abstract extensive wave measurements were collected on the north carolina virginia continental shelf in the autumn of 1999 comparisons of observations and spectral refraction computations reveal strong cross shelf decay of energetic remotely generated swell with for one particular event a maximum reduction in wave energy of 93 near the virginia coastline where the shelf is widest these dramatic energy losses were observed in light wind conditions when dissipation in the surface boundary layer caused by wave breaking whitecaps was weak and wave propagation directions were onshore with little directional spreading these observations suggest that strong dissipation of wave energy takes place in the bottom boundary layer the inferred dissipation is weaker for smaller amplitude swells for the three swell events described here observations are reproduced well by numerical model hindcasts using a parameterization of wave friction over a movable sandy bed directional spectra that are narrow off the s

water waves over a muddy bed a two layer stokes boundary layer model this study aims to develop an asymptotic theory for the flow kinematics of a thin layer of viscous mud under water surface waves the mud depth the mud stokes boundary layer thickness and the wave amplitude are assumed to be comparable with one another and much smaller than the wavelength by virtue of this sharp contrast in length scales boundary layer equations are used to describe motion of both the mud and the immediately overlying water analytical solutions are sought and explicit expressions are obtained under progressive waves for the fluid velocity fields interface wave characteristics and wave damping rate at the first order and the steady mean discharge rate and mass transport velocity of mud at the second order effects of the mud layer thickness density and viscosity ratios on these kinematic quantities are examined in detail

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

observations of nearshore wave dissipation over muddy sea beds 1 wave propagation on the louisiana inner shelf is studied using concurrent measurements at two shallow water ocean observatory sites in sedimentary environments dominated by mud at one site and sand at the other contrary to the widely accepted hypothesis that mud induced wave dissipation is important only for long waves observations show significant damping of high frequency short waves which interact weakly with the bottom the mechanism of short wave dissipation is not understood numerical simulations show that other dissipative processes such as refraction or depth limited breaking do not account for the magnitude of the observed effects independent observations of strong sediment re working during storms suggest that these effects are related to sediment resuspension processes

geologic and paleomagnetic evidence for a possible miocene collision in western panay central philippines the palawan metamorphic terrane represented by paleozoic basement on the northwestern tip of panay central philippines collided with the remainder of the island between late oligocene and middle miocene time subduction continued north and south of the collision zone along the manila and negros trenches respectively evidence for the collision is 1 the juxtaposition along a major fault of old melange terrane against highly metamorphosed paleozoic continental sediments similar to those exposed on palawan island itself and 2 miocene thrusting in the melange terrane oriented parallel to the suture paleomagnetic studies show that since the early miocene the island of panay has rotated 20 clockwise declination 20 9 inclination 22 1 for six sites whereas the northern philippines has rotated counterclockwise these contrasting rotations are consistent with a collision of the philippine arc with a northeastward converging palawan block

modelling of wave damping at guyana mud coast abstract the guyana coastal system is characterized by very thick deposits of amazon mud and high mud concentrations in its coastal waters the mud deposits can be quite soft and may liquefy under incoming waves subsequently the liquefied mud damps the incoming waves effectively this paper presents a simple model to predict wave attenuation over soft fluid mud beds this model is based on the two layer approach by gade gade h g 1958 effects of a non rigid impermeable bottom on plane surface waves in shallow water journal of marine research 16 2 61 82 which is implemented in the standard version of the state of the art wave prediction model swan input to the mud wave damping module consists of the extension thickness density and viscosity of the liquefied fluid mud layer the model is validated against small scale wave attenuation measurements carried out in a laboratory wave flume the model predictions agree favourably with the experimental data next the model is applied to predict wave height and wave attenuation in the guyana coastal system extension and thickness of the liquefiable layer could be assessed from dual frequency echo soundings in the absence of field data the density of the liquefied mud layer is obtained from literature whereas the value of the mud s viscosity had to be established by trial and error the selected value though is in the range of literature values the model predicts significant wave attenuation the computed changes in wave energy spectrum agree qualitatively with measurements in surinam whereas the decrease in significant wave height agrees more or less with historic observations along the guyana coast

plate coupling along the manila subduction zone between taiwan and northern luzon abstract we use gps data trench parallel gravity anomaly tpga and bathymetry to infer plate coupling patterns along the manila subduction zone using a block model and a fault geometry constrained by seismicity we simultaneously solve for the location of euler pole and angular velocity between the sunda and luzon blocks as well as the slip deficit rate on plate interface our estimates show that the euler pole between the sunda and luzon blocks is situated at southern palawan near 8 3 n and 119 4 e with the angular velocity of 4 6 myr 1 the estimated convergence rate along the manila trench continuously decreases southward from 91 mm yr at the northern tip of luzon to 55 mm yr north of mindoro the inversion of gps data reveals partially locked fault patches extending from the west luzon trough to the east of scarborough seamount chain the slip deficit rate in this region is in the range of 20 30 mm yr corresponding to a coupling ratio of 0 4 however the fault slip behavior is not well resolved near the north luzon trough based on a good correlation between locations of large subduction zone earthquakes and areas possessing gravity low we investigate a variety of tpga based plate coupling models assuming different scaling between tpga values and plate coupling ratios the tpga based plate coupling models offer plausible rupture scenarios which are not constrained by current gps data the partially locked fault zone near 15 16 5 n may be associated with the subducted scarborough seamount wherein oceanic floor is highly fractured the great subduction zone earthquake propagates beneath the scarborough seamount seems to be unlikely the densification of gps network in central luzon and seafloor geodetic observations close to trench axis are crucial to distinguish the detailed fault coupling patterns and the role of subducted seamounts

role of fluids in faulting inferred from stress field signatures the stress orientation signature of weak faults containing high pressure fluids has been observed for segments of the san andreas fault system in southern california the inferred lithostatic fluid pressures extend into the surrounding relatively intact rock in a zone scaling with the width of the interseismic strain accumulation repeated strain related fracturing and crack sealing may have created low permeability barriers that seal fluids into the network of currently active fractures

stress permutations three dimensional distinct element analysis accounts for a common phenomenon in brittle tectonics 1 using three dimensional 3 d distinct element modeling we explored a variety of simulations to characterize and interpret the stress permutations in brittle tectonics stress inversions of fault slip data or earthquake focal mechanisms often revealed such permutations the main aim of our study is to produce simple mechanically consistent 3 d models that account for these switches between the principal stress axes 1 2 or 2 3 even with simple boundary conditions the stress changes induced by variations in rheology are large enough to modify the local tectonic behavior and to produce permutations of principal stress axes rather than simple directional changes of stress axes which exist but often remain limited the relative variations in principal stress magnitudes are the major cause of permutations 1 2 and 2 3 in nature permutations being left apart the orientations of axes often remain tightly clustered in our experiments we adopted a ratio 2 3 1 3 of 0 5 which makes permutations difficult low favors 2 3 permutations high favors 1 2 permutations and we explored a variety of tectonic situations involving compression extension and strike slip our experiments indicate that the major causes of stress permutations are the heterogeneity of the brittle deformation e g intact rock massifs between heavily faulted deformation zones and the anisotropy of the mechanical properties that results from fracturing and faulting which concur to modify the mechanical balance inside the analyzed volume and to produce stress permutations our models show that contrasts and anisotropy in rock properties favor stress permutations of major importance is the existence of relatively resistant zones at tips of the deformed zones acting as channels where stress concentrates and switches occur because in nature such zones move in time and space it is not surprising that stress permutations are so pervasive

temporal and spatial variation of stress field in taiwan from 1991 to 2007 insights from comprehensive first motion focal mechanism catalog article i nfo the total amount of 4 761 focal mechanisms was determined based on p wave first motion polarities from 1991 to 2007 in taiwan region this dataset offers us a good opportunity to examine temporal and spatial variability of the stress field we find that the spatial variations of stress axes are mainly controlled by tectonic structures while the temporal changes are greatly influenced by the chi chi earthquake the orientation of the maximum horizontal compressive stress axes sh shows a general agreement with the direction of plate motion between a depth range of 0 30 km the 20 anticlockwise rotation of sh from the longitudinal valley lv to western taiwan is probably caused by the left lateral motion on the longitudinal valley fault lvf that has consumed part of the oblique motion of plate convergence on the other hand part of the oblique convergence is transferred into the central range and the hsuehshan range judging from counterclockwise rotation of sh from east to west and strike slip faulting in the hsuehshan range most events with a depth greater than 30 km occur offshore eastern taiwan and the azimuth of sh is close to e w directed different from nw se directed at shallow depths this may infer the existence of the transition of lithosphere rheology in offshore eastern taiwan the trends of sh in the depth of 0 10 km are strongly affected by the coseismic stress change of the chi chi earthquake in the northern half of the chi chi rupture area the trends of sh rotate 30 clockwise and the stress ratio increased by a factor of six after the mainshock the orientations of sh still differ by 30 in 2007 comparing to that in the period before the chi chi earthquake the variation of sh trend is more diverse in the southern half of the rupture area showing 20 counterclockwise rotation immediately after the 1999 mainshock followed by a clockwise rotation the trend of sh returns to the pre seismic direction of 110 in 2001 these notable changes of sh before and after the chi chi mainshock suggests that the magnitude of background stress in the rupture area is close to the coseismic stress drop we also recognize a significant 10 counterclockwise rotation of sh in the entire chi chi rupture area between 1991 and 1999 before the earthquake took place however to the south of the chi chi rupture the trends of sh remain little changed before and after the chi chi earthquake

a bingham plastic model for a muddy seabed under long waves the dynamics of cohesive mud under wave action is important to the geological processes along some coastlines because of the complex rheological properties there have been vastly different models for predicting the effect of mud on wave attenuation this paper focuses on a non newtonian behavior of highly concentrated mud whereby shearing strain occurs only if the shear stress exceeds a certain threshold called the yield stress the wave induced motion in a thin mud layer beneath a clear water layer is studied there are circumstances in which the yield stress is dominant and acts like a coulomb friction mud motion can change from continuous to intermittent as the wave attenuates the effects of the non newtonian behavior on wave damping are compared to other more common causes of dissipation

a numerical study on nonlinear surface wave evolution over viscoelastic mud abstract accurate prediction of wave energy dissipation over mud is essential for understanding nearshore circulation sediment transport processes and design of engineering projects along muddy coasts it is desirable to be able to simulate wave propagation over a wide range of mud behaviors as wave evolution is highly dependent on mud rheological characteristics in this study we develop a numerical model to simulate surface wave evolution over viscoelastic muds this new wave mud interaction model is comprised of a frequency domain phase resolving model for wave propagation that explicitly solves nonlinear wave wave interactions and a model for mud induced surface wave dissipation and modulation model results show satisfactory agreement with laboratory measurements of wave height and attenuation rates over mud and shows improvement over the model with a viscous mud mechanism the model is then used to investigate the combined effect of mud viscoelasticity and nonlinear wave wave interactions on surface wave evolution cnoidal and random wave simulations are conducted in general qualitative measures such as shape of cnoidal waves or pattern of variation in h r m s of random waves are dictated by direct mud induced damping which due to resonance effects has a substantially different structure over viscoelastic mud compared to viscous mud nonlinear interactions affect spectral shape and distribution of energy loss across the spectrum subharmonic interactions cause indirect damping of high frequencies but ameliorate damping of harmonics around mud s resonance frequency therefore neglecting mud elasticity can result in significant errors in estimation of bulk wave characteristics and spectral shape

wave dissipation by muddy seafloors 1 muddy seafloors cause tremendous dissipation of ocean waves here observations and numerical simulations of waves propagating between 5 and 2 m water depths across the muddy louisiana continental shelf are used to estimate a frequency and depth dependent dissipation rate function short period sea 4 s and swell 7 s waves are shown to transfer energy to long period 14 s infragravity waves where in contrast with theories for fluid mud the observed dissipation rates are highest the nonlinear energy transfers are most rapid in shallow water consistent with the unexpected strong increase of the dissipation rate with decreasing depth these new results may explain why the southwest coast of india offers protection for fishing and for the 15th century portuguese fleet only after large waves and strong currents at the start of the monsoon move nearshore mud banks from about 5 to 2 m water depth when used with a numerical nonlinear wave model the new dissipation rate function accurately simulates the large reduction in wave energy observed in the gulf of mexico

tertiary and quaternary magmatism in mindanao andleyte philippines geochronology geochemistry and tectonic setting abstract a set of 230 igneous rock samples from leyte and mindanao philippines has beenstudied for their major and trace element chemistry mineralogy and geochronology 40k 40ar dating several volcanic sectors have been distinguished according to their geographic position tectonic framework and geologic history a in leyte the ophiolitic basement at 50 ma is morb type but characterized by negative nb 0 anomalies typical of back arc basalts miocene 20 and 11 ma and late pliocene to quaternary magmas 2 5 0 4 show typical island are talc alkaline characteristics in eastern mindanao eocene 47 46 ma magmatic rocks are arc tholeiitic while those of miocene 18 6 ma and early pliocene 4 5 ma age are generally talc alkaline dacites of 12 ma age have adakitic compositions the young 3 ma volcanics in surigao and north davao are mostly adakitic although talc alkaline basalts and basaltic andesites dated at 0 9 ma overlie adakites in surigao b in daguma oligocene 30 ma arc tholeiitic diorites intrude older arc sequences miocene0 17 0 7 7 ma dikes and flows in both daguma and sarangani peninsula are mostly talc alkaline andesites although minor adakites dated back to 18 ma have been found overlying undated breccias in daguma in southeast daguma the tale alkaline volcanic substratum of pliocene recent volcanoes is 8 6 ma old quaternary 1 83 ma to present andesitec and dacitic volcanics in mt parker and mt matutum southwest daguma are either talc alkaline or adakitic while mt blit northwest daguma is built up by 2 5 ma old hfse enriched basalts 1 ma old normal talc alkaline andesites and 0 8 ma old adakites c in central mindanao miocene andesites yield ages of 20 and 16 ma pliocene 2 5 ma 0 volcanism is dominated by talc alkaline and hfse enriched basalts and basaltic andesites quaternary lavas i ma to present range in character from tale alkaline to shoshonitic and exhibit wide variations in their trace element behavior d miocene lavas of 19 11 ma age in zamboanga are mostly talc alkaline though a lone 0 18 ma old adakite was also encountered pliocene quaternary basalts and basaltic andesites dated back to 2 2 1 ma are hfse enriched such that some display positive nb anomalies in spidergrams adakites of 3 8 0 7 ma age are located nearer the sulu trench the compositional variety of tertiary recent island arc magmatism in leyte and mindanao couldbe linked to the tectonic settings prevailing during each magmatic phase and is here utilized to refine paleokinematic reconstruction of the tectonic history of the philippine archipelago the identification of adakites at certain time periods is used to date timings of subduction initiations and arc polarity reversals major and abrupt changes in geochemical compositions are well exhibited by the pliocene pleistocene volcanism involving the collision related shift from talc alkaline to shoshonitic volcanism in central mindanao and adakite production in all the volcanic sectors defined the origin of the hfse enrichments among lavas younger than 3 ma could be attributed to the presence in their mantle source of an oib component or alternatively to the enrichment of this mantle source through percolation of slab derived melts

geological setting of the rapu rapu gold rich volcanogenic massive sulfide deposits albay province philippines gold rich fe cu zn volcanogenic massive sulfide deposits occur within strata of probable jurassic age on rapu rapu island in albay province philippines massive sulfides at the ungay malobago and hixbar deposits are spatially associated with dacitic volcanic rocks within a highly deformed sequence of mafic volcanic and quartzofeldspathic sedimentary rocks the massive sulfide deposits formed at the stratigraphic contact between footwall dacites and hangingwall mafic volcanic and quartzofeldspathic rocks the deposits and their host strata have undergone regional metamorphism with strong penetrative deformation metamorphic mineral assemblages and textural evidence suggest that peak metamorphism was upper greenschist to lower amphibolite grade and syn d1 deformation based on the age of regional metamorphism deformation is inferred to be mid tertiary in age deformation at rapu rapu resulted in reorientation of the strata into a broad antiform with strong shallow plunging elongation fabrics overturning of the volcanic sequence that hosts the ungay malobago deposit and complex folding of the mineralized zones the present highly linear form of the ungay malobago deposit is mainly a product of this ductile strain

an improved method for determining the regional stress tensor using earthquake focal mechanism data application to the san fernando earthquake sequence the orientations of fault planes and slip directions indicated by a population of earthquake focal mechanisms can be used to determine best fit regional principal stress directions and r 2 1 3 1 a measure of relative stress magnitudes under the assumption of uniform stress in the source region this analysis allows for the possibility that failure occurs on preexisting zones of weakness of any orientation in the inversion we perform a grid search of stress models to find the one which requires the smallest total rotation of all the fault planes that is needed to match the observed and predicted slip directions the method allows for errors in orientations of both the fault planes and slip directions we have an objective means for identifying the more likely of the two possible fault planes from each focal mechanism relative to a given stress model thus we do not face the problem of ambiguity of nodal planes which plagues other analyses of this kind by using a grid search of stress models rather than a linearization scheme we are able to perform a realistic error analysis and thus establish confidence limits for the preferred regional stresses the method can be used to investigate possible stress inhomogeneities during earthquake sequences by analyzing subsets of the data population the technique has been applied to 76 events from the san fernando earthquake sequence for which we have found best fit stresses plunge and azimuth 1 7 187 2 27 281 3 62 84 and r 0 65 the average misfit between this stress model and all the data is about 8 and all but eight of the aftershocks have misfits of less than 20 these values are considerably less than the uncertainty of the focal mechanism determinations therefore significant stress inhomogeneities are not required by the data our analysis does not support the suggestion of a change in stresses during the aftershock sequence as proposed by others on the basis of an apparent change in focal mechanisms

genetic relations between the central and southern philippine trench and the sangihe trench we surveyed the junction between the central and southern philippine trench and the sangihe trench near 6on using swath bathymetry gravity and magnetics these data along with seismicity allow us to discuss the genetic relations between these trenches and the forces acting on converging plates our final model favors the northern extension of the halmahera arc up to 8on with three segments offset left laterally along nw se transform faults accretion of the northern segment to mindanao island 4 to 5 m y ago resulted in the failure within the philippine sea plate east of the arc initiation of the philippine trench between 7on and 10on agrees with the maximum recorded depth of the philippine trench floor 10 000 m below sea level and philippine sea slab 200 km south of 6on trench junction another segment of the arc is being subducted beneath the sangihe margin while south of 3on the southern segment of the halmahera arc is still active the rapid southward shallowing of the trench floor along the southern philippine trench the type of faulting affecting both sides of the trench the lack of significant interplate seismicity and the concentration of the seismicity beneath the miangas talaud ridge are interpreted as a slowing down of the subduction along this branch of the philippine trench compared with the rest of the subduction zone the sangihe deformation front has been recognized up to 7on but seems active only south of 6on

physical processes and fine grained sediment dynamics coast of surinam south america abstract the prograding holocene mud wedge between the amazon and orinoco rivers in the trade wind belt of northeastern south america provides a modern day example of muds accumulating under moderate wave energy conditions gigantic shore attached mudbanks 10 km 20 km composed partly of thixotropic fluid mud gel front this coast every 30 60 km to form a buffer to wave attack and a temporary storage for fine grained sediments this mesotidal coast tide range 2 0 m with gentle offshore slope 0 0006 allows the exposure twice a day of extensive tidal flat deposits which are backed by mangrove swamps on a well developed chenier plain complex field experiments were conducted in surinam du ing 1975 and 1977 to provide new information on process form relationships in this interesting but unusual muddy environment simultaneous measurements of waves currents tide elevation suspended sediment concentration and variations in mud density show that soft intertidal and subtidal muds are suspended at both tide and wave frequency suspended sediment concentrations typically exceed 1 000 mg l at the surface as incoming solitary like waves partially disperse fluid mud into overlying water on a falling or rising tide redeposition of mud may occur near time of high tide the strong attenuation of shallow water waves by these muds provides conditions that are favorable for further sedimentation high concentrations of suspended fluid mud together with solitary like waves from the northeast throughout the year can lead to extraordinarily high net sediment transport rates in the nearshore zone calculations based on solitary wave theory and on data obtained from this study indicate that 15 65 106 m3 of mud can move along shore each year without involving breaking waves the concept of radiation stress and a nearshore circulation cell or bedload transport farther offshore outside the zone of wave dominance wind driven currents and the guiana current combine to transport muds to the northwest consistent with the observed direction of mudflat migration

crust and upper mantle structure and its tectonic implications in the south china sea and adjacent regions abstract we present a 3d s velocity model for the crust and upper mantle of the south china sea and the surrounding regions constrained from the analysis of over 12 000 of fundamental rayleigh wave dispersion curves between 10 s and 150 s periods the lateral resolution was found to vary from 2 to 4 with the increasing period over the study region a robust scheme of debayle and sambridge allowed us to conduct the tomographic inversion efficiently for massive datasets group velocity maps varying with period show lateral heterogeneities well related to the geological and tectonic features in the study region the 3d s velocity model was constructed from the 1d structure inversion of the tomographic group velocity dispersion curves at each node the obtained average crustal structure is similar to the prem model while the average mantle velocity is typically lower than the global average the complicated 3d structures reveal three prominent features correlated with geological divisions sea basin regions island and arc regions and continental regions the derived crustal and lithospheric thicknesses range from 15 to 50 km and from 60 to 140 km respectively with the thinnest in the south china sea the thickest in eastern tibet and the yangtze block and the medium in the south china fold belt indochina and island arc regions our results further confirm that 1 a mesozoic subduction zone which is interpreted as the tectonic weak zone during the paleogene exists along the south china margin 2 the influence of the indochina extrusion along the red river fault is limited for the south china sea region 3 there is a slab remnant of the proto south china sea beneath borneo new findings suggest that the mesozoic subduction zone should be built into any evolution model for the region as well as the other two major tectonic boundaries of the red river fault and proto south china sea subduction zone

nonlinear wave dynamics in the presence of mud induced dissipation on atchafalaya shelf louisiana usa abstract the interplay between wave nonlinearity and mud induced dissipation is studied here using wave and sediment transport measurements collected in winter 2008 on the atchafalaya shelf louisiana between the 8 m and the 4 m isobaths this study focuses on the relatively energetic storm that occurred on march 4th 2 m wave height in 8 m water depth which caused significant bed reworking and left in its wake a 20 cm layer of hindered settling fluid mud while the net wave dissipation rate was maximal during the hindered settling phase after the storm consistent with previous observations significant dissipation was observed throughout the storm duration with a secondary maximum associated with the peak of the storm and the maximum bed reworking effects the effects of mud induced dissipation on the nonlinear shoaling process are investigated here using triads a newly developed spectral model for nonlinear shoaling of waves with mud parameter values estimated using a crude inverse modeling approach the model was used to separate the effects of mud induced dissipation from nonlinear transfers during both the erosion and deposition phases of bed reworking numerical simulations show that nonlinear transfers are always active even under strong dissipation regimes remarkably they also suggest that heavy mud induced dissipation causes the nonlinear energy cascade to revert direction with bulk transfers occurring from high frequency to low frequency spectral bands this supports the hypothesis that the interplay between mud induced dissipation which is much more significant than the bottom dissipation over sandy beds and wave nonlinearity can drain energy from the entire spectrum and not just from the frequency bands that interact directly with the seafloor the decay of high frequency variance induced by the reversal of spectral flux in effect reduces the nonlinearity of the wave field

wave run up on cylindrical and cone shaped foundations for offshore wind turbines abstract during the last decade several offshore wind farms were built and offshore wind energy promises to be a suitable alternative to provide green energy however there are still some engineering challenges in placing the foundations of offshore wind turbines for example wave run up and wave impacts cause unexpected damage to boat landing facilities and platforms to assess the forces due to wave run up the distribution of run up around the pile and the maximum run up height need to be known this article describes a physical model study of the run up heights and run up distribution on two shapes of foundations for offshore wind turbines including both regular and irregular waves the influence of wave steepness wave height and water depth on run up is investigated the measured run up values are compared with applicable theories and previous experimental studies predicting run up on a circular pile the results show that the shape of the foundation substantially affects the maximum run up level increasing the expected run up value a new relationship between the wave climate regular and irregular waves and the run up is suggested for this the velocity stagnation head theory is adjusted and second order stokes equations are used to calculate the wave kinematics in the crest the variation of the run up around the pile is measured and it is found that the position with the lowest run up level is located under 135 while the run up at that position amounts to approximately 40 to 50 of the maximum run up

crustal stress field in southern california and its implications for fault mechanics we present a new high spatial resolution image of stress orientation in southern california based on the inversion of earthquake focal mechanisms we use this image to study the mechanics of faulting in the plate boundary region the stress field contains significant spatial heterogeneity which in some cases appears to be a result of the complexity of faulting and in other cases appears to be a cause temporal changes in the stress field are also observed primarily related to major earthquakes the observed 15 10 rotation of the stress axes due to the 1992 m7 3 landers mainshock implies that the deviatoric stress magnitude in the crust is low of the order of 10 mpa this suggests that active faults in southern california are weak the maximum principal stress axis near the san andreas fault is often at 50 to the fault strike indicating that the shear stress on the fault is comparable to the deviatoric stress the san andreas in southern california may therefore be a weak fault in a low strength crust

philippine fault a key for philippine kinematics on the basis of new geologic data and a kinematic analysis we establish a simple kinematic model in which the motion between the philippine sea plate and eurasia is distributed on two boundaries the philippine trench and the philippine fault this model predicts a velocity of 2 to 2 5 cm yr along the fault geologic data from the visayas provide an age of 2 to 4 ma for the fault an age in good agreement with the date of the beginning of subduction in the philippine trench the origin of the philippine fault would thus be the flip of subduction from west to east after the locking of convergence to the west by the collision of the philippine mobile belt with the eurasian margin

deformation and stress states along the central segmentof the philippine fault implications to wrench fault tectonics abstract this paper discusses the space time evolution of deformation and states of stress along the central segment of the left lateral philippine fault from bondoc peninsula to leyte island this segment is characterized by deformations and stress states which do not always conform to predictions of classical wrench tectonics directions of the maximum principal stress axis 1 are in some cases oriented perpendicular to the fault trace contradicting the classical fault strike 1 angle of 30 45 in other cases it is the minimum axis 3 that orients perpendicular to the fault or in other cases parallel to it rotations of 1 directions as functions of time and space are observed all along the fault more complex stress regimes are observed in areas where the philippine fault intersects with two other major left lateral strike slip faults namely the legaspi lineament in the ragay gulf area and the sibuyan sea fault in masbate these observations strongly underline the urgence of prudence and precaution in applying wrench fault tectonic principles in natural examples

transition from oblique subduction to collision earthquakes in the southernmost ryukyu arc taiwan region tectonic characteristics of the region between taiwan and the southernmost ryukyu arc are inferred from a detailed analysis of local seismicity and source parameters of 62 recent earthquakes of 5 5 mb 6 6 five major seismogenic structures can be delineated the collision seismic zone csz the interface seismic zone isz the wadati benioff seismic zone wbsz the lateral compression seismic zone lcsz and the okinawa seismic zone osz in the csz located along the east coast of taiwan and offshore earthquake focal mechanisms show horizontal p axes distributed in two directions 287 10 and 333 16 possibly reflecting a strain partition associated with the relative plate convergence between the eurasia plate and the philippine sea plate the corresponding seismic strain tensor indicates a maximum compressive strain rate of 1 2 10 7 yr 1 along 293 and a comparable extension in vertical direction presumably resulted from plate collision in the region the geometry of the isz which is distorted significantly at its westernmost end can be approximated by a north dipping plane that is gradually pushed northward with increasing dip the seismogenic portion of the interface spans a short depth range from 10 km to 35 km a clear pattern of earthquake slip partition is observed the average slip vector residual is as large as 35 seismic strain patterns within the subducted philippine sea slab show predominantly downdip extension between 80 and 120 km and downdip compression at 270 km different from the pattern of strain segmentation observed for the rest of the ryukyu arc where the northern and southern portions are dominated by downdip extension and compression respectively owing to the large convergence obliquity the slab is descending at a rate significantly slower near taiwan than in the southern ryukyu thus we interpret the appearance of downdip extension within the subducted lithosphere as a combined result of oblique subduction and the slab s negative buoyancy a number of thrust or oblique strike slip earthquakes between isz and wbsz show a consistent pattern of lateral compression with p axes oriented roughly parallel to the local strike of the trench arc system they are probably due to the compressive strain originated from the collision and transmitted laterally within the lithosphere shallow normal faulting earthquakes show successive rotation of t axes from approximately n s in the okinawa trough to approximately e w in northeast taiwan possibly as a result of interaction between the extension from the opening of the okinawa trough and the compression from collision one normal event january 18 1991 mb 5 9 occurred in the central range with t axis roughly parallel to the structural trend of taiwan implying that the nature of the orogeny in taiwan has changed from thin skinned deformation to lithospheric collision involving the whole crust and uppermost mantle

a new insight on the geometry of subducting slabs in northern luzon philippines abstract we used hypocentral and focal mechanism data in order to characterize the tectonic configuration of northern luzon and propose a model for describing the geometry of the subducted slab of the eurasian plate beneath the northern segment of the manila trench we took into consideration some of the observed bathymetric features i e the bend in the trench line of the manila trench at 20 n lat and the collision and subduction of an extinct mid oceanic ridge at 16 17 n lat and the intraplate deformation pattern in the north luzon ridge region in depicting a model that is consistent with most of the observed features and phenomena in the study area an earthquake catalog covering the period from 1619 to 1997 was used in studying the distribution of strain energy released while earthquakes covering the instrumental years 1963 1997 were used to analyze vertical distributions of earthquakes in sections we refined the location of some of the historical earthquakes using macroseismic data the inclusion of historical records despite their inaccurate locations could prevent misidentification of some regions as aseismic zones or gaps in seismicity especially along structures with long return periods a focal mechanism database that covers the period from 1963 to 1997 was used to decipher intraplate deformation pattern and subduction process of the slab the focal mechanism database consists of first motion solutions from previous authors and cmt solutions from harvard university seismicity strain energy release and focal mechanism maps and cross sections were drawn in order to have a three dimensional visualization of the geometry and stress regimes of seismogenic zones a new model of the subducted slab of the eurasian plate beneath the manila trench is proposed the model suggests the collision and subsequent partial subduction of a buoyant plateau at around 20 n lat to explain the sharp bend in the trench line the complicated deformation pattern on the overriding plate fronting the bend and the shallow dip of the subducted slab beneath this zone a tear in the slab is also inferred to be present as evidenced by the observed gap in strain energy release and the abrupt change in dip from shallow to steep south of 18 n lat the gap in seismicity and strain energy release 65 300 km depth at around 17 n lat may be used to infer the trajectory and location of the subducted extinct mid oceanic ridge mor the aseismic behavior is probably caused by the subducted ridge which is still hot and is deforming plastically this is supported by the heatflow data which shows high values along the extinct mor the subducted part of the mor may serve as the weakest zone where this tear could be localized the tear may also explain the cause of the abrupt termination of the eastern chain of volcanoes south of 18 n lat the above model which is a refinement of the model introduced by yang et al tectonophysics 258 1996 85 is consistent with the observed seismicity and deformation pattern observed bathymetric features spatial distribution and geochemical character of volcanism in northern luzon

responses of bingham plastic muddy seabed to a surface solitary wave in this paper we investigate the dynamics of muddy seabed motions induced by a surface solitary wave the muddy seabed is characterized as bingham plastic mud we focus our attention on the situations where the horizontal scale of the wave induced mud flow is much larger than the vertical scale the thickness of the mud layer is also assumed to be much smaller than the water depth above with these simplifications the dynamic pressure in the mud column remains a constant and the vertical displacement of the water mud interface is negligible the horizontal gradient of the wave induced dynamic pressure along the water seabed interface drives the motions in the mud bed for a bingham plastic muddy seafloor the mud moves either like a solid plug flow or like a viscous fluid shear flow depending on whether the magnitude of shear stress is in excess of the yield stress velocities inside these two different flow regimes and the location s of the yield surface s vary in time as functions of water mud interfacial pressure gradient and the properties of the bingham plastic mud a semi analytical approach is developed in this paper to analyse the motions inside the mud bed under a surface solitary wave loading three possible scenarios are discussed to illustrate the complexity of the seafloor responses the formula for the damping rate caused by the energy dissipation inside the muddy seabed is also derived using realistic values of the physical parameters the present results for damping rate agree qualitatively with the available field observations

structural variety and tectonic evolution of strike slip basins related to the philippine fault system northern luzon philippines the northern part of the philippine fault zone in luzon corresponds to a complex braided system of left lateral strike slip faults the nw oriented main active branch which emerges from the philippine sea splits into an array of north striking splays responsible for the tectonic evolution of the central cordillera this complex fault pattern has favored local stress field variations strike slip basins have evolved in this framework along or in the vicinity of the main splays of the philippine fault apart from the classical pull apart tectonics on a releasing fault termination overlap or bend we describe other mechanisms such as the strike slip tilting or the warping of a strip limited by two strike slip faults the strike slip basins are good recorders of the evolution of the philippine fault system those located along the north striking cordilleran faults individualized in late early pliocene to pleistocene time when the present day philippine fault initiated but their main control is the fault shape acquired in upper middle miocene time the recent tectonic evolution of the fault system is best recorded in central luzon where the active basins trace an asymmetrical v shape with the longest branch trending ne parallel to the east luzon trough and the shortest one trending nw related to the philippine fault both the fault pattern and the basin distribution demonstrate the influence of the benham rise in the tectonic evolution of luzon the structural setting is interpreted as the result of an early miocene collision event between the benham rise and the eastern margin of luzon and subsequent inception of the nw striking strand of the philippine fault the present locations of the basins result from the interaction between the structural heritage and the present day regional plate motion