Abstract The science of the internally generated behavior and spatial organization of depositional systems has come a long way since Beerbower first coined the term “autocycles” to refer to fining-upward sequences generated by river meander migration, cutoff, and eventual return. Ongoing research has broadened the scope and scale range of known autogenic dynamics, even as a unifying theme—sediment storage and release—has emerged. Many internally generated processes do not have a single characteristic length or time scale but rather occupy a broad scale range (hence, “autocyclic” has been gradually replaced by “autogenic”). But even where they are broad, the scale ranges for autogenic processes are bounded by limiting time and length scales. The central role of sediment storage and release provides a means of estimating these limiting length and time scales based on mass balance, geometry, and mean sediment flux. Recent research has also allowed us to expand the upper limits of autogenic behavior to time scales of 10 5 to 10 6 years. Finally, we recognize that autogenic dynamics is not simply superimposed on allogenic signals but interacts strongly with, modifies, and even destroys allogenic input. That the autogenic imprint on the stratigraphic record is stronger and more complex than once thought can be seen as an opportunity to focus on using the record to learn about intrinsic surface behavior under pre-human conditions, rather than simply as an archive of externally imposed signals.

Abstract The general tendency for sediment supply to deep water to be relatively high during eustatic fall and lowstand, and relatively low during rise and highstand, is recognized in the sequence-stratigraphy literature. Much less is known about the cumulative effect of repeated eustatic cycles on net deep-water sediment delivery. Here we investigate the net effect of offshore sediment delivery during a complete eustatic cycle, which we term sediment pumping , and the possibility of cumulative sediment pumping if repeated eustatic cycles increase the net delivery of sediment to deep water averaged over several cycles. We measure sediment pumping in terms of net offshore delivery after one or more complete eustatic and associated cycles relative to delivery in the absence of cycles. Combining data from a quasi-2D laboratory experiment and a 2D geometric model, we find that net sediment pumping over isolated and superimposed base-level cycles of variable period varies from somewhat negative to strongly positive, depending on (1) time period of imposed base-level cycle, (2) sense of rotation of the spatial subsidence pattern, and (3) the phase of sediment supply relative to eustatic variation. A relatively short-period base-level cycle (i.e., period less than the basin equilibrium time) increases net pumping (relative to the constant base-level reference case) whereas a relatively longperiod cycle yields no or even negative net pumping. Short-period base-level cycles superimposed on a long-period cycle produce a strong net offshore sediment pumping. Other factors being equal, base-level cycles with basin subsidence cause substantially greater net pumping in backtilted basins than in foretilted ones. When sediment supply varies over a base-level cycle, pumping is maximized when the sediment-supply maximum occurs during eustatic falling stage or lowstand. External Controls on Deep-Water Depositional Systems SEPM Special Publication No. 92 (CD version), Copyright © 2009 SEPM (Society for Sedimentary Geology), ISBN 978-1-56576-200-8, p. 41–56.

Hillslope stability depends strongly on local conditions, such as lithology and rock strength, degree of saturation, and critical slope angle. Common triggers for slope failure include severe storms, earthquakes, and removal of material from the toe of the hillslope. In this paper, we focus on the latter, in a model in which streams incise the toe and destabilize the hillslope. We investigate possible interactions between migrating knickpoints and hillslope failures in a small-scale, steadily eroding experimental landscape that experiences steady rainfall and base-level fall conditions. We monitored knickpoint propagation and hillslope failure activity with time lapse photography over a time period in which numerous knickpoints migrated through the drainage basin. We then investigated temporal and spatial relationships between hillslope failures and knickpoints and compared these results to Monte Carlo simulations of hillslope failure distributions. When focusing along a single channel, we found that, statistically (significant at the 98% confidence level), a greater number of failures occur downstream from a migrating knickpoint. These results highlight both the organized and random nature of hillslope and knickpoint interactions.

Abstract Determination of the transport ("diffusion") coefficient, the main parameter of most forward models for generating fluvial stratigraphy, requires finding the average slope required to transport the total sediment load delivered to a given point for a given water discharge. Finding this value, in turn, requires averaging the substantial fine-scale local variability in transport capacity that one encounters in most natural rivers. The problem is especially acute for braided rivers, in which the local capacity varies strongly in time and space as channels migrate, flow shifts from one part of the channel network to another, and confluences, which account for a disproportionate share of sediment flux, form and dissipate. Here, we present a model for computing spatially averaged sediment flux in a sandy braided river system. Coupled with sediment mass balance, the sediment-flux model leads to the usual diffusion equation for surface topography. The problem of indeterminacy of channel width is dealt with by using an empirical constant value of 1.8 for the mean nondimensional (Shields) stress. We test the model by applying it to a mine-tailings fan in which all independent parameters (sediment flux, water flux, grain size, deposition pattern) are well known and constant. The statistical parameters needed to determine the transport coefficient are determined from independent measurements of the river network on the fan. Using these inputs, the model predicts the fan topography well. The model suggests that, for a highly active braided system such as this one, the effect of the fluctuations in sediment flux can increase total sediment flux by a factor of two to four relative to what would be predicted from mean values alone. The data also suggest, however, that some of the key statistical parameters vary significantly downstream along the fan. This variation may result from downstream variation in grain-size distribution, sediment flux, or both.

Abstract These notes are designed to introduce the concepts and techniques of quantitative modeling of basin subsidence histories. The notes also describe some of the methods and results of modeling the development of sedimentary sequences generated by the interaction of subsidence, sediment supply, and sea-level changes. Analyzing and modeling basin subsidence can be a powerful tool for understanding how, when, and why basins form and, thus, compliments other basin analysis techniques. Subsidence analysis is useful for understanding regional tectonic history, history of sea-level changes, and the thermal history of basins that, in turn, are essential for diagenetic and hydrocarbon studies, and it also provides a basis for interpreting lithospheric structure and mechanics. In addition, quantitative models that take into account the competing effects of varying sediment supply rates, subsidence rates, and rates of sea-level change provide insights into the relative importance of tectonics versus eustatics in the generation of the basin-filling sequences. Comparison of the observed stratigraphy to that predicted by synthetic stratigraphic models of basin sequences can be used to help determine whether or not, and to what degree, basin deposition is controlled by tectonics or eustatics. This course concentrates on the theory and application of subsidence and stratigraphic modeling by working through specific examples from real or artificial basin sequences. By the end of the course, you should be able to apply some of the basin-modeling techniques to your own studies and have a strong enough introduction to other techniques that you can understand most of the available literature.

Abstract The application of Archimedes' principle to the earth suggests that continents are buoyed up by a force equal to the weight of the displaced mantle (Turcotte and Schubert, 1982). Adjacent blocks of different thickness and/or density structure will have different relative relief (Fig. 2.1). Typical lithospheric structure beneath the continents and the oceans are shown in Figure 2.2, these values will be used in most of our discussions. Below some depth, there is no density contrast between the two adjacent columns, and asthenosphere of equal density underlies both columns (Fig. 2.2). The weight of the columns above this depth of compensation must be equal. In this model of isostasy, we can calculate the relative relief between two adjacent continental columns of differing density structure (Fig. 2.3):

Abstract The goal of geohistory analysis is to produce a graphical representation of the vertical movement of a stratigraphic horizon in a sedimentary basin as an indicator of subsidence and uplift history in the basin since the horizon was deposited (Van Hinte, 1978; Fig. 3.1). Several types of stratigraphic data are needed to do a geohistory or subsidence analysis. These data include a stratigraphic column showing the present-day thickness of the stratigraphic units, types of lithologies, ages of horizons, and estimated paleowater depths. Other types of data that are useful, although not necessary, are porosity information for the units and thermal information, if your goal is to determine thermal history of the basin. In addition, there are several assumptions and uncertainties that are built into this analysis. Most of these problems can be overcome if thick stratigraphic sections of relatively shallow-water deposits are used and only long-term, large-scale changes are studied.

Abstract In an earlier section (Chapter. 2, part C), we discussed an isostatic balance for a section of rifted lithosphere. It was noted that the long-term subsidence of the rift was associated with cooling and thickening of the lithosphere. Cooling causes the lithospheric rock to become more dense. To a first approximation, the density of the mantle lithosphere (ρL) varies with temperature ( T ) according to: If the asthenosphere is nothing more than hot (T 1 = 1333 °C) mantle lithosphere, then we can calculate a density of 3184. kg/m 3 for the asthenosphere using Equation. 4.1. This suggests a density difference of only 4.4% between the coldest mantle lithosphere and the asthenosphere, but integrated over the entire thickness of the lithosphere temperature differences make a substantial contribution to the overall isostatic balance. In the following sections we will discuss a thermal model developed by McKenzie (1978) to describe the subsidence history of rift basins.

Abstract First, we need to go back and look at simple isostatic balances (discussed in Chapter. 2, part A) in a slightly different way. Consider the two columns shown in Figure 5.1. The column on the left is a reference column, and the column on the right shows the same crustal section thickened now by a factor of MP, where 0 < p < 1. The base of the thickened crust is deflected into the mantle by an amount w due to the weight of the added crust. A local isostatic balance requires that: This result says that the weight of the mountain belt (including the portion that lies below datum) is balanced by a buoyancy force from the mantle. This is a local isostatic balance in that the deflection of the crust at any location depends only on the local amount of crustal thickening at that location. One important shortcoming of the local isostatic balance is that it neglects the lateral strength of the lithosphere. A more realistic (and highly successful) model assumes that the lithosphere responds to loads like an elastic plate overlying an inviscid fluid. The elastic plate corresponds to some poorly defined, colder portion of the thermal lithosphere, whereas the inviscid fluid corresponds to a hotter portion of the lithosphere and the asthenosphere. You should not forget that the elastic plate model is just an extension of the local isostatic compensation model.

Abstract Having discussed the primary mechanisms of subsidence we can briefly focus on the plate tectonic settings of major sedimentary basins and examine their typical subsidence histories and mechanisms (Fig. 6.1). Much has been written about the driving mechanisms of basin formation in most tectonic settings. An early overview was provided by Dickinson (1976). The purpose of this chapter is not to attempt to summarize the state of knowledge of basin evolution. Instead, we simply define each basin type, following the basin classification scheme of Dickinson (1976), and focus on a few key points regarding the mechanisms of basin evolution. Due to space limitations we do not cite all of the relevant literature but provide just a few key references. The subsidence curves (Fig. 6.1) come primarily from the published literature, augmented by analyses done by ourselves or by students in our sedimentary basins course at the University of Wyoming. All the curves have been backstripped following the local-isostatic method described above (Steckler and Watts, 1978). However, inconsistencies arise from the use of different time scales, compaction corrections and paleowater depth estimates made by the various authors. Nonetheless, the overall consistency of the subsidence curves in each of the various tectonic settings suggests that use of different scales by different workers do not generate errors large enough to mask the overall trends.

Abstract Karl Popper once wrote, “Out of theories we create a world: not the real world, but our own nets in which we try to catch the real world.” One of the more exciting applications of quantitative analysis of basin development is its use as a tool in determining the relative importance of tectonic, eustatic and climatic effects on the development of basinal stratigraphy. The goal in modeling the development of basinal stratigraphy is not to try to explain every detail observed in the real world, but instead serves as a framework from which we can view the real world. As such, simplified basin-filling models can guide us to look for those critical field relations that can be used to distinguish between the fundamental factors governing the formation of the basinal stratigraphy. Geologists tend to think of certain factors as the primary controls on the development of basin fills, including: source area uplift rate, source area lithology, rain fall, temperature regime, sea level changes, and basin subsidence rates. In contrast, basin models can address only the most basic controls, such as: basin subsidence, eustatics, volume of sediment supply (flux), sorting of the sediment supply, and rates of sediment transport. The problem is that the relationship between the geologists' factors and the modeling parameters are far from straight forward. For example, changes in rainfall might affect sorting, flux and transportability of the sediment supplied to a basin. Changes in source area lithology can affect all of the same modeling factors. Because of

Abstract One of the more long-lived debates in the earth sciences concerns the origin of sedimentary sequences observed on continental margins (Barrell, 1917; Pitman, 1978; Seuss, 1906; Sloss, 1962; Vail et al., 1977; Vella, 1965; Watts, 1982). As defined by Vail et al. (1977), these sequences are packages of conformable sediments, representing a time span of 1 to 10 My, that are bounded by unconformities or horizons that can be correlated with unconformities. Strata within individual sequences show continuous onlap onto the continental margin, whereas sequence boundaries mark abrupt seaward shifts in coastal onlap. Originally, Vail et al. (1977) argued that onlapping sediments record a marine transgression, whereas sequence boundaries signal an abrupt regression. Subsequently, more detailed studies of a small number of sequences (Vail et al., 1984) showed that onlapping sediments in the upper parts are nonmarine, suggesting that regressions are gradual rather than abrupt. Based on the apparent global synchroneity of many sequence boundaries, Vail et al. (1977) and Haq et al. (1987) propose that transgressions and regressions are caused by oscillatory eustatic variations. Some workers (Hallam, 1984; Miall, 1986) question whether or not the resolution of biostratigraphic correlations is sufficient to show that sequence boundaries are truly synchronous between basins. A fundamental difficulty in understanding the origin of sedimentary sequences is that transgressions and regressions can be caused by changes

Abstract Up until now, this course has mainly focused on subsidence in sedimentary basins—how it occurs, and how it can be determined as a function of time over the history of a basin. The pattern of subsidence in time and space largely determines the gross geometry of time-bounded sedimentary units because it controls the rate at which space is created for sedimentation (accommodation potential). On a finer scale, however, the pattern of subsidence plays a major role in determining the distribution of facies types within the sedimentary fill of the basin. The relation between patterns of subsidence and the internal character of the sediments can be understood through the use of basin-filling models, which we will study in this section. Basin-filling modeling, along with a variety of other quantitative techniques, has recently been referred to as “quantitative dynamic stratigraphy”, which can be read about in a new book by the same name (Cross, 1990). The development and application of basin-filling models is a subdiscipline that is still in its infancy, so the emphasis in this section of the course will be on the principles underlying various approaches, how the models are constructed, and the kinds of things they tell us, rather than on specific case studies and standard techniques. Because the necessary mechanics are understood better for streams than for any other sedimentary system, much of our effort will focus on sedimentation in alluvial systems. Given our present state of ignorance of many of the basic processes of basin

In this course we have summarized the techniques of quantitative basin analysis from the standpoint of modeling subsidence histories and mechanisms. We recommend becoming familiar with successful basin studies that have used these techniques to more fully understand the concepts that are only briefly covered here. The reference list accompanying these short course notes provides many, but by no means all, examples of subsidence analyses and will serve as a good starting point for more intense study. We also recommend attempting a subsidence analysis of your own on any basin of interest. The application of these techniques to basin analysis studies provides a relatively new, but very powerful, tool to understanding the role of tectonics and sea level in the formation of a sedimentary basin.

Abstract This publication is designed to introduce the concepts and techniques of quantitative modeling of basin subsidence histories. The book also describes some of the methods and results of modeling the development of sedimentary sequences generated by the interaction of subsidence, sediment supply, and sea-level changes. It concentrates on the theory and application of subsidence and stratigraphic modeling by working through specific examples from real or artificial basin sequences