The Okavango giant mafic dyke swarm (NE Botswana): its structural significance within the Karoo Large Igneous Province

The structural organization of a giant mafic dyke swarm, the Okavango complex, in the northern Karoo Large Igneous Province (LIP) of NE Botswana is detailed. This N110°E-oriented dyke swarm extends for 1500 km with a maximum width of 100 km through Archaean basement terranes and Permo-Jurassic sedimentary sequences. The cornerstone of the study is the quantitative analysis of N>170 (exposed) and N>420 (detected by ground magnetics) dykes evidenced on a ca. 80-km-long section lying in crystalline host-rocks, at high-angle to the densest zone of the swarm (Shashe area). Individual dykes are generally sub-vertical and parallel to the entire swarm. Statistical analysis of width data indicates anomalous dyke frequency (few data <5.0 m) and mean dyke thickness (high value of 17 m) with respect to values classically obtained from other giant swarms. Variations of mean dyke thicknesses from 17 (N110°E swarm) to 27 m (adjoining and coeval N70°E giant swarm) are assigned to the conditions hosting fracture networks dilated as either shear or pure extensional structures, respectively, in response to an inferred NNW–SSE extension. Both fracture patterns are regarded as inherited brittle basement fabrics associated with a previous (Proterozoic) dyking event. The Okavango N110°E dyke swarm is thus a polyphase intrusive system in which total dilation caused by Karoo dykes (estimated frequency of 87%) is 12.2% (6315 m of cumulative dyke width) throughout the 52-km-long projected Shashe section. Assuming that Karoo mafic dyke swarms in NE Botswana follow inherited Proterozoic fractures, as similarly applied for most of the nearly synchronous giant dyke complexes converging towards the Nuanetsi area, leads us to consider that the resulting triple junction-like dyke/fracture pattern is not a definitive proof for a deep mantle plume in the Karoo LIP.     

Breach erosion of earthfill dams (BEED) model

A computer model has been developed for simulation of breach erosion of earthfill dams (BEED). The model incorporates the processes of surface erosion and slope sloughing to simulate breach enlargement. Depletion of reservoir water is approximated by a volume continuity equation while broad-crested weir hydraulics is utilized to describe flow over and through the breach. Due to the implicit form of these equations, an iterative solution is proposed with convergence achieved within only a few iterations. The BEED model is written in both FORTRAN 77 and BASIC computer languages. Testing of the model using historical data of the failures of Teton and Huaccoto dams showed that timing, shape, and magnitude of the predicted outflow hydrograph were adequately simulated by this model. The same is true for the dimensions of the terminal breach. A sensitivity analysis indicated that internal friction angle and the relation for surface erosion were the major factors affecting the model results.     

Study of the air-sea interactions at the mesoscale: the SEMAPHORE experiment

The SEMAPHORE (Structure des Echanges Mer-Atmosphère, Propriétés des Hétérogénéités Océaniques: Recherche Expérimentale) experiment has been conducted from June to November 1993 in the Northeast Atlantic between the Azores and Madeira. It was centered on the study of the mesoscale ocean circulation and air-sea interactions. The experimental investigation was achieved at the mesoscale using moorings, floats, and ship hydrological survey, and at a smaller scale by one dedicated ship, two instrumented aircraft, and surface drifting buoys, for one and a half month in October-November (IOP: intense observing period). Observations from meteorological operational satellites as well as spaceborne microwave sensors were used in complement. The main studies undertaken concern the mesoscale ocean, the upper ocean, the atmospheric boundary layer, and the sea surface, and first results are presented for the various topics. From data analysis and model simulations, the main characteristics of the ocean circulation were deduced, showing the close relationship between the Azores front meander and the occurrence of Mediterranean water lenses (meddies), and the shift between the Azores current frontal signature at the surface and within the thermocline. Using drifting buoys and ship data in the upper ocean, the gap between the scales of the atmospheric forcing and the oceanic variability was made evident. A 2 °C decrease and a 40-m deepening of the mixed layer were measured within the IOP, associated with a heating loss of about 100 W m⁻². This evolution was shown to be strongly connected to the occurrence of storms at the beginning and the end of October. Above the surface, turbulent measurements from ship and aircraft were analyzed across the surface thermal front, showing a 30% difference in heat fluxes between both sides during a 4-day period, and the respective contributions of the wind and the surface temperature were evaluated. The classical momentum flux bulk parameterization was found to fail in low wind and unstable conditions. Finally, the sea surface was investigated using airborne and satellite radars and wave buoys. A wave model, operationally used, was found to get better results compared with radar and wave-buoy measurements, when initialized using an improved wind field, obtained by assimilating satellite and buoy wind data in a meteorological model. A detailed analysis of a 2-day period showed that the swell component, propagating from a far source area, is underestimated in the wave model. A data base has been created, containing all experimental measurements. It will allow us to pursue the interpretation of observations and to test model simulations in the ocean, at the surface and in the atmospheric boundary layer, and to investigate the ocean-atmosphere coupling at the local and mesoscales.     

Distinct brief major events in the Karoo large igneous province clarified by new Ar/Ar ages on the Lesotho basalts

Recent mineral separate ages obtained on the Karoo large igneous province (southern Africa) suggest that the province was built by several distinct magmatic pulses over a rather long period on the order of 5–6 Ma concerning the main erupted volume [Jourdan, F., Féraud, G., Bertrand, H., Kampunzu, A.B., Tshoso, G., Watkeys, M.K., Le Gall., B., 2005. The Karoo large igneous province: Brevity, origin, and relation with mass extinction questioned by new ⁴⁰Ar/³⁹Ar age data, Geology 33, 745–748]. Although this apparently atypical province is dated in more detail compared to many other large igneous provinces, volumetrically important areas still lack sufficient high-quality data. The timing of the Karoo province is crucial as this event is correlated with the breakup activity of the Gondwana supercontinent. The Lesotho basalts represent a major lava sequence of the province, but have not yet been precisely dated by systematic analysis of mineral separates. We analyzed plagioclase separates from five lava flows encompassing the complete 1.4-km-thick Lesotho sequence from top to bottom using the ⁴⁰Ar/³⁹Ar method. We obtained five plateau and mini-plateau ages statistically indistinguishable and ranging from 182.3 ± 1.6 to 181.0 ± 2.0 Ma (2σ). We derived an apparent maximum duration for this event of  0.8 Ma by neglecting correlated errors embedded in the age uncertainties.

A critical review of previous ages obtained on the Lesotho sequence [Duncan R.A., Hooper, P.R., Rehacek, J., Marsh, J.S., Duncan, A.R., 1997. The timing and duration of the Karoo igneous event, southern Gondwana. Journal of Geophysical Research 102, 18127–18138] shows that groundmass analyses are unreliable for high-resolution geochronology, due to alteration and ³⁹Ar recoil effects. Discrepancy between our ages and a previous plagioclase age at  184 Ma obtained by the later workers is tentatively attributed to the heterogeneity of the monitor used and/or cryptic excess ⁴⁰Ar. The current age database suggests that at least three temporally and spatially distinct brief major events (the Lesotho and southern Botswana lava piles and the Okavango dyke swarm) are so far recognized in the Karoo province. Identification of brief and volumetrically important Karoo magmatic events allows detecting the migration of the Karoo magmatism and potentially the stress regime that affected the southern African lithosphere at this time. A filtered compilation of 60 ages obtained with homogeneous intercalibrated standards suggests a shorter duration for the main pulses of the magmatism between 3 and 4.5 Ma, compared to a whole province duration of  10 Ma, between  182 and  172 Ma.     

Isotopic fractionation of Cu in tektites

Tektites are terrestrial natural glasses of up to a few centimeters in size that were produced during hypervelocity impacts on the Earth’s surface. It is well established that the chemical and isotopic composition of tektites is generally identical to that of the upper terrestrial continental crust. Tektites typically have very low water content, which has generally been explained by volatilization at high temperature; however, the exact mechanism is still debated. Because volatilization can fractionate isotopes, comparing the isotopic composition of volatile elements in tektites with those of their source rocks may help to understand the physical conditions during tektite formation.Interestingly, volatile chalcophile elements (e.g., Cd and Zn) seem to be the only elements for which isotopic fractionation is known so far in tektites. Here, we extend this study to Cu, another volatile chalcophile element. We have measured the Cu isotopic composition for 20 tektite samples from the four known different strewn fields. All of the tektites (except the Muong Nong-types) are enriched in the heavy isotopes of Cu (1.98 < δ⁶⁵Cu < 6.99) in comparison to the terrestrial crust (δ⁶⁵Cu ≈ 0) with no clear distinction between the different groups. The Muong Nong-type tektites and a Libyan Desert Glass sample are not fractionated (δ⁶⁵Cu ≈ 0) in comparison to the terrestrial crust. To refine the Cu isotopic composition of the terrestrial crust, we also present data for three geological reference materials (δ⁶⁵Cu ≈ 0).An increase of δ⁶⁵Cu with decreasing Cu abundance probably reflects that the isotopic fractionation occurred by evaporation during heating. A simple Rayleigh distillation cannot explain the Cu isotopic data and we suggest that the isotopic fractionation is governed by a diffusion-limited regime. Copper is isotopically more fractionated than the more volatile element Zn (δ^66/64Zn up to 2.49‰). This difference of behavior between Cu and Zn is predicted in a diffusion-limited regime, where the magnitude of the isotopic fractionation is regulated by the competition between the evaporative flux and the diffusive flux at the diffusion boundary layer. Due to the difference of ionic charge in silicates (Zn²⁺ vs. Cu⁺), Cu has a diffusion coefficient that is larger than that of Zn by at least two orders of magnitude. Therefore, the larger isotopic fractionation in Cu than in Zn in tektites is due to the significant difference in their respective chemical diffusivity.     

40Ar/39Ar and (U‐Th)/He model age signatures of elusive Mercurian and Venusian meteorites

No meteorites from Mercury and Venus have been conclusively identified so far. In this study, we develop an original approach based on extensive Monte Carlo simulations and diffusion models to explore the radiogenic argon (⁴⁰Ar*) and helium (⁴He*) loss behavior and the range of ⁴⁰Ar/³⁹Ar and (U‐Th)/He age signatures expected for a range of crystals if meteorites from these planets were ever to be found. We show that we can accurately date the crystallization age of a meteorite from both Mercury and Venus using the ⁴⁰Ar/³⁹Ar technique on clinopyroxene (± orthopyroxene) and that its ⁴⁰Ar/³⁹Ar age should match the Pb‐Pb age. At the surface of Mercury, phases like albite and anorthite will exhibit a complete range of ⁴⁰Ar* loss ranging from 0% to 100%, whereas merrillite and apatite will show 100% ⁴He* loss. By measuring the crystal size and diffusion parameters of a series of plagioclase crystals, one can inverse the ⁴⁰Ar* loss value to estimate the maximum temperature experienced by a rock, and narrow down the possible pre‐ejection location of the meteorite at the surface of Mercury. At the surface of Venus, plagioclase and phosphate phases will only record the age of ejection. The (U‐Th)/He systematics of merrillite and apatite will be, respectively, moderately and strongly affected by ⁴He* loss during the transit of the meteorite from its host planet to Earth. Finally, meteorites from Mercury or Venus will each have their own ⁴⁰Ar/³⁹Ar and (U‐Th)/He isotopic age and ³⁸Ar_c cosmic ray exposure age signatures over a series of different crystal types, allowing to unambiguously recognize a meteorite for any of these two planets using radiogenic and cosmogenic noble gases.     

The problem of inherited Ar* in dating impact glass by the Ar/Ar method: Evidence from the Tswaing impact crater (South Africa)

The Tswaing meteorite impact crater is a 1.13 km diameter structure located in the 2.05 Ga Nebo granite of the Bushveld Complex. The impact age had previously been determined by fission track dating to 220 ± 104 ka. ⁴⁰Ar/³⁹Ar step-heating and total fusion experiments performed on single- and multi-grain impact glass aliquots gave apparent ages ranging from 1.0 ± 0.3 Ma to 204 ± 6 Ma. These “ages” indicate that the radiogenic Ar derived from the target rocks has not been completely degassed as a result of the impact process, despite fusion of the target material. Results of step-heating experiments imply that the trapped within the glass is located in two distinct reservoirs thought to be the glass matrix and fluid/vapor inclusions (or un-melted residual clasts). Calculations assuming an age of 0.2 ± 0.1 Ma for Tswaing (fission track data) reveal that the amount of inherited ⁴⁰Ar*() relative to the pre-impact concentration varies from 0.015% to 4.15%. The spread defined by likely reflects the various quench rates experienced by the glass, most certainly due to the pre-impact position of the sample relative to the center of the crater. We compare the influence of on the apparent ⁴⁰Ar/³⁹Ar age determination of five impact structures. Our calculations show that the main characteristic controlling the age offset (for a given proportion of ) is the age difference between the impact and the target rocks (i.e., the ⁴⁰Ar* concentration in the target rock). The buffer effect for a given crater structure can be predicted knowing the age of the basement and having a rough estimation of the age of the crater structure itself. The occurrence of is likely influenced by (1) the degree of polymerization (i.e., silicate structure complexity) of the target rock and presumably related to the diffusivity of Ar in the melt and glass, (2) the Ar partial pressure at the grain boundary, (3) the quantity of energy involved in the impact, and (4) the porosity of the target rocks. For glass that inevitably suffers inherited and/or excess ⁴⁰Ar*, the use of the inverse isochron technique can be appropriate but should be applied with careful statistical treatment.