A detailed analysis of the February 1996 aftershock sequence in the eastern Pyrenees, France

Following the 1996 February 18 M _L = 5.2 earthquake in the Agly massif in the eastern French Pyrenees, we installed a temporary network of seismometers around the epicentre. In this paper, we analyse 336 well-located aftershocks recorded from February 19 to February 23 by 18 temporary stations and two permanent stations located less than 35  km from the epicentre. Most aftershocks have been located with an accuracy better than 1.5  km in both horizontal and vertical positions. Their spatial distribution suggests the reactivation of a known fault system. We determined 39 fault-plane solutions using P -wave first motions. Despite their diversity, the focal mechanisms yield an E–W subhorizontal T-axis. We also determined fault-plane solutions and principal stress axes using the method developed by Rivera & Cisternas (1990 ) for the 15 best-recorded events. We obtain a pure-shear-rupture tectonic regime under N–S subhorizontal compression and E–W subhorizontal extension. These principal stress axes, which explain the focal mechanisms for at least 75 per cent of the 39 aftershocks, are different from the axes deduced from the main shock. The post-earthquake stress field caused by the main-shock rupture, modelled as sinistral strike slip on three vertical fault segments, is computed for various orientations and magnitudes of the regional stress field, assumed to be horizontal. The aftershock distribution is best explained for a compressive stress field oriented N30°E. Most aftershocks concentrate where the Coulomb failure stress change increases by more than 0.2  MPa. The diversity of aftershock focal mechanisms, poorly explained by this model, may reflect the great diversity in the orientations of pre-existing fractures in the Agly massif.     

Moment tensors of microearthquakes from the Eyjafjallajökull volcano in South Iceland

Analyses of relative P - and S -wave amplitudes of 15 selected earthquakes ( M _L <2.3) from a seismic swarm, which occurred in May and June 1994 at the Eyjafjallajökull volcano in South Iceland, reveal similar radiation patterns, a thrust-type double-couple with an additional source component. All focal solutions have nearly vertical T -axes and horizontally oriented P -axes, with E-W-oriented nodal planes. The volume increase corresponding to an isotropic source component is estimated to be in the range of 24 m³. The temporal and spatial seismic pattern, small magnitude range, focal mechanisms and depth range of the Eyjafjallaökull earthquakes indicate vertical intrusion of magma into a confined region at the northern flank of the volcano.     

Palaeointensity results from Ethiopian basalts: implications for the Oligocene geomagnetic field strength

From a large collection of Ethiopian flood basalts (~30  Myr old) sampled for magnetostratigraphy, ⁴⁰Ar/³⁹Ar geochronology and geochemical investigations, 47 samples were selected in order to test their suitability for Thellier palaeointensity experiments. Only 17 samples from eight individual flows yielded reliable palaeointensity estimates, with flow-mean virtual dipole moments ranging from 3.0 to 10.5 × 10²²  A  m² .
  A critical review of the Oligocene palaeointensity data set, including these new Ethiopian data, indicates an Oligocene mean virtual dipole moment of 5.1 ± 2.5 × 10²²  A  m² for the complete data set. After applying mild selection criteria, the reduced data set yields a mean value of only 4.6 ± 1.9 × 10²²  A  m² . This value is significantly lower than the present-day field strength but is higher than the Mesozoic dipole low mean field. This low Oligocene field might be in agreement with the high palaeosecular variation and rather high non-dipole field invoked around 30  Ma. However, the Oligocene data set is largely dependent on the palaeointensity determinations from Armenia, obtained mainly from baked contacts, which show amazingly low dispersion at both flow and between-flow levels. More data are needed to reduce the weight of these determinations on the mean value and avoid a possible bias.     

Temporal changes in shear-wave splitting at an isolated swarm of small earthquakes in 1992 near Dongfang, Hainan Island, southern China

An isolated swarm of small earthquakes occurred in 1992, near Dongfang on Hainan Island, southern China. The Institute of Geophysics, State Seismological Bureau of China, monitored the swarm with five DCS-302 digital accelerometers for three months from 1992 June 1. 18 earthquakes, with magnitudes M _L ranging from 1.8 to 3.6, were well located by five stations, and shear-wave splitting varying azimuthally was analysed on 27 seismic records from these events. The mean polarization azimuth of the faster shear wave was WNW. Time delays between the split shear waves at two stations varied with time and space. The time delays at one station fell abruptly after earthquakes of magnitudes 3.1 and 3.6, but did not change significantly at the second station. This behaviour is consistent with the delay-time changes being caused by changes in the aspect ratio of vertical liquid-filled (EDA) cracks. Thus, the variation in shear-wave-splitting time delay could be due to changes in crustal stress related to nearby small-magnitude earthquake activity. The connection between earthquake activity and crustal stress variation measured by shear-wave splitting leaves the door open for possible observations of crustal stress transients related to the onset of an earthquake; however, our data cannot be considered as definite evidence for such precursors.     

Subsurface mass redistributions at Mount Etna (Italy) during the 1995–1996 explosive activity detected by microgravity studies

Gravity changes are presented from a series of field microgravity surveys conducted at Mt Etna between August 1994 and November 1996, a period including the 1995–1996 explosive summit activity. Data were collected along a microgravity network of 69 stations at a monthly to annual sampling rate, depending on each subarray of the network.
  Results show that seasonal changes in water level within the volcano may induce gravity changes of up to 20  μgal on Etna's southern slope, and indicate that significant magma movement occurred within and below Etna's edifice between 1994 and 1996. In particular, between September 1994 and October 1995, a mass increase of 2 × 10¹⁰  kg occurred 2000  m beneath the summit craters. Between October 1995 and July 1996 this mass was lost, while another 2 × 10¹⁰  kg was injected at about 1000  m  a.s.l. into the 1989 fracture system. From the gravity data alone, it is not possible to distinguish whether the first shallow intrusion (1994–1995) was then injected laterally into the 1989 fracture, or summit activity was fed by the first shallow intrusion, while new magma entered the 1989 fracture system.
  While magma was being redistributed within the volcanic edifice, measurements along an E–W-trending profile on the southern slope of the volcano detected some 1.5 × 10¹¹  kg of magma accumulating 2–3  km below sea level between October 1995 and November 1996.     

Relocation of M ≥ 2 events of the 1989 Dôbi seismic sequence in Afar: evidence for earthquake migration

The Ethiopian side of central Afar was struck in August 1989 by the largest seismic sequence (three 6.1 ≤ M _s ≤ 6.3 events, 15 with M _s or m _b ≥ 5.0) since that of Serdo in 1969. Using the Djibouti seismological network, we relocated 297 of the events of that sequence. As most of the large events took place outside the network, we assessed the accuracy and stability of earthquake relocations by using three different velocity models and two relocation codes to try to relate individual shocks to distinct faults and surface breaks. A majority of the events apparently occurred underneath the floor of the Dôbi graben, an area about 45  km long and 15  km wide, rupturing boundary and inner floor faults, in agreement with the surface cracks and scarps that we mapped in the area. The relocation shows that the principal events propagated about 50  km northwestwards along the graben in the first 40  hr. A day and a half after the beginning of the sequence, smaller events ( M ≤ 4) started to propagate more than 55  km eastwards, towards Asal Lake. Using two three-component stations installed near the Ethiopian border, we could determine reliable depths for 21 events. The depths are compatible with a seismogenic crust about 14  km thick in the Dôbi and Hanle graben area. Although the Dôbi sequence ruptured about 50  km of the fault array extending from Serdo to Asal, where the regional stress was released by earthquakes in 1969 and 1978, respectively, a seismic gap about 50  km long still subsists along the northern part of the Gaggade region (Der'êla half-graben).     

Seismic quiescence before the M 7, 1988, Spitak earthquake, Armenia

A detailed analysis of the 35  yr of seismicity between 1962 and 1997 using a gridding technique shows that the M 7, Spitak earthquake of 1988 December 7 was preceded by a quiescence anomaly that started at approximately 1984±0.5, and lasted about 5±0.5  yr, up to the main shock. This quiescence anomaly had a radius of about 20±3  km, estimated from circular areas with 75 per cent rate decrease, centred at the point of maximum significance of the anomaly. The quiescence was clearly present in the aftershock volume during the 5  yr before the 1988 main shock, but its statistically strongest expression was located 30  km NW of the epicentre. This anomaly fulfills the association rules between precursory quiescence anomalies and main shocks, even for a tight definition, and is therefore proposed as a case of precursory quiescence. The largest value of the standard deviate Z , found by random selection of samples by gridding, was Z =14 for a time window of T _w=3  yr, using a sample size of N =300 events. This makes this anomaly the strongest observed so far, and it is the first documented in an environment of continental collision. There are no false alarms exceeding in significance the precursor. The Armenian earthquake catalogue used for this study had 4600 earthquakes with M ≥ M _min=2.2 in the area bounded by 39.5° to 42°N/42.5° to 47°E. From the point of view of homogeneous reporting this is the best catalogue we have analysed so far. The limits of the data used and the density of the grid are dictated by the data, and have no influence on the results. The choice of free parameters does not influence the results significantly within the following limits: 100≤ N ≤500, 2≤ T _w≤7, 2.2≤ M _min≤2.8.     

Constraints on the frequency–magnitude relation and maximum magnitudes in the UK from observed seismicity and glacio-isostatic recovery rates

Earthquake populations have recently been shown to have many similarities with critical-point phenomena, with fractal scaling of source sizes (energy or seismic moment) corresponding to the observed Gutenberg–Richter (G–R) frequency–magnitude law holding at low magnitudes. At high magnitudes, the form of the distribution depends on the seismic moment release rate M˙ and the maximum magnitude m _max . The G–R law requires a sharp truncation at an absolute maximum magnitude for finite M˙ . In contrast, the gamma distribution has an exponential tail which allows a soft or 'credible' maximum to be determined by negligible contribution to the total seismic moment release. Here we apply both distributions to seismic hazard in the mainland UK and its immediate continental shelf, constrained by a mixture of instrumental, historical and neotectonic data. Tectonic moment release rates for the seismogenic part of the lithosphere are calculated from a flexural-plate model for glacio-isostatic recovery, constrained by vertical deformation rates from tide-gauge and geomorphological data. Earthquake focal mechanisms in the UK show near-vertical strike-slip faulting, with implied directions of maximum compressive stress approximately in the NNW–SSE direction, consistent with the tectonic model. Maximum magnitudes are found to be in the range 6.3–7.5 for the G–R law, or 7.0–8.2 m _L for the gamma distribution, which compare with a maximum observed in the time period of interest of 6.1 m _L . The upper bounds are conservative estimates, based on 100 per cent seismic release of the observed vertical neotectonic deformation. Glacio-isostatic recovery is predominantly an elastic rather than a seismic process, so the true value of m _max is likely to be nearer the lower end of the quoted range.     

Thermochronometry and palaeomagnetism of the Archaean Nelshoogte Pluton, South Africa

²⁰⁷Pb/²⁰⁶Pb single-grain zircon, ⁴⁰Ar/³⁹Ar single-grain hornblende and biotite, and ⁴⁰Ar/³⁹Ar bulk-sample muscovite and biotite ages from the Nelshoogte trondhjemite pluton located in eastern Transvaal, South Africa, show that this granitoid had a protracted thermal history spanning 3213±4  Ma to about 3000  Ma. Whole-rock ⁴⁰Ar/³⁹Ar ages from cross-cutting dolerite dykes indicate that these were intruded at about 1900  Ma. There is no evidence of this or other, later events significantly affecting the argon systematics of the minerals from the pluton dated by the ⁴⁰Ar/³⁹Ar method.
  The pluton has a well-defined palaeomagnetic pole which is dated at 3179±18 (2 σ ) Ma by ⁴⁰Ar/³⁹Ar dating of hornblende. This pole (18°N, 310°E, A ₉₅=9°) yields a palaeolatitude of 0°, significantly different from other Archaean poles from the Kaapvaal Craton. The palaeolatitude difference implies that there was significant apparent polar wander during the Archaean. A second, overprinting magnetization seen in the pluton is also seen in the lower-Proterozoic dolerite dykes, and is consistent with other lower-Proterozoic (2150–1950  Ma) poles for southern Africa.     

The Gulf of Aqaba earthquake swarm of 1983 January-April

Summary. In the period 1983 January 21 -April 20, more than 500 local earthquakes ( M _L≤ 4.85) occurred in the Gulf of Aqaba area between latitudes 29°00'and 29°25'and longitudes 34°30'and 34°45'. Most of the activity including the largest shocks was restricted to the area between latitudes 29°07'and 29° 15'and longitudes 34°33'and 34°42'where the NW Atiya regional dyke crosses the area and is horizontally displaced by NE strike-slip faults. The first-motion directions of four large shocks, including the largest, at both UNJ and HLW stations are in agreement with a strike-slip mechanism at a NE-trending fault in this area. The b value showed a temporal increase with time from 0.43 to 0.69. This, together with other geological and geophysical observations may indicate that subsurface magmatic activity has affected the stressed crustal rocks, thus triggering earthquake activity.
This swarm and historical information indicate that the Gulf of Aqaba-Dead Sea Jordan transform is characterized by both swarm and foreshock-aftershock types of seismic activity and therefore the relatively large proportion of non-seismic slip along the southern part of this transform may actually be higher if swarm-type activities are considered.     

Aseismic fault movement before the 1995 Kobe earthquake detected by a GPS survey: implication for preseismic stress localization?

By inversion analysis of the baseline changes and horizontal displacements observed with GPS (Global Positioning System) during 1990–1994, a high-angle reverse fault was detected in the Shikoku-Kinki region, southwest Japan. The active blind fault is characterized by reverse dip-slip (0.7±0.2  m yr⁻¹ within a layer 17–26  km deep) with a length of 208±5  km, a (down-dip) width of 9±2  km, a dip-angle of 51°±2° and a strike direction of 40°±2° (NE). Evidence from the geological investigation of subfaults close to the southwestern portion of the fault, two historical earthquakes ( M _L=7.0, 1789 and 6.4, 1955) near the centre of the fault, and an additional inversion analysis of the baseline changes recorded by the nationwide permanent GPS array from 18 January to 31 December 1995 partially demonstrates the existence of the fault, and suggests that it might be a reactivation of a pre-existing fault in this region. The fact that hardly any earthquakes ( M _L>2.0) occurred at depth on the inferred fault plane suggests that the fault activity was largely aseismic. Based on the parameters of the blind fault estimated in this study, we evaluated stress changes in this region. It is found that shear stress concentrated and increased by up to 2.1 bar yr⁻¹ at a depth of about 20  km around the epicentral area of the 1995 January 17  Kobe earthquake ( M _L=7.2, Japan), and that the earthquake hypocentre received a Coulomb failure stress of about 5.6 bar yr⁻¹ during 1990–1994. The results suggest that the 1995  Kobe earthquake could have been induced or triggered by aseismic fault movement.     

Inference of mantle viscosity from GRACE and relative sea level data

Gravity Recovery And Climate Experiment (GRACE) satellite observations of secular changes in gravity near Hudson Bay, and geological measurements of relative sea level (RSL) changes over the last 10 000 yr in the same region, are used in a Monte Carlo inversion to infer-mantle viscosity structure. The GRACE secular change in gravity shows a significant positive anomaly over a broad region (>3000 km) near Hudson Bay with a maximum of ∼2.5 μGal yr⁻¹ slightly west of Hudson Bay. The pattern of this anomaly is remarkably consistent with that predicted for postglacial rebound using the ICE-5G deglaciation history, strongly suggesting a postglacial rebound origin for the gravity change. We find that the GRACE and RSL data are insensitive to mantle viscosity below 1800 km depth, a conclusion similar to that from previous studies that used only RSL data. For a mantle with homogeneous viscosity, the GRACE and RSL data require a viscosity between  1.4 × 10²¹  and  2.3 × 10²¹  Pa s. An inversion for two mantle viscosity layers separated at a depth of 670 km, shows an ensemble of viscosity structures compatible with the data. While the lowest misfit occurs for upper- and lower-mantle viscosities of  5.3 × 10²⁰  and  2.3 × 10²¹  Pa s, respectively, a weaker upper mantle may be compensated by a stronger lower mantle, such that there exist other models that also provide a reasonable fit to the data. We find that the GRACE and RSL data used in this study cannot resolve more than two layers in the upper 1800 km of the mantle.     

Permo-Carboniferous (Kiaman) palaeointensity results from the western Bohemian Massif, Germany

Thellier palaeointensity analyses were carried out on samples from three quartz porphyry intrusions located within the western part of the Bohemian Massif close to the German-Czech Republic border. These intrusives have been dated at 280-275 Ma. which places them in the main part of the Permo-Carboniferous (Kiaman) Reversed Superchron (P-CRS). This age is consistent with the mean characteristic remanence direction calculated at D/ I = 204/- 1, giving a VGP at 36N, 341E. Palaeointensity values range from 6.6 μT to 22.1 μT, with a mean value of 12.0 ± 1.3 μT, which has a corresponding VDM value of (3.0 ± 1.2) × 10²² A m². This suggests that the geomagnetic field strength recorded by the quartz porphyries is only 37 per cent of the current value, supporting recent studies that report a field strength of between 25 and 40 per cent of the present-day value during the P-CRS.     

New techniques for the analysis of earthquake sources from local array data with an application to the 1993 Scotts Mills, Oregon, aftershock sequence

We analysed aftershocks recorded by a temporary digital seismic network following the moderate M _w = 5.5 1993, Scotts Mills, Oregon, earthquake. A technique to retrieve source moment tensors from local waveforms was developed, tested, and applied to 41 small earthquakes ( M _w ranging from 1.6 to 3.2). The derived focal mechanisms, although well resolved, are highly variable and do not share a common nodal plane. In contrast, the majority of the events, relocated with a joint hypocentre determination algorithm, collapse to a well-focused plane. The incompatibility of the nodal planes of most events with the plane defined by their locations suggests that the aftershocks did not occur on the fault plane, but tightly around it, outlining the rupture area rather than defining it. Furthermore, the moment tensors reveal stable P -axes, whereas T  -axes plunges are highly dispersed. We detect a rotation of average T  -axis plunge with depth, indicating a change from shallower, predominantly dip-slip mechanisms to deeper strike-slip mechanisms. These characteristics are difficult to explain by remnant stress concentrations on the main-shock rupture plane or asperity- and barrier-type models. We suggest that the aftershocks occurred under the ambient regional stress, triggered by a sudden weakening of the region surrounding the main-shock slip, rather than from a shear stress increase due to the main shock.     

Thrust faulting and crust—upper mantle structure in East Australia

Summary. The mid-crustal earthquake of 1973 March 9 (m_b= 5.5, h ≤ 20 km) located 60 km south-west of Sydney, Australia, provides unambiguous evidence of contemporary thrust faulting in South-eastern Australia — a region of high heat flow and Cenozoic basaltic volcanism. Aftershock locations suggest a steeply dipping fault in the depth range from 8 to 24 km with a lateral extent of about 8 km. The mechanism solution is consistent with a tectonic stress field that is dominated by east—west horizontal compression. A seismic moment of 5.7 ± 10²³± 20 per cent dyne-cm was computed from surface-wave amplitudes. Minimum values of slip and stress drop, 2 cm and 1 bar respectively, were estimated from the moment and a fault size taken' from aftershock locations.
Refinement modelling by a controlled Monte Carlo technique was used to provide unbiased models directly from multimode group velocities. The dispersion of fundamental and higher mode surface waves recorded at the digital high-gain station at Charters Towers, Queensland, and the WWSSN station at Adelaide, South Australia, is satisfied by crust- and upper-mantle models which have neither pronounced S-wave low-velocity zones nor thick high-velocity lids within 140 km of the Earth's surface. These models have subcrustal shear velocities of 4.20–4.32 km/s which are 0.4–0.5 km/s slower than Canadian shield shear velocities (CANSD).     

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.     

Temporal and spatial variations of post-seismic deformation following the 1999 Chi-Chi, Taiwan earthquake

We use GPS displacements collected in the 15 months after the 1999 Chi-Chi, Taiwan earthquake  ( M _w 7.6)  to evaluate whether post-seismic deformation is better explained by afterslip or viscoelastic relaxation of the lower crust and upper mantle. We find that all viscoelastic models tested fail to fit the general features in the post-seismic GPS displacements, in contrast to the satisfactory fit obtained with afterslip models. We conclude that afterslip is the dominant mechanism in the 15-month period, and invert for the space–time distribution of afterslip, using the Extended Network Inversion Filter. Our results show high slip rates surrounding the region of greatest coseismic slip. The slip-rate distribution remains roughly stationary over the 15-month period. In contrast to the limited coseismic slip on the décollement, afterslip is prominent there. Maximum afterslip of 0.57 m occurs downdip and to the east of the hypocentral region. Afterslip at hypocentral depths is limited to the southern part of the main shock rupture, with little or no slip on the northern section where coseismic slip was greatest. Whether this results from along strike variations in frictional properties or dynamic conditions that locally favour stable sliding is not clear. In general, afterslip surrounds the area of greatest coseismic slip, consistent with post-seismic slip driven by the main shock stress change. The total accumulated geodetic afterslip moment is  3.8 × 10¹⁹ N m  , significantly more than the seismic moment released by aftershocks,  6.6 × 10¹⁸ N m  . Afterslip and aftershocks appear to have different temporal evolutions and some spatial correlations, suggesting that aftershock rates may not be completely controlled by the rate of afterslip.     

Strain accumulation along an oblique plate boundary: the Reykjanes Peninsula, southwest Iceland

We use annual GPS observations on the Reykjanes Peninsula (RP) from 2000 to 2006 to generate maps of surface velocities and strain rates across the active plate boundary. We find that the surface deformation on the RP is consistent with oblique plate boundary motion on a regional scale, although considerable temporal and spatial strain rate variations are observed within the plate boundary zone. A small, but consistent increase in eastward velocity is observed at several stations on the southern part of the peninsula, compared to the 1993–1998 time period. The 2000–2006 velocities can be modelled by approximating the plate boundary as a series of vertical dislocations with left-lateral motion and opening. For the RP plate boundary we estimate left-lateral motion  18⁺⁴₋₃ mm yr⁻¹  and opening of  7⁺³₋₂ mm yr⁻¹  below a locking depth of  7⁺¹₋₂ km  . The resulting deep motion of  20⁺⁴₋₃ mm yr⁻¹  in the direction of  N(100⁺⁸₋₆)°E  agrees well with the predicted relative North America–Eurasia rate. We calculate the areal and shear strain rates using velocities from two periods: 1993–1998 and 2000–2006. The deep motion along the plate boundary results in left-lateral shear strain rates, which are perturbed by shallow deformation due to the 1994–1998 inflation and elevated seismicity in the Hengill–Hrómundartindur volcanic system, geothermal fluid extraction at the Svartsengi power plant, and possibly earthquake activity on the central part of the peninsula.     

Use of 40Ar/39Ar K-feldspar thermochronology in basin thermal history reconstruction: an example from the Big Lake Suite granites, Warburton Basin, South Australia

The potential use of ⁴⁰Ar/³⁹Ar thermochronologic data from K-feldspars in reconstructing basin thermal history has been evaluated using the example of the Warburton/Cooper/Eromanga Basin, Australia's largest onshore oil- and gas-producing basin. Results from ⁴⁰Ar/³⁹Ar step-heating experiments reveal details of the evolution of the basin system, including the following: (1) the operation of high geothermal gradient regimes during the earliest basin evolution, suggesting that basin formation was active rather than passive; (2) slow cooling from a Permo-Triassic temperature peak of at least 250–300°C; (3) a rise in thermal gradients to contemporary bottom hole temperatures in the last 5–10 Myr; and (4) spatially variable recrystallization events between 100 and 50 Ma and at around 20 Ma. Initial microstructural observations serve as a useful predictor of the quality and nature of the obtainable age information. Data from 'pristine' K-feldspars may constrain the peak temperature conditions experienced in the basin, the basin's early thermal history and also any recent changes in thermal gradient. Contrasting data from texturally modified K-feldspars may constrain times of thermal transients and/or fluid flow, with the preferred interpretation that K-feldspars recrystallize in response to such events. The Warburton/Cooper/Eromanga Basin example suggests that the ⁴⁰Ar/³⁹Ar technique may serve as a useful adjunct to apatite and zircon fission track analysis and conventional organic maturation indices in basin thermal history analysis.     

Slip rate on the Dead Sea transform fault in northern Araba valley (Jordan)

The Araba valley lies between the southern tip of the Dead Sea and the Gulf of Aqaba. This depression, blanketed with alluvial and lacustrine deposits, is cut along its entire length by the Dead Sea fault. In many places the fault is well defined by scarps, and evidence for left-lateral strike-slip faulting is abundant. The slip rate on the fault can be constrained from dated geomorphic features displaced by the fault. A large fan at the mouth of Wadi Dahal has been displaced by about 500 m since the bulk of the fanglomerates were deposited 77–140 kyr ago, as dated from cosmogenic isotope analysis (¹⁰Be in chert) of pebbles collected on the fan surface and from the age of transgressive lacustrine sediments capping the fan. Holocene alluvial surfaces are also clearly offset. By correlation with similar surfaces along the Dead Sea lake margin, we propose a chronology for their emplacement. Taken together, our observations suggest an average slip rate over the Late Pleistocene of between 2 and 6 mm yr⁻¹, with a preferred value of 4 mm yr⁻¹. This slip rate is shown to be consistent with other constraints on the kinematics of the Arabian plate, assuming a rotation rate of about 0.396° Myr⁻¹ around a pole at 31.1°N, 26.7°E relative to Africa.