modeling precipitation over ice sheets an assessment using greenland the interaction between ice sheets and the rest of the climate system at long time scales is not well understood and studies of the ice ages typically employ simplified parameterizations of the climate forcing on an ice sheet it is important therefore to understand how an ice sheet responds to climate forcing and whether the reduced approaches used in modeling studies are capable of providing robust and realistic answers this work focuses on the accumulation distribution and in particular considers what features of the accumulation pattern are necessary to model the steady state response of an ice sheet we examine the response of a model of the greenland ice sheet to a variety of accumulation distributions both observational datasets and simplified parameterizations the predicted shape of the ice sheet is found to be quite insensitive to changes in the accumulation the model only differs significantly from the observed ice sheet for a spatially uniform accumulation rate and the most important factor for the successful simulation of the ice sheet s shape is that the accumulation decreases with height according to the ability of the atmosphere to hold moisture however the internal ice dynamics strongly reflects the influence of the atmospheric circulation on the accumulation distribution

controls on the strength of coupling among climate erosion and deformation in two sided frictional orogenic wedges at steady state 1 many important insights regarding the coupling among climate erosion and tectonics have come from numerical simulations using coupled tectonic and surface process models however analyses to date have left the strength of the coupling between climate and tectonics uncertain and many questions unanswered we present an approximate analytical solution for two sided orogenic wedges obeying a frictional rheology and in a condition of flux steady state that makes explicit the nature and sensitivity of the coupling between climate and deformation we make the simplifying assumption that the wedge grows in a self similar fashion consistent with airy isostasy such that topographic taper is invariant with orogen width tectonic influx rate and climate we illustrate first how and why the form of the erosion rule matters to orogen evolution and then derive a physically based orogen scale erosion rule we show that steady state orogen width crest elevation and crustal thickness are controlled by the ratio of accretionary flux to erosional efficiency to a power dictated by the erosion process remarkably we show that for most combinations of parameters in the erosion law rock uplift rate is more strongly controlled by erosional efficiency than it is by the accretionary flux further assuming frontal accretion with no underplating the spatial distribution of erosional efficiency dictates the relative rock uplift rates on the pro wedge and retro wedge and the time averaged trajectories of rocks through the orogen the restriction to invariant frictional properties is conservative in these respects systems subject to positive feedback between erosion and rheology will exhibit even stronger coupling among climate erosion and deformation than shown here

mechanics of fold and thrust belts and accretionary wedges cohesive coulomb theory a critically tapered fold and thrust belt or submarine accretionary wedge is one that is on the verge of coulomb failure everywhere including its base where frictional sliding along a decollement is assumed to be occurring cohesion within a wedge can add significantly to the overall strength near the toe the effect of this is to decrease the near toe taper leading to a critical topographic profile that is concave upward if the decollement is planar we obtain an approximate self consistent solution for the state of stress within a thin skinned cohesive critical coulomb wedge and determine the relationship between the wedge taper and its strength and basal friction the theory is then applied to the presently deforming fold and thrust belt of western taiwan fitting of theoretical critical wedge shapes to topographic profiles and measurements of the step up angles of thrust faults from the basal decollement are used to constrain the taiwan wedge strength parameters an attractive assertion fully consistent with all the observations is that the mechanics of fold and thrust belts and accretionary wedges is governed by normal frictional and fracture strengths of rocks measured in the laboratory in particular if byerlee s law b 0 85 is adopted as the coefficient of sliding friction on the base we find a coefficient of internal friction 0 9 1 0 in the wedge and a wedge cohesion so 5 20 mpa other solutions having strengths and ambient stresses up to 4 times lower than this can also however satisfy the data

growth and erosion of fold and thrust belts with an application to the aconcagua fold and thrust belt argentina 1 the development of topography within and erosional removal of material from an orogen exerts a primary control on its structure we develop a model that describes the temporal development of a frontally accreting critically growing coulomb wedge whose topography is largely limited by bedrock fluvial incision we present general results for arbitrary initial critical wedge geometries and investigate the temporal development of a critical wedge with no initial topography increasing rock erodibility and or precipitation decreasing mass flux accreting to the wedge front increasing wedge sole out depth decreasing wedge and basal decollement overpressure and increasing basal decollement friction lead to narrow wedges large power law exponent values cause the wedge geometry to quickly reach a condition in which all material accreted to the front of the wedge is removed by erosion we apply our model to the aconcagua fold and thrust belt in the central andes of argentina where wedge development over time is well constrained we solve for the erosional coefficient k that is required to recreate the field constrained wedge growth history and these values are within the range of independently determined values in analogous rock types using qualitative observations of rock erodibilities within the wedge we speculate that power law exponents of 1 3 m 0 4 and 2 3 n 1 characterize the erosional growth of the aconcagua fold and thrust belt this general model may be used to understand the development of mountain belts where orogenic wedges grow as they deform at their coulomb failure limit

improving groundwater project analysis with geophysical data when interpreting the earth s subsurface for groundwater contamination we face a fundamental problem we typically sample one millionth of the relevant material if only these data are used to determine the material between boreholes and the continuity of high hydraulic conductivity units the answers will be ambiguous for example consider a site where three boreholes intersect two hydrofacies figure 1 if no other information is available all of the interpretations shown and many more not shown are equally likely

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

orographic precipitation and the relief of mountain ranges 1 the influence of climate on the relief of mountain ranges has long been supposed but has received little quantitative investigation the simplified system of a longitudinal river profile in a bedrock channel coupled to a physically based model of orographic precipitation is used to explore how orographically influenced patterns of precipitation affect the relief of fluvially sculpted mountain ranges it is shown that the impact on the orography depends on the nature of the precipitation regime significant changes in relief also occur for imposed temperature changes of realistic magnitude although the ability of the relief to respond to any given change depends sensitively on its response time and the timescale of the forcing

on steady states in mountain belts the dynamic system of tectonics and erosion contains important feedback mechanisms such that orogenic systems tend toward a steady state this concept is often invoked but the nature of the steady state is commonly not specified we identify four types of steady state that characterize the orogenic system and illustrate these cases by using numerical model results and natural examples these types are 1 flux steady state 2 topographic steady state 3 thermal steady state and 4 exhumational steady state they refer to the erosional flux the topography the subsurface temperature field and the spatial pattern of cooling ages respectively models suggest that the topography will reach a steady mean form at the scale of an orogenic belt but perfect topographic steady state is unlikely to be achieved at shorter length scales thermal steady state is a precondition for exhumational steady state and in the case of temperature dependent deformation topographic steady state exhumational steady state is characterized by reset age zones spatially nested according to closure temperature as illustrated in natural systems from new zealand the cascadia accretionary margin and taiwan

a fluid model for the shape of accretionary wedges in this paper the accretionary wedges landward of most oceanic trenches are modelled assuming a newtonian fluid rheology it is shown that lubrication theory is applicable and that the addition of material to the wedges can be neglected in determining the shape of the wedges the predicted shape of the wedges is in good agreement with bathymetric profiles across the kurile ryukyu and aleutian wedges from the shape of the above wedges it is deduced that accretionary sediments have viscosities between 1018 and 1019 poise

geochemical characterization of the depositional environment of paleocene and eocene sediments of the tertiary central basin of svalbard abstract the tertiary central basin ctb of svalbard has been studied extensively in the past because it contains sediments from the paleocene eocene thermal maximum petm previous studies have found that a negative carbon isotopic extrusion cie is presents without showing elevated concentration of organic carbon in this study we extend the record of the ctb with a comprehensive inorganic and organic geochemical data set based on the sediments encountered in the research well bh 10 2008 comprising almost 1100 m of cored sediments major and trace elements analyses in combination with mineralogical carbon isotopic total organic carbon and total sulphur measurements were used to reconstruct the depositional environments of the lower paleocene basilika fm up to the eocene aspelintoppen fm several depositional environments have been suggested which are consistent with the sedimentological interpretation of the respective well section back calculation of the initial organic carbon concentration shows that bacterial sulphate reduction had a major impact on the persevered organic carbon inventory of the lower basilika fm and of the paleocene eocene transition section based on organic carbon loss introduced by bacterial decay as well as measured isotopic composition we estimate an initial isotopic composition of the evaluated sediments which helps to elucidate on the fate of labile organic matter and explains the unusually low organic carbon inventory of the paleocene eocene transition sediments

feedbacks among climate erosion and tectonics in a critical wedge orogen the interactions among climate erosion and tectonics have long been of interest to geologists but have not yet been united into a single theoretical framework in this study representations of orographic precipitation and fluvial erosion are combined with the concept of a critical wedge orogen the idealized framework captures the basic system dynamics it also allows for a formal analysis of the precipitation and tectonic interactions in terms of feedback factors and gains and so permits a quantitative comparison of their relative strengths the constraint of self similar growth in a critical wedge orogen acts as a tectonic governor whereby changes in orogen size are strongly damped we determine that this negative tectonic feedback is stronger than the precipitation feedback which may be negative or positive depending on whether precipitation increases or decreases with orogen size for an extreme positive feedback unconstrained runaway growth is possible ultimately leading to plateau formation when orographic precipitation leads to a significant rain shadow there is a strong partitioning of rock uplift rates that favors the wet windward flank of the orogen this flank also dominates the response of the whole orogen to changes in climate or tectonic forcing except in the case where a large fraction of the eroded material is recycled into the wedge finally it is demonstrated that the response time of the orogen depends on the feedbacks and is proportional to the gain of the system

climate and tectonic controls on glaciated critical taper orogens abstract in this work we present the results of a new analytical model that examines the coupling between glacial erosion and orogen development surface processes are assumed to be glacially dominated and tectonic activity is controlled by critical wedge mechanics in these circumstances we find that orogen width is strongly dependent on both the rate of accretion and on the rate of precipitation the orogen size is linked to tectonic and climate changes via proportionality constants the orogen width scales with the rate of accretion to between the 2 3rd and 2nd powers and with the rate of precipitation to the 1 3rd and 5 4th powers the value of the proportionality constants varies with the taper angle and the relative rates of ice deformation and sliding in all cases the sensitivities are higher than those calculated for fluvially eroding critical wedge orogens analytical solutions are supported by the results of a numerical flow line model the flow line model further predicts that uplift rates will be highest at and just below the equilibrium line altitude if glacial ablation is largely a calving process rock uplift in the wedge will be strongly focused towards the toe the predicted response time of glaciated orogens to changes in climate and tectonic forcing is dependent upon the constants of erosion the rheology and the rate of precipitation these analyses predict that actual orogens have variable e folding response times for example approximately 1 5 myr for the southern alps of new zealand and the olympic mountains of washington state and approximately 5 5 myr for the european alps

rheological dependence of extension in wedge models of convergent orogens abstract although several mechanical models of extension in contractional orogens have been proposed in many cases these models have been based on viscous constitutive laws lack of more complicated rheological models has primarily been due to the technical difficulty of solving for deformation with complex rheological behavior and boundary conditions a finite element model is presented that allows for solution of the wedge mechanical problem with viscous non linear viscous and coulomb plastic rheologies this model is used to investigate extension in a doubly vergent wedge model in which deformation is primarily contractional driven by a reversal of basal traction proposed models showing extension in the crust at shallow crustal levels in spite of overall contraction and contraction at depth are confirmed but only for high argand number low crustal viscosity at low argand number little or no extension is observed non linear stress dependence power law flow also suppresses extension but not entirely some extension is observed at high effective argand number however with a coulomb plastic rheology no extension is observed for a wide range of parameter values no extension is observed for models with viscous deformation at depth and plastic deformation near the surface this suggests that the extension is a function of the rheologic model not the effective strength of the crust coeval extension and contraction in convergent orogens with steady external forcing therefore requires viscous deformation to very shallow depths

spatial correlation structure estimation using geophysical and hydrogeological data the large spatial variability of hydraulic properties in natural geologic systems over a wide range of scales and the difficulty of collecting representative and sufficient hydraulic property measurements using conventional sampling techniques render estimation of spatial correlation parameters difficult further compounding the estimation problem is the observation that the integral scale estimate is a function of the measurement support scale to mitigate these problems we investigate the use of tomographic geophysical data in combination with hydrogeological data in the spatial correlation estimation procedure two synthetic case studies were investigated where the scale of the geophysical measurements were varied relative to the scale of the hydrogeological properties the spatial correlation structure parameter estimation procedure was performed in the spectral domain where analysis of data having different support scales and spatial sampling windows was facilitated comparison of the spatial correlation structure parameters estimated from measured data with those of the synthetic aquifers revealed which type of data tomographic hydrogeological or a combination of both was most effective for recovering spatial correlation statistics under different sampling heterogeneity conditions these synthetic case studies suggest that collection of a few tomographic profiles and interpretation of these profiles together with limited well bore data can yield correlation structure information that is otherwise obtainable only from extensive hydrological sampling

geomorphic limits to climate induced increases in topographic relief recognition of the potential for strong dynamic coupling between atmospheric and tectonic processes has sparked intense cross disciplinary investigation and debate on the question of whether tectonics have driven long term climate change or vice versa it has been proposed that climate change might have driven the uplift of mountain summits through an isostatic response to valley incision because isostasy acts to compensate mean elevations the debate hinges on the question of whether climate change can significantly increase topographic relief or more precisely increase the volume of missing mass between summits and ridges here we show that in tectonically active mountain ranges geomorphic constraints allow only a relatively small increase in topographic relief in response to climate change thus although climate change may cause significant increases in denudation rates potentially establishing an important feedback between surficial and crustal processes neither fluvial nor glacial erosion is likely to induce significant isostatic peak uplift

response time of glaciers as a function of size and mass balance 1 theory a simple interpretation of the traditional definitions of glacier and ice sheet response time e g thickness divided by mass balance rate h b suggests that larger glaciers respond more slowly than small glaciers to a perturbation in climate however with reasonable choices for mass balance behavior a scaling analysis shows that the response time of valley glaciers decreases as a function of increasing size when other variables are held constant in essence this is because larger valley glaciers push further into the ablation zone and ablation increases more rapidly than the thickness so that h b gets smaller ice sheets have different mass balance regimes than valley glaciers and as they grow larger the ablation does not increase faster than the thickness therefore as ice sheets grow in surface area the response time increases for both ice sheets and valley glaciers the response time also depends on a mass balance index which is defined as the slope of the balance curve as a function of horizontal distance along the glacier surface rather than elevation the response time decreases as the balance index increases so for valley glaciers an increase in the balance index and an increase in glacier size have opposite effects on the response time the balance index is typically larger for smaller valley glaciers therefore a small glacier will typically respond faster than a large glacier but this quicker response is because of mass balance considerations and not because of the dynamic characteristics of the glacier arising from its small size as often assumed

erosion and exhumation in accretionary orogens experimental and geological approaches 1 this study examines the role of erosion on exhumation and deformation history of accretionary orogens through experimental and geological methods fault propagation kinematics and exhumation rates in an eroded thrust wedge are analyzed through scaled sandbox simulations by testing different parameters friction rheologic layering and sediment input output ratio a model coulomb wedge is submitted to erosion under flux steady state conditions that is the volume of eroded material remains equal to the volume of newly accreted material maintaining a constant surface slope during shortening results show that the way of exhumation depends on the internal dynamics of thrust wedges and conversely on how this dynamics is modified by erosion in the eroded thrust wedges the diversity of exhumation patterns is controlled by the mode of fault propagation depending on the basal friction high or low the vertical component of exhumation is generally higher for the wedges with high basal friction than for low friction wedges the uplift of material occurs along a cluster of subvertical thrusts in the middle part of the eroded thrust wedge with low basal friction the material is exhumed along a series of inclined 20 50 thrusts in the rear of the high friction wedge the zone of maximum exhumation is generally localized in the middle part of the thrust wedge and migrates toward the backstop with continued shortening the vertical exhumation rate increases with time and the material accreted later is rapidly transferred to the main exhumation zone compared to the material accreted during the early stages the exhumation occurs at rates 2 times faster in the wedges with half the thickness than in full thickness wedges the total quantity of eroded material at the end of convergence constitutes 36 50 of initial model area for the thrust wedges with low basal friction and 20 40 for the high basal friction wedges the extent of basal underplating increases with total shortening the area of basal underplated material constitutes up to 30 and 40 of the eroded thrust wedge area for the models with low and high basal friction respectively presence of decollements in the accreted series allows underplating of thrust units developing an anticlinal stack whose growth and location is favored by erosion our results are compared to different present day active and ancient convergent orogens himalaya taiwan southern alps of new zealand and cascadia orogen allowing the different specific exhumation processes to be characterized

controls on the patterns of topography and erosion rate in a critical orogen 1 we solve for the steady state patterns of erosion rate and topography in a critical wedge to understand the mutual adjustment of tectonics erosion and topography in the absence and presence of spatial variations in precipitation rate we consider steady state systems in which tectonics favors a critical topographic form assumed to be the mean elevation across the range and in which surface erosion by rivers and hillslopes operates on a two dimensional landscape we find that 1 a nonuniform critical topographic form implies a nonuniform pattern of ridge valley relief and hence a nonuniform pattern of erosion rate and 2 when the system is forced by local variations in precipitation rate maintenance of the critical topographic form requires a local response of rock uplift that greatly dampens changes in topography we apply these concepts to the western side of the olympic mountains of washington state where mean elevation ridge valley relief and precipitation rates increase from the coast to the topographic crest of the range we find that the main control on the erosion rate pattern is the pattern in mean elevation and the amount of precipitation in contrast the pattern of precipitation is only a minor control as a whole our work demonstrates an approach for developing the theoretical context that is necessary for interpreting spatial associations between patterns in topography precipitation and erosion in natural orogens

critical form and feedbacks in mountain belt dynamics role of rheology as a tectonic governor 1 the rheology that governs deformation within a convergent orogen also controls its topographic form or in other words the covaration of its height or thickness and its width under conditions pertaining to small orogens and for rheologies ranging from linear viscous to coulomb plastic we show that this topographic form is insensitive to the distributions of fluxes into and out of the wedge and is thus a fundamental property of the system it can therefore be thought of as a critical topographic form directly analogous to the well studied case of a coulomb plastic rheology which predicts a constant critical taper angle independent of its size the tendency of the system to evolve toward a critical topographic form can be regarded as a tectonic governor that strongly damps the response of the orogen to changes in climate or tectonic forcing scaling relationships can be derived for the variation of orogen width height and exhumation rate as a function of accretionary flux and precipitation rate this study explores how sensitive these scaling relationships are to the assumed rheology it is found that the scaling relationships vary by less than a factor of 2 across the range of typical geological rheologies from coulomb friction to power law viscous to linear viscous these scaling relationships provide a first order representation of the behavior of convergent wedges and are relatively insensitive to the underlying assumptions

constraints on extension of continental lithosphere the extension of the continental lithosphere is considered in terms of the deformation of a thin viscous sheet with a temperature dependent rheology which is governed primarily by the strength of the upper mantle the lithosphere is subjected either to a constant strain rate or to a strain rate decreasing with time corresponding to the stretching of a region of fixed volume with a constant velocity boundary condition the progress of the extension is influenced by the decrease in strength of the lithosphere owing to its attenuation and by an increase in strength resulting from cooling the relative importance of these two effects depends on two parameters a dimensionless strain rate and the total strain for a given strain rate the extending lithosphere goes through two phases as the strain increases in the first the effect of thinning the lithosphere is predominant and the average strength of the lithosphere remains close to and perhaps slightly below its original strength in the second the effect of cooling of the lithosphere leads to very rapid increase in strength as a function of strain resulting presumably in an effective limit on the degree of stretching that the lithosphere experiences for the case of strain rate decreasing with time when the initial strain rates are of order 10 15s 1 this maximum extension is between 100 and 150 for a strain rate of 10 l4 s 1 the maximum lies around 400 to 500 at higher strain rates than this very large extension may occur before cooling significantly affects the strength of the lithosphere this model accords well with the constraints available from the continental margins of the north atlantic

analogue model of gravitational collapse and surface extension during continental convergence most mountain belts occur where continents collide so it is not surprising that the dominant form of surface deformation is shortening but there is also abundant geological evidence for extension in the central parts of many mountain belts at the same time as shortening occurs elsewhere previous models for extension require temporal changes in the thermal structure of the lithosphere1 2 the rate of convergence2 3 the strength of the crust3 or the geometry of accretion4 here we present a simple model in which no such changes are required for surface extension during the convergent thickening of a viscous crusta layer convergence is driven by the motion of two plates at the base of the layer in the area where one plate subducts under the other the surface begins to extend soon after the start of convergence eventually stretching by more than 60 per cent extension occurs because gravity drives horizontal flow faster at the free surface than in the centre of a viscous layer that is fixed at its base in real mountain belts mid crustal weaknesses may allow the depth dependent motion required for surface extension

length scales for continental deformation in convergent divergent and strike slip environments analytical and approximate solutions for a thin viscous sheet model the deformation of a thin viscous layer that has a moving boundary is investigated for comparison with zones of deformation in the continental lithosphere exact analytical solutions for the case of a newtonian fluid and approximate solutions for the case of fluid with power law rheology show that when the imposed velocity vector is normal to the boundary congressional or extensional regime the deformation field decays away from the boundary with a characteristic length scale 1 3 to 1 10 the wavelength of the imposed boundary velocity distribution for n between 1 and 10 where n is the stress strain exponent in the rheology in contrast when the imposed velocity vector is parallel to the boundary transcurrent regime the length scale of the deformation field is approximately 4 times smaller in each case these length scales decrease approximately as n the difference in length scales arises even in the absence of any buoyancy forces acting on thickened or thinned crust such forces would modify the ratio of length scales but not sufficiently to affect this result

messinian climate change and erosional destruction of the central european alps at the end of the miocene the european alps ceased outward expansion and tectonic uplift and exhumation shifted into the orogen interior this shift is consistent with a change from orogenic construction to orogenic destruction reflecting an increase in the ratio of erosional flux to accretionary flux the coincidence of this change with an increase in sediment yield from the alps suggests a climate driven increase in erosional flux the timing of deformation and sediment release from the southern alps indicates that the tectonic change occurred synchronous with the last phase of the messinian salinity crisis we attribute the increase in erosional flux to a climatic shift to wetter conditions throughout europe likely augmented by the base level fall that occurred during the mediterranean dessication this climate change is represented in the stratigraphic record by the lago mare deposits of the mediterranean salinity crisis

hydrogeological characterization of the south oyster bacterial transport site using geophysical data a multidisciplinary research team has conducted a field scale bacterial transport study within an uncontaminated sandy pleistocene aquifer near oyster virginia the overall goal of the project was to evaluate the importance of heterogeneities in controlling the field scale transport of bacteria that are injected into the ground for remediation purposes geochemical hydrological geological and geophysical data were collected to characterize the site prior to conducting chemical and bacterial injection experiments in this paper we focus on results of a hydrogeological characterization effort using geophysical data collected across a range of spatial scales the geophysical data employed include surface ground penetrating radar radar cross hole tomography seismic cross hole tomography cone penetrometer and borehole electromagnetic flowmeter these data were used to interpret the subregional and local stratigraphy to provide high resolution hydraulic conductivity estimates and to provide information about the log conductivity spatial correlation function the information from geophysical data was used to guide and assist the field operations and to constrain the numerical bacterial transport model although more field work of this nature is necessary to validate the usefulness and cost effectiveness of including geophysical data in the characterization effort qualitative and quantitative comparisons between tomographically obtained flow and transport parameter estimates with hydraulic well bore and bromide breakthrough measurements suggest that geophysical data can provide valuable high resolution information this information traditionally only partially obtainable by performing extensive and intrusive well bore sampling may help to reduce the ambiguity associated with hydrogeological heterogeneity that is often encountered when interpreting field scale bacterial transport data

steady state erosion of critical coulomb wedges with applications to taiwan and the himalaya 1 orogenic structure appears to be partially controlled by the addition to and removal of material from the mountain belt by tectonic accretion and geomorphic erosion respectively we developed a coupled erosion deformation model for orogenic wedges that are in erosional steady state and deform at their coulomb failure limit erosional steady state is reached when all material introduced into the wedge is removed by erosion that is limited by the rate at which rivers erode through bedrock we found that the ultimate form of a wedge is controlled by the wedge mechanical properties sole out depth of the basal decollement erosional exponents basin geometry and the ratio of the added material flux to the erosional constant as this latter ratio is increased wedge width and surface slopes increase we applied these models to the taiwan and himalayan orogenic wedges and found that despite a higher flux of material entering the former the inferred ratio was larger for the latter calculated values for the erodibility of each wedge showed at least an order of magnitude lower value for the himalaya relative to taiwan these values are consistent with the lower precipitation regime in the himalaya relative to taiwan and the exposure of crystalline rocks within the himalayan orogenic wedge independently determined rock erodibility estimates are consistent with the accretionary wedge sediments and metasediments and the crystalline and high grade metamorphic rocks exposed within taiwan and the himalaya respectively therefore differences in rock type and climate apparently lead to key differences in the erosion and hence orogenic structure of these two mountain belts

asymmetric growth of the pyrenees revealed through measurement and modeling of orogenic fluxes orogenic asymmetry is ubiquitous in terms of topography structure and metamorphism but the relative influences of climate and tectonics in determining such asymmetry remain contentious this study demonstrates that for the case of the european pyrenees orogenic asymmetry is not defined by inherited structural asymmetry nor by orographic asymmetry of precipitation but is the simple consequence of the inherent tectonic growth of a doubly vergent orogen the pyrenees record approximately 60 myr of upper crustal accretion erosion and sedimentation that developed from late cretaceous to early miocene times approximately 165 km of plate convergence was accommodated by four phases of orogenic growth that can be summarized in terms of their generic orogenic development as follows 1 early inversion of passive margin extensional faults initiating the pro wedge 2 growth by a combination of frontal accretion and underplating to the pro wedge symmetrical and steady record of erosion and foreland basin sedimentation 3 cessation of pro wedge frontal accretion with propagation of the retro wedge deformation front 4 accelerated growth of pro wedge antiformal stack induced by underplating accelerated erosion over antiformal stack and increase in sediment discharge to pro foreland basin stabilization of outer pro wedge by sediment blanketing discrete element model experiments are used to explore the natural tendency of the system towards asymmetry model output demonstrates that during active frontal accretion to the pro wedge of an orogen the maxima in erosional denudation is offset towards the pro wedge relative to the maximum depths of exhumation this asymmetry is fundamentally driven by the progressive reduction in advection velocities of material within the orogen from the pro wedge into the retro wedge hence for an efficient erosional regime a dynamic coupling between erosion and structural uplift of the pro wedge will reduce activity of the retro wedge and hence drive asymmetry in terms of low temperature thermochronometry erosion and sedimentation the models also demonstrate the tendency for pro wedge growth through frontal accretion to slow however the cessation of frontal accretion as recorded by the south pyrenean thrust belt phases 3 and 4 requires additional forcing it is proposed that the thick wedge top sedimentation recorded at this time increased taper angles and hence stabilized the wedge accretionary influx into the orogen at this time was accomplished through underplating these processes governing asymmetry do not require a spatially variable climate or a laterally heterogeneous underthrust plate

ground penetrating radar for environmental applications abstract ground penetrating radar gpr is a near surface geophysical technique that can provide high resolution images of the dielectric properties of the top few tens of meters of the earth in applications in contaminant hydrology radar data can be used to detect the presence of liquid organic contaminants many of which have dielectric properties distinctly different from those of the other solid and fluid components in the subsurface the resolution approximately meter scale of the radar imaging method is such that it can also be used in the development of hydrogeologic models of the subsurface required to predict the fate and transport of contaminants gpr images are interpreted to obtain models of the large scale architecture of the subsurface and to assist in estimating hydrogeologic properties such as water content porosity and permeability its noninvasive capabilities make gpr an attractive alternative to the traditional methods used for subsurface characterization

implications of sediment flux dependent river incision models for landscape evolution 1 developing a quantitative understanding of the factors that control the rate of river incision into bedrock is critical to studies of landscape evolution and the linkages between climate erosion and tectonics current models of long term river network incision differ significantly in their treatment of the role of sediment flux we analyze the implications of various sediment flux dependent incision models for large scale topography in an attempt 1 to identify quantifiable and diagnostic differences between models that could be detected from topographic data or from the transient responses of perturbed systems and 2 to explain the apparent ubiquity of mixed bedrock alluvial channels in active orogens although certain forms of the various models can be discarded as inconsistent with morphological data we find that the relative intrinsic concavity indices of detachment and transport limited systems defined herein largely dictate whether the various models can be tied to distinctive steady state morphologies preliminary data suggest that no such diagnostic differences may exist and other methods must be developed to test models accordingly we develop and explore differences in the scaling behavior of topographic relief and the extent of detachment versus transport limited channels as a function of rock uplift rate that may allow discrimination among various models further we explore potentially diagnostic differences in the rates and patterns of transient channel response to changes in rock uplift rate in addition to general differences between detachment and transport limited systems our analysis identifies an interesting hysteresis in landscape evolution hybrid channels at the threshold between detachment and transport limited conditions are expected to act as detachment limited systems in response to an increase in rock uplift rate or base level fall and as transport limited systems in response to a decrease in rock uplift rate especially during postorogenic topographic decline the analyses presented set the stage for field studies designed to test quantitatively the various river incision models that have been proposed

time lapse inversion of crosswell radar data the combination of differential radar tomography with conventional tracer and or hydraulic tests facilitates high resolution characterization of subsurface heterogeneity and enables the identification of preferential flow paths in dynamic imaging each tomogram is typically inverted independently under the assumption that data sets are collected quickly relative to changes in the imaged property e g attenuation or velocity however such snapshot tomograms may contain large errors if the imaged property changes significantly during data collection acquisition of less data over a shorter time interval could ameliorate the problem but the resulting decrease in ray density and angular coverage could degrade model resolution to address these problems we propose a new sequential approach for time lapse tomographic inversion the method uses space time parameterization and regularization to combine data collected at multiple times and to account for temporal variation the inverse algorithm minimizes the sum of weighted squared residuals and a measure of solution complexity based on an a priori space time covariance function and a spatiotemporally variable mean we demonstrate our approach using a synthetic 2 d time lapse x z t data set based loosely on a field experiment in which difference attenuation radar tomography was used to monitor the migration of a saline tracer in fractured rock we quantitatively show the benefits of space time inversion by comparing results for snapshot and time lapse inversion schemes inversion over both space and time results in superior estimation error model resolution and data reproduction compared to conventional snapshot inversion finally we suggest strategies to improve time lapse cross hole inversions using ray based inversion constraints and a modified survey design in which different sets of rays are collected in alternating time steps