crustal structure and surface zonation of the canadian appalachians implications of deep seismic reflection data in 1986 1181 km of marine seismic reflection data was collected to 18 20 s of two way traveltime in the gulf of st lawrence area the seismic profiles sample all major surface tectono stratigraphic zones of the canadian appalachians they complement the 1984 deep reflection survey northeast of newfoundland together the seismic profiles reveal the regional three dimensional geometry of the orogen three lower crustal blocks are distinguished on the seismic data they are referred to as the grenville central and avalon blocks from west to east the grenville block is wedge shaped in section and its subsurface edge follows the form of the appalachian structural front the grenville block abuts the central block at mid crustal to mantle depths the avalon block meets the central block at a steep junction that penetrates the entire crust consistent differences in the seismic character of the moho help identify boundaries of the deep crustal blocks the moho signature varies from uniform over extended di

stress field constraints on intraplate seismicity in eastern north america focal mechanisms of 32 north american midplate earthquakes mo 3 8 6 5 were evaluated to determine if slip is compatible with a broad scale regional stress field derived from plate driving forces and if so under what conditions stress regime pore pressure and frictional coefficient using independent information on in situ stress orientations from well bore breakout and hydraulic fracturing data and assuming that the regional principal stresses are in approximately horizontal and vertical planes 10o the constraint that the slip vector represents the direction of maximum resolved shear stress on the fault plane was used to calculate relative stress magnitudes defined by the parameterb s2 s3 s s3 from the fault stress geometry as long as the focal mechanism has a component of oblique slip i e the b axis does not coincide with the intermediate principal stress direction this calculation identifies which of the two nodal planes is a geometrically possible slip plane gephart 1985 slip in a majority of the earthquakes 25 of 32 was found to be geometrically compatible with reactivation of favorably oriented preexisting fault planes in response to the broad scale uniform regional stress field slip in five events was clearly inconsistent with the regional stress field and appears to require a localized stress anomaly to explain the seismicity significantly all five of these events occurred prior to 1970 when many regional networks were installed and their focal mechanisms are inconsistent with more recent solutions of nearby smaller events the frictional likelihood of the geometrically possible slip on the selected fault planes was evaluated in the context of conventional frictional faulting theory the ratio of shear to normal stress on the fault planes at hypocentral depth was calculated relative to an assumed regional stress field regional stress magnitudes were determined from 1 s s3 ratios based on the frictional strength of optimally oriented faults the basis for the linear brittle portion of lithospheric strength profiles 2 the computed relative stress magnitude b values and 3 a vertical principal stress assumed equal to the lithostat two end member possibilities were examined to explain the observed slip in these less than optimally oriented fault planes first the frictional coefficient was held constant on all faults hydrostatic pore pressure was assumed regionally and the fault zone pore pressure was determined since pore pressure is a measurable quantity with real limits in the crust p0 s3 this end member case was used to determine which of the geometrically possible slip planes were frictionally likely slip planes alternately pore pressure was fixed at hydrostatic everywhere and the required relative lowered frictional coefficient of the fault zone was computed slip in 23 of the 25 geometrically compatible earthquakes was determined to also be frictionally likely in response to an approximately horizontal and vertical regional stress field derived from plate driving forces whose magnitudes are constrained by the frictional strength of optimally oriented faults assuming hydrostatic pore pressure regionally the conditions for slip on these 23 relatively well oriented earthquake faults were determined relative to this regional crustal strength model and require only moderate increases in pore pressure between about 0 4 0 8 of lithostatic hydrostatic is about 0 37 of lithostatic or alternately moderate lowering 50 of the frictional coefficient on the faults which slipped superlithostatic pore pressures are not required focal mechanisms for the two other earthquakes with slip vectors geometrically consistent with the regional stress field however did require pore pressures far exceeding the least principal stress or extremely low coefficients of friction these events may reflect either local stress rotations undetected with current sampling or poorly constrained focal mechanisms the analysis also confirmed a roughly north to south contrast in stress regime between the central eastern united states and southeastern canada previously inferred from a contrast in focal mecha nisms between the two areas most central eastern united states earthquakes occur in response to a strike slip stress regime whereas the southeastern canadian events require a thrust faulting stress regime this contrast in stress regime with a constant maximum horizontal stress orientation determined by far field plate driving forces requires a systematic lateral variation in relative stress magnitudes superposition of stresses due to simple flexural models of glacial rebound stresses are of the correct sense to explain the observed lateral variation but maximum computed rebound related stress magnitude changes are quite small about 10 mpa and do not appear large enough to account for the stress regime change if commonly assumed stress magnitudes determined from frictional strength apply to the crust at seismogenic depths

an improved estimate of lithospheric thickness based on a reinterpretation of tilt data from pleistocene lake algonquin a reinterpretation of relative uplift data associated with the deformed strandline of pleistocene lake algonquin is presented the new analysis is based on a dynamic model that allows for the finite rate of retreat of the laurentide ice sheet and the viscous relaxation of the earth s mantle attention is also paid to the effects of the shape of the load profile on the interpretation in the context of this model the thickness of the lithosphere is constrained to be not more than 85 km this value is based on the assumption of a model ice sheet with a parabolic cross section if the load has a sharper edge a thickness of 110 km is consistent with the data these revised values slightly exceed walcott s original estimate they are on the other hand distinctly lower than the value of 200 km which has recently been inferred by peltier for the thickness of the north american lithosphere the relations between the different estimates are briefly discussed it is suggested that a reconciliation of the estima

crust and upper mantle velocity structure in the appalachian orogenic belt implications for tectonic evolution the crust and upper mantle structure of the northeastern united states neus was studied by combining teleseismic and regional body wave observations with surface wave dispersion measurements the velocity models suggest that structures as much as 200 km deep and even deeper can be correlated with surficial geologic and tectonic features and that the grenville and appalachian orogenic belts show marked differences in crustal structure this has the important implication that major orogenic belts have effects that reach well into the lithosphere and are stable for extended periods of time perhaps as long as 1 b y regional travel times recorded across the neus seismic network indicate that the northern appalachians are characterized by a well defined two layer crust with a relatively high velocity lower layer the crust of the grenville province in new york state is slightly thinner than that of the new england appalachians and it is vertically homogeneous with nearly constant p and s velocities lateral variations in structure were studied using time term analysis and teleseismic p wave residuals the suture between the grenville and appalachian provinces in the northeastern united states probably occurs along a north northeast trending belt extending from northwestern vermont to southwestern connecticut high crustal velocities and or crustal thinning a linear gravity high a serpentinite belt precambrian uplifts and the taconic thrusts are found along much of this belt comparison of seismic velocities with resistivity measurements suggests that the lower crust of the grenville province may be composed of rocks with hydrous mineral phases resulting in lowered velocities and a higher poisson s ratio alternatively the rocks of the lower crust beneath the grenville province may be similar to those found on the surface while higher velocity mafic mineralogies are prevalent in the lower crust of the appalachians this is consistent with the hypothesis that the grenville crust underwent substantial reactivation thickened and became vertically uniform during the grenville orogeny in contrast the rocks of the appalachian belt probably were associated with a cycle of oceanic opening and closure which suggests an ensimatic origin of the lower crust in this region recent cocorp results indicate that the orogenic belt of the southern appalachians is allochthonous and is underlain by a basal decollement that extends eastward to the coastal plain however contrasts in structural style and tectonic evolution between the northern and southern appalachians suggest that it may not be possible to extrapolate the cocorp results to the northern appalachians north of about 41 latitude and that the suture between the grenville province and the appalachian orogenic belt in new england is located in the vicinity of or just east of the serpentinite belt

control of tectonic inheritance on continental intraplate strain rate and seismicity abstract present day deformation and seismicity of continental lithosphere are characterized by a first order dichotomy between plate boundary zones pbz and stable continental regions scr whereas the former are associated with high strain and seismicity rates the latter tend to remain un deformed except in localized regions of higher strain and seismicity commonly related to fossilized paleo pbz acting as locally weaker domains because of their low amplitudes these intraplate strain and seismicity rates are particularly difficult to measure and characterize in this study we propose a simple model to explain and quantify first order continental strain rate variations focusing on intraplate regions assuming near failure equilibrium on 1d lithosphere profiles we derive steady state strain rates driven by tectonic forces as a function of rheological models that include new strain weakening rheologies in order to simulate tectonic inheritance within this framework inherited strain weakening plays a fundamental role in allowing for and explaining strain and seismicity concentration in intraplate weak zones our model predicts strain rates in intraplate weak zones up to two to three orders of magnitude higher compared to stable cratons due to the effect of tectonic inheritance on rheology weakening these model predictions are in agreement with empirical estimations in intraplate regions and can provide a conceptual framework for characterization of scr strain and seismicity rates

seismotectonics of northeastern united states and adjacent canada data for local earthquakes recorded by a network of stations in northeastern united states and adjacent canada were analyzed to study the seismicity the relationship between earthquakes and known faults the state of stress and crustal and upper mantle velocity structure in addition portable seismographs were deployed in the field to study aftershocks as a result accurate locations for about 364 local earthquakes 2 mb 5 and 22 focal mechanism solutions were determined a comparison of the spatial distribution of these events 1970 1979 with historical earthquakes 1534 1959 reveals that seismic activity in the northeast is relatively stationary in space those areas that have had little or no seismicity historically are relatively aseismic today whereas the historically active areas are also active today the instrumental locations historical seismicity and focal mechanism solutions show an internal consistency that help us distinguish two distinct seismogenic provinces 1 the adirondack western quebec province is a northwesterly trending zone of seismic activity about 200 km wide and at least 500 km long extending from the se adirondacks into western quebec canada thrust faulting on planes striking nnw to nw appears to predominate and the inferred axis of maximum horizontal compression is largely uniform and trends wsw nearly parallel to the calculated absolute plate motion of north america little or no seismicity is found where anorthosite outcrops at the surface correlations between gravity anomalies and earthquake locations suggest that seismic activity in this zone is localized to regions of steep ne or sw gradient in bouguer anomalies this zone does not appear to extend southeastward to boston as proposed by some workers 2 the appalachian province is a northeasterly trending zone of seismic activity extending from northern virginia to new brunswick canada highangle reverse or thrust faulting on n to ne trending planes appears to predominate the western margin of this province however appears to be relatively aseismic we attribute this relative lack of activity to one or more of the following the presence of igneous activity postdating rifting of africa from north america the occurrence of intense metamorphism during the acadian orogeny which may have annealed preexisting faults and the predominance of ductile as opposed to brittle deformation in the geologic past the inferred axis of maximum horizontal compression along the eastern margin of the appalachians is rather uniform and trends w wnw almost perpendicular to the magnetic lineations offshore we suggest that this w wnw compression reflects the gravitational force arising from horizontal density variations in the oceanic lithosphere as it cools and moves away from spreading centers

seismic activity and structural features in the charlevoix region quebec the charlevoix region is historically the most active earthquake zone in eastern canada understanding the links between its seismicity and the faults of the region is important for the assessment

the thickness of the continental lithosphere two independent sets of geophysical data are shown to be particularly sensitive to the thickness of the continental lithosphere both are connected with the response of the planet to the last deglaciation event of the current ice age the first set of data consists of relative sea level histories in the age range 0 10 kyr b p from sites along the eastern seaboard of north america such data appear to require a continental lithospheric thickness in excess of 200 km the second set of data consists of the ils record of the motion of the rotation pole over the time range 0 75 years b p which reveals a secular drift at the rate of 0 95 0 15 deg 106 years toward hudson bay when the viscosity profile of the mantle is fixed to that required to explain free air gravity and relative sea level data from within the ice margin and the observed nontidal acceleration of planetary rotation then these observations also suggest a continental lithospheric thickness of the same order

three dimensional interpretation of gravity data from sedimentary basins using an exponential density depth function in mapping the topography of the basement of deep sedimentary basins by gravity modelling the accuracy can be improved by incorporating an exponential increase in density with depth for calculating the gravity effect of a three dimensional 3d structure with such an exponential density depth relation a frequency domain forward algorithm based on series expansion is presented the numerical evaluation of which can be performed efficiently by fast fourier transform the algorithm can be applied in a recursive procedure to give the inverse solution in terms of basement relief the inversion procedure is satisfactorily tested on a 2d synthetic example and a 3d field example of gravity data from the western margin of the pannonian basin in eastern austria where up to 2 2 km of tertiary sediments overlie an igneous or metamorphic basement the results are confirmed by basement intersections in several wells

mechanisms of crustal deformation thirty sixthwilliam smith lecture the energy source for crustal deformation is isotopic heating and secular cooling of the mantle in a true solid heat would be lost to the surface by conduction however solid state creep processes allow the earth s solid mantle to exhibit a fluid behaviour thus thermal convection can convert heating into directed motion variations in temperature lead to variations in densitythrough thermal expansion and contraction although the resultant body forces are vertical horizontal variations in temperature lead to horizontal body forces at an ocean ridge the variations in the thickness of the lithosphere lead to the elevation of the ridge the result is a substantial horizontal ridge push force this force along with the trench pull force on the descending lithosphere is the primary force driving plate tectonics lithospheric thinning must occur beneath continental rifts and plateau uplifts i suggest that lithospheric thinning is caused by the diapiric penetration of hot asthenospheric rock to the crust mantle boundary this penetration may also cause lithospheric stoping and plateau uplifts lithospheric thinning leads to large tensional stresses above the zone of thinning these forces are generally much larger than the compressional forces generated by crustal thinning the forces generated by lithospheric thinning are quite large and may be responsible for the propagation of continental rifts in island arcs and andean type orogenic belts both lithospheric thinning and crustal thickening result in large tensional forces probably the most complex zone of crustal deformation is the continental collision variations in types of continental collisions are discussed

maximum horizontal stress orientations in the cooper basin australia implications for plate scale tectonics and local stress sources summary borehole breakouts and drilling induced tensile fractures ditfs were interpreted in 61 wells in the cooper basin indicating an average maximum horizontal stress orientation of 101 n a total of 890 borehole breakouts and 608 ditfs were interpreted in the cooper basin the approximately east west maximum horizontal stress orientation is consistent over much of the basin except in the patchawarra trough where maximum horizontal stress rotates to a northwest southeast orientation this rotation in maximum horizontal stress orientation is consistent with in situ stress data to the northwest of the cooper basin the stress field in the cooper basin appears to mark the apex of a major horseshoe shaped rotation in maximum horizontal stress direction across central eastern australia finite element modelling of the in situ stress field of the indo australian plate iap using a range of plate scale tectonic forces is able to match the regional maximum horizontal stress orientation over most of australia reasonably well including the mean east west maximum horizontal stress orientation in the cooper basin however plate boundary scale modelling does not adequately match the horseshoe shaped stress rotation across central eastern australia the average east west maximum horizontal stress orientation in the cooper basin indicates that stresses from tensional forces acting along the tonga kermadec subduction zone are not transmitted into the interior of the australian plate the majority of the tensional forces associated with the tonga kermadec subduction zone are most likely accommodated along the numerous spreading centres within the lau havre backarc basin a number of more localized stress anomalies have also been identified these cannot be explained by plate scale tectonic forces and are possibly a result of geological structure and or density contrasts locally perturbing the stress field

intraplate stresses and seismicity the role of topography and density heterogeneities abstract the stress induced by the topography and by density heterogeneities in the lithosphere has been computed for two different tectonic regions of the north american plate the lithosphere was assumed to consist of a two dimensional layered elastic slab over an inviscid fluid the calculations show that these stresses are locally comparable in magnitude to the regional tectonic stress 20 50 mpa two areas have been analyzed in greater detail the southeastern united states and the rio grande rift the stress in the southeastern united states is computed along a cross section perpendicular to the strike of the appalachian mountains the computed stress appears largest beneath the blue ridge and valley and ridge provinces where the appalachian seismicity is located the stress in the western united states is computed along a cross section perpendicular to the rio grande rift large stress is shown beneath the rift and the adjacent regions the stress induced by surface and internal loading is extensional in the rio grande rift and compressional in the colorado plateau and in the great plains

seismicity and faulting in the charlevoix zone of the st lawrence valley abstract continuous monitoring of the seismically active charlevoix zone since late 1977 has delineated the boundaries of the earthquake zone under the north shore and has enabled an association to be made with surface mapped faults that strike along the st lawrence valley and dip to the southeast within the active zone aseismic volumes are found under parts of the st lawrence river and another elongated zone of activity is found to the southeast on unmapped faults dipping under the river the previous suggestions that the activity is related to old rift faults which have been later weakened by a meteor impact would seem to be upheld

heat flow and structure of the lithosphere in the eastern canadian shield twenty two new determinations of heat flow and radiogenic heat production in the superior and grenville provinces of the canadian shield are presented the new data and previously published measurements strongly constrain the thermal structure of the eastern canadian shield in the abitibi greenstone belt heat flow gradually increases from 29 mw m 2 near the grenville front to 44 mw m 2 east of the kapuskasing uplift this heat flow variation is interpreted in terms of crustal thickening and increased thickness of a tonalitic layer with average heat production of about 1 1 w m 3 this interpretation based on estimated heat production of major rock types in the region is consistent with crustal models derived from recent seismic reflection and refraction studies it also leads to an estimate of about 12 mw m 2 for the mantle heat flow beneath the area the average heat flow in the grenville province 41 10 mw m 2 is the same as that of the superior province this similarity and the lack of significant variation of heat flow across the grenville front indicate that the crust on both sides of the front has similar heat production and thus composition in the western part of the grenville province heat flow reaches high values in the vicinity of the boundary between the allochtonous polycyclic and monocyclic belts where enriched granitic plutons are found in the crystalline terranes in the central part of the grenville province heat flow and heat production are related to each other the parameters of the linear heat flow heat production relationship qr 30 2 in mw m 2 and d 7 1 1 7 km are close to those of the much younger appalachian province implying that the higher appalachian heat flow is due solely to higher heat production in the upper crust the data provide no evidence for variation of mantle heat flow between the superior grenville and appalachian provinces whose tectonic ages range between 2700 and 400 ma the small value of the mantle heat flow about 12 mw m 2 implies that the depth to the 450 c isotherm which controls the effective elastic thickness of the lithosphere is very sensitive to crustal heat production

anisotropy of the flexural response of the lithosphere in the canadian shield 1 we show that the two dimensional coherence between topography and bouguer gravity is anisotropic throughout most of the canadian shield excepted hudson bay basin we interpret the axis of maximum coherence as the direction where the lithosphere is weakest in the southeastern shield and the appalachians the azimuthally averaged coherence cannot be fitted by the thin elastic plate model in the other regions it yields the effective elastic thickness the average coherence increases in the direction perpendicular to the main tectonic boundaries the axis of high electrical conductivity and the seismic fast axis this consistency suggests that the same stress field affected the entire lithosphere during the last tectonic event

focal mechanisms of earthquakes in southern quebec southeastern ontario and northeastern new york with implications for regional seismotectonics and stress field characteristics focal mechanisms have been derived for 30 earthquakes in the adirondack mountains western quebec seismic zone from a modified version of the computer program focmec using first motions of pn and pg waves and sv p amplitude ratios measured from sg and pg phases reverse faulting on medium to steeply dipping planes is prevalent focal depth and velocity structure can vary within realistic limits without much influence on the focal mechanism solution deviatoric compressive stress axes have a strong horizontal component but the trend varies most align approximately northeast southwest in accordance with both the orientation predicted by models based on plate tectonics and that derived from most other studied earthquakes in the region however for several earthquakes mechanisms have been obtained that deviate from this pattern notably for events near 45 5 n to 46 n 75 w and the north gower main shock aftershock series of october 1983 possible explanations for these anomalies in terms of local structural disturbances or reactivation of preexisting zones of weakness with different structural orientation than those created during the current tectonic regime are presented different hypotheses for the seismotectonic evolution and processes of the area are presented and discussed in addition to orientations of nodal planes and deviatoric stress axes there are other indicators such as the seismicity distribution and patterns of earthquake sequences of inhomogeneity in the proposed adirondack mountains western quebec seismic zone

intraplate seismicity and stress in the southeastern united states abstract the stress induced by topography and by density heterogeneities in the lithosphere has been computed for two seismic regions in the southeastern united states the southern appalachian mountains and the south carolina coastal plain the lithosphere was assumed to be a three dimensional layered elastic slab overlying an inviscid fluid the calculations indicate that the local stress is of the same magnitude as the tectonic stress tens of mpa but contradictory conclusions are suggested by the comparison of stress differences and principal directions with earthquake locations and focal mechanisms in the southern appalachians the seismicity correlates well with the stress maximum and focal mechanisms agree with the calculations where regional stress is combined with locally induced stress in south carolina a combination of regional and local stress can explain the orientation of focal mechanisms but stress does not concentrate in the region of charleston and thus other factors must play a role in the charleston seismicity

ridge forces absolute plate motions and the intraplate stress field torque poles are calculated for a variety of possible forces acting on the plates including ridge push slab pull and collisional resistance these torque poles are then compared to the directions of absolute plate motions there is a strong correlation between ridge torque poles and the azimuth of absolute plate motions for the north american south american pacific cocos and eurasian plates simple slab pull torques correlate well with absolute motion azimuths for the pacific nazca and cocos plates and moderately well with the absolute motion azimuth of the indoaustralian plate collisional resistance torque poles correlate with the absolute motion azimuth of the eurasian plate only the correlations are presented as further evidence that the absolute reference frame for plate motion is determined by the surface plates themselves torque poles for various forces are also compared with several long wavelength features of the global intraplate stress field that also tend to be aligned with absolute motion directions in general ridge torque directions agree well with the orientations of maximum horizontal stresses for stable north america western europe and south america and provide an alternative explanation for the alignment in terms of ridge push forces rather than basal drag collisional resistance forces can also explain the alignment of stresses in western europe for the indoaustralian plate the torque pole for collisional resistance forces is consistent with the general pattern of stresses in at least the western half of the plate but is not a good predictor of the entire data set for the plate other processes in addition to collisional resistance must be important for the indo australian plate ridge push forces may account for a significant portion of the long wavelength features of the intraplate stress field especially away from continental collisions such a conclusion is consistent with negative buoyancy of the slab being an important component of the driving mechanism as previously suggested slab forces may be largely balanced within the subducted slab itself and thus have limited effect on deformation of the surface plates

some seismological and geometric features of intraplate earthquakes to understand the nature and cause of intraplate earthquakes we carried out a systematic search for common features associated with them various geological geophysical and seismic data were analyzed and geometric and seismic parameters were compared for 29 intraplate events m 4 5 occurring worldwide the results suggest that intraplate earthquakes occur by the reactivation of pre existing zones of weakness in a compressional stress regime which is generally oriented parallel to the absolute direction of plate motion two styles of faulting were observed for type a earthquakes observed in the eastern united states eastern china western europe and west africa deformation occurs by strike slip motion on steeply dipping faults type b events are usually associated with thrust or normal faulting and were observed in eastern canada fennoscandia australia and peninsular india focal depths did not explain the differences between type a and type b earthquakes rather the differences are due to the perturbation of a regional stress field by a local stress field the regional stress field for type a events is primarily due to ridge push forces the method of stress perturbation varied from region to region the earthquakes in the type b category were further divided into three subcategories type b1 occurring in eastern canada and fennoscandia type b2 occurring in australia peninsular india and central brazil and type b3 occurring in the gulf coast of the u s the inferred causes of the perturbing stresses are deglaciation intracontinental resistive forces resulting from the collision of plates and the deposition of a large thickness of sediments respectively type a events were found to be associated with intersecting faults and generally supported the intersection model of talwani 1988 and had analogs in plate boundary events

stress induced by topography and crustal density heterogeneities implication for the seismicity of southeastern canada abstract the stress induced by the topography and density heterogeneities in the lithosphere has been calculated for a large region covering most of southeastern canada the calculations assume that the lithosphere consists of a three dimensional layered elastic slab overlaying an inviscid fluid the calculations show that the deviatoric stress difference induced by density heterogeneities and the topographic load is of the same order as the regional tectonic stress tens of mpa stresses appear largest along the saint lawrence valley where most of the seismicity is located e g the charlevoix and the lower saint lawrence areas the charlevoix seismic zone was analyzed in greater detail and theoretical focal mechanisms were computed this analysis shows that the stress induced by surface and internal loading is compressional in this area and the principal stress orientation is consistent with earthquake focal mechanisms and in situ measurements

numerical modeling of flexural stresses in the upper crust of ne brazil implications for stress regimes in intraplate settings abstract we model the contribution of the stresses generated by the topography and density heterogeneities in the stress field of ne brazil and the interplay between these and regional stresses our model consists of a horizontally layered elastic slab the lithosphere that floats above an inviscid fluid the asthenosphere where the surface load is generated by the topography and bathymetry and the internal loads are calculated from the bouguer anomalies the deviatoric flexural stresses are calculated in the fourier domain our results showed local flexural stresses with magnitudes comparable to those of the plate wide stresses and mainly controlled by the long wavelength bouguer anomalies scenarios for the total stress field were calculated by superposing different values for an e w oriented regional compression we predict the stress axes orientations and stress regimes for the equatorial margin and the pernambuco shear zone pesz region however we highlight the importance of spreading stresses not considered here in rotating the maximum horizontal compression shmax predictions the prediction of stress regime and stress axes orientation in pesz using low to zero regional compression suggest that a decoupling of the shallow stresses from the plate wide stresses might occur in this region another possible explanation found was that the spreading stresses might counterbalance the regional stresses we conclude that the superposition of principally flexural and plate wide stresses explains reasonably the observed stress regimes and shmax orientations in a great part of ne brazil furthermore flexural stresses are playing an important role in the reactivation of the pesz possibly controlling the stress field

aftershock sequences of the 1982 miramichi new brunswick earthquakes the mb 5 7 miramichi new brunswick earthquake of 9 january 1982 was located by comparison to the aftershock distribution in central new brunswick canada at 47 00 n 66 60 w with a focal depth of 7 km and an uncertainty of 3 km on each of the spatial coordinates the main shock was followed by a long and complex aftershock sequence with four events m 4 5 and more than 800 events m 1 0 until 30 june 1982 the aftershocks were studied by three temporary field networks set up in the epicentral areas from 10 to 20 january 2 to 7 april and 17 to 23 june 1982 and by the permanent canadian seismograph networks the aftershock activity was concentrated in a volume approximately 6 km ns 6 km ew 8 km deep and arranged in a north trending conjugate v pattern the miramichi earthquakes were predominatly thrust earthquakes occurring on north trending mid to high angle faults in the upper 8 km of the crust in the presence of an east west compressive horizontal regional stress field detailed field surveys have not yet been able to associate the earthquakes with any known preexisting faults the main shock had a moment of 2 2 0 7 1024 dyne cm and a stress drop of 35 to 70 bars the rupture started at 7 km depth on a plane striking 195 dipping 50 down to the west and propagated up to the east with a rake angle of 120 aftershocks deeper than 3 5 km had composite mechanisms very similar to the main shock s but shallower events took place on more steeply dipping planes up to 78 dip the rupture planes of the main shock and its principal aftershocks may have been similarly curved the magnitude frequency relationship b 0 8 of the aftershocks was comparable to the regional one the aftershock activity decayed as t 0 8 and had the highest rate of activity following the principal aftershock on 11 january that aftershock mb 5 4 appears to have ruptured the east dipping plane conjugate to the main shock s plane two other principal aftershocks mb 5 1 on 9 january and mb 5 0 on 31 march were associated with the main shock rupture plane while the fourth large event mb 4 7 on 16 june occurred 30 km west of the other activity and was seemingly unrelated however all the earthquakes were apparently associated with the same granitic terrain the miramichi earthquakes did not produce any significant surface faulting one small 25 mm surface displacement was found but it was identified as a secondary stress induced movement along a preexisting joint

stress channelling and partitioning of seismicity in the charlevoix seismic zone quebec canada summary the charlevoix seismic zone csz in the st lawrence valley of quebec is historically the most active in eastern canada the structurally complex region comprises rift faults formed during the opening of the iapetus ocean superimposed by a 350 ma meteorite impact structure resulting in a circular highly fractured zone although seismicity is localized along two steeply dipping planar rift parallel zones previous work indicates that most of the large scale rift faults bound the low magnitude background seismicity rather than generate earthquakes themselves in order to gain insight into the mechanics of the partitioning of this seismicity a 2 d model of the csz was built using the stress analysis code flac the rift faults are represented by frictional discontinuities the heavily fractured impact structure is represented by an elastic continuum of reduced modulus boundary displacements are used to generate a regional stress field with the major horizontal component in the direction of tectonic loading given a high strength the rift faults have little effect on the stress patterns stress trajectories naturally flow around the crater of reduced elastic modulus leaving the fractured area with lower stresses than the background level however when the rift faults have low strength they are unable to support stress trajectories inclined to them due to the resolved shear stress exceeding their strength this prevents trajectories from flowing out of the rift effectively channelling higher magnitude stresses into the region of the impact structure between the faults low strength bounding faults can thus explain the localization of seismicity into linear bands rather than distributed seismicity throughout the impact structure it also explains how the rift faults act as boundaries to regions of low magnitude seismicity these results indicate that the interplay between faults of varying strength and zones of differing elastic modulus can give rise to complicated stress patterns and can explain many of the seismicity patterns observed in the csz this has implications for other intraplate seismic zones as it shows an example of how regional weak faults can modify stress conditions around local structures and drive seismicity the results are particularly relevant for other regions located within rifted crust such as the new madrid seismic zone which possibly display evidence of stress channelling

characteristics of the western quebec seismic zone a comparison of geological data shows that most of the earthquakes of western quebec lie within or near the eastern and western boundaries of a grenville ca 1000 ma central metasedimentary belt and in the area of the younger fault zone along the ottawa river although most of the earthquakes have occurred within terrain of grenville age or older the younger rift zone is characterized by the larger earthquakes temiskaming m6 2 1935 cornwall massena m5 7 1944 the earthquakes also correlate with gradient zones around two prominent anomalies in the smoothed bouguer gravity field the pattern of aeromagnetic anomalies and drainage channels reflects preeambrian structural trends and indicates the central metasedimentary belt continues beneath paleozoic sediments to lake ontario the strikes of magnetic gravity and drainage features north of lake ontario are parallel to features mapped to the south of the lake geological data and geophysical trends suggest that a cataclastic zone is continuous alon

characterisation of the seismological pattern in a slowly deforming intraplate region central and western france abstract we analyzed the seismicity of central and western france using historical data a compilation of all recorded earthquakes from 1962 to 2002 4574 events relocated and all published focal mechanisms 119 focal solutions the aim is to understand what are the causes of earthquakes and stress accumulation in a slowly deforming intraplate region the distribution of earthquakes and focal mechanisms is first correlated with recognized faults geological structures and tomographic images then in order to better understand the distribution of hypocenters and seek deeper crustal sources for stress accumulation euler solutions are computed from the available bouguer anomaly data the analysis of the obtained pattern for heat flow values provides a better understanding of the concentration of seismicity in some particular zones two different behaviors of this slowly deforming intraplate region are evidenced one is linked to the presence of a hot spot under the massif central the other to reactivation of the hercynian structural heritage our results highlight that several possible sources for earthquake clustering can be invoked in central western france

the 1982 cocrust seismic experiment across the ottawa bonnechere graben and grenville front in ontario and quebec summary during the summer of 1982 the canadian consortium for crustal reconnaissance using seismic techniques cocrust conducted a major long range seismic refraction and wide angle reflection experiment across the grenville province of the canadian shield three seismic lines each approximately 300 km in length were located i along the ottawa bonnechere graben ii perpendicular to the graben and iii perpendicular to the crenville front geological evidence indicates that the graben may have originated from a failed arm of the st lawrence rift system and the grenville front marks the boundary between the grenville province and the much older superior province other geological features of the grenville province that were traversed by the profiles were the central metasedimentary belt and the central gneiss belt analysis of the data involved conventional traveltime procedures coupled with the use of synthetic seismogram analysis using programs that were written to handle laterally heterogeneous structures results from the survey indicate a variation in near surface seismic velocity

thick and deformed antarctic sea ice mapped with autonomous underwater vehicles satellites have documented trends in antarctic sea ice extent and its variability for decades but estimating sea ice thickness in the antarctic from remote sensing data remains challenging in situ observations needed for validation of remote sensing data and sea ice models are limited most have been restricted to a few point measurements on selected ice floes or to visual shipboard estimates here we present three dimensional 3d floe scale maps of sea ice draft for ten floes compiled from two springtime expeditions by an autonomous underwater vehicle to the near coastal regions of the weddell bellingshausen and wilkes land sectors of antarctica mean drafts range from 1 4 to 5 5 m with maxima up to 16 m we also find that on average 76 of the ice volume is deformed ice our surveys indicate that the floes are much thicker and more deformed than reported by most drilling and ship based measurements of antarctic sea ice we suggest that thick ice in the near coastal and interior pack may be under represented in existing in situ assessments of antarctic sea ice and hence on average antarctic sea ice may be thicker than previously thought

spatial variation in the frequency magnitude distribution of earthquakes under the tectonic framework in the middle east abstract spatial variations of seismic energy released and b value over the middle east region are investigated based on a seismicity catalog from 1995 to 2007 the goal is to use these seismic parameters and based on other geodetic and geophysical observations such as gps measurements strain rate model fault distribution focal mechanism crustal model q model and gravity measurements etc to uncover spatial patterns that seem anomalous areas of high energy released cumulative seem to correspond to the areas of relatively high b values areas of high energy released and high b values seem to correspond very well with the location of continental collision where earthquake activities are high the divergent boundary between arabia and african plates and subduction zone at makran seem to correspond to low to moderate energy release northern pamir azerbaijan caucasus the lower part of zagros mountains eastern turkey owen fracture zone strait of bob el mandeb and south of the sulaiman shear zone seem to correspond to high cumulative energy released high strain rate high b values and high fault density while the central and eastern iran southern zagros northern pakistan gulf of aden alborz southwest of the caspian sea western caucasus and kopeh dagh seem to correspond with lower b values the cross section map for hindu kush shows general decreasing of the b values with depth however a region of high b value is observed underneath pamir at depths from 170 to 230 km this anomaly region can be due to dehydration of pamir crustal slab at depth