topography and tectonics of the central new madrid seismic zone results of numerical experiments using a three dimensional boundary element program we present results of a series of numerical experiments designed to test hypothetical mechanisms that drive deformation in the new madrid seismic zone experiments are constrained by subtle topography and the distribution of seismicity in the region we use a new boundary element algorithm that permits calculation of the three dimensional deformation field surface displacement fields are calculated for the new madrid zone under both far field plate tectonic scale and locally derived driving strains results demonstrate that surface displacement fields cannot distinguish between either a far field simple or pure shear strain field or one that involves a deep shear zone beneath the upper crustal faults thus neither geomorphic nor geodetic studies alone are expected to reveal the ultimate driving mechanism behind the present day deformation we have also tested hypotheses about strain accommodation within the new madrid contractional step over by including linking faults two southwest dipping and one vertical recently inferred from microearthquake data only those models with step over faults are able to predict the observed topography surface displacement fields for long term relaxed deformation predict the distribution of uplift and subsidence in the contractional step over remarkably well generation of these displacement fields appear to require slip on both the two northeast trending vertical faults and the two dipping faults in the step over region with very minor displacements occurring during the interseismic period when the northeast trending vertical faults are locked these models suggest that the gently dipping central step over fault is a reverse fault and that the steeper fault extending to the southeast of the step over acts as a normal fault over the long term

rise and fall of the southern santa cruz mountains california from fission tracks geomorphology and geodesy the santa cruz mountains are closely associated with a left bend along the right lateral san andreas fault the loma prieta area on the northeast side of the san andreas consists of fault bounded blocks that rise along active deeply rooted reverse and oblique slip faults six samples from a transect across this area yield concordant apatite fission track ages averaging 4 6 0 5 ma these ages date the time of cooling below 110 c and suggest that about 3 km of unroofing has occurred over the last 4 6 m y allowing for current elevations of about 1 km this suggests an average uplift rate of the order of 0 8 mm yr over the last 4 6 m y to further define the extent and distribution of this young uplift we used morphometric analyses of the youthful topography of the area steep drainage slopes and high local relief indicate that the area northeast of the san andreas forms a well defined zone of high uplift in contrast the region on the southwest side of the san andreas is characterized by broad upwarping and folding more subdued topography old fission track ages and mean quaternary uplift rates of 0 1 0 4 mm yr geodetic data 1906 san francisco earthquake and subsequent strain transients 20 years of interseismic deformation 1989 loma prieta earthquake and 2 years of post loma prieta earthquake strain show that the southern santa cruz mountains repeatedly rise and subside through a complex sequence of bay area deformation events an additional deformation element that involves reverse slip averaging 2 3 mm yr along the foothills thrust system must occur to explain the longer term uplift pattern in the loma prieta area since the late pliocene the asymmetry in deformation on opposite sides of the san andreas probably reflects the contrasting rock types on opposite sides of the fault the influence of preexisting structures and the interaction with neighboring faults of the san andreas system

paleoclimatic signature in terrestrial flood deposits large scale process simulation was used to reconstruct the geologic evolution during the past 600 000 years of an alluvial fan in northern california in order to reproduce the sedimentary record the simulation accounted for the dynamics of river flooding sedimentation subsidence land movement that resulted from faulting and sea level changes paleoclimatic trends induced fluctuations in stream flows and dominated the development of the sedimentary deposits the process simulation approach serves as a quantitative means to explore the genesis of sedimentary architecture and its link to past climatic conditions and fault motion

fitting degradation of shoreline scarps by a nonlinear diffusion model the diffusion model of degradation of topographic features is a promising means by which vertical offsets on holocene faults might be dated in order to calibrate the method we have examined present day profiles of wave cut shoreline scarps of late pleistocene lakes bonneville and lahontan it may be assumed that these scarps were initially at least as steep as the angle of repose offsets range from 1 to 12 m and present slope angles range from 9 to 29 a parameter called apparent diffusion age defined as half the mean square horizontal extent of the slope function of each profile is plotted as a function of scarp offset the points show a clear trend of apparent age increasing nearly linearly with offset if linear diffusion held and scarps were initially vertical apparent diffusion age would be the same for all the shoreline profiles the increasing trend can only partly be explained by nonvertical initial scarp slope and therefore the rate of transport of material downslope must increase significantly faster than a linear law in the range of slopes spanned by the data the transport law must become linear at small slope to reduce scatter between profiles with varying ambient slopes the transport law adopted for the purpose of dating is k0 1 5s2 where s is local slope the transport coefficient k0 is correlated inversely with fan slope suggesting that there is a dependence on the particle size distribution a table is included that allows easy application of the model to scarps with simple initial shape

late cenozoic uplift of mountain ranges and global climate change chicken or egg the high altitudes of most mountain ranges have commonly been ascribed to late cenozoic uplift without reference to when crustal thickening and other tectonic processes occurred deep incision and recent denudation of these mountain ranges abundant late cenozoic coarse sediment near them and palaeobotanical evidence for warmer climates where high mountain climates today are relatively cold have traditionally been interpreted as evidence for recent uplift an alternative cause of these phenomena is late cenozoic global climate change towards lower temperatures increased alpine glaciation a stormier climate and perturbations to humidity vegetative cover and precipitation

influence of late cenozoic mountain building on ocean geochemical cycles in a steady state ocean input fluxes of dissolved salts to the sea must be balanced in mass and isotopic value by output fluxes for the elements strontium calcium and carbon rivers provide the primary input whereas marine biogenic sedimentation dominates removal dissolved fluxes in rivers are related to rates of continental weathering which in turn are strongly dependent on rates of uplift the largest dissolved fluxes today arise in the himalayan and andean mountain ranges and the tibetan plateau during the past 5 m y uplift rates in these areas have increased significantly this suggests that weathering rates and river fluxes may have increased also the oceanic records of carbonate sedimentation level of the calcite compensation depth and delta sup 13 c and delta sup 87 sr in biogenic sediments are consistent with a global increase in river fluxes since the late miocene the cooling of global climate over the past few million years may be linked to a decrease in atmospheric co sub 2 driven by enhanced continental weathering in these tectonically active regions

fault zone connectivity slip rates on faults in the san francisco bay area california the slip rate of a fault segment is related to the length of the fault zone of which it is part in turn the slip rate of a fault zone is related to its connectivity with adjoining or contiguous fault zones the observed variation in slip rate on fault segments in the san francisco bay area in california is consistent with connectivity between the hayward calaveras and san andreas fault zones slip rates on the southern hayward fault taper northward from a maximum of more than 10 millimeters per year and are sensitive to the active length of the maacama fault

fluvial landsculpting and the fractal dimension of topography abstract quantitative models of landform development can help us to understand the evolution of mountains and regional topography and the effects of tectonic motions and climate on landscape including its fractal geometry this paper presents a general and powerful three dimensional model of fluvial erosion and deposition at hill to mountain range scale the model works by accumulating the effects of randomly seeded storms or floods precipitons that cause diffusional smoothing then move downslope on digital topography grids that erode portions of elevation differences that transport a slope limited amount of eroded material and that deposit alluvium when their sediment carrying capacity is exceeded the iteration of these simple and almost linear rules produce very complicated simulated landscapes demonstrating that complex landscapes do not require complex laws each process implemented in the model is affected differently by changes in horizontal scale erosion a scale free process roughens topography at all wavelengths this roughening is balanced by diffusive processes scaling as 1 l 2 at short wavelengths and deposition scaling scaling as 1 l 2 at long wavelenghts such a mixture of scale free and scale dependent processes can produce multifractal behavior in the models the fractal dimension of the model topography is much more sensitive to climatic variables than to tectonic uplift landscape evolution may be fractal but it does not seem to be chaotic analysis of topography of areas in southern arizona using variograms shows approximately fractal behavior with mean fractal dimension around 2 2 2 3 departures from an exact fractal relatioship imply that the topography is in detail multifractal the fractal dimension at short wavelengths is less than that at long wavelengths this variation could either be caused by the relative strengths of diffusive and erosional processes shaping the topography or a result of changes in climatic or tectonic conditions still preserved in the landscape

modification of wave cut and faulting controlled landforms from a casual observation that the form of degraded fault scarps resembles the error function this investigation proceeds through an elementary diffusion equation representation of landform evolution to the application of the resulting equations to the modern topography of scarplike landforms the morphologic observations can be analyzed either in the form of one or more cross strike elevation profiles or in the form of the slope offset plot a point plot of maximum scarp slope versus scarp offset working with either or both of these data representations for nine geologic structures which range in age from 3 to 400 ka b p and in offset from 1 to 50 m we apply analytical solutions for the vertical initial value scarp the vertical continuous offset scarp and the finite slope initial value scarp the model calculations are intrinsically ambiguous yielding as the final answer only the product t in the case of the initial value problem or the product a 1 in the case of the repeated faulting problem here t is the age of a single scarp forming event 2a is the vertical slip rate and is the mass diffusivity a single profile across three sea cliffs along the santa cruz california coast is analyzed as three separate initial value problems a reasonably constrained age for the sea cliff standing above the highway 1 platform returns 11 gkg 1 gkg 1 m2 ka with this we can date the two older sea cliffs in fact we do the converse age estimates for these two older sea cliffs based on a uniform rate of uplift both yield the same as for the lower sea cliff we treat a single profile of the raymond fault in pasadena san marino in terms of the repeated faulting problem for it the uplift rate of r crook and others yields 16 gkg the very substantial preexisting offset across the raymond fault must have been buried leveled some 230 ka b p when the modern topography began to form our analysis of the lake bonneville shoreline scarps reveals a dependence of t on 2a suggestive of nonlinear modification processes this appearance is treated with the finite slope initial value scarp model to determine 1 1 gkg for the lake bonneville shoreline scarps the suggestion of m n machette that approximately 100 000 year old meter high scarps are unobservable in weakly consolidated alluvial terranes of the basin and range and rio grande rift valley provinces can be formulated as 1 gkg the coincidence between this inequality and the lake bonneville shoreline is striking and it suggests that the value of 1 gkg may be generally applicable as a good first approximation to the modification of alluvial terranes within the semiarid regions of the western united states the lake bonneville shoreline is the basis for dating four sets of fault scarps in west central utah the drum mountains fault scarps can be modeled in several different circumstances but the most likely interpretation is that these fault scarps formed as the result of a single episode of normal faulting 3 6 to 5 7 ka b p the younger age is associated with quite low initial slope angles 25 the other three sets of fault scarps show no evidence for finite initial value slopes fault scarps along the eastern base of the fish springs range are very young 3 ka b p we estimate the age of fault scarps along the western flank of the oquirrh mountains to be 32 ka b p which meets the weak geologic constraint that they be older than the lake bonneville shoreline fault scarps along the northeastern margin of the sheeprock mountains are even older 53 ka b p an intriguing consequence of our single event analysis of these scarps is that an 11 5 m offset occurred in a single earthquake

channel changes in badlands stripping of the vegetation and soil from a 13 hectare site in virginia underlain by coastal plain sediments created a rapidly evolving badland topography two types of channels developed 1 sand bed alluvial channels were graded to transport the bed material load supplied from slope erosion with available runoff but they also generally eroded their beds slowly and 2 steeper bedrock floored channels incised rapidly in bedrock channels the erosion rate was proportional to the 4 9ths power of drainage area and the 2 3rds power of gradient these exponents are consistent with a model in which the erosion rate is proportional to the bed shear during high flows due to rapid mass wasting and reduced runoff the alluvial channels became as much as 50 steeper during the winter than the summer with an attendant yearly cycle of winter aggradation and summer entrenchment the gradients their seasonal variability and their downstream hydraulic geometry were consistent with the predictions of total load transport formulas for sand beds and high loads the hydraulic geometry of alluvial channels in the virginia badlands were similar to that on the morrison formation in the western united states

quantification of soil production and downslope creep rates from cosmogenic 10be accumulations on a hillslope profile average soil transport rates over a period of 3500 yr on a convex soil mantled hillslope have been quantified using a mass balance model that incorporates the soil concentration of the cosmogenic isotope 10 be the 10 be model results support the assumption used in most geomorphic models that the soil creep rate is proportional to surface gradient the predicted diffusion coefficient is 360 55 cm 3 yr 1 cm 1 contour length and the average rate of soil production is 0 026 0 007 cm yr within the uncertainty of this technique the data do not reject g k gilbert9s hypothesis that some hillslopes may exist in a condition of dynamic equilibrium with a uniform soil production rate however the model does not require an assumption of dynamic equilibrium and may be an approach that uniquely allows the quantification of a local soil production rate law

complex fault interactions in a restraining bend on the san andreas fault southern santa cruz mountains california the unusual oblique thrust mechanism of the october 18 1989 loma prieta earthquake focused attention on the complex tectonic setting of this segment of the san andreas fault near the mainshock epicenter the san andreas fault curves to the left defining a restraining bend the large thrust component of the mainshock focal mechanism is consistent with the horizontal compression expected across restraining bends however repeated loma prieta type earthquakes cannot exclusively produce the observed topography of the southern santa cruz mountains the highest point of which experienced subsidence during the 1989 earthquake in this paper the authors integrate seismic geomorphic and tectonic data to investigate the possibility that motions on faults adjacent to the san andreas fault play an important role in producing the observed topography the three dimensional geometry of active faults in this region is imaged using the loma prieta preshock and aftershock sequences the most conspicuous features of the fault geometries at depth are 1 the presence of two distinct zones of seismicity corresponding to the san andreas and the sargent berrocal fault zones 2 the concave upward shape of the loma prieta rupture surface 3 the reduction in dip of the deepest portions of the rupture plane asmore the mainshock hypocenter is approached 4 the apparent transfer of shallow slip in some areas from faults in the san andreas fault zone to those in the sargent berrocal fault zone and 5 the presence of a deep northeasterly dipping plane associated with the sargent berrocal fault zone the authors find that a model of fault interactions which involves displacement on faults in both the san andreas and the sargent berrocal fault zones is consistent with loma prieta coseismic displacements preshock and aftershock seismicity and observed topography less

evolution of the northern santa cruz mountains by advection of crust past a san andreas fault bend the late quaternary marine terraces near santa cruz california reflect uplift associated with the nearby restraining bend on the san andreas fault excellent correspondence of the coseismic vertical displacement field caused by the 17 october 1989 magnitude 7 1 loma prieta earthquake and the present elevations of these terraces allows calculation of maximum long term uplift rates 1 to 2 kilometers west of the san andreas fault of 0 8 millimeters per year over several million years this uplift in concert with the right lateral translation of the resulting topography and with continual attack by geomorphic processes can account for the general topography of the northern santa cruz mountains

forcing of late cenozoic northern hemisphere climate by plateau uplift in southern asia and the american west geologic evidence indicates that net vertical uplift occurred on a large kilometer scale and at accelerating rates during the middle and late cenozoic in plateaus of southern asia and the american west based on this evidence general circulation model sensitivity tests were run to isolate the unique effects of plateau uplift on climate the experiments simulated significant climatic changes in many places some far from the uplifted regions the basic direction of most of these simulated responses to progressive uplift is borne out by changes found in the geologic record winter cooling of north america northern europe northern asia and the arctic ocean summer drying of the north american west coast the eurasian interior and the mediterranean winter drying of the north american northern plains and the interior of asia and changes over the north atlantic ocean conducive to increased formation of deep water the modeled changes result from increased orographic diversion of westerly winds from cyclonic and anticyclonic surface flow induced by summer heating and winter cooling of the uplifted plateaus and from the intensification of vertical circulation cells in the atmosphere caused by exchanges of mass between the summer heated and winter cooled plateaus and the mid latitude oceans disagreements between the geologic record and the model simulations in alaska and the southern rockies and plains may be related mainly to the lack of narrow mountain barriers in the model orography taken together the observed regional trends comprise much of the pattern of late cenozoic climatic deterioration in the northern hemisphere that culminated in the plio pleistocene ice ages the success of the uplift sensitivity experiment in simulating the correct pattern and sign of most of the observed regional climatic trends points to uplift as an important forcing function of late cenozoic climatic change in the northern hemisphere at time scales longer than orbital variations however the modest amplitude of the uplift induced cooling simulated at high latitudes indicates a probable need for additional climatic forcing

strike slip faults on mars observations and implications for global tectonics and geodynamics 1 the tectonic evolution of mars has been driven primarily by the interaction of tharsis induced loading stresses with a uniform contractional stress field leading to global assemblages of graben and wrinkle ridges until recently strike slip faults have appeared to be largely absent from the tectonic record we here present evidence for a new set of noachian to early hesperian strike slip faults southwest of tharsis the best preserved fault in the newly mapped set has a length of 200 km and measured lateral offset of 5 9 km we use the locations of both previously mapped and newly identified strike slip faults to investigate the tectonic and geodynamic evolution of mars we model the tectonic evolution of the planet in response to the evolving balance between the loading and contractional stresses finding widespread regions in which strike slip faulting was favored the observed tectonic history is consistent with a scenario in which loading induced stresses peaked early in mars history followed by the growth of contractional stresses leading to a shift in the preferred locus of strike slip faulting from regions southwest of tharsis in the noachian to early hesperian to northwest of tharsis in the early amazonian the contractional strain history inferred from the ages and locations of the strike slip faults is consistent with the strains calculated from models of the thermal evolution of mars in which the gradual secular cooling of the planet over time was augmented by plume induced contraction during tharsis formation in the noachian

correlation between the geologic record and computed plate motions for central california we have attempted to correlate geologic events in part of central california with computed motions of oceanic plates relative to cratonic north america three composite terranes are considered 1 the sierra nevada 2 the coast ranges northeast of the san andreas fault and 3 the salinian block in the sierra nevada jurassic plutonism ending about 147 m y b p and cretaceous plutonism 120 80 m y b p correlate with deduced farallon north america fa na convergence as does the nevadan orogeny 158 153 m y b p however the gap in plutonism 146 121 m y b p outlasted by far an apparent minimum in convergence 145 135 m y b p magmatism ceased in the sierra and uplift slowed during the laramide orogeny 75 45 m y b p when tectonic activity had moved east to the site of the rockies presumably because of a low angle of subduction but volcanisin resumed in the sierra about 33 m y b p indicating that a steeper angle was reestablished apparently east west compression slackened when the pacific north america pa na transform regime developed basalt was erupted in the southern sierra 12 3 m y b p and during this same time the range started its most rapid rise as a normal fault bounded block the coast ranges north east of the san andreas fault are underlain by the franciscan complex the older parts of the complex were probably assembled in a subduction zone far south of their present location but we believe that the same fa na convergence that was cited for sierran plutonism was responsible despite the disparity in latitude as the farallon plate was probably in contact with most of the western margin of north and south america a perceived strong pulse in convergence 100 85 m y b p evidently produced the coast range thrust which on geologic grounds alone most likely originated between 96 and 88 m y b p some of the northwesterly transport of franciscan rocks relative to north america is ascribed to oblique kula north america ku na convergence 85 43 m y b p during which period the coast range thrust was revived about 60 m y b p the farallon plate again influenced the california margin after the kuia plate had moved northward and fa na convergence caused the accretion of the coastal belt and king range terranes of the neo franciscan the advent of pacific north america pa na interaction was marked by local volcanism which advanced from the southeast toward the northwest in consonance with the passage of the mendocino triple junction dickinson and snyder 1978 in the salinian block plutons 107 82 m y old can be ascribed to the same fa na convergence as the contemporary cretaceous plutons of the sierra nevada although the salinian intrusives originated much farther south paleomagnetic results champion et al 1980 from the salinian related pigeon point formation 75 71 m y in age suggest about 2500 km of poleward essentially coastwise movement this postulated large movement during part or all of the last 75 m y poses the greatest dilemma in our study the kuia is the only known plate whose computed obliquely convergent motion is appropriate for the late cretaceous paleogene part of the migration and even this motion is inadequate for the estimated transport within the most likely time window furthermore the farallon plate seems ideally suited for the laramide orogeny 70 45 m y b p both plates could not have been in contact with the same part of the margin at the same time oblique convergence of the kula plate could possibly have caused large coastwise transport and also the laramide orogeny the coast ranges northeast of the san andreas fault and the salinian block both experienced strike slip faulting and en echelon folding correlative with pa na interaction indicating that they were close together by neogene time faulting and folding have been particularly marked since the inception of spreading in the gulf of california and attachment of westernmost california to the pacific plate about 4 m y b p the rise of individual coastal ranges commencing 3 1 m y b p may possibly have been caused by slight convergence which appears in computed pa na relative motions from 5 m y b p to the present

the growth of geological structures by repeated earthquakes 2 field examples of continental dip slip faults a strong test of our understanding of the earthquake cycle is the ability to reproduce extant fault bounded geological structures such as basins and ranges which are built by repeated cycles of deformation along strike slip faults the coseismic and interseismic deformation can be nearly equal in magnitude and opposite in sign resulting in little permanent deformation except for the fault offset for dip slip faults portions of the crust are lifted and dropped and so buoyancy forces are exerted the seismic and interseismic deformations do not balance and structures grow and become subject to erosion and deposition we consider three examples for which the structure and fault geometry are well known the white wolf reverse fault in california site of the 1952 kern county m 7 3 earthquake the lost river normal fault in idaho site of the 1983 borah peak m 7 0 earthquake and the cricket mountain normal fault in utah site of quaternary slip events basin stratigraphy and seismic reflection records are used to profile the structure and coseismic deformation measured by leveling surveys is used to estimate the fault geometry to reproduce these structures we add the deformation associated with the earthquake cycle the coseismic slip and postseismic relaxation to the flexure caused by the observed sediment load treating the crust as a thin elastic plate overlying a fluid substrate the cumulative deformation is principally dependent on the elastic plate thickness modestly sensitive to the sediment substrate density difference and insensitive to the fluid viscosity for the 4 to 8 ma structures we deduce a longterm flexural rigidity of 2 15 1019 nm this is equivalent to an elastic plate thickness of 2 4 km for a young s modulus of 2 5 1010 nm 2 this value is found where independent estimates of the elastic thickness from the coherence between surface topography and gravity yield values of about 4 km but where coseismic fault slip extends to a depth of 10 15 km thus much of the seismogenic crust must weaken substantially during the life of active faults causing the fault bounded basins to narrow over time

geodetic estimate of coseismic slip during the 1989 loma prieta california earthquake offsets in the relative positions of geodetic stations resulting from the loma prieta earthquake can be explained with a dislocation model that includes buried oblique slip on a rupture surface extending 37 km along the strike of the san andreas fault dipping 70 degree to the sw and extending from a depth of about 5 to 17 5 km assuming uniform slip on a rectangular surface the mean values for a range of reasonable fault geometries are 1 6 plus minus 0 3 m right lateral strike slip and 1 2 plus minus 0 3 m reverse slip slip on an adjacent extension of the rupture to the southeast recorded in the aftershock sequence is not well constrained by the geodetic data the geodetic data clearly preclude rupture extending near the surface

hillslope and channel evolution in a marine terraced landscape santa cruz california a flight of marine terraces along the central california coastline provides a unique setting for the study of topographic evolution wavecut platforms mantled by 2 6 m of marine terrace cover deposits are separated by 10 50 m tall decaying sea cliffs paleocliff edges become more rounded with age yet the details of the profiles and frequent bedrock exposure on the upper slopes imply weathering limited transport five bedrock stream channels etched through the marine terrace sequence display one to three distinct convexities in their longitudinal profiles detailed hand level surveys of the hillslopes and of the stream channel longitudinal profiles constrain hillslope evolution and channel incision components of a numerical model of landscape evolution we account for regolith production as a function of regolith depth in accord with the field observation that hillslope processes are presently dominated by the activities of burrowing rodents the transport process is taken to be diffusive stream incision is assumed to be controlled by stream power for which we use the surrogate of local drainage area slope product best fits of the numerical model to field data imply hillslope diffusivity is 10 m2 kyr 1 regolith production rate on bare bedrock is 0 3 m kyr 1 and falls off rapidly with regolith cover and the constant controlling the efficiency of stream incision is 5 to 7 10 7m 1 kyr1

erosional dynamics flexural isostasy and long lived escarpments a numerical modeling study erosional escarpments are common features of high elevation rifted continents fission track data suggest that these escarpments form by base level lowering and or marginal uplift during rifting followed by lateral retreat of an erosion front across tens to hundreds of kilometers previous modeling studies have shown that this characteristic pattern of denudation can have a profound impact upon marginal isostatic uplift and the evolution of offshore sedimentary basins yet at present there is only a rudimentary understanding of the geomorphic mechanisms capable of driving such prolonged escarpment retreat in this study we present a nonlinear two dimensional landscape evolution model that is used to assess the necessary and sufficient conditions for long term retreat of a rift generated escarpment the model represents topography as a grid of cells with drainage networks evolving as water flows across the grid in the direction of steepest descent the model accounts for sediment production by weathering fluvial sediment transport bedrock channel erosion and hillslope sediment transport by diffusive mechanisms and by mass failure numerical experiments presented explore the effects of different combinations of erosion processes and of dynamic coupling between denudation and flexural isostatic uplift model results suggest that the necessary and sufficient conditions for long term escarpment retreat are 1 incising bedrock channels in which the erosion rate increases with increasing drainage area so that the channels steepen and propagate headward 2 a low rate of sediment production relative to sediment transport efficiency which promotes relief generating processes over diffusive ones 3 high continental elevation which allows greater freedom for fluvial dissection and 4 any process including flexural isostatic uplift that helps to maintain a drainage divide near an escarpment crest flexural isostatic uplift also facilitates escarpment retreat by elevating topography in the vicinity of an eroding escarpment thereby increasing channel gradients and accelerating erosion which in turn generates additional isostatic uplift of all the above conditions high continental elevation is common to most rift margin escarpments and may ultimately be the most important factor

evolution of the santa cruz mountains california through tectonic growth and geomorphic decay i describe a strategy for modeling a specific large scale topographic feature that recognizes the discrepancy between tectonic and geomorphic scales and allows for important feedbacks between these processes in the geomorphic model the largest channel in each 4 4 km cell is explicitly treated the channel incision being driven by a stream power rule hillslopes internal to each cell respond rapidly to incision of the local channel prompting their treatment in the model as steady state forms i present analytic expressions for the local relief which is controlled by the local channel incision rate and rate constants associated with the dominant hillslope process either diffusive or landsliding initial conditions include any preexisting topography and a lightly etched a channel network the model is illustrated using the santa cruz mountains california i show that an earlier hypothesis that the santa cruz mountains are due to advection of topography past an uplift source related to a restraining bend in the san andreas fault remains viable when more specific geomorphic processes are considered the local tectonics are represented by crustal thickening necessitated by conservation of mass in the bend right lateral translation of crust with respect to the bend and flexural response to the distributed load the model produces reasonable geographic and statistical distributions of topography using rates of tectonic and geomorphic processes that are within the range of those measured locally this argues that this restraining bend is a long lived feature of the san andreas fault the slightly higher crest and shorter length of the southern santa cruz mountains may potentially be explained by slightly faster slip of the crust west of the fault with respect to the bend than of that to the east deflection is minimal beneath these relatively small mountains disallowing the interesting feedback of erosionally driven uplift the predicted pattern of exhumation could be used to guide the search for fission track sites that could aid in the testing of the model

the growth of geological structures by repeated earthquakes 1 conceptual framework in many places earthquakes with similar characteristics have been shown to recur if this is common then relatively small deformations associated with individual earthquake cycles should accumulate over time to create geological structures following this paradigm we show that existing models developed to describe leveling line changes associated with the seismic cycle can be adapted to explain geological features associated with a fault in these models an elastic layer containing the fault overlies a viscous half space with a different density fault motion associated with an earthquake results in immediate deformation followed by a long period of readjustment as stresses relax in the viscous layer and isostatic equilibrium is restored deformation is also caused as a result of the loading and unloading due to sediment deposition and erosion in this paper the parameters that control the growth of dip slip structures are identified we find that the flexural rigidity of the crust or the apparent elastic thickness provides the main control of the width of a structure the loading due to erosion and deposition of sediment determines the ratio of uplift to subsidence between the two sides of the fault the flexure due to sediment load is much more important in this respect than whether the fault is normal or reverse in character we find that in general real structures are associated with apparent elastic thicknesses of 4 km or less and thus with very low flexural rigidities

effect of height and orientation microclimate on geomorphic degradation rates and processes late glacial terrace scarps in central idaho terrace scarps can serve as a nearly ideal natural laboratory for the study of the evolution of slopes this paper examines the effects of scarp size height and orientation microclimate by keeping constant variables such as age lithology and regional climate if a scarp degrades as a closed system and downslope movement is directly proportional to surface gradient the evolution of the scarp is modeled by the diffusion equation for a group of scarps of same age and known starting angle the diffusion equation model predicts the relation between maximum scarp angle and scarp height h late pleistocene terrace scarps now as steep as 33 25 as well as measured angles of repose for sand and gravel require a starting angle as steep as 33 5 for latest pleistocene idaho and utah scarps as h increases is gentler more degraded than modeled by the diffusion equation with a constant rate coefficient the degradation rate coefficient c increases tenfold with scarp height it should not change with scarp height if downslope movement is solely determined by surface gradient to the first power soil wash appears to be responsible for this departure from the diffusion equation model for transport rate by soil wash is a function of scarp size height south facing scarps are less vegetated and more degraded than north facing scarps for scarps 2 m high the degradation rate c on s facing scarps is 2 times that on n facing scarps for 10 m scarps it is 5 times the observed dependence of the rate coefficient c on scarp height can be removed by normalizing c to values for west facing scarps of the same height the residual c values calculated by this method correlate well with differences in incident solar radiation resulting from the different scarp orientations and maximum gradients this correlation demonstrates the importance of orientation on slope processes and their rates through the differences in freeze thaw cycles soil moisture and vegetative cover scarp morphology may be used to estimate age if one accounts for the effects of climate and for scarp height orientation and lithology for example using the dated bonneville shoreline scarps for calibration and comparing only scarps of equal height we estimate the drum mountains fault scarps to be 9 000 yr old this age is about twice that produced by previous diffusion equation calculations that have not accounted for the height as we have here but it is the same as independent geologic estimates of their age

10be chronometry of bedrock to soil conversion rates abstract we report concentrations of cosmogenic 10 be t 1 2 1 5 10 6 yrs in soil excavated from a soil mantled hillslope in black diamond mines regional park contra costa county california the most striking features of the data are 1 the similarity in the downward decreasing trends of 10 be concentrations in two soil profiles collected 75 m apart 2 the coincidence in each soil profile of the soil bedrock interface as defined by visual inspection of soil pits and the level at which 10 be concentrations attain very low values 4 10 6 atoms g and 3 the extremely low 10 be concentrations in the underlying regolith 0 5 10 6 atoms gram the inventory of 10 be in these soils is low equivalent to about 6000 yrs of 10 be accumulation in a soil initially containing no 10 be on the basis of these measurements and with the aid of simple models of soil 10 be motions on the hillslope we conclude that 10 be loss from the surface is dominated by its removal in soil by creep we calculate local rates of bedrock to soil conversion of between 0 15 and 0 27 km 10 6 yrs comparing these with uplift rates determined for coastal regions of california indicates that soil creep alone is capable of removing soil from the local geomorphic system at a rate equivalent to the rate of uplift of much of the coast

erosion thresholds and land surface morphology we propose a graphical technique to analyze the entirety of landforms in a catchment to define quantitatively the spatial variation in the dominance of different erosion processes high resolution digital elevation data of a 1 2 km2 hilly area where the channel network had been mapped in the field were used in the digital terrain model topog to test threshold theories for erosion the land surface was divided into 20 m2 elements whose shapes were then classified as convergent planar or divergent the entire landscape plotted on a graph of area per unit contour length against surface gradient shows each planform plotting as a separate field a simple steady state hydrologic model was used to predict zones of saturation and areas of high pore pressure to mimic the extreme hydrologic events responsible for erosive instability of the land surface the field observation that saturation overland flow is rare outside convergent zones provided a significant constraint on the hydrologic parameter in the model this model was used in threshold theories to predict areas of slope instability and areas subject to erosion by saturation overland flow both of which can contribute to channel initiation the proportion of convergent elements predicted to exceed the threshold varies greatly with relatively small changes in surface resistance demonstrating a high sensitivity to land use such as cattle grazing overall the landscape can be divided using erosion threshold lines into areas prone to channel instability due to runoff and stable areas where diffusive transport predominates

the quaternary marine terraces of santa cruz california evidence for coseismic uplift on two faults the coast parallel elevation pattern of marine terraces along 50 km of coast centered about santa cruz california provides an opportunity to assess the relative roles of candidate uplift mechanisms we assemble 130 existing latitude longitude elevation inner edge locations from three of these marine terraces to constrain elastic dislocation models of coseismic uplift uplift resulting from repeated oblique slip seismic events similar to the loma prieta earthquake m s 6 9 of october 17 1989 appears to explain well the portion of the terrace uplift pattern south of santa cruz depending on whether the lowest terrace is assumed to correspond to isotope stage 5c or 5e the calculated mean repeat time for loma prieta events is either 490 or 820 yr corresponding to either 4 1 or 2 4mm yr slip rates on the loma prieta fault residual uplift not explained by such events increases monotonically to the north toward ano nuevo where the san gregorio fault zone intersects the coastline we assess whether this residual uplift may result from repeated events on the san gregorio fault given only the strike of the san gregorio fault we use an extensive parameter search to determine the geometry of and the style of motion on the implied rupture plane that best explains this residual uplift the fault parameters that minimize the misfit between data and model imply activity of a moderately steep 70 north east dipping san gregorio fault that ruptures in oblique thrust slip from 15 km to the surface the thrust slip component dominates the uplift pattern resulting in uplift of the coastline from santa cruz northward these parameters match available geological geodetic and seismic data on the dip and sense of motion on strictures within the san gregorio fault zone the san gregorio fault has therefore been an active late pleistocene participant in north american pacific plate boundary motion we conclude that the complex pattern of uplift along this 50 km reach of coastline can be explained by coseismic mechanisms tied to san andmas bend related fault structures onshore and to the san gregorio fault zone offshore as both of these mechanisms result locally in considerable shore normal uplift gradients caution is warranted in applying common correlation techniques based solely on uplift rates in the assignment of terrace ages in this and other such settings in addition while it is tempting to use terrace culmination offsets as a means of estimating san andreas fault dextral slip rates we provide a formalism to determine under what conditions such calculations are meaningful we conclude that the terrace offsets along the santa cruz coastline result not from dextral slip on the san andreas but from differential indentation of the loma prieta uplift pattern by the various sea level highstand coastlines

application of a diffusion erosion model to alluvial channels which degrade due to base level lowering the applicability of a diffusion erosion model that predicts transient longitudinal profiles of degrading alluvial channels has already been demonstrated in laboratory experiments and is here applied to some field examples by fitting simulated longitudinal profiles to measured ones for several streams with discharge 0 04 q5 680 m3s 1 diffusion erosion coefficients k were determined which are well correlated to discharge and drainage area headward knickpoint migration which results from base level lowering was also treated through the diffusion erosion model theoretically the distance of a knickpoint from the mouth of a tributary should be proportional to the square root of the time that elapsed since the base level was lowered field observations are in accordance with these theoretical considerations use of the model is discussed for prediction of erosion and for assessment of the age of alluvial fault scarps crossed by streams

scarp degraded by linear diffusion inverse solution for age under the assumption that landforms unaffected by drainage channels are degraded according to the linear diffusion equation a procedure is developed to invert a scarp profile to find its diffusion age diffusion age having dimension length 2 is the product of diffusivity times chronological age the second moment of scarp slope grows linearly with age this fact together with an assumption about initial scarp shape allows the inverse determination of diffusion age age found assuming a vertical initial scarp is termed apparent age any nonvertical initial scarp profile has a nonzero initial apparent age true diffusion age differs from apparent age by a fraction of scarp offset squared the inverse procedure applied to synthetic data yields the following rules of thumb evidence of initial scarp shape has been lost when apparent age reaches twice its initial value if a scarp is formed by two events the inversion gives their offset weighted mean age with an error that is a fraction of offset squared a scarp that appears to have been formed by one event may have been formed by two with an interval between them as large as apparent age after scarps of two fault traces have diffused to appear as one the error in inferred age may be as large as half the apparent age variation of apparent age along strike would indicate multiple fault traces the simplicity of scarp profile measurement and this inversion makes profile analysis attractive if linearity of the flow law time for a free face to be reduced to the angle of repose and variation of diffusivity with climate and material could be established profile analysis would become a reliable dating technique

coseismic deformation and dislocation models of the 1989 loma prieta earthquake derived from global positioning system measurements global positioning system gps surveys were used to determine coseismic displacements within the epicentral region of the 1989 loma prieta earthquake near santa cruz the california highway department caltrans conducted the preearthquake survey in february and march 1989 using single frequency gps receivers in march 1990 we reoccupied eight of the caltrans stations using dual frequency receivers relative displacements were computed by differencing the station coordinates determined from these two surveys displacements were determined relative to sites remote from the earthquake by incorporating gps measurements from loma prieta station lp1 to the very long baseline interferometry vlbi site at fort ord examination of repeated measurements in the single frequency preearthquake survey indicates precision is 10 mm in the horizontal although it may be as low as 40 mm in the weakly connected part of the network the precision in the vertical is generally 30 40 mm in contrast the largest observed coseismic displacement is 410 mm horizontal and 341 mm vertical we use a quasi newton method to estimate the parameters of the uniform rectangular dislocation source that best fits the observed displacements the best fitting dislocation dips 68 southwest extending from 3 to 16 km depth strikes n66 w and has 1 8 m of dip slip and 1 2 m of strike slip the best fitting dislocation overlaps the aftershocks but is rotated 15 counterclockwise from the trend of the afterhock zone we find that the data can be fit nearly as well with a dislocation that is constrained to lie in the plane of the aftershocks in this case the fault extends from 5 to 15 km depth the constrained model explains 96 of the data both models have a moment of m0 2 9 1019 n m we also use the quasi newton method to estimate rectangular dislocation models from the leveling data of marshall et al 1991 and the trilateration and gps data of it lisowski et al 1990 in contrast to the previous studies we find that both of these data sets are independently best fit with dislocations embedded within the aftershock zone

results from a new model of river basin evolution this paper briefly describes a model of the erosional development of catchments and their channel networks the model differentiates between the dominant transport processes in hillslope and channels the development of channels and hillslopes occurs in an integrated manner as a function of physically observable mechanisms the growth of a river basin is qualitatively described the model concepts are used to study the basin during periods of growth transient periods as well as during dynamic equilibrium this leads to hypotheses about the relationship between slopes relief tectonic uplift erodability runoff and catchment area it is shown that the model leads to very reasonable and desirable behaviour of hillslopes during retreat and degradation

escarpment evolution on high elevation rifted margins insights derived from a surface processes model that combines diffusion advection and reaction experiments with a surface processes model of large scale 1 1000 km long term 1 100 m y erosional denudation are used to establish the controls on the evolution of a model escarpment that is related to the rifting of a continent the model describes changes in topographic form as a result of simultaneous short and long range mass transport representing hillslope diffusive processes and fluvial transport advection respectively fluvial entrainment is modeled as a first order kinetic reaction which reflects the credibility of the substrate and therefore the fluvial system is not necessarily carrying at capacity one dimensional and planform models demonstrate that the principal controls on the evolution of an initially steep model escarpment are 1 antecedent topography drainage 2 the timescale or equivalently a length scale in the fluvial entrainment reaction 3 the flexural response of the lithosphere to denudation and 4 the relative efficiencies of the short and long range transport processes when rainfall and substrate lithology are uniform a significant amount of discharge draining over the escarpment top causes it to degrade only when the top of the model escarpment coincides with a drainage divide can escarpment retreat occur for these conditions an additional requirement for retreat of a model escarpment without decline is a long reaction time scale for fluvial entrainment this corresponds to a substrate that is hard to detach by fluvial erosion and therefore to fluvial erosion that is not transport limited continuous backtilting of an escarpment due to flexural isostatic uplift in response to denudational unloading helps maintain the scarp top as a divide it is essential if the escarpment gradient is to be preserved during retreat in a uniform lithology low flexural rigidities promote steep and slowly retreating escarpments for given rainfall and substrate conditions the morphology of a retreating model escarpment is determined by the ratio of the short range diffusive and long range advective transport efficiencies a low ratio which is interpreted to correspond to a relatively arid climate and weathering limited conditions promotes steep sharp topped escarpments with straight main slopes and escarpment retreat occurs over a wide range of height scales a high ratio interpreted to correspond to a more humid temperate climate produces a convex upper slope and concave lower slope morphology and only major escarpments are predicted to preserve a high scarp gradient lithological contrasts in the model produce more complex morphologies and predict the formation of scarps crowned by an erosionally resistant caprock however resistant caprocks are not an essential requirement for model scarps to retreat we conclude that the inferred controls and model behavior are both consistent with the present day morphology of rifted continental margins and with modern conceptual models of landscape evolution

the morphology of strike slip faults examples from the san andreas fault california fault zones are seldom straight nor are fractures within them continuous permanent strain fields are imprinted near a fault of complex geometry as a result of continued slip these strain fields may be manifest as vertical deformation near the fault strain fields that increase area result in subsidence and strain fields that reduce area result in uplift although these effects are modified by erosion and deposition the general features of fault zone morphology are a consequence of the geometry and slip distribution along strike graphical views of strain fields associated with simple vertical strike slip faults are derived using boundary element methods if the fault is straight and parallel to the applied shear strain the hills and valleys resulting from finite fault rupture are separated by the fault trace this is rarely observed in nature instead strike slip faults cut through hills and valleys apparently of their own making we show this to be a consequence of segments of the fault being nonparallel to each other and to the overall applied shear strain coseismic strains are large but transient since they are substantially annulled by rupture on contiguous fault segments we test our theory with examples from parts of the san andreas fault system in california the morphologies of the parkfield region and the southern san andreas fault in the coachella valley are shown to be consistent with plausible distributions of slip on mapped surface faults the location of the salton sea appears to be a consequence of a reduction of slip on the san andreas fault toward the brawley seismic zone

seismic moment frequency relation for shallow earthquakes regional comparison we determine the parameter values for the seismic moment frequency relation using flinn engdahl s regionalization of global seismicity and the harvard centroid moment tensor data the earthquake size distribution is approximated by the gamma law a version of the gutenberg richter distribution with an exponential taper at the maximum moment there is no statistically significant variation of the fj value the analog of the b value for all seismic regions except for the midocean ridge systems for the latter regions fj 0 93 whereas for all other zones fj 0 63 the maximum moment mxg can be statistically evaluated only for subduction zones treated as a whole mxg 10 21 2 x 10 22 newton m which corresponds to the worldwide mxg value for other regions as well as for single subduction zones mxg is determined by comparing the number of events in each zone with the seismic moment rate calculated on the basis of the nuvel 1 model of plate motion for subduction zones we obtain the estimate of mxg which agrees with the statistical value giving evidence that most tectonic deformation is released by earthquakes we test the hypothesis that no statistically significant variation in mxg occurs in subduction and continent collision zones the hypothesis cannot be rejected with available data for the midocean ridges the mxg estimate cannot be unambiguously determined it is quite possible that the estimate is significantly biased because only a small part of tectonic deformation at the ridges is released by earthquakes these results have importance in evaluating seismic risk and may also lead to developing a physical theory for earthquake generation