Evolution of the Kivu Rift,East Africa: interplay among tectonics,sedimentation and magmatism

The seismically and volcanically active Kivu Rift, in the western branch of the East African Rift System, is a type locale for studies of high‐elevation, humid‐climate rift basins, as well as magmatic basin development. Interpretations of offshore multi‐channel seismic (MCS) reflection data, terrestrial radar topography, lake bathymetry and seismicity data recorded on a temporary array provide new insights into the structure, stratigraphy and evolution of the Kivu rift. The Kivu rift is an asymmetric graben controlled on its west side by a ca. 110 km‐long, N‐S striking border fault. The southern basins of the lake and the upper Rusizi river basin are an accommodation zone effectively linking 1470 m‐high Lake Kivu to 770 m‐high Lake Tanganyika. MCS data in the eastern Kivu lake basin reveal a west‐dipping half graben with at least 1.5 km of sedimentary section; most of the ca. 2 km of extension in this sub‐basin is accommodated by the east‐dipping Iwawa normal fault, which bounds an intrabasinal horst. Lake Kivu experienced at least three periods of near desiccation. The two most recent of these approximately correlate to the African Megadrought and Last Glacial Maximum. There was a rapid lake level transgression of at least 400 m in the early Holocene. The line load of the Virunga volcanic chain enhances the fault‐controlled basin subsidence; simple elastic plate models suggest that the line load of the Virunga volcanic chain depresses the basin by more than 1 km, reduces flank uplift locally and broadens the depocentre. Not only do the voluminous magmatism and degassing to the lake pose a hazard to the riparian population, but our studies demonstrate that magmatism has important implications for short‐term processes such as lake levels, inflow and outlets, as well as long term modification of classic half‐graben basin morphology.     

The Pan-African continental margin in northeastern Africa: evidence from a geochronological study of granulites at Sabaloka, Sudan

Ion microprobe zircon ages, a Nd model age and RbSr whole-rock dates are reported from the high-grade gneiss terrain at Sabaloka on the River Nile north of Khartoum, formally considered to be part of the Archaean/early Proterozoic Nile craton. The granulites, which are of both sedimentary and igneous derivation, occur as remnants in migmatites. Detrital zircon ages range from ≈ 1000 to ≈ 2650 Ma and prove the existence of Archaean to late Proterozoic continental crust in the sedimentary source region. The Nd model age for one sedimentary granulite is between 1.26 (T_CHUR) and 1.70 (T_DM) Ga and provides a mean crustal residence age for the sedimentary precursor. Igneous zircons in enderbitic gneiss crystallized at 719 ± 81 Ma ago, an age that also corresponds to severe Pb loss in the detrital zircons and which probably reflects the granulite event at Sabaloka. The RbSr data indicate isotopic homogenization at about 700 Ma ago in the granulites and severe post-granulite disturbance at ≈ 570 Ma in the migmatites. We associate this disturbance with hydration, retrograde metamorphism and anatexis that produced undeformed granites ≈ 540 Ma ago. The ≈ 700 Ma granulite event at Sabaloka suggests that this part of the Sudan belongs to the Pan-African Mozambique belt while the ancient Nile craton lay farther west. The gneisses studied here may represent the infrastructure of the ancient African continental margin onto which the juvenile arc assemblage of the Arabian-Nubian shield was accreted during intense horizontal shortening and crustal interstacking of a major collision event.     

Characteristics of the seismic attenuation in two tectonically active zones of Southern Europe

The seismic energy attenuation in the frequency range of 1–18 Hz was studied in the two tectonically active zones of Irno Valley (Southern Italy) and Granada Basin (South-East Spain). Data were recorded by short period vertical components seismographs for low-magnitude local earthquakes. The method of coda waves, assuming singleS toS scattering approximation, was used to calculate the quality factorQ from the two data set. Results show a quality factor increasing with frequency, following the empirical lawQ=Q _o f ⁿ .Q _o andn are lower for the Irno Valley than for Granada. This result is interpreted in terms of different scattering environments present in the two investigated areas.     

Iron porphyrins in the Permian Kupferschiefer of the lower Rhine Basin, N.W. Germany

Immature samples of the Permian Kupferschiefer from the Lower Rhine Basin in N.W. Germany were analysed for tetrapyrrole pigment type and abundance. The sediment, thought to have been deposited in a marine regime with enhanced salinity, was found to contain high concentrations of metalloporphyrins. The porphyrins are complexed to nickel (Ni) and oxovanadium (V=0), but high abundances of iron (Fe) porphyrins were also detected using UV/visible spectroscopy and mass spectrometry. The presence in the latter of series of aetioporphyrins, cycloalkanoporphyrins, di-cycloalkanoporphyrins and benz-cycloalkanoporphyrins was confirmed by accurate mass measurements; HPLC co-injection of deoxophylloerythroetioporphyrin (C₃₂ DPEP) with the demetallated iron porphyrins indicated its presence in the sediment as an iron complex. The study provides the first evidence for the occurrence of Fe porphyrins in geological samples other than coals and lignites, and reports the highest concentrations in sedimentary organic matter to date.     

Assessing possible changes in flood frequency due to climate change in mid‐sized watersheds in New York State,USA

The frequency of flooding is often presumed to increase with climate change because of projected increases in rainfall intensities. In this paper, using 50‐plus years of historical discharge and meteorological data from three watersheds in different physiographic regions of New York State, USA, we find that annual maximum stream discharges are associated with 20% or less of the annual maximum rainfall events. Instead of rainfall events, approximately 20% of annual maximum stream discharges are associated with annual maximum snowmelt events while 60% of annual maximum discharges are associated with moderate rainfall amounts and very wet soil conditions. To explore the potential for changes in future flood risk, we employed a compound frequency distribution that assumes annual maximum discharges can be modelled by combining the cumulative distribution functions of discharges resulting from annual maximum rainfall, annual maximum snowmelt, and occurrences of moderate rain on wet soils. Basing on a compound frequency distribution comprised of univariate general extreme value (GEV) and gamma distributions, we found that a hypothetical 20% increase in the magnitude of rainfall‐related stream discharge results in little change in 96th percentile annual maximum discharge. For the 99th percentile discharge, two waterbodies in our study had a 10% or less increase in annual maximum discharge when annual maximum rainfall‐related discharges increased 20% while the third waterbody had a 16% increase in annual maximum discharges. Additionally, in some cases, annual maximum discharges could be offset by a reduction in the discharge resulting from annual maximum snowmelt events. While only intended as a heuristic tool to explore the interaction among different flood‐causing mechanisms, use of a compound flood frequency distribution suggests a case can be made that not all waterbodies in humid, cold regions will see extensive changes in flooding due to increased rainfall intensities. Copyright © 2011 John Wiley & Sons, Ltd.     

ENSO,the IOD and the intraseasonal prediction of heat extremes across Australia using POAMA-2

The simulation and prediction of extreme heat over Australia on intraseasonal timescales in association with the El Niño–Southern Oscillation (ENSO) and the Indian Ocean Dipole (IOD) is assessed using the Bureau of Meteorology’s Predictive Ocean Atmosphere Model for Australia (POAMA). The analysis is based on hindcasts over 1981–2010 and focuses on weeks 2 and 3 of the forecasts, i.e. beyond a typical weather forecast. POAMA simulates the observed increased probabilities of extreme heat during El Niño events, focussed over south eastern and southern Australia in SON and over northern Australia in DJF, and the decreased probabilities of extreme heat during La Niña events, although the magnitude of these relationships is smaller than observed. POAMA also captures the signal of increased probabilities of extreme heat during positive phases of the IOD across southern Australia in SON and over Western Australia in JJA, but again underestimates the strength of the relationship. Shortcomings in the simulation of extreme heat in association with ENSO and the IOD over southern Australia may be linked to deficiencies in the teleconnection with Indian Ocean SSTs. Forecast skill for intraseasonal episodes of extreme heat is assessed using the Symmetric Extremal Dependence Index. Skill is highest over northern Australia in MAM and JJA and over south-eastern and eastern Australia in JJA and SON, whereas skill is generally poor over south-west Western Australia. Results show there are windows of forecast opportunity related to the state of ENSO and the IOD, where the skill in predicting extreme temperatures over certain regions is increased.     

Testing and modeling of the mechanical behavior of dolomite in the Wudongde hydropower plant

In order to understand the mechanical behaviours of the surrounding rocks in the underground caverns of the Wudongde hydropower plant, triaxial tests are performed on a type of dolomite. It is revealed that damage induced by crack development is the main factor controlling the nonlinear plastic deformation and failure behaviour of the dolomite in both pre- and post-peak regimes. Based on this understanding, a coupled elastoplastic damage model is developed for capturing the dolomite’s mechanical behaviours. In the model, the effects of plasticity and damage on rocks is described by introducing plastic hardening and damage softening commonly in the plastic yield surface. Which are both derived from a suitable Helmholtz free energy function. The model is used to simulate the triaxial tests. Comparisons between test results and the numerical modelling show that the developed model is capable of describing the macro mechanical behaviours of the Wudongde dolomite.     

Molecular dynamics simulations of pressure and temperature effects on MgSiO3 and Mg2SiO4 melts and glasses

Ab-initio interionic potentials for Mg²⁺, Si⁴⁺, and O^2– have been used in molecular dynamics (MD) simulations to investigate diffusivity changes, pressure-induced structural transitions, and temperature effects on polymerization in MgSiO₃ and Mg₂SiO₄ melts and glasses. The potential gives reasonable agreement with the 0.1 MPa radial distribution function of MgSiO₃ glass. Maxima in the diffusion coefficients of Si⁴⁺ and O^2– occur as pressure is increased on the MgSiO₃ melt. The controlling structural mechanism for this behavior is the Q¹ species of SiO₄ tetrahedra. Mg²⁺ diffusion coefficients decrease monotonically with pressure in both melt compositions. Increasing Mg²⁺ coordination number and population of 3- and 4-membered SiO₄ rings with pressure combine to hinder translation of the Mg²⁺ ions. The dominant changes in structure with pressure are a decrease in the intertetrahedral (Si-O--Si) angle up to approximately 4 g/cm³ and coordination changes of the ions above this density. Temperature effects on viscosity in these simulated melts are indirectly studied by analyzing polymerization changes with temperature. Polymerization and coordination numbers increase with decreasing temperature and a small quench rate effect is observed. Fair agreement is found between the MD simulations and experimental equation of state for Mg₂SiO₄, but the equation of state predictions for MgSiO₃ melts are much less accurate. The zero pressure volume, V ₀, is significantly higher and K ₀ is lower in the simulations than empirical values. The inadequacies reflect error in using the ionic approximation for polymerized systems and a need to collect more data for a variety of molecular configurations in the development of ab-initio potentials.