stuck in the mud earthquake nucleation and propagation through accretionary forearcs 1 subduction zone earthquakes can propagate to the surface causing large seafloor displacements resulting in tsunamis this requires the earthquake to rupture through clay rich sediments of the accretionary wedge which are largely aseismic as found previously the frictional properties of a range of wet clays at low slip velocity are velocity strengthening thus inhibiting earthquake nucleation however at high slip velocity the same materials weaken almost immediately resulting in a negligible critical slip weakening distance and fracture energy we interpret this behaviour as rapid thermal pressurization of the pore fluid within the clay gouge the lack of fracture energy can explain how a large rupture propagating from depth might not be arrested by clay rich velocity strengthening sediments as is commonly seen the results suggest that generally earthquakes may be difficult to nucleate on mature faults dominated by clay but the propagation of earthquakes through these zones is energetically very favourable

comparison of smectite and illite rich gouge frictional properties application to the updip limit of the seismogenic zone along subduction megathrusts along plate boundary subduction thrusts the transformation of smectite to illite within fault gouge at temperatures of 150 c is one of the key mineralogical changes thought to control the updip limit of seismicity if correct this hypothesis requires illite rich gouges to exhibit frictionally unstable velocity weakening behavior here we report on laboratory experiments designed to investigate the frictional behavior of natural and synthetic clay rich gouges we sheared 5 mm thick layers of commercially obtained pure ca smectite a suite of smectite quartz mixtures and natural illite shale grain size ranging from 2 to 500 m in the double direct shear geometry to shear strains of 7 30 at room humidity and temperature xrd analyses show that the illite shale contains dominantly clay minerals and quartz within the clay sized fraction 40 mpa our data specifically the velocity strengthening behavior of illite shale under a wide range of conditions do not support the hypothesis that the smectite illite transition is responsible for the seismic aseismic transition in subduction zones we suggest that other depth and temperature dependent processes such as cementation consolidation and slip localization with increased shearing may play an important role in changing the frictional properties of subduction zone faults and that these processes in addition to clay mineralogy should be the focus of future investigation

new insights into deformation and fluid flow processes in the nankai trough accretionary prism results of ocean drilling program leg 190 moore g f taira a klaus a becker l boeckel b cragg b a dean a fergusson c l henry p hirano s hisamitsu t et al 2001 new insights into deformation and fluid flow processes in the nankai trough accretionary prism results of ocean drilling program leg 190 geochemistry geophysics geosystems 2 article no 2001gc000166

sediment deformation and hydrogeology of the nankai trough accretionary prism synthesis of shipboard results of odp leg 131 the main objective of leg 131 was to provide data on the deformational processes and associated hydrogeology of the nankai prism toe drilling succeeded for the first time in the history of ocean drilling in penetrating the complete sedimentary sequence to basaltic basement reaching 1327 mbsf metres below seafloor with good core recovery 55 excellent correlation of the lithology and structure including the frontal thrust and the decollement with seismic reflection images was also determined bedding dips faults and shear bands analyzed in the cores confirm the pattern of deformation to be mainly due to nw se shortening as expected from the plate tectonic convergence vector below the decollement no significant deformation features were observed indicating that the decollement is a sharp discontinuity in stress transmission physical properties data show major discontinuities at the decollement notably an increase in porosity below the later this may indicate excess pore pressure in the subducted section and decollement zone a less marked increase in porosity below the frontal thrust may reflect the youthfulness of this feature attempts to make downhole measurements were severely hampered by unstable hole conditions but useful constraints have been placed on the thermal regime and some calibration of laboratory physical properties toin situ conditions has been provided andin situ stress and pore pressure were measured in the uppermost sediments evidence of channelized fluid flows is inconclusive no sharp geochemical signatures or unequivocal geochemical anomalies indicative of channelized fluid flow were found thermal measurements are not significantly different from those predicted by a purely conductive heat flow model a signature of low chloride pore water near the decollement may partly be related to smectite diagenesis but may also be due to episodic fluid flow events we conclude that dewatering probably occurred dominantly through diffuse flow throughout the accreted sediments at this site

spectrum of pore types and networks in mudrocks and a descriptive classification for matrix related mudrock pores matrix related pore networks in mudrocks are composed of nanometer to micrometer size pores in shale gas systems these pores along with natural fractures form the flow path permeability network that allows flow of gas from the mudrock to induced fractures during production a pore classification consisting of three major matrix related pore types is presented that can be used to quantify matrix related pores and relate them to pore networks two pore types are associated with the mineral matrix the third pore type is associated with organic matter om fracture pores are not controlled by individual matrix particles and are not part of this classification pores associated with mineral particles can be subdivided into interparticle interp pores that are found between particles and crystals and intraparticle intrap pores that are located within particles organic matter pores are intrap pores located within om interparticle mineral pores have a higher probability of being part of an effective pore network than intrap mineral pores because they are more likely to be interconnected although they are intrap om pores are also likely to be part of an interconnected network because of the interconnectivity of om particles in unlithifed near surface muds pores consist of interp and intrap pores and as the muds are buried they compact and lithify during the compaction process a large number of interp and intrap pores are destroyed especially in ductile grain rich muds compaction can decrease the pore volume up to 88 by several kilometers of burial at the onset of hydrocarbon thermal maturation om pores are created in kerogen at depth dissolution of chemically unstable particles can create additional moldic intrap pores

elevated pore pressure and anomalously low stress in regions of low frequency earthquakes along the nankai trough subduction megathrust 1 recent seismic reflection and ocean bottom seismometer obs studies reveal broad regions of low seismic velocity along the nankai subduction plate boundary megathrust offshore sw japan these low velocity zones lvz s extend 55 km landward from the trench corresponding to depths of 10 km below sea floor here we estimate the in situ pore pressure and stress state within these lvz s by combining p wave velocities obtained from the geophysical surveys with new well constrained empirical relations between p wave velocity porosity and effective mean stress defined by laboratory deformation tests on drill core samples of the incoming oceanic sediment we document excess pore pressures of 17 87 mpa that extend 55 km into the subduction zone indicating that trapped pore fluids support 45 91 of the overburden stress along the base of the upper plate and surrounding major fault zones the resulting effective stresses in the lvz are limited to 1 3 of the values expected for non overpressured conditions these low effective stresses should lead to a mechanically weak and predominantly aseismic plate boundary the region of lowest effective stress coincides with precisely located very low frequency earthquakes providing the first quantitative evidence linking these anomalous slip events to low stress and high pore pressure

strength characteristics of japan trench borehole samples in the high slip region of the 2011 tohoku oki earthquake abstract the 2011 tohoku oki earthquake demonstrated that the shallowest reaches of plate boundary subduction megathrusts can host substantial coseismic slip that generates large and destructive tsunamis contrary to the common assumption that the frictional properties of unconsolidated clay rich sediments at depths less than 5 km should inhibit rupture we report on laboratory shearing experiments at low sliding velocities 1 mm s using borehole samples recovered during iodp expedition 343 jfast spanning the plate boundary decollement within the region of large coseismic slip during the tohoku earthquake we show that at sub seismic slip rates the fault is weak sliding friction s 0 2 0 26 in contrast to the much stronger wall rocks s 0 5 the fault is weak due to elevated smectite clay content and is frictionally similar to a pelagic clay layer of similar composition the higher cohesion of intact wall rock samples coupled with their higher amorphous silica content suggests that the wall rock is stronger due to diagenetic cementation and low clay content our measurements also show that the strongly developed in situ fabric in the fault zone does not contribute to its frictional weakness but does lead to a near cohesionless fault zone which may facilitate rupture propagation by reducing shear strength and surface energy at the tip of the rupture front we suggest that the shallow rupture and large coseismic slip during the 2011 tohoku earthquake was facilitated by a weak and cohesionless fault combined with strong wall rocks that drive localized deformation within a narrow zone

porosity of sediments in accretionary prisms and some implications for dewatering processes rapid dewatering largely through loss of porosity is one of the most striking changes affecting deformed sediments in accretionary prisms analysis of porosity from sample and seismic data confirms that porosity in accretionary prisms is significantly less than that in basinal strata of similar lithology more specifically the porosity depth gradient increases arcward across the prism until porosities less than 10 occur within a few kilometers of the surface in deeper parts of the prism dewatering and volume reduction occur not only during loss of porosity but also during dehydration and metamorphic reactions calculations involving volume change lead us to define a function which reflects total volume reduction the distribution of the relative intensity of dewatering within the prism can be calculated from the porosity contours and trajectories of sediment elements an initial attempt indicates that in the toe region the rate of dewatering is accelerated at depth relative to the basinal section and may be greatest near the base of the prism toward the rear of the prism the maximum rate of dewatering occurs at progressively shallower depths several lines of evidence suggest that in addition to gradual diffuse dewatering localized zones of strong dewatering are associated with zones of intense shear this model of shear dewatering based on critical state soil mechanics suggests that sediments adjacent to the basal decollement should be strongly dewatered

elevated fluid pressure and extreme mechanical weakness of a plate boundary thrust nankai trough subduction zone pore fluid pressure is a key parameter governing both the shear strength of faults and the conditions for seismic and tsunamigenic slip on subduction zone megathrusts however quantitative constraints on this parameter based on in situ data in these or any faults have proven elusive using seismic interval velocities derived from a three dimensional seismic reflection survey we compute porosity fluid pressure and effective stress at the plate interface in the nankai trough subduction zone offshore southwestern japan to 20 km down the fault dip using well constrained locally calibrated empirical relationships between velocity porosity and effective stress we show that the fault and immediately subjacent footwall are nearly undrained implying that the subduction megathrust slips under a shear stress of

on the mechanical role of smectite in subduction zones smectites a group of hydrated phyllosilicate minerals are often a main component of sediments entering a subduction zone they are common in hemipelagic sediments and consist of an authigenetic component that arises from the alteration of volcanic ash intercalated in the sediment and a detrital component smectite rich strata are structurally weak retain water and are highly porous a review of sediment lithology in subduction zones around the world suggests that decollement zones may preferentially form in smectite rich horizons one consequence of this observation is that the assignment of high frictional resistance to materials in a decollement may be unwarranted at deeper structural levels smectite is transformed to illite a mechanically stronger mineral estimates of the depth of this transition coincide with the depth at which earthquakes are first detected along plate interfaces i speculate that the smectite illite transition in subduction zones may be an important component of sediment consolidation processes that trigger interplate seismicity

geological controls on matrix permeability of devonian gas shales in the horn river and liard basins northeastern british columbia canada abstract controls of matrix permeability are investigated for devonian gas shales from the horn river and liard basins in northeastern british columbia canada mineralogy is varied with high carbonate high quartz and moderate quartz carbonate and clay rich strata quartz content varies between 2 and 73 carbonate varies between 1 and 93 and clay varies between 3 and 33 the toc content ranges between 0 3 and 6 wt and porosity varies between about 1 and 7 for horn river basin samples quartz is mainly biogenic in origin derived from radiolarians toc content increases with the quartz content suggesting the toc and quartz both are derived from siliceous phytoplankton a positive relationship between porosity and quartz content is due to the positive relationship between quartz and toc matrix permeability parallel to bedding varies between 7 5e 02 and 7 1e 07 md at an effective stress of 15 mpa variation in permeability is due to a complex combination of factors that includes origin and distribution of minerals pore size distribution and fabric mercury intrusion capillary curves indicate that the higher matrix permeability values 2e 03 md occurs in samples that contain interconnected pore apertures greater than 16 m even when these samples may contain less macropores than low permeability samples the fabric of high permeability samples can be either isotropic or anisotropic however permeability of anisotropic samples is more sensitive to changes in effective stress than isotropic samples more highly anisotropic samples contain moderate amounts of quartz carbonate and in some clay high permeability samples that contain a more balanced ratio between micro meso and macroporosity would not only have faster flow rates but also greater access to sorbed gas within the microporosity compared to samples that lack mesopores several muskwa samples compared to evie and besa river samples contain higher quartz moderate clay and high toc content coupled with high permeability less sensitivity to effective stress and balanced ratios between micro meso and macroporosity would be a lower exploration risk due a greater propensity to fracture the ability to produce and store hydrocarbons due to higher toc contents and greater communication between macropores and micropores in the organic and clay fractions

critical taper model of fold and thrust belts and accretionary wedges the fold and thrust belts and submarine accretionary wedges that lie along compressive plate boundaries are one of the best understood deformational features of the earth s upper crust although there is considerable natural variation among the many fold and thrust belts and accretionary wedges that have been recognized and explored several features appear to be universal in cross section fold and thrust belts and accretionary wedges occupy a wedge shaped deformed region overlying a basal detachment or decollement fault the rocks or sediments beneath this fault show very little deformation the decollement fault characteristically dips toward the interior of the mountain belt or in the case of a submarine wedge toward the island arc the topography in contrast slopes toward the toe or deformation front of the wedge deformation within the wedge is generally dominated by imbricate thrust faults verging toward the toe and related fault bend folding two north american fold and thrust belts that exhibit these features are shown in figure 1 neither of these two examples is tectonically active today the southern canadian fold and thrust belt was active during the late jurassic and cretaceous 150 100 ma whereas the southern appa lachians were deformed during the late carboniferous to permian alle ghenian orogeny 300 250 ma figure 2 shows two examples that are currently active the taiwan fold and thrust belt produced by the sub duction of the eurasian plate beneath the philippine sea plate suppe 1981 1987 and the barbados accretionary wedge produced by the sub

effects of smectite to illite transformation on the frictional strength and sliding stability of intact marine mudstones 1 at subduction zones earthquake nucleation and coseismic slip occur only within a limited depth range known as the seismogenic zone one leading hypothesis for the upper aseismic seismic transition is that transformation of smectite to illite at 100 150 c triggers a change from rate strengthening frictional behavior that allows only stable sliding to rate weakening behavior considered a prerequisite for unstable slip previous studies on powdered gouges have shown that changes in clay mineralogy alone are unlikely to control this transition but associated fabric and cementation developed during diagenesis remain possible candidates we conducted shearing experiments designed specifically to evaluate this hypothesis by using intact wafers of mudstone from ocean drilling program site 1174 offshore sw japan which have undergone progressive smectite transformation in situ we sheared specimens along a sawcut in a triaxial configuration oriented parallel to bedding at normal stresses of 20 150 mpa and a pore pressure of 1 mpa during shearing we conducted velocity stepping tests to measure the friction rate parameter a b friction coefficient ranges from 0 28 0 40 and values of a b are uniformly positive both are independent of clay transformation progress our work represents the most direct and comprehensive test of the clay transformation hypothesis to date and suggests that neither illitization nor accompanying fabric development and cementation trigger a transition to unstable frictional behavior we suggest that strain localization in combination with precipitation of calcite and quartz is a viable alternative that is consistent with both field observations and recent conceptual models of a heterogeneous seismogenic zone

depth distribution of coseismic slip along the nankai trough japan from joint inversion of geodetic and tsunami data coseismic slip distribution on the fault plane particularly in the downdip direction associated with large subduction earthquakes can be estimated by joint inversion of geodetic and tsunami data two large earthquakes the 1944 tonankai earthquake mw 8 1 and the 1946 nankaido earthquake mw 8 3 occurred on the nankai trough southwestern japan where the philippine sea plate is subducting beneath the eurasian plate the source areas of these events extended over both land and ocean coseismic crustal movements on land were measured by leveling while those in ocean were recorded as tsunami waveforms on tide gauges the coseismic slip distribution inverted from these data shows that the slip on the shallower part of the fault plane is comparable to that on the deeper parts this indicates that large coseismic slip can occur beneath an accretionary wedge where current seismicity is low the result has implications for other subduction zones having a similar tectonic environment such as for the pacific northwest region of the united states although the possibility of a large earthquake there is still debated should a large subduction earthquake occur in this region coseismic slip on the shallow part could be large and the potential for large tsunamis is high

coupled fluid flow and deformation modeling of the frontal thrust region of the kumano basin transect japan implications for fluid pressures and decollement downstepping 1 results of coupled fluid flow and deformation modeling illustrate the influence of frontal thrust slip on excess pore pressures and decollement downstepping in the kumano basin transect offshore japan although permeable trench sediments and a buried sandstone layer provide potential fluid escape pathways the numerical modeling suggests that excess pore pressure could build beneath the frontal thrust within hemipelagic clay of the shikoku basin facies simulated excess pore pressures range between 45 and 77 of the overburden weight 3 km landward of the deformation front and are very sensitive to the permeability of the shikoku basin clay and the coefficient of friction of the fault zone excess pore pressures increase as the fault coefficient of friction is increased suggesting that volume change due to shear within the footwall sediments is an important contributor to excess pore pressures raising the coefficient of friction also increases shear within the sediments beneath the fault zone suggesting downward migration of the decollement zone this shearing focuses in the weaker shikoku basin clays rather than an underlying turbidite unit suggesting that lithologic differences play an important role in controlling deformation location

shear strength of sediments approaching subduction in the nankai trough japan as constraints on forearc mechanics 1 the mechanical behavior of the plate boundary fault zone is of paramount importance in subduction zones because it controls megathrust earthquake nucleation and propagation as well as the structural style of the forearc in the nankai area along the nantroseize kumano drilling transect offshore sw japan a heterogeneous sedimentary sequence overlying the oceanic crust enters the subduction zone in order to predict how variations in lithology and thus mechanical properties affect the formation and evolution of the plate boundary fault we conducted laboratory tests measuring the shear strengths of sediments approaching the trench covering each major lithological sedimentary unit we observe that shear strength increases nonlinearly with depth such that the apparent coefficient of friction decreases in combination with a critical taper analysis the results imply that the plate boundary position is located on the main frontal thrust further landward the plate boundary is expected to step down into progressively lower stratigraphic units assisted by moderately elevated pore pressures as seismogenic depths are approached the decollement may further step down to lower volcaniclastic or pelagic strata but this requires specific overpressure conditions high taper angle and elevated strengths in the toe region may be local features restricted to the kumano transect

unconventional shale gas systems the mississippian barnett shale of north central texas as one model for thermogenic shale gas assessment shale gas resource plays can be distinguished by gas type and system characteristics the newark east gas field located in the fort worth basin texas is defined by thermogenic gas production from low porosity and low permeability barnett shale the barnett shale gas system a self contained source reservoir system has generated large amounts of gas in the key productive areas because of various characteristics and processes including 1 excellent original organic richness and generation potential 2 primary and secondary cracking of kerogen and retained oil respectively 3 retention of oil for cracking to gas by adsorption 4 porosity resulting from organic matter decomposition and 5 brittle mineralogical composition the calculated total gas in place gip based on estimated ultimate recovery that is based on production profiles and operator estimates is about 204 bcf section 5 78 10 9 m 3 1 73 10 4 m 3 we estimate that the barnett shale has a total generation potential of about 609 bbl of oil equivalent ac ft or the equivalent of 3657 mcf ac ft 84 0 m 3 m 3 assuming a thickness of 350 ft 107 m and only sufficient hydrogen for partial cracking of retained oil to gas a total generation potential of 820 bcf section is estimated of this potential approximately 60 was expelled and the balance was retained for secondary cracking of oil to gas if sufficient thermal maturity was reached gas storage capacity of the barnett shale at typical reservoir pressure volume and temperature conditions and 6 porosity shows a maximum storage capacity of 540 mcf ac ft or 159 scf ton

geochemical and geological characteristics of the lower cambrian shales in the middle upper yangtze area of south china and their implication for the shale gas exploration abstract during the early cambrian the middle upper yangtze area of south china evolved to an open shelf with anoxic bottom waters prevailing in most places accumulated black shales in this area named niutitang formation qiongzhusi formation in some places thus become the prospective resources for shale gas exploration and development based on compiled petrological and organic geochemistry data and other information obtained by relevant laboratory tests and field investigation we documented regional distributional patterns of lower cambrian organic rich shales in the middle upper yangtze area according to their lateral extent thickness total organic carbon thermal maturity and burial depth this study also features the general reservoir properties of dark shales in this area together these results highlight several most favorable areas including youyang xiushan huayuan suiyang daozhen dejiang shiqian sinan yinjiang jiangkou tongren songtao and huangping shibing for shale gas exploration

geological and geochemical characteristics of marine continental transitional shale from the lower permian taiyuan formation taikang uplift southern north china basin abstract in recent years extensive studies have been conducted on conventional shale gas resources but little work has been performed on the resource potential of marine continental shales this paper takes the lower permian taiyuan shale in the taikang uplift of the southern north china basin as an example to characterize a transitional shale reservoir the results indicate that the taiyuan shale has high toc values averaging 2 67 and contains type iii gas prone organic matter the shale has high ro values averaging 3 44 and tmax ranges from 424 c to 601 c with respect to fracability or brittleness the content of brittle minerals in the taiyuan shale ranges from 21 to 71 averaging 45 55 and the average content of clay minerals is 44 1 based on sem and sem with ar ion milling observation interparticle pores intraparticle pores and microfractures are developed in the taiyuan shale which can provide reservoir space for the accumulation of shale gas organic matter pores rarely appear due to the organic type and high maturity the results of this study can guide the shale gas exploration of taiyuan formation in the taikang uplift of the southern north china basin

accretionary prisms in subduction earthquake cycles the theory of dynamic coulomb wedge 1 we expand the theory of critically tapered coulomb wedge for accretionary prisms by considering stress changes in subduction earthquake cycles building on the coulomb plasticity of the classical theory we assume an elastic perfectly coulomb plastic rheology and derive exact stress solutions for stable and critical wedges the new theory postulates that the actively deforming most seaward part of an accretionary prism the outer wedge overlies the updip velocity strengthening part of the subduction fault and the less deformed inner wedge overlies the velocity weakening part the seismogenic zone during great earthquakes the outer wedge is pushed into a compressively critical state with an increase in basal and internal stresses and pore fluid pressure after the earthquake the outer wedge returns to a stable state the outer wedge geometry is controlled by the peak stress of the updip velocity strengthening part of the subduction fault achieved in largest earthquakes the inner wedge generally stays in the stable regime throughout earthquake cycles acting as an apparent backstop and providing a stable environment for the formation of forearc basins the new theory has important implications for the studies of the updip limit of the seismogenic zone the evolution of accretionary prisms and forearc basins activation of splay faults and tsunami generation evolution of the fluid regime and mechanics of frontal prisms at margins dominated by tectonic erosion

heat flow and fluid flow regime in the western nankai accretionary prism the thermal structure of the western nankai subduction zone has been studied through surface heat flow measurements with conventional geothermal probes and estimation of heat flow from depths of gas hydrate bottom simulating reflectors temperature and thermal conductivity measurements at odp site 808 located at the toe of the nankai accretionary prism added important information on any variation in heat flow with depth and thus on thermal effects of pore fluid flow in the prism surface heat flow is generally higher than 130 mw m 2 on the floor of the nankai trough such high heat flow cannot be explained simply by conductive cooling of the subducting shikoku basin lithosphere recent closely spaced probe measurements revealed that the heat flow on the trough floor shows local variation with particularly high values near the deformation front the heat flow distribution suggests that warm fluid flows seaward possibly along the decollement and locally rises toward the seafloor along some channels resulting in both the regional and local high heat flow anomalies contrary to this inference no active localized fluid flow was detected at site 808 the temperatures measured at six depths down to 347 m below seafloor mbsf just below the frontal thrust can be well explained by conductive heat flow of 125 130 mw m2 no geochemical anomaly associated with present day pore fluid movement was observed in cores recovered from the frontal thrust or the decollement on the other hand the high heat flow at site 808 requires some heat source below 350 mbsf there is a possibility that the fluid flow in the nankai accretionary prism is localized and transient and thus was not detected at this single drill site

fracture development in paleozoic shale of chongqing area south china part one fracture characteristics and comparative analysis of main controlling factors abstract natural macro and micro fractures in the longmaxi shale lower silurian from yuye 1 well and the niutitang shale lower cambrian from yuke 1 well are characterized regional geological structural evolution and shale geochemistry including toc and mineral composition etc are analyzed the controlling factors of natural fractures in palaeozoic shale in the southeast of chongqing have been confirmed which have been analyzed qualitatively or semi quantitatively the results show that fractures are abundant in both the longmaxi shale and niutitang shale due to the multiphase strong tectonic movement in the southeast of chongqing therefore tectonic factors are the principal factors that control fracture development in the paleozoic shale of this region followed by the mineral composition of the shale vertically the mineral composition of shale in yuke 1 well has a high dispersion degree and a strong heterogeneity the dispersion coefficients of quartz and dolomite were 0 48 and 1 respectively therefore the development and distribution of fractures are controlled by quartz and dolomite content in the niutitang shale moreover the quartz content shows positive correlation with the fractures quantity and fractures are most abundant where dolomite content up to 64 7 while the mineral composition of shale in yuye 1 well has a low dispersion degree and a strong homogeneity in a vertical position the dispersion coefficients of quartz and plagioclase were 0 15 and 0 26 respectively which shows unconspicuous control effect for fractures therefore the development and distribution of fractures are mainly controlled by the toc content moreover the toc content shows positive correlation with the fractures quantity fractures are sealed by calcite quartz and pyrite although sealed fractures do not contribute to the permeability of the reservoir they are important planes of weakness that tend to be reactivated by hydraulically induced fractures sem images show that intergranular pore is connected by microfractures thereby the validity of porosity has been improved the mixed network of microfractures and micropore provide enhanced permeability or storage capacity for the reservoir

mineralogy and trace element geochemistry of gas shales in the united states environmental implications abstract this paper presents a compilation of published mineralogic and trace element data from nine gas shales in the united states formations analyzed include the antrim bakken barnett eagle ford haynesville marcellus new albany utica and woodford these mineralogic and trace element data can be used to assess the potential for environmental impacts during hydraulic fracturing impacts addressed in this study include 1 the potential for acid rock drainage generation during gas shale weathering 2 the distribution of trace elements in gas shales and comparison with regulatory guidelines and 3 the implications for environmental management of well cuttings the use of the mineralogic data to assess the fracability of the gas shales is also considered compilations of the mineralogy and geochemistry of gas shales can be a valuable resource for managing real and perceived environmental problems associated with their exploitation comprehensive environmental assessment to fully address these issues in addition to other potential environmental impacts will require collection and collation of additional data on the mineralogy and trace element geochemistry of gas and other hydrocarbon producing shales

a note on the frictional strength of laumontite from cajon pass california laumontite mineralization is pervasive in joints and shear zones encountered in the cajon pass drillhole in southern california in order to determine whether a gouge composed of this hydrated mineral affects shear strength in a manner similar to low strength day rich fault gouges frictional sliding experiments were performed under dry saturated and high pore pressure conditions at effective pressures up to 450 mpa coefficients of friction ranged between 0 66 and 0 84 consistent with most crustal rocks and well above the values typical of clay rich san andreas fault gouges saturation state had no effect on strength or sliding stability these results suggest that the presence of laumontite in shear zones at cajon pass will not affect the shear strength of the rock in a way that can account for the inferred low ambient shear stresses

evolution of the nankai trough decollement from the trench into the seismogenic zone inferences from three dimensional seismic reflection imaging we mapped the amplitude of the nankai trough subduction thrust seismic reflection from the trench into the seismogenic zone with three dimensional seismic reflection data the decollement thrust forms within the lithologically homogeneous lower shikoku basin facies along an initially nonreflective interface the reflection develops from a porosity contrast between accreted and underthrust sedimentary material because of accretionary wedge consolidation and rapid loading and delayed consolidation of the underthrust section a decollement amplitude map shows a significant decline from high amplitudes at the trench to barely detectable levels 25 30 km landward three other observations coincide with the amplitude decline 1 the decollement initially steps down to deeper stratigraphic levels 2 the wedge taper increases dramatically and 3 the thrust becomes seismogenic the amplitude decline and the coincident decollement and accretionary wedge tectonic and seismogenic behavior are attributed to the loss of fluids and potentially loss of excess fluid pressures downdip along the subduction thrust

laboratory results indicating complex and potentially unstable frictional behavior of smectite clay a central problem in explaining the apparent weakness of the san andreas and other plate boundary faults has been identifying candidate fault zone materials that are both weak and capable of hosting earthquake like unstable rupture our results demonstrate that smectite clay can be both weak and velocity weakening at low normal stress 30 mpa our data are consistent with previous work which has focused on higher normal stress conditions 50 mpa and greater and found only velocity strengthening if natural fault zones contain significant smectite one key implication of our results is that localized zones of high pore pressure which reduce effective normal stress could be important in controlling potential sites of earthquake nucleation our experiments indicate that friction of smectite is complex and depends upon both sliding velocity and normal stress this complexity highlights the need for detailed experiments that reflect in situ conditions for fault gouges

structure of the nankai trough accretionary zone from multichannel seismic reflection data new multichannel seismic reflection data collected over the nankai trough image the accretionary complex in two areas the international program of ocean drilling leg 87 transect area western area and the region of upcoming ocean drilling program leg 131 eastern area the incoming shikoku basin sedimentary section consists of hemipelagic clays and thin terrigenous turbidites the basin section is overlain by a trench wedge that is 12 16 km wide and 350 750 m thick at the thrust front accretionary deformation begins in a protothrust zone that is characterized by thickening and seaward tilting of the trench wedge the zone in the western area is 4 5 km wide and is characterized by kink folds the zone in the eastern area is only 2 5 km wide and does not exhibit such folds the frontal thrusts in each area are imaged as fault plane reflections and ramp upward from within the basin hemipelagic section the overthrusting sediments form fault bend folds over these ramps thrust spacing at the toe of the slope is 1 5 2 5 km the second thrust cuts up from an inferred decollement within the shikoku basin sedimentary section in the eastern area a reflection marking the top of the basin pelagic sediment section changes from normal to reversed polarity about 6 3 km seaward of the thrust front and underlies the entire protothrust zone this reflector continues with reversed polarity under the accretionary complex and is at the level of the basal decollement the underlying basin pelagic section is apparently thrust undisturbed beneath the accretionary prism the reversal of polarity indicates a change in reflection coefficient that is due to a combination of decreasing seismic velocity and density across the interface this decrease in velocity and density may indicate that the decollement is a zone of high porosity due to fluid expulsion from deeper within the accretionary prism the reflections from the first and second thrusts are also reversed polarity possibly indicating that they also are pathways of fluid expulsion the critical wedge taper of the western area is greater than that of the eastern area an observation that is consistent with the existence of an overpressured decollement in the eastern area

accretion and erosion in subduction zones the role of fluids convergence occurs today along about 60 000 km of plate boundaries parsons 1981 in this paper we consider frontal subduction which occurs along about 40 000 km of convergent plate boundaries using the kine matic solution of de mets et al 1990 and including 0 5 km 2 yr1 of marginal basin openings otsuki 1989 we find that 2 95 km 2 per year disappear in frontal subduction at an average rate of 7 3 cm yr the upper overriding plate is significantly modified not only through the addition of a large amount of magmatic material markhinin 1968 but also through tectonic processes which result either in outgrowth of the margin through addit ion of sediments silver 1969 or in recession through frontal and subcrustal erosion which appears to be linked to sediment subduction coats 1962 gilluly 1963 the various phenomena observed must in some way be related to variations in mechanical coupling between the two plates uyeda kanamori 1979 had proposed that the transition from erosion to accretion results from increasing mechanical coupling the marianas trench was presented as typical of an crosional margin and the chile trench of an accreting margin ruff kanamori 1980 and kana mori 1986 showed that the magnitude of the subduction earthquakes increases with the convergence velocity and the inverse of the age of the subducted lithosphere this correlation was qualitatively explained by the control of age on the d ip of the platc and of convergence velocity on the

abundance of smectite and the location of a plate boundary fault barbados accretionary prism isolating the respective roles of factors responsible for the initiation and localization of fault zones remains one of the more important goals of research in neotectonics the frontal decollement zone of the northern barbados accretionary prism provides several important clues as to how the interwoven variables of clay mineralogy fluid flow chemical interactions and sediment physical properties affect strain localization this plate boundary fault is centered at a lithologic contact between miocene claystone and fine grained oligocene turbidites the fault zone is nearly 40 m thick at ocean drilling program site 948 and its upper part passes through smectite rich deposits a sharp minimum in percent smectite and a maximum in percent illite mark the base of the decollement there is a consistent increase in percent smectite with distance above the base of the decollement but the top of the decollement is poorly defined by clay mineralogy the intrinsic mechanical weakness of strata with abundant smectite group clays probably influences where the fault tip propagates into the undeformed stratigraphy of the atlantic abyssal plain a second inherited parameter is the local abundance of radiolarians which contribute to higher than normal porosities sediment shear strength also decreases because pore pressure within the fault zone is significantly greater than hydrostatic the principal cause of excess pore pressure seems to be updip fluid advection in theory however decreases in pore fluid salinity and porosity collapse should increase the amount of physicochemical stress generated by expandable clay minerals the imported fluid is unusually low in salinity because it has migrated from zones of deeper seated dehydration reactions if fresher pore water migrates to the propagating tip of the decollement its arrival should increase smectite swelling and reduce the shear strength of the mudstone even more the location and evolution of the decollement therefore are controlled by a complicated interplay of static factors inherited from the abyssal atlantic stratigraphy and dynamic factors associated with episodic fluid flow and changing fluid chemistry

quantifying the mediterranean freshwater budget throughout the late miocene new implications for sapropel formation and the messinian salinity crisis the cyclic sedimentary record of the late miocene mediterranean shows a clear transition from open marine to restricted conditions and finally to evaporitic environments associated with the messinian salinity crisis this evolution has been attributed to changes in mediterranean atlantic connectivity and regional climate which has a strong precessional pulse 31 coupled climate simulations with different orbital configurations have been combined in a regression model that estimates the evolution of the freshwater budget of the mediterranean throughout the late miocene the study suggests that wetter conditions occur at precession minima and are enhanced at eccentricity maxima we use the wetter peaks to predict synthetic sapropel records using these to retune two mediterranean sediment successions indicates that the overall net freshwater budget is the most likely mechanism driving sapropel formation in the late miocene our sapropel timing is offset from precession minima and boreal summer insolation maxima during low eccentricity if the present day drainage configuration across north africa is used this phase offset is removed if at least 50 more water drained into the mediterranean during the late miocene capturing additional north african monsoon precipitation for example via the chad eosahabi catchment in libya in contrast with the clear expression of precession and eccentricity in the model results obliquity which is visible in the sapropel record during minimum eccentricity does not have a strong signal in our model by exploring the freshwater evolution curve in a box model that also includes mediterranean atlantic exchange we are able for the first time to estimate the mediterranean s salinity evolution which is quantitatively consistent with precessional control additionally we separate and quantify the distinct contributions regional climate and tectonic restriction make to the lithological changes associated with the messinian salinity crisis the novel methodology and results of this study have numerous potential applications to other regions and geological scenarios as well as to astronomical tuning

intraoceanic thrusts in the nankai trough off the kii peninsula implications for intraplate earthquakes 1 we identified intraoceanic thrusts developed as imbricate structures within the subducting philippine sea plate off the kii peninsula in central japan manifesting as strong amplitude reflections observed in an industry standard three dimensional 3d seismic reflection data set these imbricate intraoceanic thrusts cut through the oceanic crust as a discontinuous thrust plane striking approximately parallel to the trench in our survey area large intraplate earthquakes with moment magnitudes mw over 7 occurred on 5 september 2004 causing strong ground motions on the islands of japan and tsunami waves the locations of the intraoceanic thrusts recognized in the seismic data are distributed around the estimated hypocenters of the mainshocks and aftershocks of the 2004 earthquakes furthermore their geometry extracted from the 3d seismic data could explain the kind of complex rupture pattern observed during the 2004 events therefore we propose that the intraoceanic thrusts are seismogenically active

updip limit of the seismogenic zone beneath the accretionary prism of southwest japan an effect of diagenetic to low grade metamorphic processes and increasing effective stress off southwest japan the seaward limit of coseismic displacement or updip limit of the seismogenic zone of the 1946 mw 8 3 thrust earthquake reaches to 4 km depth and 40 km landward of the trench this limit coincides with the estimated location of the 150 c isotherm and has been linked to changes in physical properties associated with the smectite to illite clay mineral transition here we show that this limit correlates with a suite of diagenetic to low grade metamorphic processes characterized by 1 declining fluid production and decreasing fluid pressure ratio and 2 active clay carbonate and zeolite cementation and the transition to pressure solution and quartz cementation these diagenetic to low grade metamorphic changes cause the onset of velocity weakening during thrust faulting an increase in effective stress and strengthening of the hanging wall which together combine to produce recordable earthquakes

effect of hydration state on the frictional properties of montmorillonite based fault gouge 1 we report on laboratory experiments examining the effect of hydration state on the frictional properties of simulated clay and quartz fault gouge we tested four mixtures of ca montmorillonite and quartz 100 70 50 and 30 montmorillonite at four hydration states dry 16 0 wt water we controlled the hydration state using either oven drying for 13 wt h2o under conditions of controlled relative humidity for each clay quartz mixture and hydration state we measured frictional properties over a range of normal stresses 5 100 mpa and sliding velocities 1 300 m s we observe a systematic decrease in the coefficient of friction with increasing water content normal stress and clay content values of for 50 50 mixtures range from 0 57 to 0 64 dry and decrease to 0 21 0 55 for the most hydrated cases wet for layers of 100 montmorillonite ranges from 0 41 0 62 dry to 0 03 0 29 wet as water content is increased from 0 to 20 0 wt the friction rate parameter a b becomes increasingly positive variation in a b values decreases dramatically as normal stress increases if our experimental results can be applied to natural fault gouge the combination of stress state hydration state and quartz content that facilitates unstable fault behavior implies that the onset of shallow seismicity in subduction zones is more complicated than a simple transition from smectite to illite

controls on methane sorption capacity of indian coals details of methane sorption properties of some indian coals ranging in rank from high volatile bituminous c 0 62 ro max to medium volatile bituminous coal 1 46 ro max were investigated to determine the major factors influencing gas adsorption capacity and desorption rate variables studied included moisture content mineral matter content rank and coal type maceral composition adsorption isotherm analysis of dry coals showed that adsorption capacity followed a second order polynomial trend with rank equilibrium moist samples showed a linear increase in adsorption capacity with rank and had a significantly reduced adsorption capacity compared to the dry coals with the reduction being related to the moisture content mineral matter acted as a simple diluent to the gas adsorption capacity of the coals and was found to be nonadsorbent adsorption capacity moist was reduced by 0 32 cm3 g 10 ft3 gas t for every 1 increase in the ash yield bright and dull coal lithotypes showed strong separation in their adsorption capacities on an as analyzed basis with the bright coals adsorbing greater quantities of gas on a dmmf basis however no relationship was observed between coal type and gas storage capacity comparison of maceral composition with adsorption capacity dry confirmed this observation effective diffusivity de of methane through the coal was seen to be affected by coal type and rank effective diffusivity de decreased as rank increased which is related to the increasing microporous nature of the coal bulk coals tested had 2 3 times larger effective diffusivities than bright coals and dull coals had intermediate rates the larger de values for the bulk coals may be related to the presence of collodetrinite coupled with mineral matter which acts as higher permeability pathways for the gas compared with the high ash dull coals dominated by inertinite macerals differences were noted with similar age permian australian coals where methane adsorption capacities were larger and coal type had a stronger influence on adsorption capacity begin page 202

friction experiments on saturated sediments and their implications for the stress state of the nankai and barbados subduction thrusts the nankai and barbados forearcs have low stress subduction thrusts the sediments entering the subduction zone and namely the material in the decollement zones have been well characterized by numerous deep sea drilling legs and studies of the recovered cores nankai has high heat flow and significant amounts of illite while barbados is a smectite dominated system based on results from ring shear 2 mpa normal stress and direct shear 30 mpa tests on marine sediments and mineral standards this translates into a residual frictional resistance of r 0 25 and r 0 11 in clay horizons respectively such values agree with theoretical estimates from critical wedge theory nankai b 0 16 0 26 and barbados b 0 06 0 09 and fault spacing geometries from seismic profiles nankai b 0 12 0 23 and barbados b 0 11 0 19 maximum pore pressure ratios of 0 85 and 0 73 for nankai and barbados respectively allow us to estimate effective shear stresses as a function of friction coefficient and density of the sediment gouge to reach only 10 mpa or less in the frontal 50 km of the decollement zone respectively our data support the contention that fluid pressure transients and sediment composition contribute equally to the weakness along plate boundary faults down to the seismogenic zone with the first probably dominating the shallow decollement shear velocity stepping tests show that the clay dominated gouges strengthen velocity irrespective of the clay mineralogy and hence suggest that clay transformation does not affect the onset of seismogenesis

transition of accretionary wedge structures around the up dip limit of the seismogenic subduction zone the nankai accretionary prism is divided into three segments outer and inner wedges and their transition zone these wedges reflect different aspects of wedge taper internal deformation and basal plate boundary fault the outer wedge is characterized by narrow critical taper internal deformation by in sequence fold and thrust and aseismic decollement the inner wedge represents a stable narrow taper weakly deformed internal structure and seismogenic plate boundary fault along its base the transition zone between the two wedges shows large critical taper with steep surface slope internal structure of out of sequence thrust and step down of decollement onto the sediment oceanic basement interface the trench slope break and oceanward margin of forearc basin is located around the landward edge of this transition zone these common aspects might be related to the lithification of both accreted and underthrust sediments and the resultant switch of the plate boundary fault deformation and lithification process recorded in exhumed on land melange of accretionary complexes suggest that the step down of the plate boundary decollement occurs around the up dip limit of seismogenic subduction zone

organic rich marcellus shale lithofacies modeling and distribution pattern analysis in the appalachian basin the marcellus shale is considered to be the largest unconventional shale gas resource in the united states two critical factors for unconventional shale reservoirs are the response of a unit to hydraulic fracture stimulation and gas content the fracture attributes reflect the geomechanical properties of the rocks which are partly related to rock mineralogy the natural gas content of a shale reservoir rock is strongly linked to organic matter content measured by total organic carbon toc a mudstone lithofacies is a vertically and laterally continuous zone with similar mineral composition rock geomechanical properties and toc content core log and seismic data were used to build a three dimensional 3 d mudrock lithofacies model from core to wells and finally to regional scale an artificial neural network was used for lithofacies prediction eight petrophysical parameters derived from conventional logs were determined as critical inputs advanced logs such as pulsed neutron spectroscopy with log determined mineral composition and toc data were used to improve and confirm the quantitative relationship between conventional logs and lithofacies sequential indicator simulation performed well for 3 d modeling of marcellus shale lithofacies the interplay of dilution by terrigenous detritus organic matter productivity and organic matter preservation and decomposition affected the distribution of marcellus shale lithofacies distribution which may be attributed to water depth and the distance to shoreline the trend of normalized average gas production rate from horizontal wells supported our approach to modeling marcellus shale lithofacies the proposed 3 d modeling approach may be helpful for optimizing the design of horizontal well trajectories and hydraulic fracture stimulation strategies

the effect of mineralogy and effective normal stress on frictional strength of sheet silicates abstract phyllosilicates are common constituents of upper crustal faults and subduction forearcs we studied the effect of mineralogy and controlled effective normal stress between 5 mpa and 100 mpa on frictional strength although the effect of mineralogy on frictional strength of single phase phyllosilicate gouges has been previously studied the influence of effective normal stress has not we conducted water saturated and vacuum dry frictional tests on single phase phyllosilicate gouges using a triaxial apparatus minerals included talc pyrophyllite kaolinite lizardite illite montmorillonite chlorite muscovite phlogopite and biotite particle size 20 mpa dry friction coefficients are constant and wet friction coefficients show a small increase this is attributed to loss of water and increased contact area the results indicate that frictional strength of clay rich faults increases at depths less than 1 km under hydrostatic pore fluid pressures

do subducting seamounts generate or stop large earthquakes seamount subduction is a common process in subduction zone tectonics contradicting a widely held expectation that subducting seamounts generate large earthquakes seamounts subduct largely aseismically producing numerous small earthquakes on rare occasions when they do produce relatively large events the ruptures tend to be complex suggesting multiple rupture patches or faults we explain that the seismogenic behavior of these seamounts is controlled by the development and evolution of an adjacent fracture network during subduction and cannot be described using the frictional behavior of a single fault the complex structure and heterogeneous stresses of this network provide a favorable condition for aseismic creep and small earthquakes but an unfavorable condition for the generation and propagation of large ruptures

frontal wedge deformation near the source region of the 2011 tohoku oki earthquake 1 we report an uplift of 5 m with a horizontal displacement of more than 60 m due to the 2011 tohoku oki earthquake the uplift was measured by an ocean bottom pressure gauge installed before the earthquake on a frontal wedge which formed an uplift system near the japan trench horizontal displacements of the frontal wedge were measured using local benchmark displacements obtained by acoustic ranging before and after the earthquake the average displacements at the frontal wedge were 58 m east and 74 m east southeast these results strongly suggest a huge coseismic slip beneath the frontal wedge on the plate boundary the estimated magnitude of the slip along the main fault was 80 m near the trench our results suggest that the horizontal and vertical deformations of the frontal wedge due to the slip generated the tremendous tsunami that struck the coastal area of northeastern japan

mechanisms of shale gas storage implications for shale gas exploration in china this article reviews the mechanisms of shale gas storage and discusses the major risks or uncertainties for shale gas exploration in china at a given temperature and pressure the gas sorption capacities of organic rich shales are primarily controlled by the organic matter richness but may be significantly influenced by the type and maturity of the organic matter mineral composition especially clay content moisture content pore volume and structure resulting in different ratios of gas sorption capacity gsc to total organic carbon content for different shales in laboratory experiments the gsc of organic rich shales increases with increasing pressure and decreases with increasing temperature under geologic conditions assuming hydrostatic pressure gradient and constant thermal gradient the gsc increases initially with depth due to the predominating effect of pressure passes through a maximum and then decreases because of the influence of increasing temperature at greater depth this pattern of variation is quite similar to that observed for coals and is of great significance for understanding the changes in gsc of organic rich shales over geologic time as a function of burial history at an elevated temperature and pressure and with the presence of moisture the gas sorption capacities of organic rich shales are quite low as a result adsorption alone cannot protect sufficient gas for high maturity organic rich shales to be commercial gas reservoirs two models are proposed to predict the variation of gsc and total gas content over geologic time as a function of burial history high contents of free gas in organic rich shales can be preserved in relatively closed systems loss of free gas during postgeneration uplift and erosion may result in undersaturation the total gas contents lower than the sorption capacity and is the major risk for gas exploration in marine organic rich shales in china

comparison of frictional strength and velocity dependence between fault zones in the nankai accretionary complex accretionary complexes host a variety of fault zones that accommodate plate convergence and internal prism deformation including the decollement imbricate thrusts and out of sequence thrusts or splays these faults especially the decollement and major splay faults are considered to be candidates for hosting slow slip events and large magnitude earthquakes but it is not clear what modes of slip should be expected at shallow levels or how they are related to fault rock frictional properties we conducted laboratory experiments to measure the frictional properties of fault and wall rock from three distinct fault zone systems sampled during integrated ocean drilling program expedition 316 and ocean drilling program leg 190 to the nankai trough offshore japan these are 1 a major out of sequence thrust fault termed the megasplay site c0004 2 the frontal thrust zone a region of diffuse thrust faulting near the trench site c0007 and 3 the decollement zone sampled 2 km from the trench site 1174 at 25 mpa effective normal stress at slip rates of 0 03 100 m s and in the presence of brine as a pore fluid we observe low friction 0 46 for all of the materials we tested however the weakest samples 0 30 are from the decollement zone material from the megasplay fault is significantly weaker than the surrounding wall rocks a pattern not observed in the frontal thrust and decollement all samples exhibit primarily velocity strengthening frictional behavior suggesting that earthquakes should not nucleate at these depths a consistent minimum in the friction rate parameter a b at sliding velocities of 1 3 m s 0 1 0 3 m d is observed at all three sites suggesting that these shallow fault zones may be likely to host slow slip events

origin of a zone of anomalously high porosity in the subduction inputs to nankai trough abstract one poorly understood feature of the subduction inputs to the nankai trough subduction zone sw japan is a stratigraphic interval with an anomalously high porosity zone hpz which is up to 240 m thick and located within the clay and volcanic ash rich shikoku basin facies to investigate the origin of their peculiar physical properties we integrated logging while drilling lwd data shipboard density measurements and visual descriptions of core samples recovered from four drill sites of the ocean drilling program and integrated ocean drilling program we combined those observations with scanning electron microscopy sem and laboratory consolidation tests on both hpz samples and artificial mixtures of ash glass shards smectite and vesicular pumice smectite lwd data indicate that the hpz mudstones have a large proportion of dispersed volcanic ash 20 30 the consolidation tests show that the rate of porosity loss with increasing effective stress consolidation behavior is consistent among hpz specimens and matches artificial mixtures containing up to 60 volcanic material however absolute values of porosity remain higher for hpz samples compared to artificial mixtures so processes in addition to the mechanical effects of volcanic clasts must be contributing we suggest that hydration and partial dissolution turns clusters of volcanic glass into aggregates with distinct microfabric sem images confirm the presence of strengthened grain to grain contacts which probably inhibit the collapse of the intra aggregate pore space the aggregates behave like grains so that cohesive strength of the bulk sediment and rate of porosity loss remain nearly unchanged during burial the two step diagenetic process of dissolution and precipitation depends critically on a threshold abundance of fine grained dispersed volcanic ash pumice older units in the shikoku basin with only traces of dispersed ash show no such effects hpzs should be expected in other subduction zones with similar compositional and diagenetic prerequisites

crustal velocity field of southwest japan subduction and arc arc collision we investigate crustal deformation in southwest japan over a 3 year period revealed by a permanent dense global positioning system gps array southwest japan is a part of the amurian plate a microplate moving about 10 mm yr toward the east with respect to the eurasian plate it overrides the philippine sea plate at the nankai trough and collides with the northeast japan arc in the central part of japan in this paper we first derive gps site velocities relative to the stable part of the amurian plate in order to isolate signals of crustal deformation caused by the subduction and or the collision the velocity field has a conspicuous feature indicating the interseismic elastic loading by the philippine sea plate slab at the nankai trough characterized by the northwestward movements of points throughout the studied area their amplitudes are the largest at the pacific coast and decay toward the japan sea coast with a subtle systematic shift of azimuths a model assuming an elastic half space the convergence rate at the nankai trough based on a refined euler vector and the strength of the coupling inferred from a thermal model could explain the velocity field in the western part of the studied area to a large extent those in the eastern part systematically deviate from them and the residual components there show east west shortening and north south extension this may represent crustal thickening and trenchward extrusion of crustal blocks caused by the collision between southwest and northeast japan this suggests that the collision between southwest and northeast japan gives rise to not only crustal thickening but also trenchward extrusion of crustal block a velocity contrast was found across the median tectonic line the largest inland active fault in japan but the current permanent gps network is not dense enough for us to discuss its coupling depth

splay fault slip in a subduction margin a new model of evolution in subduction zones major thrusts called splay faults are thought to slip coseismically during large earthquakes affecting the main plate interface we propose an analytical condition for the activation of a splay fault based on force balance calculations and suggest thrusting along the splay fault is generally conditioned by the growth of the accretionary wedge or by the erosion of the hanging wall in theory normal slip on the splay fault may occur when the decollement has a very low friction coefficient seaward such a low friction also implies an unstable extensional state within the outer wedge finite element elasto plastic calculations with a geometry based on the nankai kumano section were performed and confirm that this analytical condition is a valid approximation furthermore localized extension at a shallow level in the splay hanging wall is observed in models for a wide range of friction coefficients from 0 to the value of internal friction coefficient of the rock here equals to 0 4 the timing of slip established for the splay fault branch drilled on nankai kumano transect suggests a phase of concurrent splay and accretionary wedge growth 2 ma to 1 5 ma followed by a locking of the splay 1 3 ma active extension is observed in the hanging wall this evolution can be explained by the activation of a deeper and weaker decollement followed by an interruption of accretion activation of a splay as a normal fault as hypothesized in the case of the tohoku 2011 earthquake can be achieved only if the friction coefficient on the decollement drops to near zero we conclude that the tectonic stress state largely determines long term variations of tightly related splay fault and outer decollement activity and thus influences where and how coseismic rupture ends but that occurrence of normal slip on a splay fault requires coseismic friction reduction

a deep strong reflector in the nankai accretionary wedge from multichannel seismic data implications for underplating and interseismic shear stress release 1 we present the results of structural and stratigraphic interpretations for the multichannel seismic profiles to reveal the structure of the nankai accretionary wedge off shikoku island on the basis of reflection characteristics we identify three major seismic units that correspond to a unit consisting of both oligocene miocene accretionary prism and the sedimentary cover unit a a unit of underthrust miocene pliocene hemipelagic sediments unit b and a unit of subducting pre miocene oceanic crustal rocks unit c according to relative strength of the topmost reflector of unit b the nankai accretionary wedge can be divided into three different structural zones along the cross strike direction apparent decollement reflector zone low amplitude reflector lar zone and deep strong reflector dsr zone despite the velocity uncertainty involved in the depth conversion and therefore thickness estimation error the decollement step down and the possible duplex structures in the lar and dsr zones support that the landward thickness variation of unit b is caused by the structural thickening such as the underplating process structural features and reverse polarity suggest that the dsr as an upper interface of the underplated layer unit b corresponds to a shear plane or a fault zone with high pore fluid pressures we found that most of the dsr was observed in a region which has not suffered from the 1946 coseismic dislocation even though it was locked during the interseismic period the dsr as a shear plane or a fault zone constrains a regional variety of the coseismic dislocation around the updip limit of the nankai seismogenic zone considering all of the dsr characteristics we propose a hypothesis for the dsr s development and role the dsr may play an important role as a possible releaser of the shear stress energy stored at the plate boundary around the updip limit of the nankai seismogenic zone a good correlation between the theoretical vertical displacement and seafloor topographic feature support the possible slip along the dsr as a shear plane

shallow subduction zone earthquakes and their tsunamigenic potential summary we have examined the source spectra of all shallow subduction zone earthquakes from 1992 to 1996 with moment magnitude 7 0 or greater as well as some other interesting events in the period range 1 20 s by computing moment rate functions of teleseismic p waves after comparing the source spectral characteristics of tsunami earthquakes earthquakes that are followed by tsunamis greater than would be expected from their moment magnitude with regular events we identified a subclass of this group slow tsunami earthquakes this subclass consists of the 1992 nicaragua the 1994 java and the february 1996 peru earthquakes we found that these events have an anomalously low energy release in the 1 20 s frequency band with respect to their moment magnitude although their spectral drop ov is comparable to those of the other earthquakes from an investigation of the centroid and body wave locations it appears that most earthquakes in this study conformed to a simple model in which the earthquake nucleates in a zone of compacted and dehydrated sediments and ruptures up dip until the stable sliding friction regime of unconsolidated sediments stops the propagation sedimentstarved trenches e g near jalisco can produce very shallow slip because the fault material supports unstable sliding the slow tsunami earthquakes also ruptured up dip however their centroid is located unusually close to the trench axis the subduction zones in which these events occurred all have a small accretionary prism and a thin layer of subducting sediment ocean surveys show that in these regions the ocean floor close to the trench is highly faulted we suggest that the horst and graben structure of a rough subducting oceanic plate will cause contact zones with the overriding plate making shallow earthquake nucleation and up dip propagation to the ocean floor possible the rupture partly propagates in sediments making the earthquake source process slow two factors have to be considered in the high tsunami generating potential of these events first the slip propagates to shallow depths in low rigidity material causing great deformation and displacement of a large volume of water second the measured seismic moment may not represent the true earthquake displacement because the elastic constants of the source region are not taken into account in the standard cmt determination

broad weak regions of the nankai megathrust and implications for shallow coseismic slip abstract deep within the nankai trough subduction zone the plate boundary thrust slips along a well imaged megasplay fault system during the megathrust earthquakes that regularly strike southwest japan the routing of the active plate boundary thrust along an upward branching splay fault causes deep underthrusting of an unusually thick section of material attached to the subducting ocean crust here we present three dimensional seismic reflection data that shows this unusually thick section is fluid rich sediment that results in broad weakly coupled regions of the megathrust down into the updip end of the seismogenic zone the weakly coupled regions lie above an underthrust low seismic impedance presumed to be low velocity and low density sediment section that is between one and two kilometers thick and covers 3300 km2 at least one eighth of the total rupture area of the 1944 mw 8 1 tonankai earthquake this underthrust sediment section is broader and much deeper than inferred at other margins sediment underthrusting into the normally well coupled seismogenic zone likely releases fluids and elevates fluid pressure which reduces inter plate coupling along portions of the megasplay fault and allows coseismic rupture to propagate to unusually shallow depths and generate large tsunami as inferred for the 1944 tonankai event splay faults may be a common yet transient mechanism for developing weak subduction zone thrusts

interactions between deformation and fluids in the frontal thrust region of the nantroseize transect offshore the kii peninsula japan results from iodp expedition 316 sites c0006 and c0007 integrated ocean drilling program iodp expedition 316 sites c0006 and c0007 examined the deformation front of the nankai accretionary prism offshore the kii peninsula japan in the drilling area the frontal thrust shows unusual behavior as compared to other regions of the nankai trough drilling results integrated with observations from seismic reflection profiles suggest that the frontal thrust has been active since 0 78 0 436 ma and accommodated 13 to 34 of the estimated plate convergence during that time the remainder has likely been distributed among out of sequence thrusts further landward and or accommodated through diffuse shortening unlike results of previous drilling on the nankai margin porosity data provide no indication of undercompaction beneath thrust faults furthermore pore water geochemistry data lack clear indicators of fluid flow from depth these differences may be related to coarser material with higher permeability or more complex patterns of faulting that could potentially provide more avenues for fluid escape in turn fluid pressures may affect deformation well drained sand rich material under the frontal thrust could have increased fault strength and helped to maintain a large taper angle near the toe recent resumption of normal frontal imbrication is inferred from seismic reflection data associated decollement propagation into weaker sediments at depth may help explain evidence for recent slope failures within the frontal thrust region this evidence consists of seafloor bathymetry normal faults documented in cores and low porosities in near surface sediments that suggest removal of overlying material overall results provide insight into the complex interactions between incoming materials deformation and fluids in the frontal thrust region

compositional and fluid pressure controls on the state of stress on the nankai subduction thrust a weak plate boundary abstract we show that both fault mineralogy and regional excess fluid pressure contribute to low resolved shear stresses on the nankai subduction plate boundary off southwest japan ring and direct shear tests indicate that saturated clay minerals in the fault possess intrinsically low residual friction coefficients r at stress levels between 1 0 and 40 mpa the direct shear r values for purified smectite are 0 14 0 02 for illite 0 25 0 01 and for chlorite 0 26 0 02 for point load velocities of 0 0001 mm s these clay minerals dominate the nankai subduction decollement zone illite plus quartz is mechanically important in the altered incoming muroto section and the predicted decollement r should lie between 0 2 and 0 32 this low residual strength together with elevated fluid pressure limits shear stresses to below 4 mpa within the frontal 50 km of the subduction system consistent with the low wedge taper in this region a higher wedge taper off the ashizuri peninsula indicates basal shear stresses rise slightly along strike towards this region our analysis indicates lower fluid pressures must predominantly be responsible because only small second order along strike variations in r are predicted to occur as a result of variations in smectite and total clay content these variations should be further reduced at depth under the wedge as smectite is diagenetically altered to illite however our data suggest the low r values of the clay rich decollement still limit shear stresses to between 17 and 29 mpa within the frontal 50 km of the wedge consistent with other estimates of plate boundary weakness indeed we propose that it should be expected that subduction plate boundaries like nankai will be weak because of the intrinsic presence of clay rich faults and moderate fluid overpressures our data do not support the hypothesis that the smectite to illite reaction directly controls the onset of seismogenic behavior deep in the nankai system because there is already a mechanical dominance of illite rather than smectite in the shallow decollement zone and we find all the clay phases tend to velocity strengthen however temperature activated clay diagenesis and dehydration may cause secondary changes in the fault properties and state of stress across the up dip limit of the seismogenic zone

microphysical model predictions of fault restrengthening under room humidity and hydrothermal conditions from logarithmic to power law healing the maximum fault strength and rate of interseismic fault strengthening healing are of great interest to earthquake hazard assessment studies as they directly relate to event magnitude and recurrence time previous laboratory studies have revealed two distinct frictional healing behaviors referred to as dieterich type and non dieterich type healing these are characterized by respectively log linear and power law increase in the strength change with time to date there is no physical explanation for the frictional behavior of fault gouges that unifies these seemingly inconsistent observations using a microphysical friction model previously developed for granular fault gouges we investigate fault strengthening analytically and numerically under boundary conditions corresponding to laboratory slide hold slide tests we find that both types of healing can be explained by considering the difference in grain contact creep rheology at short and long time scales under hydrothermal conditions favorable for pressure solution creep healing exhibits a power law evolution with hold time with an exponent of 1 3 and an apparent cutoff time of hundreds of seconds under room humidity conditions where grain contact deformation exhibits only a weak strain rate dependence the predicted healing also exhibits a power law dependence on hold time but it can be approximated by a log linear relation with of a few seconds we derive analytical expressions for frictional healing parameters i e healing rate cutoff time and maximum healing of which the predictions are consistent with numerical implementation of the model finally we apply the microphysical model to small fault patches on a natural carbonate fault and interpret the restrengthening during seismic cycles

fractured shale gas systems the first commercial united states natural gas production 1821 came from an organic rich devonian shale in the appalachian basin understanding the geological and geochemical nature of organic shale formations and improving their gas producibility have subsequently been the challenge of millions of dollars worth of research since the 1970s shale gas systems essentially are continuous type biogenic predominant thermogenic or combined biogenic thermogenic gas accumulations characterized by widespread gas saturation subtle trapping mechanisms seals of variable lithology and relatively short hydrocarbon migration distances shale gas may be stored as free gas in natural fractures and intergranular porosity as gas sorbed onto kerogen and clay particle surfaces or as gas dissolved in kerogen and bitumen five united states shale formations that presently produce gas commercially exhibit an unexpectedly wide variation in the values of five key parameters thermal maturity expressed as vitrinite reflectance sorbed gas fraction reservoir thickness total organic carbon content and volume of gas in place the degree of natural fracture development in an otherwise low matrix permeability shale reservoir is a controlling factor in gas producibility to date unstimulated commercial production has been achievable in only a small proportion of shale wells those that intercept natural fracture networks in most other cases a successful shale gas well requires hydraulic stimulation together the devonian antrim shale of the michigan basin and devonian ohio shale of the appalachian basin accounted for about 84 of the total 380 bcf of shale gas produced in 1999 however annual gas production is steadily increasing from three other major organic shale formations that subsequently have been explored and developed the devonian new albany shale in the illinois basin the mississippian barnett shale in the fort worth basin and the cretaceous lewis shale in the san juan basin in the basins for which estimates have been made shale gas resources are substantial with in place volumes of 497 783 tcf the estimated technically recoverable resource exclusive of the lewis shale ranges from 31 to 76 tcf in both cases the ohio shale accounts for the largest share

spectrum of slip behaviour in tohoku fault zone samples at plate tectonic slip rates the tohoku earthquake ruptured a fault that also generates slow slip events laboratory experiments on rock samples from the fault show that this spectrum of slip behaviours is a natural consequence of shearing at slow plate convergence rates

detailed coseismic slip distribution of the 1944 tonankai earthquake estimated from tsunami waveforms coseismic slip distribution on the fault plane of the 1944 tonankai earthquake is estimated from inversion of tsunami waveforms three improvements from a previous study satake 1993 are made these are 1 smaller subfaults are used to resolve detailed slip distribution 2 the sub faults fit better to the plate interface geometry and 3 finer and more accurate bathymetry data is used the inversion result shows that a maximum slip of about 3 m occurred on the plate interface off shima peninsula the total seismic moment is estimated to be 2 0 1021 nm mw 8 2 the result confirms that the 1944 tonankai earthquake did not rupture the plate interface beneath the tokai region and supports the existence of the seismic gap in the tokai region suggested by ishibashi 1981 the slip of about 1 5 m on the plate interface beneath atsumi peninsula northeast of the large slip area is necessary to explain the observed tsunami waveforms although no seismic moment release was estimated from strong motion data by kikuchi et al 1999 this may suggest that the rupture beneath atsumi peninsula was slow

reservoir evaluation of the lower silurian longmaxi formation shale gas in the southern sichuan basin of china abstract the lower silurian longmaxi formation the palaeozoic is organic rich black shale in the southern sichuan basin the yangtze plate of china this study analyses the lateral extent and thickness burial depth total organic carbon content and thermal maturity of the longmaxi formation black shale as the key features of the shale gas reservoir the thickness of the black shale ranges from 10 to 170 m the thickest reservoir is located in changning naxi yongchuan region the toc of the shale at the bottom of the formation 50 m thickness is above 2 0 the lateral distribution of toc varies with the lateral distribution of thickness with the maximum toc in the gongxian luzhou yongchuan region the burial depth ranges from 2000 to 4500 m the shale is thermally over mature the evaluation of reservoir characteristics indicates that the longmaxi formation has conditions appropriate for shale gas accumulation and thus resource potential in the southern sichuan basin of china the objective of this preliminary evaluation of the reservoir characteristics is to locate potential areas favourable for exploration the most favourable areas are defined here as those where the thickness of black shale is more than 100 m and the burial depth is less than 3000 m these cover approximately 12 600 km 2 the most favourable areas which cover an area of approximately 5100 km 2 are located in the northeast luzhou region