Coseismic and post-seismic signatures of the Sumatra 2004 December and 2005 March earthquakes in GRACE satellite gravity

The GRACE satellite mission has been measuring the Earth's gravity field and its temporal variations since 2002 April. Although these variations are mainly due to mass transfer within the geofluid envelops, they also result from mass displacements associated with phenomena including glacial isostatic adjustment and earthquakes. However, these last contributions are difficult to isolate because of the presence of noise and of geofluid signals, and because of GRACE's coarse spatial resolution (>400 km half-wavelength). In this paper, we show that a wavelet analysis on the sphere helps to retrieve earthquake signatures from GRACE geoid products. Using a wavelet analysis of GRACE geoids products, we show that the geoid variations caused by the 2004 December ( M _w= 9.2) and 2005 March ( M _w= 8.7) Sumatra earthquakes can be detected. At GRACE resolution, the 2004 December earthquake produced a strong coseismic decrease of the gravity field in the Andaman Sea, followed by relaxation in the area affected by both the Andaman 2004 and the Nias 2005 earthquakes. We find two characteristic timescales for the relaxation, with a fast variation occurring in the vicinity of the Central Andaman ridge. We discuss our coseismic observations in terms of density changes of crustal and upper-mantle rocks, and of the vertical displacements in the Andaman Sea. We interpret the post-seismic signal in terms of the viscoelastic response of the Earth's mantle. The transient component of the relaxation may indicate the presence of hot, viscous material beneath the active Central Andaman Basin.     

Tectonic significance of the present relief of the Betic Cordillera

Tortonian calcarenites of the Betic Cordillera were deposited in coastal or very shallow marine environments and represent an ideal marker for estimating vertical movements from the late Miocene to the Present. A map showing the heights at which these Tortonian marine rocks are situated has a clear correlation with the present relief, indicating that today's relief has been formed since the Tortonian. There is also a good correlation between present relief and the Bouguer anomaly distribution in the Betic Cordillera, as well as with crustal thickness. Likewise, the present relief is directly related to the geodynamic setting of a horizontal N–S to NNW–SSE compression and an almost perpendicular extension, along with isostatic readjustment, existing in the Betic Cordillera from the Tortonian. As a result of these regional stresses, faults and folds have produced notable vertical movements. The highest rates of uplift of the Betic Cordillera coincide with large antiforms, in particular those of the Sierra Nevada and the Sierra Filabres. Several subsiding sectors also exist (for example, the Granada Basin or the Guadalquivir Basin). The foreland Guadalquivir Basin has a complex history because the uplift in its eastern sector and subsidence in the western sector coexisted during the late Tortonian. Today the whole Betic Cordillera is characterized by differential regional uplift, even in the aforementioned subsiding sectors.     

The wavelet method as an alternative for reducing ionospheric effects from single-frequency GPS receivers

The ionospheric effect is one of the major errors in GPS data processing over long baselines. As a dispersive medium, it is possible to compute its influence on the GPS signal with the ionosphere-free linear combination of L1 and L2 observables, requiring dual-frequency receivers. In the case of single-frequency receivers, ionospheric effects are either neglected or reduced by using a model. In this paper, an alternative for single-frequency users is proposed. It involves multiresolution analysis (MRA) using a wavelet analysis of the double-difference observations to remove the short- and medium-scale ionosphere variations and disturbances, as well as some minor tropospheric effects. Experiments were carried out over three baseline lengths from 50 to 450 km, and the results provided by the proposed method were better than those from dual-frequency receivers. The horizontal root mean square was of about 0.28 m (1σ).     

Impact of channel bar form and grain size variability on estimated stranding risk of juvenile brown trout during hydropeaking

Hydropeaking leads to artificial fluctuations in discharge and corresponding water levels with pronounced dewatering areas between base and peak flow along gravel bars and channel banks. In the present study, 16 hydropeaking reaches in Austria were investigated to assess possible differences in the estimated stranding risk for young of the year brown trout according to different gravel bar types and differences in microtopography roughness. Based on hydrodynamic‐numerical modelling, a predictive habitat modelling approach was implemented in the study design. Accompanied by grain size sampling along the various channel bars, a conceptual stranding risk model (SRM) was developed. The results showed that a high variability in estimated stranding risk exists for the tested sites considering discharge ratios of 1:3, 1:5 and 1:10. With respect to the discussion of establishing legal thresholds for ramping ratios in discharge, it was documented that, exemplarily, a discharge ratio base flow/peak flow of 1:5 (winter base flow conditions) could cause minor differences in the spatial extent of dewatering areas and the related estimated stranding risk for juvenile brown trout compared to a ratio of 1:2 for summer base flow conditions. Microtopographic roughness was addressed due to sampling and analysis of grain size distributions. Statistical testing of grain size distributions revealed significant differences between the surface material compositions of the investigated gravel bars. Those differences are evident, particularly for the coarser fraction (d₉₀), which is important as cover for young of the year brown trout. These aspects of grain size in habitat use and hydraulics have been addressed in the conceptual SRM. The results showed that point bar morphology, in particular, was less sensitive to the risk of stranding compared to, for example, alternating gravel bars. Considering the multiple pressures for alpine rivers, the improvement of structural features due to bar formation and related self‐forming processes is discussed as a possible alternative for future mitigation measures to reduce the negative impacts of hydropeaking. Copyright © 2014 John Wiley & Sons, Ltd.     

Stress Changes Modelled for the Sequence of Strong Earthquakes in the South Iceland Seismic Zone Since 1706

The South Iceland seismic zone is, roughly speaking, situated between two sections of the mid-Atlantic ridge, i.e., the Reykjanes Ridge southwest of Iceland and the Eastern Volcanic Zone on the island. It is a transform zone, where earthquakes are expected to occur on E-W-trending left-lateral shear faults, equivalent to conjugate, N-S-oriented right-lateral, rupture planes. In fact, earthquakes take place on en-échelon N-S-oriented faults, which is indicated by the distribution of main shock intensities, aftershocks as well as by surface fault traces. The stress field continuously generated in the fault zone by opening of the adjacent ridges is computed and superimposed on the stress field changes induced by a series of 13 earthquakes (M 6) between 1706 and 2000. The level of the pre-seismic stress field is analysed as well as the size of the area under high stress. Finally, the post-seismic stress field of June 2000 is analysed, to see where high stresses might have accumulated. The modelling indicates that the rupture planes located on separated parallel N-S-striking zones are dense enough to lead to an area-wide stress release by the series of events. The obtained pre-seismic stress level for most events is high and stable with the exception of situations when several strong shocks occur over a time span of several days, i.e., display typical main shock-aftershock patterns. The size of areas under high stress aside from of the rupture plane, i.e., where no event occurs at the specific time, is of medium to small size.     

A review of new interpretations of the tectonostratigraphy,geochemistry and evolution of the Onverwacht Suite,Barberton Greenstone Belt,South Africa

The Paleoarchean (ca. 3.5–3.3 Ga) Onverwacht Suite (OS) of the Barberton Greenstone Belt consists of a 15‐km thick imbricate tectonic stack of seven complexes consisting predominantly of volcanic rocks and intrusions. Tectonostratigraphically from base to top they are the Sandspruit, Theespruit, Komati, Hooggenoeg, Noisy, Kromberg and Mendon Complexes. The Hooggenoeg and Noisy Complexes in the middle of the OS are separated by a significant unconformity resulting from the uplift of the submarine lavas and deep erosion, demonstrating the onset of tectonic accretion prior to 3455 Ma. The basic lavas of the tectonostratigraphic lower (Theespruit, Sandspruit and Komati) and upper (Mendon) complexes are composed of komatiite, komatiitic basalt and high-MgO basalt, whereas those in the middle part (Hooggenoeg and Kromberg) are predominantly high- to low-MgO tholeiitic basalts. Felsic volcanic rocks and intrusions are important in two of the complexes (Theespruit and Noisy). The ultramafic to basaltic lavas show REE patterns that are almost flat and resemble those of modern MORB, whereas those of the felsic rocks are flat from Lu to Gd and moderately to strongly enriched in LREE, similar to modern arcs. Average ε_{Nd (T)} values are close to depleted mantle growth curves. In MORB-normalised multi-element diagrams, the komatiitic to basaltic rocks exhibit flat patterns from Lu through La and consistent relative enrichment in the elements Pb, U, Th, Ba and Cs. Apart from the Komati Complex, the majority of the lavas show significant negative Nb and Ta anomalies. Enrichment in non-conservative incompatible elements (Cs, Ba, Th, LREE) relative to conservative elements (Ta, Nb, Zr, Hf, Ti, Y, HREE) shows that the komatiitic to basaltic magmas were generated from metasomatised mantle above subducting altered oceanic crust. The geochemistry of the felsic rocks indicates an origin by melting of subducted amphibolite and eclogite. The tectonostratigraphy and the geochemical characteristics of the lavas and intrusions are consistent with successive obduction and accretion of segments of oceanic crust formed in back-arc basins and volcanic arcs.     

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.     

Thresholds of intrabed flow and other interactions of turbidity currents with soft muddy substrates

Controlled laboratory experiments reveal that the lower part of turbidity currents has the ability to enter fluid mud substrates, if the bed shear stress is higher than the yield stress of the fluid mud and the density of the turbidity current is higher than the density of the substrate. Upon entering the substrate, the turbidity current either induces mixing between flow‐derived sediment and substrate sediment, or it forms a stable horizontal flow front inside the fluid mud. Such ‘intrabed’ flow is surrounded by plastically deformed mud; otherwise it resembles the front of a ‘bottom‐hugging’ turbidity current. The ‘suprabed’ portion of the turbidity current, i.e. the upper part of the flow that does not enter the substrate, is typically separated from the intrabed flow by a long horizontal layer of mud which originates from the mud that is swept over the top of the intrabed flow and then incorporated into the flow. The intrabed flow and the mixing mechanism are specific types of interaction between turbidity currents and muddy substrates that are part of a larger group of interactions, which also include bypass, deposition, erosion and soft sediment deformation. A classification scheme for these types of interactions is proposed, based on an excess bed shear stress parameter, which includes the difference in the bed shear stress imposed by the flow and the yield stress of the substrate and an excess density parameter, which relies on the density difference between the flow and the substrate. Based on this classification scheme, as well as on the sedimentological properties of the laboratory deposits, an existing facies model for intrabed turbidites is extended to the other types of interaction involving soft muddy substrates. The physical threshold of flow‐substrate mixing versus stable intrabed flow is defined using the gradient Richardson number, and this method is validated successfully with the laboratory data. The gradient Richardson number is also used to verify that stable intrabed flow is possible in natural turbidity currents, and to determine under which conditions intrabed flow is likely to be unstable. It appears that intrabed flow is likely only in natural turbidity currents with flow velocities well below ca 3·5 m s^?1, although a wider range of flows is capable of entering fluid muds. Below this threshold velocity, intrabed flow is stable only at high‐density gradients and low‐velocity gradients across the upper boundary of the turbidity current. Finally, the gradient Richardson number is used as a scaling parameter to set the flow velocity limits of a natural turbidity current that formed an inferred intrabed turbidite in the deep‐marine Aberystwyth Grits Group, West Wales, United Kingdom.