Linking source and sink: Evaluating the balance between onshore erosion and offshore sediment accumulation since Gondwana break-up, South Africa

The geomorphic origin and evolution of the tectonically unique interior highland of southern Africa, the Kalahari Plateau, and its flanking low-lying coastal planes, remain largely unresolved because of a lack of regional quantitative analyses of its uplift and erosion history. Here we focus on the southern Cape, South Africa and link onshore denudation, based on new apatite fission track thermochronology results, to offshore sediment accumulation, using abundant well data and a seismic reflection profile. We attempt to relate source and sink in order to resolve some first order issues concerning timing of the exhumation and development of the topographic features of southern Africa. The volume of sediment accumulated off South Africa's south coast is calculated using 173 wells and a seismic reflection profile. A total, uncompacted, sediment volume of 268,500 km³ accumulated off South Africa's south coasts since  136 Ma, in the Outeniqua and Southern Outeniqua Basins. Accumulation volumes and rates were highest in the early Cretaceous (48,800 × 10⁴ km³ at  8150 km³/Ma from  136 to 130 Ma, and 57,500 × 10⁴ km³ at 5750 km³/Ma from  130 to 120 Ma) and mid–late Cretaceous (83,700 × 10⁴ km³ at 3200 km³/Ma from  93 to 67 Ma). Volumes and accumulation rates were lowest for the early–mid-Cretaceous (47,400 × 10⁴ km³ at 1750 km³/Ma from  120 to 93 Ma) and the Cenozoic (31,200 × 10⁴ km³ at 450 km³/Ma from  67 to 0 Ma). Although our analysis shows that the accumulated volume of offshore sediments does not match the calculated volume of onshore erosion, as quantified through apatite fission track thermochronology (e.g. Tinker, J.H., de Wit, M.J., Brown, R., 2008. Mesozoic exhumation of the 439 southern Cape, South Africa, quantified using apatite fission track thermochronology. Tectonophysics, doi: 10.1016/j.tecto.2007.10.009), the timing of increased sediment accumulation closely matches the timing of increased onshore denudation. This suggests that the greatest volumes of material were transported from source to sink during two distinct Cretaceous episodes, and that the processes driving onshore denudation decreased by an order of magnitude during the Cenozoic.     

Geochemistry of the Jurassic Mirdita Ophiolite (Albania) and the MORB to SSZ evolution of a marginal basin oceanic crust

The Middle Jurassic Mirdita Ophiolite in northern Albania is part of an ophiolite belt occurring between the Apulian and Pelagonian subcontinents in the Balkan Peninsula. The upper mantle and crustal units of the Mirdita Ophiolite show major changes in thickness, rock types, and chemical compositions from west to east as a result of its complex evolution in a suprasubduction zone (SSZ) environment. The  3–4-km-thick Western Mirdita Ophiolite (WMO) includes lherzolite–harzburgite, plagioclase–lherzolite, plagioclase–dunite in its upper mantle units and a plutonic complex composed of olivine gabbro, troctolite, ferrogabbro, and gabbro. These peridotites and gabbroic rocks are overlain directly by a  600-m-thick extrusive sequence containing basaltic pillow lavas and hyaloclastites. Sheeted dikes are rare in the WMO. The  12-km-thick Eastern Mirdita Ophiolite (EMO) includes tectonized harzburgite and dunite with extensive chromite deposits, as well as ultramafic cumulates including olivine clinopyroxenite, wehrlite, olivine websterite, and dunite forming a transitional Moho with the overlying lower crustal section. The plutonic rocks are made of pyroxenite, gabbronorite, gabbro, amphibole gabbro, diorite, quartz diorite, and plagiogranite. A well-developed sheeted dike complex has mutually intrusive relations with the underlying isotropic gabbros and plagiogranites and feeds into the overlying pillow lavas. Dike compositions change from older basalt to basaltic andesite, andesite, dacite, quartz diorite, to late-stage andesitic and boninitic dikes as constrained by crosscutting relations. The  1.1-km-thick extrusive sequence comprises basaltic and basaltic andesitic pillow lavas in the lower 700 m, and andesitic, dacitic and rhyodacitic massive sheet flows in the upper 400 m. Rare boninitic dikes and lavas occur as the youngest igneous products within the EMO. The basaltic and basaltic andesitic rocks of the WMO extrusive sequence display MORB affinities with Ti and Zr contents decreasing upsection (TiO₂ = 3.5–0.5%, Zr = 300–50 ppm), while _Nd(T) (+ 8 to + 6.5) varies little. These magmas were derived from partial melting of fertile MORB-type mantle. Fractional crystallization was important in the evolution of WMO magmas. The low Ti and HREE abundances and Cs and Ba enrichments in the uppermost basaltic andesites may indicate an increased subduction influence in the evolution of the late-stage WMO magmas. Basaltic andesites in the lower 700 m of the EMO volcanic sequence have lower TiO₂ ( 0.5%) and Zr ( 50 ppm) contents but _Nd(T) values (+ 7 to + 6.5) are similar to those of the WMO lavas. These rocks show variable enrichment in subduction-enriched incompatible elements (Cs, Ba, Th, U, LREE). The basaltic andesites through dacites and boninites within the upper 400 meters of EMO lavas show low TiO₂ ( 0.8–0.3%) and _Nd(T) (+ 6.5 to + 3.0). The mantle source of these rocks was variably enriched in Th by melts derived from subducted sediments as indicated by the large variations in Ba, K, and Pb contents. EMO boninitic dikes and lavas and some gabbroic intrusions with negative _{Nd (T)} values (− 1.4 and − 4.0, respectively) suggest that these magmas were produced from partial melting of previously depleted, ultra-refractory mantle. The MORB to SSZ transition (from west to east and stratigraphically upwards in the Mirdita Ophiolite and the progression of the _Nd(T) values from + 8.0 to − 4.0 towards the east resulted from an eastward shift in protoarc–forearc magmatism, keeping pace with slab rollback in this direction. The mantle flow above the retreating slab and in the arc-wedge corner played a major role in the evolution of the melting column, in which melt generation, aggregation/mixing and differentiation occurred at all levels of the sub-arc/forearc mantle. The SSZ Mirdita Ophiolite evolved during the intra-oceanic collapse and closure of the Pindos marginal basin, which had a protracted tectonic history involving seafloor spreading, protoarc rifting, and trench-continent collision.     

The Cerro Morado Andesites: Volcanic history and eruptive styles of a mafic volcanic field from northern Puna, Argentina

The Upper Miocene Cerro Morado Andesites constitutes a mafic volcanic field (100 km²) composed of andesite to basaltic andesite rocks that crop out 75 km to the east from the current arc, in the northern Puna of Argentina. The volcanic field comprises lavas and scoria cones resulting from three different eruptive phases developed without long interruptions between each other. Lavas and pyroclastic rocks are thought to be sourced from the same vents, located where orogen-parallel north-south faults crosscut transverse structures.The first eruptive phase involved the effusion of extensive andesitic flows, and minor Hawaiian-style fountaining which formed subordinate clastogenic lavas. The second phase represents the eruption of slightly less evolved andesite lavas and pyroclastic deposits, only distributed to the north and central sectors of the volcanic field. The third phase represents the discharge of basaltic andesite magmas which occurred as both pyroclastic eruptions and lava effusion from scattered vents distributed throughout the entire volcanic field. The interpreted facies model for scoria cones fits well with products of typical Strombolian-type activity, with minor fountaining episodes to the final stages of eruptions.Petrographic and chemical features suggest that the andesitic units (SiO₂ > 57%) evolved by crystal fractionation. In contrast, characteristics of basaltic andesite rocks are inconsistent with residence in upper-crustal chambers, suggesting that batches of magmas with different origins or evolutive histories arrived at the surface and erupted coevally.Based on the eruptive styles and lack of volcanic quiescence gaps between eruptions, the Cerro Morado Andesites can be classified as a mafic volcanic field constructed from the concurrent activity of several small, probably short-lived, monogenetic centers.     

Experimental frictional heating of coal gouge at seismic slip rates: Evidence for devolatilization and thermal pressurization of gouge fluids

High velocity (1 m/s) friction experiments on bituminous coal gouge display several earthquake-related phenomena, including devolatilization by frictional heating, gas pressurization, and slip weakening. Stage I is characterized by sample shortening and reduction in the coefficient of friction (μ) from  1 to 0.6. Stage II is characterized by high frequency ( 5 Hz) oscillations in stress and strain records and by gas emissions. Stage III is marked by rapid weakening (μ  0.1 to 0.35) and sample shortening, together with continued gas emissions. Stage IV produces stable stress records and continued weakness (μ  0.2), but without gas emission. Stage I shortening is due to compaction of the gouge and the weakening is attributed to mechanical or thermal effects. Stage II behavior is interpreted as due to coal gasification and fluctuations in fluid pressure, resulting in high frequency stick-slip type behavior. Dramatic reduction in shear stress in stage III is attributed to gas pressurization by pore collapse and corresponds to a frictional instability, analogous to nucleation of an earthquake. Microstructural observations indicate the deformation was brittle during stages I and II but ductile during stages III and IV. Time dependent finite element frictional heat models indicate the center of the samples became hot ( 900 °C) during stage II, whereas the edge of samples remained relatively cold (< 300 °C). Vitrinite reflectance of coal samples shows an increase in reflectance from  0.5 to  0.8% over the displacement interval 20–40 m (20–40 s), indicating that the reflectance responds to frictional heating on a short time scale. The energy expended per unit area in these low stress, large displacement experiments is similar to that of higher stress ( 50 MPa), short displacement ( 1 m) earthquakes ( 10⁷ J/m²).     

Mesozoic plutons of the Yidun Arc, SW China: U/Pb geochronology and Hf isotopic signature

The Yidun Arc is a Triassic volcanic arc located between the Songpan Garzê Fold Belt and the Qiangtang Block, southwest China. To constrain the age of a number of the major granitic plutons from the Yidun Arc, laser ablation ICP-MS U/Pb analysis of zircon was conducted. Hafnium isotope data was also acquired through laser-ablation multicollector ICPMS analysis of zircon, with the aim of gaining insight into the age and nature of the source region of the plutons. Three age groups have been identified from seven granite samples: Early–Middle Triassic ( 245 to 229 Ma), Late Triassic ( 219 to 216 Ma) and Cretaceous ( 105 to 95 Ma). Hafnium analysis shows the Triassic granites to have negative and variable ε_Hf values and Mesoproterozoic ( 1.6 Ga) depleted-mantle model ages, which is interpreted to reflect derivation from an isotopically heterogeneous, largely crustal source. The Cretaceous granite shows higher and less variable ε_Hf values and slightly younger model ages ( 1.3 Ga), and is interpreted to be derived from melting of a more homogeneous crustal source. A depleted-mantle model age of  1.5 Ga is calculated from the pooled Triassic and Cretaceous samples. The source region for these magmas may be tentatively correlated with Mesoproterozoic material of the Yangtze Craton, which has been suggested to underlie the Yidun Arc; however, further work is necessary to demonstrate this suggestion.     

A late Quaternary paleotemperature record from Hanging Lake, northern Yukon Territory, eastern Beringia

The late Quaternary paleoclimate of eastern Beringia has primarily been studied by drawing qualitative inferences from vegetation shifts. To quantitatively reconstruct summer temperatures, we analyzed lake sediments for fossil chironomids, and additionally we analyzed the sediments for fossil pollen and organic carbon content. A comparison with the δ¹⁸O record from Greenland indicates that the general climatic development of the region throughout the last glaciation–Holocene transition differed from that of the North Atlantic region. Between  17 and 15 ka, mean July air temperature was on average 5°C colder than modern, albeit a period of near-modern temperature at  16.5 ka. Total pollen accumulation rates ranged between  180 and 1200 grains cm^− 2 yr^− 1. At  15 ka, approximately coeval with the Bølling interstadial, temperatures again reached modern values. At  14 ka, nearly 1000 yr after warming began, Betula pollen percentages increased substantially and mark the transition to shrub-dominated pollen contributors. Chironomid-based inferences suggest no evidence of the Younger Dryas stade and only subtle evidence of an early Holocene thermal maximum, as temperatures from  15 ka to the late Holocene were relatively stable. The most recognizable climatic oscillation of the Holocene occurred from  4.5 to 2 ka.     

Brittle–ductile–brittle deformation during cooling of tonalite (Adamello, Southern Italian Alps)

Late- to post-magmatic deformation in slightly diachronous contiguous intrusions of the north-western Adamello batholith (Southern Alps, Italy) is recorded as, from oldest to youngest: (i) joints, (ii) solid-state ductile shear zones, (iii) faults associated with epidote-K-feldspar veins and (iv) zeolite veins and faults. Structures (ii) to (iv) are localized on the pervasive precursory network of joints (i), which developed during the earliest stages of pluton cooling. High temperature ( 500 °C), ductile overprinting of joints produced lineations, defined by aligned biotite and hornblende, on the joint surfaces and highly localized mylonites. The main phase of faulting, producing cataclasites and pseudotachylytes, occurred at  250 °C and was associated with extensive fluid infiltration. Cataclasites and pseudotachylytes are clustered along different E–W-striking dextral strike-slip fault zones correlated with the activity of the Tonale fault, a major tectonic structure that bounds the Adamello batholith to the north. Ductile deformation and cataclastic/veining episodes occurred at P = 0.25–0.3 GPa during rapid cooling of the batholith to the ambient temperatures ( 250 °C) that preceded the exhumation of the batholith. Timing of the sequence of deformation can be constrained by ³⁹Ar–⁴⁰Ar ages of  30 Ma on pseudotachylytes and various existing mineral ages. In the whole composite Adamello batholith, multiple magma pulses were intruded over the time span 42–30 Ma and each intrusive body shows the same ductile-to-brittle structural sequence localized on the early joint sets. This deformation sequence of the Adamello might be typical of intrusions undergoing cooling at depths close to the brittle–ductile transition.     

Neoproterozoic adakitic plutons in the northern margin of the Yangtze Block, China: Partial melting of a thickened lower crust and implications for secular crustal evolution

Numerous Neoproterozoic felsic and mafic–ultramafic intrusions occur in the Hannan region at the northern margin of the Yangtze Block. Among these, the Wudumen and Erliba plutons consist of granodiorites and have SHRIMP zircon U–Pb ages of  735 Ma. The rocks have high K₂O (0.8–3.6 wt.%) and Na₂O (4.4–6.4 wt.%) and low MgO (0.4–1.7 wt.%). They also have high Sr/Y (32–209) and (La/Yb)_n ratios (4.4–38.6). Their εNd values range from − 0.41 to − 0.92 and zircon initial ¹⁷⁶Hf/¹⁷⁷Hf ratios from 0.282353 to 0.282581. These geochemical features are similar to those of adakitic rocks produced by partial melting of a thickened lower crust. Our new analytical results, combined with the occurrence of voluminous arc-related mafic–ultramafic intrusions emplaced before 740 Ma, lead us to propose that the crustal evolution in the northern margin of the Yangtze Block during Neoproterozoic involved: (1) rapid crustal growth and thickening by underplating of mafic magmas from the mantle which was modified by materials coming from the subducting oceanic slab from  1.0 to  0.74 Ga, and (2) partial melting of the thickened lower crust due to a thermal anomaly induced by upwelling of asthenosphere through an oceanic slab window, producing the  735 Ma adakitic Wudumen and Erliba plutons. Our model suggests that the crustal thickness was more than 50 km at the northern margin of the Yangtze Block at  735 Ma, and rule out the possibility of a mantle plume impact causing the > 735 Ma magmatism in the region.     

Cretaceous and Cenozoic cooling history across the ultrahigh pressure Tongbai–Dabie belt, central China, from apatite fission-track thermochronology

The crystalline terrane of the Tongbai–Dabie region, central China, comprising the Earth's largest ultrahigh-pressure (UHP) exposure was formed during Triassic collision between the Sino–Korean and Yangtze cratons. New apatite fission-track (AFT) data presented here from the UHP terrane, extends over a significantly greater area than reported in previous studies, and includes the (eastern) Dabie, the Hong'an (northwestern Dabie) and Tongbai regions. The new data yield ages ranging from 44 ± 3 to 142 ± 36 Ma and mean track lengths between 10 and 14.4 μm. Thermal history models based on the AFT data taken together with published ⁴⁰Ar/³⁹Ar, K–Ar, apatite and zircon (U–Th)/He and U–Pb data, exhibit a three-stage cooling pattern that is similar across the study region, commencing with an Early Cretaceous rapid cooling event, followed by a period of relative thermal stability during which rocks remained at temperatures within the AFT partial annealing zone (60–110 °C) and ending with a possible renewed phase of accelerated cooling during Pliocene to Recent time. The first cooling phase followed large-scale transtensional deformation between 140 and 110 Ma and is related to Early Cretaceous eastward tectonic escape and Pacific back arc extension. Between this phase and the subsequent slow cooling phase, a transition period from 120 to 80 Ma (to 70 to 45 Ma along the Tan–Lu fault) was characterised by a relatively low cooling rate (3–5 °C/Ma). This transition is likely related to a tectonic response associated with the mid-Cretaceous subduction of the Izanagi–Pacific plate as well as lithospheric extension and thinning in eastern Asia. The present regional AFT age pattern is therefore basically controlled by the Early Cretaceous rapid cooling event, but finally shaped through active Cenozoic faulting. Following the transition phase the subsequent slow cooling phase pattern implies a net reduction in horizontal compressional stress corresponding to increased extension rates along the continental margin due to the decrease in plate convergence. Modelling of the AFT data suggests a possible Pliocene–Recent cooling episode, which may be supported by increased rates of sedimentation observed in adjacent basins. This cooling phase may be interpreted as a response to the far-field effects of the frontal India–Eurasia collision to the west. Approximate estimates suggest that the total amount of post 120 Ma denudation across the UHP orogen ranged from 2.4 to 13.2 km for different tectonic blocks and ranged from 0.8 to 9.7 km during the Cretaceous to between 1.7 and 3.8 km during the Cenozoic.     

Anatolian surface wave evaluated at GEOFON Station ISP Isparta, Turkey

We have studied seismic surface waves of 255 shallow regional earthquakes recently recorded at GEOFON station ISP (Isparta, Turkey) and have selected these 52 recordings with high signal-to-noise ratio for further analysis. An attempt was made by the simultaneous use of the Rayleigh and Love surface wave data to interpret the planar crust and uppermost mantle velocity structure beneath the Anatolian plate using a differential least-square inversion technique. The shear-wave velocities near the surface show a gradational change from approximately 2.2 to 3.6 km s^− 1 in the depth range 0–10 km. The mid-crustal depth range indicating a weakly developed low velocity zone has shear-wave velocities around 3.55 km s^− 1. The Moho discontinuity characterizing the crust–mantle velocity transition appears somewhat gradual between the depth range  25–45 km. The surface waves approaching from the northern Anatolia are estimated to travel a crustal thickness of  33 km whilst those from the southwestern Anatolia and part of east Mediterranean Sea indicate a thicker crust at  37 km. The eastern Anatolia events traveled even thicker crust at  41 km. A low sub-Moho velocity is estimated at  4.27 km s^− 1, although consistent with other similar studies in the region. The current velocities are considerably slower than indicated by the Preliminary Reference Earth Model (PREM) in almost all depth ranges.     

Insights into Wasatch fault vertical slip rates using the age of sediments in Timpanogos Cave, Utah

Timpanogos Cave, located near the Wasatch fault, is about 357 m above the American Fork River. Fluvial cave sediments and an interbedded carbonate flowstone yield a paleomagnetic and U–Th depositional age of 350 to 780 ka. Fault vertical slip rates, inferred from calculated river downcutting rates, range between 1.02 and 0.46 mm yr^− 1. These slip rates are in the range of the 0–12 Ma Wasatch Range exhumation rate ( 0.5–0.7 mm yr^− 1), suggesting that the long-term vertical slip rate remained stable through mid-Pleistocene time. However, the late Pleistocene (0–250 ka) decelerated slip rate ( 0.2–0.3 mm yr^− 1) and the accelerated Holocene slip rate ( 1.2 mm yr^− 1) are consistent with episodic fault activity. Assuming that the late Pleistocene vertical slip rate represents an episodic slowing of fault movement and the long-term (0–12 Ma) average vertical slip rate, including the late Pleistocene and Holocene, should be  0.6 mm yr^− 1, there is a net late Pleistocene vertical slip deficit of  50–75 m. The Holocene and late Pleistocene slip rates may be typical for episodes of accelerated and slowed fault movement, respectively. The calculated late Pleistocene slip deficit may mean that the current accelerated Wasatch fault slip rate will extend well into the future.     

Paleomagnetic constraints on the paleogeography and oroclinal bending of the Devonian volcanic arc in Kazakhstan

Numerical modelling, incorporating coupling between surface processes and induced flow in the lower continental crust, is used to address the Quaternary evolution of the Gulf of Corinth region in central Greece. The post-Early Pleistocene marine depocentre beneath this Gulf overlies the northern margin of an older (Early Pleistocene and earlier) lacustrine basin, the Proto Gulf of Corinth Basin or PGCB. In the late Early Pleistocene, relief in this region was minimal but, subsequently, dramatic relief has developed, involving the creation of  900 m of bathymetry within the Gulf and the uplift by many hundreds of metres of the part of the PGCB, south of the modern Gulf, which forms the Gulf's main sediment supply. It is assumed that, as a result of climate change around 0.9 Ma, erosion of this sediment source region and re-deposition of this material within the Gulf began, both processes occurring at spatial average rates of  0.2 mm a^− 1. Modelling of the resulting isostatic response indicates that the local effective viscosity of the lower crust is  4 × 10¹⁹ Pa s, indicating a Moho temperature of  560 °C. It predicts that the  10 mm a^− 1 of extension across this  70 km wide model region, at an extensional strain rate of  0.15 Ma^− 1, is partitioned with  3 mm a^− 1 across the sediment source,  2 mm a^− 1 across the depocentre, and  5 mm a^− 1 across the ‘hinge zone’ in between, the latter value being an estimate of the extension rate on normal faults forming the major topographic escarpment at the southern margin of the Gulf. This modelling confirms the view, suggested previously, that coupling between this depocentre and sediment source by lower-crustal flow can explain the dramatic development in local relief since the late Early Pleistocene. The effective viscosity of the lower crust in this region is not particularly low; the strong coupling interpreted between the sediment source and depocentre results instead from their close proximity. In detail, the effective viscosity of the lower crust is expected to decrease northward across this model region, due to the northward increase in exposure of the base of the continental lithosphere to the asthenosphere; in the south the two are separated by the subducting Hellenic slab. The isostatic consequences of such a lateral variation in viscosity provide a natural explanation for why, since  0.9 Ma, the modern Gulf has developed asymmetrically over the northern part of the PGCB, leaving the rest of the PGCB to act as its sediment source.     

Cretaceous–Tertiary geology of the Gangdese Arc in the Linzhou area, southern Tibet

Knowledge of the Cretaceous–Tertiary history of upper crustal shortening and magmatism in Tibet is fundamental to placing constraints on when and how the Tibetan plateau formed. In the Lhasa terrane of southern Tibet, the widely exposed angular unconformity beneath uppermost Cretaceous–lower Tertiary volcanic-bearing strata of the Linzizong Formation provides an excellent geologic and time marker to distinguish between deformation that occurred before vs. during the Indo-Asian collision. In the Linzhou area, located  30 km north of the city of Lhasa, a > 3-km-thick section of the Linzizong Formation lies unconformably on Cretaceous and older rocks that were shortened by both northward- and southward-verging structures during the Late Cretaceous. The Linzizong Formation dips northward in the footwall of a north-dipping thrust system that involves Triassic–Jurassic strata and a granite intrusion in the hanging wall. U–Pb zircon geochronologic studies show that the Linzizong Formation ranges in age from 69 Ma to at least 47 Ma and that the hanging wall granite intrusion crystallized at  52 Ma, coeval with dike emplacement into footwall Cretaceous strata. ⁴⁰Ar/³⁹Ar thermochronologic studies suggest slow cooling of the granite between 49 and 42 Ma, followed by an episode of accelerated cooling to upper crustal levels beginning at  42 Ma. The onset of rapid cooling was coeval with the cessation of voluminous arc magmatism in southern Tibet and is interpreted be a consequence of either (1) Tertiary thrusting in this region or (2) regional rock uplift and erosion following removal of overthickened Gangdese arc lower crust and upper mantle or break-off of the Neo-Tethyan oceanic slab.     

Constraining models of the tectonic setting of the giant Olympic Dam iron oxide–copper–gold deposit, South Australia, using deep seismic reflection data

In the Gawler Craton, the completeness of cover concealing the crystalline basement in the region of the giant Olympic Dam Cu–Au deposit has impeded any sufficient understanding of the crustal architecture and tectonic setting of its IOCG mineral-system. To circumvent this problem, deep seismic reflection data were recently acquired from  250 line-km of two intersecting traverses, centered on the Olympic Dam deposit. The data were recorded to 18 s TWT ( 55 km). The crust consists of Neoproterozoic cover, in places more than 5 km thick, over crystalline basement with the Moho at depths of 13–14 s TWT ( 40–42 km). The Olympic Dam deposit lies on the boundary between two distinct pieces of crust, one interpreted as the Archean–Paleoproterozoic core to the craton, the other as a Meso–Neoproterozoic mobile belt. The host to the deposit, a member of the  1590 Ma Hiltaba Suite of granites, is situated above a zone of reduced impedance contrast in the lower crust, which we interpret to be source-region for its  1000 °C magma. The crystalline basement is dominated by thrusts. This contrasts with widely held models for the tectonic setting of Olympic Dam, which predict extension associated with heat from the mantle producing the high temperatures required to generate the Hiltaba Suite granites implicated in mineralization. We use the seismic data to test four hypotheses for this heat-source: mantle underplating, a mantle-plume, lithospheric extension, and radioactive heating in the lower crust. We reject the first three hypotheses. The data cannot be used to reject or confirm the fourth hypothesis.     

Analysis of crustal deformation in Luzon, Philippines using geodetic observations and earthquake focal mechanisms

We utilize regional GPS velocities from Luzon, Philippines, with focal mechanism data from the Harvard Centroid Moment Tensor (CMT) Catalog, to constrain tectonic deformation in the complex plate boundary zone between the Philippine Sea Plate and Eurasia (the Sundaland block). Processed satellite imagery and digital elevation models are used with existing gravity anomaly, seismicity, and geologic maps to define a suite of six elastic blocks. Geodetic and focal mechanism data are inverted simultaneously to estimate plate rotations and fault-locking parameters for each of the tectonic blocks and faults comprising Luzon. Major tectonic structures that were found to absorb the plate convergence include the Manila Trench (20–100 mm yr^− 1) and East Luzon Trough ( 9–15 mm yr^− 1)/Philippine Trench ( 29–34 mm yr^− 1), which accommodate eastward and westward subduction beneath Luzon, respectively; the left-lateral strike-slip Philippine Fault ( 20–40 mm yr^− 1), and its northward extensions, the Northern Cordillera Fault ( 17–37 mm yr^− 1 transtension), and the Digdig Fault ( 17–27 mm yr^− 1 transpression). The Macolod Corridor, a zone of active volcanism, crustal thinning, extension, and extensive normal and strike-slip faulting in southwestern Luzon, is associated with left-lateral, transtensional slip of  5–10 mm yr^− 1. The Marikina Fault, which separates the Central Luzon block from the Southwestern Luzon block, reveals  10–12 mm yr^− 1 of left-lateral transpression. Our analysis suggests that much of the Philippine Fault and associated splays are locked to partly coupled, while the Manila and Philippine trenches appear to be poorly coupled. Luzon is best characterized as a tectonically active plate boundary zone, comprising six mobile elastic tectonic blocks between two active subduction zones. The Philippine Fault and associated intra-arc faults accommodate much of the trench-parallel component of relative plate motion.     

Evolution of mineralizing fluids in the Zn–Pb–Cu(–Ag ± Au) skarn and epithermal deposits of the world-class San Martín district, Zacatecas, Mexico

Mineralizing fluids at the San Martín skarn show an evolution characterized by prograde and retrograde associations. The prograde mineral associations consist of (1) a massive garnet zone, (2) a tremolite ± garnet zone, and (3) a late association of quartz, sphalerite, calcite and fluorite lining the vugs in the garnet zone. The fluids of the prograde associations exhibit decreasing temperatures of homogenization (Th) and variable salinities. The fluids of the massive garnet zone have salinities of 36 wt.% NaCl equiv. and Th of 645 to 570 °C, corresponding to pressures of 1055 bar. At the tremolite ± garnet zone, Th range from 438 to 354 °C. In the late association at the endoskarn, the following evolution can be drawn: (a) salinities of 50 to 42 wt.% NaCl equiv., and Th of 455 to 346 °C in quartz, (b) salinities of 46 wt.% NaCl equiv., and Th of 415 to 410 °C in sphalerite, (c) salinities of 50 to 37 wt.% NaCl equiv., and Th of 479 to 310 °C in calcite, (d) salinities of 33 to 28 wt.% NaCl equiv. and of 24 to 22 wt.% KCl in fluorite, and (e) two types of fluids with salinities of 2 and 39 wt.% NaCl equiv. and Th 344 and 300 °C, respectively, in later saccharoidal quartz segregations. The retrograde mineral associations comprise pervasive propylitic alteration to carbonization, and mantos with sulfides. Fluids in epidote have salinities of 7.6 wt.% NaCl equiv. and Th of 287 to 252 °C, and in calcite have salinities of 9.2 to 1 wt.% NaCl equiv. and Th of 188 to 112 °C. Fluids in the sulfide assemblages in the mantos have salinities of 8 to 3 wt.% NaCl equiv. and Th 300 °C, with corresponding pressures of 94 bar. Fluids in late epithermal veins close to the intrusive body have salinities of 10 to 5 wt.% NaCl equiv. and Th of 275 to 200 °C, and distal veins show salinities of 2 to 1 wt.% NaCl equiv. and Th of 160 °C.     

Finite-fault slip models for the 15 April 1979 (Mw 7.1) Montenegro earthquake and its strongest aftershock of 24 May 1979 (Mw 6.2)

We revisit the April 1979 Montenegro earthquake sequence to invert for finite-fault slip models for the mainshock of 15 April 1979 (M_w 7.1) and of the strongest aftershock of 24 May 1979 (M_w 6.2) using P, SH and SV waveforms, retrieved from IRIS data center. We also used body waveform modelling inversion to confirm the focal mechanism of the mainshock as a pure thrust mechanism and rule out the existence of considerable strike slip component in the motion. The mainshock occurred along a shallow (depth 7 km), low angle (14°) thrust fault, parallel to the coastline and dipping to the NE. Our preferred slip distribution model for the mainshock indicates that rupture initiated from SE and propagated towards NW, with a speed of 2.0 km/s. Moment was released in a main slip patch, confined in an area of L  50 km × W  23 km. The maximum slip ( 2.7 m) occurred  30 km to the NW of the hypocenter (location of rupture initiation). The average slip is 49 cm and the total moment release over the fault is 4.38e19 Nm. The slip model adequately fits the distribution of the M_w ≥ 4.3 aftershocks, as most of them are located in the regions of the fault plane that did not slip during the mainshock. The 24 May 1979 (M_w 6.2) strongest aftershock occurred  40 km NW of the mainshock. Our preferred slip model for this event showed a characteristic two-lobe pattern, where each lobe is  7.5 × 7.5 km². Rupture initiated in the NW lobe, where the slip obtained its maximum value of 45 cm, very close to the hypocenter, and propagated towards the south-eastern lobe where it reached another maximum value — for this lobe — of 30 cm, approximately 10 km away from the hypocenter. To indirectly validate our slip models we produced synthetic PGV maps (Shake maps) and we compared our predictions with observations of ground shaking from strong motion records. All comparisons were made for rock soil conditions and in general our slip models adequately fit the observations especially at the closest stations where the shaking was considerably stronger. Through the search of the parameter space for our inversions we obtained an optimum location for the mainshock at 42.04°N and 19.21° E and we also observed that better fit to the observations was obtained when the fault was modeled as a blind thrust fault.     

Seismic properties of the upper crust in the central Friuli area (northeastern Italy) based on petrophysical data

The compressional and shear wave velocities have been measured at room temperature and pressure up to 450 MPa on 5 sedimentary rock samples, representative of the most common lithologies of the upper crust in the central Friuli area (northeastern Italy). At 400 MPa confining pressure the Triassic dolomitic rock shows the highest velocities (Vp  7 km/s, Vs  3.6 km/s), the Jurassic and Triassic limestones samples intermediate velocities (Vp  6.3 /s, Vs  3.5 km/s) and the Cenozoic and Paleozoic sandstones the lowest velocities (Vp  6.15 km/s, Vs  3.35 km/s). The Paleozoic sandstone sample is characterized by the strongest anisotropy (10%) and significant birefringence (0.2 km/s) is found only on the Cenozoic sandstone sample. We elaborated the synthetic profiles of seismic velocities, density, elastic parameters and reflection coefficient, related to 4 one-dimensional geological models extended up to 22 km depth. The synthetic profiles evidence high rheological contrasts between Triassic dolomitic rocks and the soft sandstones and the Jurassic limestones. The Vp profiles obtained from laboratory measurements match very well the in-situ Vp profile measured by sonic log for the limestones and dolomitic rocks, supporting our one-dimensional modelling of the calcareous-carbonatic stratigraphic series. The Vp and Vs values of the synthetic profiles are compared with the corresponding ones obtained from the 3-D tomographic inversion of local earthquakes. The laboratory Vp are generally higher than the tomographic ones with major discrepancies for the dolomitic lithology. The comparison with the depth location of seismicity reveals that the seismic energy is mainly released in correspondence of high-contrast rheological boundaries.     

The age and petrology of the Chimbadzi Hill Intrusion, NW Zimbabwe: first evidence for early Paleoproterozoic magmatism in Zimbabwe

The mafic-ultramafic Chimbadzi Hill intrusion in the NW of the Zimbabwe craton is a dyke with inward-dipping margins comprising magnetite peridotite, troctolite and magnetite melatroctolite. The magnetite peridotite is composed of about equal amounts of V- and Ti-bearing magnetite and olivine (Fo₆₀). The troctolite is composed of about 50% olivine (Fo_50-54), 40% plagioclase (An_53-58), 7% clinopyroxene and minor apatite and magnetite with ilmenite lamellae. Geochemical trends suggest that the Chimbadzi Hill Intrusion formed by fractional crystallisation from a single initial magma. However, the more primitive magnetite peridotite overlies the more evolved troctolite in the intrusion. This ‘apparent’ inverted stratigraphy may be due to emptying of a fractionated magma chamber from the top, or to floor subsidence during intrusion.U–Pb dating on baddeleyite reveals that the age of the Chimbadzi Hill Intrusion is 2262 ± 2 Ma. This age does not correspond to any known tectono-thermal event in the Zimbabwe Craton or adjacent metamorphic belts. It is 300 Ma younger than the late Archean Great Dyke, and 230 Ma older than other Paleoproterozoic events in and around the craton. Therefore, it may represent a so far undocumented very early Proterozoic igneous event in the Zimbabwe Craton. The intrusion represents a vanadium resource for Zimbabwe, with titanium potentially being mined as by-product.     

Age and origin of middle Neoproterozoic mafic magmatism in southern Yangtze Block and relevance to the break-up of Rodinia

SHRIMP U–Pb zircon age, geochemical and Sm–Nd isotopic data are reported for mid-Neoproterozoic volcanic rocks and mafic intrusions in northern Guangxi (Guibei) and western Hunan (Xiangxi) Provinces along the southern margin of the Yangtze Block. The mafic igneous rocks studied are generally synchronous, dated at  765 Ma. The least-contaminated dolerite samples from Xiangxi are characterized by high εNd(T) value of 3.3 to 5.3 and OIB-type geochemical features, indicating that they were derived from an OIB-like mantle source in a continental rift setting. The spilites and gabbros in Guibei show basaltic compositions transitional between the tholeiitic and calc-alkaline series. Despite depletion in Nb and Ta relative to La and Th, they have Zr/Sm = 27–35 and Ti/V = 30–40, affinitive to intraplate basalts. Their εNd(T) values are variable, ranging from − 1.2 to 3.2 for the spilites and from − 1.7 to 2.9 for the gabbros, suggesting that these spilites and gabbros crystallized from crustal-contaminated mafic magmas derived from a metasomatised subcontinental lithospheric mantle source. We conclude that the  765 Ma mafic magmatic rocks in Guibei and Xiangxi were formed in a single continental rift setting as part of the broadly concurrent  780–750 Ma rift magmatism over much of South China, which may be related to the plume activities during the breakup of Rodinia.