Numerical modelling study on the flexural uplift of the Transantarctic Mountains

In this study, based on a 2-D thermomechanical finite element model, the uplift of the Transantarctic Mountains (TAM) is discussed in relation to the flexural uplift of a rheologically layered lithosphere, which is described by Vening-Meinesz's cantilever kinematics. The general model behaviour shows that the thickness of the crust and the geothermal gradient in the lithosphere are the principal factors in controlling the effective elastic thickness ( T _e). Although T _e is also significantly dependent on the magnitude of the uplift and the wet or dry rheological condition of rocks, these two factors do not have a dominant influence on the half-wavelength of the TAM. The model with a plausible crustal structure beneath Antarctica shows that the thermal structure beneath East Antarctica is the critical factor, controlling the half-wavelength of the TAM. If there is a significant radiogenic heat source in the Antarctic lithosphere, T _e beneath East Antarctica is estimated to be 50 km, at most, and the lithosphere has no potential to explain an exceptionally large-scale half-wavelength of the TAM. Even for the model without a heat source, if East Antarctica is significantly thermally influenced by West Antarctica, T _e is estimated to be about 60 km, and it is difficult to reproduce the half-wavelength of the TAM. Contrarily, when a radiogenic heat source is absent and the thermal structure beneath East Antarctica is not significantly affected by that beneath West Antarctica, the rheological structure beneath East Antarctica has the potential to reproduce the half-wavelength of the TAM ( T _e∼ 100 km). Thus, the presence of a radiogenic heat source in the crust and mantle and the thermal influence of West Antarctica on East Antarctica are crucial factors in the reproduction of the half-wavelength found in the TAM.     

Uplift and heat flow following the injection of magmas into the lithosphere

Summary. The thermal effect of a rapid injection of hot magmas into the lower part of the lithosphere is modelled as an increase in heat production through the invaded region. The change in surface heat flow and the uplift resulting from the thermal expansion are determined in three-dimensional axially symmetric geometry: they are expressed as the space time convolutions of a Green's function with the anomalous heat production.
The anomalies with shorter wavelength (compared to the lithospheric thickness) are attenuated. This filtering affects the surface uplift more than the heat flow anomaly; the attenuation effect is larger when only the lower part of the lithosphere is invaded.
The uplift time constant is of the same order as the heat conduction time if the lower lithosphere is invaded by magmas at a moderate rate (i.e. the rate of injection does not exceed the equivalent of 0.1 per cent of the lithospheric volume in 10⁶yr). Fifty per cent of the total uplift takes place in about 80 × 10⁶yr for a lithosphere 100 km thick. The uplift is slightly faster when the whole lithosphere is invaded. The heat flow anomaly is delayed when the lower part of the lithosphere is invaded.
The spatial extent and the timing of the uplift and heat flow anomalies are critical in determining the mechanism's feasibility. Magma injections explain rapid uplifts [> 100 m (10⁶ yr)⁻¹] only if the magma is supplied at a very high rate (i.e. at least 10 per cent of the lithosphere volume per 10⁶yr). It is a feasible mechanism for uplifts that occur over longer periods of time (≊ 30 × 10⁶yr) such as those that seem to have occurred when the African plate came to rest with respect to the mantle.     

Modelling sedimentation in ocean trenches: the Nankai Trough from 1 Ma to the present

Abstract Sediment accumulation within ocean trenches located at actively accreting convergent margins is determined by an interplay between sediment supply, sediment subduction/accretion at the toe of the overriding accretionary complex and the rate of subduction. Modelling trench sedimentation provides insight into the principal controlling factors, and a means of deriving, from the pattern of sedimentation, how factors, such as the rates of sediment supply and subduction, have varied over the period of accumulation of the trench sediments. Two DSDP-ODP drill sites within the Nankai Trough reveal a coarsening-upward megasequence, indicating a progressive facies transition from abyssal muds to outer-trench silts to inner-trench sands. The changing geometry of the trench-wedge over the past 1 Myr has been determined by modelling variations in net sediment flux for two trench-perpendicular profiles. The models were constrained to fit the stratigraphy at the drill sites, and the simulated present-day geometries of the trench were matched with those shown on seismic reflection profiles by successive adjustment of the model. Results from both sites confirm a ‘slow’ subduction rate of <20 km Myr^-1. At the south-western site (582), the width of the trench-wedge has ranged from 13 to 21 km over the past 1 Myr. To the north-east, at Site 808, the width has ranged from 7 to 13 km over the past 0.5 Myr. These changes in trench-wedge width are primarily the result of large changes in sediment supply rate. The subduction of the Shikoku Ridge, a fossil spreading centre adjacent to Site 808, has had a major influence on the style of sedimentation within the trench. The style of accretion from the trench to the toe of the accretionary complex has important implications for geometrical adjustment of the trench-wedge. Thrust displacement lifts the protothrust region out of the trench, resulting in a decreased width. This is followed by a phase of increasing width as the trench-wedge adjusts towards a new equilibrium. The cyclical, episodic accretion process results in a periodic second-order variation in trench-fill size that is superimposed on primary trends determined by variations in sediment supply rates and subduction rates over time.     

Regional variations of heat flow differences with depth in Alberta, Canada

Summary. Differences between estimated average heat flow values for the Mesozoic and Cenozoic formations ( Q ₁) and estimated average heat flow values for the Palaeozoic formations below the erosional unconformity ( Q ₂) are calculated for the Alberta part of the western Canadian sedimentary basin. Significant heat flow differences exist for these two intervals and the map of Δ Q = Q ₁– Q ₂ shows that Q ₂ is generally greater than Q ₁ in the western and south-western part of Alberta, while in the northern part of the province Q ₂ is generally less than Q ₁. The regional variations of Δ Q are large, with standard deviation of 26 mW m⁻² and average value –13.5 mW m⁻². A regional trend of Δ Q correlates with topographic relief and the hydraulic head variations in the basin. It is shown that there is a heat flow increase with depth in water recharge areas and a decrease in heat flow with depth in the low topographic elevation water discharge areas when comparing the average heat flow in Mesozoic + Cenozoic and Palaeozoic formations.     

Thermal evolution of the oceanic crust; its dependence on spreading rate and effect on crustal structure

Summary. The temperature field and rates of cooling and solidification of the oceanic crust and upper mantle at an ocean ridge have been calculated as a function of spreading rate. The thermal model of the accretion process incorporates latent heat release associated with solidification of the basalt. liquid forming the ocean crust and uses a heat supply boundary condition on the vertical ridge axis model boundary. It is assumed that while oceanic layer 2 cools rapidly by hydrothermal circulation, oceanic layer 3 cools predominantly by conduction. Basalt liquid injection into the upper part of oceanic layer 3 is shown to solidify instantaneously while that injected into lower crustal levels takes up to 0.4 Myr to solidify. Material solidifying instantaneously is interpreted as corresponding to the dolerite unit of the ocean crust while that taking a finite time to cool is interpreted as corresponding to the gabbroic unit. The rate of cooling of the crust is shown to be faster for slower spreading rates and consequently the thicknesses of the dolerite and gabbro units are predicted to thin and thicken respectively with increase in spreading rate. The width of the molten region, or magma chamber, within the crust at the ridge axis is shown to be approximately proportional to spreading rate with chamber half widths of 1.5 and 10.0 km for half spreading rate of 1.0 and 6.0 cm yr⁻¹. Below a critical half spreading rate of about 0.65 cm yr⁻¹ no molten region exists and the crust is entirely doleritic.     

Simple equations of sedimentation: applications to sequence stratigraphy

Abstract An equation to relate the thickness of sediment deposited (ΔSed), eustatic sea-level change (ΔE), and subsidence (ΔSub), to changes in depth of water (ΔD) is: ΔSub +ΔE-ΔSed =ΔD.
Using existing sea-level curves, the equation shows that some transgressive-regressive sequences in a foreland basin and a composite seismic facies sequence on a passive margin cannot result solely from eustatic variation. In each case, the space created by subsidence is greater than that provided by eustatic rise. However, eustatic variation could have triggered sequence development if superimposed on a basin with relatively constant values of the other parameters. Short-period sea-level fluctuations with high rates of change, exceeding 70–100 m Myr^-1 for periods less than 2–3 Myr, affect the stratigraphy and sedimentology more than longer period, higher amplitude variations.
Clinoforms are generated because of lateral variations in sedimentation rate compared to the rate of creation of accommodation space. These variations may result from differing sedimentation rates, subsidence rates, or rates of eustatic change, superimposed on a basin with lateral sediment supply. Clinoform slopes and curvatures are interpre table in terms of these variables as well as the type of sediment supplied and the energy distribution in the basin.
These equations put some well-known geological principles on a simple quantitative basis. They force precision in definition of variables, and may lead to further development of quantitative techniques in stratigraphy and sedimentology.     

Stress drops, wave spectra and recurrence intervals of great earthquakes — implications of the Etorofu earthquake of 1958 November 6

Summary . The great Etorofu earthquake of 1958 November 6 is characterized by a relatively small aftershock area (70 × 150 km²) and an extremely large felt area. The felt area is more extensive than those of any other large earthquakes which have occurred in the southern Kurile to northern Japan arc since the beginning of this century. The mechanism is a pure thrust fault typical of most great earthquakes in island arcs. A body wave magnitude of m _b = 8.2 is obtained at periods around 6 s using more than 40 observations, although an m _b value of only 7.6–7.7 would be expected empirically from the observed surface wave magnitude of M _s= 8.1–8.2. Both an unusually large felt area and a high m _b indicate a dominance of high-frequency components in the seismic waves. A seismic moment of M _o= 4.4 × 10²⁸ dyne cm is determined from long-period surface waves from which a high stress drop of Δσ = 78 bar is obtained using a relatively small aftershock area. Historic data indicate an anomalously long time interval between the 1958 event and any earlier great earthquake from the same source region. The observed high stress drop can be interpreted as a consequence of this long intervening period through which strain built up. The dominance of the high-frequency seismic waves can then be interpreted as a result of this high stress drop. Stress drops, seismic wave spectra and recurrence intervals of great earthquakes are in this way closely related to each other. The 1958 event may represent a high strength extreme of stochastic fluctuation of fracture strength relevant to great earthquakes.     

Estimating palaeorelief from detrital mineral age ranges

We propose a method that uses the increase in mineral age with elevation in some bedrock landscapes to quantify palaeotopographic relief from the age range of detrital minerals in coeval sediment. We use the rate at which mineral age changes with elevation (its age-gradient, d t /d z ) and its age range (Δ t ) in the sediment to invert for relief: Δ z =Δ t /(d t /d z ). Relief inversion requires a single-grain dating precision high enough that detrital grains originate from resolvably different elevations (e.g. laser microprobe ⁴⁰Ar/³⁹Ar fusion). The technique assumes that there is no change in mineral age during erosion and transport, that sediment is mixed well enough and (or) sampled sufficiently to capture the extrema of mineral ages, and that isochrons were horizontal during erosion. Subject to these constraints, inversion of the age range of individual grains in synorogenic sedimentary sequences allows quantitative estimation of relief development for eroded mountain ranges. This method provides the only direct quantitative measure of palaeorelief, a poorly constrained, but important aspect of many geological, geomorphological and geodynamic models.     

THERMAL CONVECTION IN THE INTERIOR OF THE EARTH

The paper deals with thermal convection in the shell of the earth, caused by various assumed zonal temperature perturbations. One temperature perturbation here treated is that due to the difference in the temperature distribution under a continental crust made up of 10 km. of granite on top of 20 km. of basaltic material and a sub-oceanic crust consisting of 25 km. of basalt. The kinematic viscosity v was assumed to be 3 × 1021, as estimated recently by N. A. Haskell from a study of the uplift of Fennoscandia after the ice load. It is found that when g, ν, p and a, the coefficient of volume expansion, are constant, the velocities and their gradients are proportional to the amplitude of the temperature perturbation and to (ag/ν), while the stresses are independent of the viscosity. The velocities are found to be of the order of I cm./year. The shearing stress exerted by the convective substratum on the crust is of the order of 10⁷ dyn./cm.², while the normal stresses are about 10 times larger. The crust is pushed upwards under the warmer (continental) regions and pulled downwards under colder (oceanic) regions. The maximum stress-difference occurs at the bottom of the crust over the centre of the oceans or continents. The surface inequalities are nearly compensated.     

How does alluvial sedimentation at range fronts modify the erosional dynamics of mountain catchments?

At the geological time scale, the way in which the erosion of drainage catchments responds to tectonic uplift and climate changes depends on boundary conditions. In particular, sediment accumulation and erosion occurring at the edge of mountain ranges should influence the base level of mountain catchments, as well as sediment and water discharges. In this paper, we use a landform evolution model (LEM) to investigate how the presence of alluvial sedimentation at range fronts affects catchment responses to climatic or tectonic changes. This approach is applied to a 25 km × 50 km domain, in which the central part is uplifted progressively to simulate the growth of a small mountain range. The LEM includes different slope and river processes that can compete with each other. This competition leads to ‘transport‐limited’, ‘detachment‐limited’ or ‘mixed’ transport conditions in mountains at dynamic equilibrium. In addition, two end‐member algorithms (the channellized‐flow and the sheet‐flow regimes) have been included for the alluvial fan‐flow regime. The three transport conditions and the two flow algorithms represent six different models for which the responses to increase of rock uplift rate and/or cyclic variation of the precipitation rate are investigated. Our results indicate that addition of an alluvial apron increases the long‐term mountain denudation. In response to uplift, mountain rivers adapt their profile in two successive stages; first by propagation of an erosion wave and then by slowly increasing their channel gradients. During the second stage, the erosion rate is almost uniform across the catchment area at any one time, which suggests that dynamic equilibrium has been reached, although the balance between erosion and rock uplift rates has not yet been achieved. This second stage is initiated by the uplift of the mountain river outlets because of sedimentation aggradation at the mountain front. The response time depends on the type of water flow imposed on the alluvial fans domains (× by 1.5 for channelized flow regime and by 10 for the sheet flow one). Cyclic variations of precipitation rate generate cyclic incisions in the alluvial apron. These incision pulses create knick‐points in the river profile in the case of ‘detachment‐limited’ and ‘mixed’ river conditions, which could be mistaken for tectonically induced knick‐points. ‘Transport‐limited’ conditions do not create such knick‐points, but nevertheless trigger erosion in catchments. The feedbacks linked to sedimentation and erosion at range front can therefore control catchment incision or aggradation. In addition, random river captures in the range front trigger auto‐cyclic erosion pulses in the catchment, capable of generating incision–aggradation cycles.     

The dynamical origin of subduction zone topography

Summary. Subduction zones are expressed topographically by long linear oceanic trenches flanked by a low outer rise on the seaward side and an island arc on the landward side. This topographic structure is reflected in free air gravity anomalies, suggesting that much of the topography originates from dynamical forces applied at the base of the crust. We have successfully reproduced the general topographic features of subduction zones by supposing that the stresses generated by the bending of the viscous lower lithosphere as it subducts are transmitted through the thin elastic upper portion of the lithosphere. The trench is due to a zone of extensional flow (associated with low pressure) in the upper part of the viscous lithosphere.
The stresses in the subducting slab are computed using a finite element technique, assuming a Maxwell viscoelastic constitutive relation. Various dips (10 to 90°) were investigated, as well as depth dependent and non-Newtonian (power law, n = 3) viscosities. Observed subduction zone dimensions are well reproduced by these models. The effective viscosity required at mid-depth in the lithosphere is about 6 × 10²² P. This low value is probably due to the stress dependence of the effective viscosity. However, these models also show that the topography of the subduction zone depends primarily upon the geometry of the subducting slab (dip, radius of curvature of the bend) rather than upon its rheology. Shear stresses beneath the trench reach maxima of approximately 50 MPa. An interesting feature of some solutions is a dynamically supported bench or platform between the trench and island arc.     

关于海洋岛屿成因分类的新意见

岛屿对于发展海洋经济,确定海洋权属,以及国家安全等方面都具有重要地位。因此,岛屿的成因分类研究具有重要理论意义和应用价值。传统上将岛屿分为大陆岛、海洋岛(火山岛与珊瑚岛)和冲积岛。这种分类已经不能适应现代地球学科的新进展。本文根据板块构造理论和大洋地貌体系,提出新的岛屿分类意见,即分为内力和外力两个成因系列,包括近岸大陆岛、隆起大陆岛、大陆火山岛、岛弧陆块岛、岛弧火山岛、俯冲增生岛、无震海岭火山岛、微型陆块岛、海山火山岛、中脊火山岛、构造断层岛、河口沙岛、障壁岛、侵蚀沙岛、珊瑚岛15个类型。     

Tectonic controls on Miocene sedimentation in the Southern Taranaki Basin and implications for New Zealand plate boundary deformation

Miocene strata in the southern Taranaki Basin (STB), up to 3 km thick, provide a distal record of erosion associated with plate boundary deformation in New Zealand. 2D and 3D seismic reflection data tied to drillhole stratigraphy have been used to constrain four main phases of basin development. These are: (a) Early Miocene (22–19 Ma) subsidence, dominantly bathyal water depths and deposition of minor submarine fans along the eastern basin margin. (b) Middle Miocene (19–14 Ma) widespread submarine fan deposition on a bathyal basin floor in the central STB. (c) Rapid Middle–Late Miocene (14–7 Ma) progradation of the shelf break northwards across the STB. (d) Widespread uplift and erosion of the STB during the latest Miocene–Pliocene (7–4.5 Ma). Bathyal water depths and fan deposition in the Early Miocene were influenced by vertical motions on major reverse faults and regional subsidence produced by subduction of the Pacific plate beneath northern New Zealand. Subsequent submarine fan deposition and northward shelf‐break progradation reflect increasing input of terrigenous material, primarily eroded from an uplifting region to the south of the STB. Sedimentation patterns in the STB are consistent with the age and locations of conglomerates deposited in onshore West Coast basins, related to this uplift and erosion. Sediment transport in the West Coast region was mainly parallel to NNE trending active reverse faults, and in the STB was perpendicular to the NE‐SW orientated shelf break, especially from ca. 14–7 Ma, when sedimentation rates exceeded fault‐displacement rates. Increases in sedimentation rates in the STB coincide with regional increases in the rates of shortening that appear to reflect plate boundary‐wide events and have been attributed to, or correlated with, increases in the plate convergence rate. Miocene sedimentation patterns in the STB thus reflect both intra‐basinal deformation and tectonic signals from the wider developing New Zealand plate boundary.     

A model of relief forming by tectonic uplift and valley incision in orogenesis

A model for the change in shape of interfluves by concurrent tectonics and denudation was developed based on the morphometric attributes of landforms. The tributaries flowing down valley side slopes and dissecting low-relief surfaces on the interfluves are one of the most important elements for relief-forming processes. They were named as β-tributaries. The valleyhead altitude ( H ), the junction altitude ( L ) and the valley length ( l ) of the β-tributary were measured. The altitudinal difference ( h=H−L ), which indicates the local relief, and the average slope of tributary (tan θ= h / L ) were calculated. The regression analyses among H , L , l , h and tan θ indicate that the valley length, relief and slope increase with an increase in valleyhead altitude.
Based on the functional relations above, the shape of interfluve is a function of uplift, altitude and erosion. This model is used to illustrate the change in cross-section of interfluves during a period of sustained rock uplift.
Successive changes in shape of interfluve can be divided into two substages: (1) the early substage, characterized by trapezoidal cross-sections with the original low-relief surfaces and residual, shallow stream networks on the ridges; and (2) the later substage, characterized by a triangular cross-section, with the original low-relief surfaces removed and with the interfluves lowered by headward erosion of β-tributaries. In central Japan, the transitional relief from the trapezoidal to the triangular cross-section appears when the ridges of interfluves attain elevations about 1600–2000 m above sea level.     

Assessing the utility of 210Pb geochronology for estimating sediment accumulation rates on river floodplains in Fiji

Low‐energy gamma ray spectroscopy has been employed to estimate floodplain sedimentation rates using measurements of ²¹⁰Pb in floodplain alluvium. The utility of the technique is assessed through the analysis of excess (unsupported) ²¹⁰Pb profiles in three sediment cores taken from the floodplain of the Labasa River on Vanua Levu in northern Fiji. A low‐energy germanium spectrometer (LEGe) was used for the nondestructive determination of excess ²¹⁰Pb in a region cultivated intensively with sugarcane. Measured average historical (c. 25 years) vertical accretion rates are between 2.2 and 4.4 cm yr^?1. The findings are broadly comparable with published sedimentation rates from analyses of radionuclide profiles elsewhere in the tropical South Pacific Islands, but the rates are higher than those measured previously at the same Labasa River sites using ¹³⁷Cs profiles. Accelerated soil erosion owing to cane burning and land tillage seems to be largely responsible for sediment production, although flood‐related effects such as channel accretion by coarse bedload and the emplacement of large organic debris also influence floodplain sedimentation. However, application of the ²¹⁰Pb technique in Fiji (and perhaps neighbouring island countries) is found to have serious drawbacks compared to the more robust ¹³⁷Cs method, owing principally to the low ²¹⁰Pb concentrations in the sandy alluvial sediment tested.     

Joint inversion of temperature,vitrinite reflectance and fission tracks in apatite with examples from the eastern North Sea area

As sediment accumulation indicates basin subsidence, erosion often is understood as tectonic uplift, but the amplitude and timing may be difficult to determine because the sedimentary record is missing. Quantification of erosion therefore requires indirect evidence, for example thermal indicators such as temperature, vitrinite reflectance and fission tracks in apatite. However, as always, the types and quality of data and the choice of models are important to the results. For example, considering only the thermal evolution of the sedimentary section discards the thermal time constant of the lithosphere and essentially ignores the temporal continuity of the thermal structure. Furthermore, the types and density of thermal indicators determine the solution space of deposition and erosion, the quantification of which calls for the use of inverse methods, which can only be successful when all models are mutually consistent. Here, we use integrated basin modelling and Markov Chain Monte Carlo inversion of four deep boreholes to show that the erosional pattern along the Sorgenfrei–Tornquist Zone (STZ) in the eastern North Sea is consistent with a tectonic model of tectonic inversion based on compression and relaxation of an elastic plate. Three wells in close proximity SW of the STZ have different data and exhibit characteristic differences in erosion estimates but are consistent with the formation of a thick chalk sequence, followed by minor Cenozoic erosion during relaxation inversion. The well on the inversion ridge requires ca. 1.7 km Jurassic-Early Cretaceous sedimentation followed by Late Cretaceous–Palaeocene erosion during inversion. No well demands thick Cenozoic sedimentation followed by equivalent significant Neogene exhumation. When data are of high quality and models are consistent, the thermal indicator method yields significant results with important tectonic and geodynamic implications.     

Lithospheric flexure due to prograding sediment loads: implications for the origin of offlap/onlap patterns in sedimentary basins

Abstract Simple elastic plate models have been used to determine the stratigraphic patterns that result from prograding sediment loads. The predicted patterns, which include coastal offlap/onlap and downlap in a basinward direction, are generally similar to observations of stratal geometry from Cenozoic sequences of the U.S. Atlantic and Gulf Coast margins. Coastal offlap is a feature of all models in which the water depth and elastic thickness of the lithosphere, T _e (which is a measure of the long-term strength of the lithosphere), are held constant, and is caused by a seaward shift in the sediment load and its compensation as progradation proceeds. The coastal offlap pattern is reduced if sediments prograde into a subsiding basin, since subsidence causes an increase in the accommodation space and loading landward of a prograding wedge. The stratal geometry that results is complex, however, and depends on the sediment supply, the amount of subsidence, and T _e. If the sediment supply to a subsiding basin proceeds in distinct 'pulses' (due, say, to different tectonic events in a source region) then it is possible to determine the relationship between stratal geometry and T _e. Coastal offlap and downlap are features of most models where the lithosphere either has a constant T _e slowly increases T_e with time, or changes T _e laterally; however, in the case where sediments prograde onto lithosphere that rapidly increases T _e with rime, the offlap can be replaced by onlap. Lithospheric flexure due to prograding sediment loads is capable of producing a wide variety of stratal geometries and may therefore be an important factor to take into account when evaluating the relative role of tectonics and eustatic sea-level changes in controlling the stratigraphic record.     

Jurassic synorogenic basin filling in western Korea: sedimentary response to inception of the western Circum-Pacific orogeny

This is the first sedimentologic and stratigraphic attempt to demonstrate Jurassic subduction-induced basin-filling processes in the early stage of the western Circum-Pacific orogeny. The Chungnam Basin in western Korea was filled with a Lower to Middle Jurassic nonmarine succession, the Nampo Group, whose deposition postdated the Triassic final assembly of Chinese continental blocks. The Nampo Group consists of two repeated, fining- to coarsening-upward alluvio-lacustrine sequences, separated by an interval of thick breccia–gravel progradation deposits and its related strong proximal unconformities. No temporal variation in the degree of chemical weathering, along with the predominance of coals and a tropic to subtropic paleoflora, reveals little or no climate fluctuations during deposition of the Nampo Group. The observed relationships provide a record of sedimentation most likely controlled by temporal variations of tectonically driven sediment flux. Such syntectonic sedimentation of the Chungnam Basin occurred at a convergent margin of continental-arc setting during the Daebo orogeny, synchronous with the early subduction of the western paleo-Pacific ocean that resulted in formation of an accretionary complex along the East Asian continental margin during Jurassic time. Hence, synorogenic deposition in the Chungnam Basin is interpreted as sedimentary response to subduction–accretion of the western paleo-Pacific plate.     

Simple analytic model for subduction zone thermal structure

A new analytic model is presented for the thermal structure of subduction zones. It applies to the deeper regions of a subduction zone, where the overriding mantle is no longer rigid but flows parallel to the slab surface. The model captures the development of one thermal boundary layer out into the mantle wedge, and another into the subducting slab. By combining this model with the analytic model of Royden (1993a , b ), which applies to regions in which the overriding plate is rigid, a nearly complete analytic model for the thermal structure of a steady-state subduction zone can be achieved. A good agreement is demonstrated between the output of the combined analytic model and a numerical finite element calculation. The advantages of this analytic approach include (1) efficiency (only limited computing resources are needed); (2) flexibility (non-linear slab shape, and processes such as erosion, and shear heating are easily incorporated); and (3) transparency (the effect of changes in input variables can be seen directly).     

The thermal field in a basin after a sudden passive pure shear lithospheric extension and sublithospheric mechanical erosion:the case of the Tuscan Basin (Italy)

A simple new model for sudden lithospheric thinning that considers the crust to be stretched and the lower layer of the lithosphere to be partially stretched and partially mechanically eroded is proposed. This model allows calculation of the thermal field of the lithosphere during the initial warming phase and the surface uplift.
Application of this model to the Tuscan Basin explains the high regional heat flux density values (>100 mW m⁻² ), the tectonic subsidence (about 1 km) and the average uplift (>400 m) observed in this region well.