Chemical and textural equilibration of garnet during amphibolite facies metamorphism: The influence of coupled dissolution–reprecipitation

Metamorphic equilibration requires chemical communication between minerals and may be inhibited through sluggish volume diffusion and or slow rates of dissolution in a fluid phase. Relatively slow diffusion and the perceived robust nature of chemical growth zoning may preclude garnet porphyroblasts from readily participating in low‐temperature amphibolite facies metamorphic reactions. Garnet is widely assumed to be a reactant in staurolite‐isograd reactions, and the evidence for this has been assessed in the Late Proterozoic Dalradian pelitic schists of the Scottish Highlands. The 3D imaging of garnet porphyroblasts in staurolite‐bearing schists reveals a good crystal shape and little evidence of marginal dissolution; however, there is also lack of evidence for the involvement of either chlorite or chloritoid in the reaction. Staurolite forms directly adjacent to the garnet, and its nucleation is strongly associated with deformation of the muscovite‐rich fabrics around the porphyroblasts. “Cloudy” fluid inclusion‐rich garnet forms in both marginal and internal parts of the garnet porphyroblast and is linked both to the production of staurolite and to the introduction of abundant quartz inclusions within the garnet. Such cloudy garnet typically has a Mg‐rich, Mn‐poor composition and is interpreted to have formed during a coupled dissolution–reprecipitation process, triggered by a local influx of fluid. All garnet in the muscovite‐bearing schists present in this area is potentially reactive, irrespective of the garnet composition, but very few of the schists contain staurolite. The staurolite‐producing reaction appears to be substantially overstepped during the relatively high‐pressure Barrovian regional metamorphism reflecting the limited permeability of the schists in peak metamorphic conditions. Fluid influx and hence reaction progress appear to be strongly controlled by subtle differences in deformation history. The remaining garnet fails to achieve chemical equilibrium during the reaction creating distinctive patchy compositional zoning. Such zoning in metamorphic garnet created during coupled dissolution–reprecipitation reactions may be difficult to recognize in higher grade pelites due to subsequent diffusive re‐equilibration. Fundamental assumptions about metamorphic processes are questioned by the lack of chemical equilibrium during this reaction and the restricted permeability of the regional metamorphic pelitic schists. In addition, the partial loss of prograde chemical and textural information from the garnet porphyroblasts cautions against their routine use as a reliable monitor of metamorphic history. However, the partial re‐equilibration of the porphyroblasts during coupled dissolution–reprecipitation opens possibilities of mapping reaction progress in garnet as a means of assessing fluid access during peak metamorphic conditions.     

An investigation of upper mantle heterogeneity beneath the Archaean and Proterozoic crust of western Canada from Lithoprobe controlled-source seismic experiments

Observations of upper mantle reflectivity at numerous locations around the world have been linked to the presence of a heterogeneous distribution of rock types within a broad layer of the upper mantle. This phenomenon is observed in wide-angle reflection data from Lithoprobe's Alberta Basement Transect [the SAREX and Deep Probe experiments of 1995] and Trans-Hudson Orogen Transect [the THoRE experiment of 1993]. SAREX and Deep Probe image the Archaean lithosphere of the Hearne and Wyoming Provinces, whereas THoRE images the Archaean and Proterozoic lithosphere of the Trans-Hudson Orogen and neighbouring areas.Finite-difference synthetic seismograms are used to constrain the position and physical properties of the reflective layer. SAREX/Deep Probe modelling uses a 2-D visco-elastic finite-difference routine; THoRE modelling uses a pseudospectral algorithm. In both cases, the upper mantle is parameterized in terms of two media. One medium is the background matrix; the other is statistically distributed within the first as a series of elliptical bodies. Such a scheme is suitable for modelling: (1) variations in lithology (e.g., a peridotite matrix with eclogite lenses) or (2) variations in rheology (e.g., lenses of increased strain within a less strained background).The synthetic seismograms show that the properties of heterogeneities in the upper mantle do not change significantly between the two Lithoprobe transects. Beneath the Trans-Hudson Orogen in Saskatchewan, the layer is best modelled to lie at depths between 80 and 150 km. Based on observations from perpendicular profiles, anisotropy of the heterogeneities is inferred. Beneath the Precambrian domains of Alberta, 400 km to the west, upper mantle heterogeneities are modelled to occur between depths of 90 and 140 km. In both cases the heterogeneous bodies within the model have cross-sectional lengths of tens of kilometers, vertical thicknesses less than 1 km, and velocity contrasts from the background of − 0.3 to − 0.4 km/s. Based on consistency with complementary data and other results, the heterogeneous layer is inferred to be part of the continental lithosphere and may have formed through lateral flow or deformation within the upper mantle.     

Attenuation relations for strong seismic ground motion in the Himalayan region

Strong motion data from various regions of India have been used to study attenuation characteristics of horizontal peak acceleration and velocity. The strong ground motion data base considered in the present work consists of various earthquakes recorded in the northern part of India since 1986 with magnitudes 5.7 to 7.2. Using these data, relations for horizontal peak acceleration and velocity, which are $$\begin{gathered} log_{10} a = 1.14 + 0.31M + 0.65log_{10} R \hfill \\ log_{10} v = 0.571 + 0.41M + 0.768log_{10} R \hfill \\ \end{gathered} $$ have been proposed wherea is the peak horizontal acceleration in cm/sec²,v is the peak horizontal velocity in mm/sec,M is body wave magnitude, andR is the hypocentral distance in km. The proposed relations are in reasonable agreement with the small amount of strong ground motion data available for the northern part of India. The present results will be useful in estimating strong ground motion parameters and in the earthquake resistant design in the Himalayan region.     

The New Multichannel Radiospectrograph ARTEMIS-IV/HECATE, of the University of Athens

We present the new solar radiospectrograph of the University of Athensoperating at the Thermopylae Station since 1996. Observations cover thefrequency range from 110 to 688 MHz. The radiospectrograph has a 7-meterparabolic antenna and two receivers operating in parallel. One is a sweepfrequency receiver and the other a multichannel acousto-optical receiver.The data acquisition system consists of a front-end VME based subsystem anda Sun Sparc-5 workstation connected through Ethernet. The two subsystems areoperated using the VxWorks real-time package. The daily operation is fullyautomated: pointing of the antenna to the sun, starting and stopping theobservations at pre-set times, data acquisition, data compression by`silence suppression', and archiving on DAT tapes. The instrument can beused either by itself to study the onset and evolution of solar radio bursts or in conjunction with other instruments including theNançay Decametric Array and the WIND/WAVES RAD1 and RAD2 low frequencyreceivers to study associated interplanetary phenomena.     

Flexure of Europa's lithosphere due to ridge-loading

Using photoclinometry, topographic profiles across europan ridges have been produced. These profiles allow the identification of bulges in the terrain adjacent to the ridges. The bulges are assumed to have been produced by flexure of the elastic lithosphere due to the load of the ridges, which lie along cracks in the crust. The distance from the crack to these “fore-bulges” depends on the thickness of the elastic plate being flexed. Based on a survey of ridges in Galileo images with resolution <300 m/pixel, the thickness of the elastic lithosphere has been determined by this method at a wide variety of sites along the leading and trailing hemispheres of Europa. The average thickness is about 200 m. The elastic lithosphere underneath smooth dilational bands tends to be thicker than plains morphology, an effect that is pronounced at Thynia Linea and Astypalaea Linea. Among the ridges investigated here, more recent loading correlates with a thicker elastic lithosphere, which may either reflect an intrinsically thicker layer, or less viscous relaxation over the shorter time period.     

Evaluating the effect of land development on sediment transport using a probability density function

An important problem in sedimentation analysis is the development of a channel section that preserves, as best as possible, the current sedimentation regime even though the flood frequency tendencies have been altered due to land development within the catchment. In order to accomplish this task, a methodology is needed that estimates sediment transport capacity for various channel configurations. Such a procedure is described which allows the computation of the total sediment transport capacity for each of several T-year return frequency runoff hydrographs. This information is used to obtain an approximate probability distribution for the total sediment transport capacity, and the mean and standard deviation of this distribution are computed.Comparing the results for the catchment in its present state with a future developed state, using a selection of new channel parameters, indicates how to improve the channel to control changes in sedimentation due to development. The analysis procedure provides a basis for estimating a new channel configuration such that the new flow conditions retain, as best as possible, the existing condition sedimentation effects, and hence retain the natural sediment supply and transport trends even though runoff flow rates have changed due to land development within the catchment.The results of Wilson Creek are typical of the several sites examined, see Table 3 below. The T=2, T=5, T=25, and T=100 year values for total sediment transport capacity, in kilotons, are 6.9, 39.4, 61.3, and 96.7 with a mean of 17.1 and standard deviation of 19.3. After development with no change in the channel the respective values increase to: 17.9, 84.6, 128.1, and 258.0 with a mean of 39.1 and standard deviation of 44.3. A new channel can be constructed which will reduce these sediment transport capacity values, after development, to 5.2, 41.0, 62.0, and 124.8 with a mean of 17.4 and standard deviation of 22.0.     

Changes in current regime and turbidity in response to a freshwater pulse in the eastmain estuary

In September 1984, the freshwater input to the Eastmain River (James Bay, Canada) was increased by a factor of 50 over a 6-d period during a controlled reservoir discharge. Changes to the current, salinity and turbidity regimes were monitored during the peak runoff. Estuarine salinity values fell rapidly with increasing mean flow, as did the amplitude of the semi-diurnal tidal currents. A large increase in bottom shear stress dispersed the settled suspension layer into the water column, raising concentrations of suspended matter in the estuary by a factor of 4 in 3 d. The peak values exceeded 150 mg I^?1. This led to erosion of the river silt deposits, with the export of an estimated 6 × 10⁴ metric tons of sediments. After the reduction of discharge, current values returned to their normal range within a day, whereas upstream salinity intrusion occurred at a slower rate. Horizontal diffusivities of about 100 m² s^?1 were required to match the observed lag.     

Interannual variability of Kelvin wave propagation in the wave guides of the equatorial Indian Ocean,the coastal Bay of Bengal and the southeastern Arabian Sea during 1993–2006

The observed variability of the Kelvin waves and their propagation in the equatorial wave guide of the Indian Ocean and in the coastal wave guides of the Bay of Bengal (BoB) and the southeastern Arabian Sea (AS) on seasonal to interannual time scales during years 1993–2006 is examined utilizing all the available satellite and in-situ measurements. The Kelvin wave regime inferred from the satellite-derived sea surface height anomalies (SSHA) shows a distinct annual cycle composed of two pairs of alternate upwelling (first one occurring during January–March and the second one occurring during August–September) and downwelling (first one occurring during April–June and the second one occurring during October–December) Kelvin waves that propagate eastward along the equator and hit the Sumatra coast and bifurcate. The northern branches propagate counterclockwise over varied distances along the coastal wave guide of the BoB. The potential mechanisms that contribute to the mid-way termination of the first upwelling and the first downwelling Kelvin waves in the wave guide of the BoB are hypothesized. The second downwelling Kelvin wave alone reaches the southeastern AS, and it shows large interannual variability caused primarily by similar variability in the equatorial westerly winds during boreal fall. The westward propagating downwelling Rossby waves triggered by the second downwelling Kelvin wave off the eastern rim of the BoB also shows large interannual variability in the near surface thermal structure derived from SODA analysis. The strength of the equatorial westerlies driven by the east–west gradient of the heat sources in the troposphere appears to be a critical factor in determining the observed interannual variability of the second downwelling Kelvin wave in the wave guides of the equatorial Indian Ocean, the coastal BoB, and the southeastern AS.     

Three-dimensional visualization of mission planning and control forthe NPS autonomous underwater vehicle

The Naval Postgraduate School (NPS) is constructing a small autonomous underwater vehicle (AUV) with an onboard mission control computer. The mission controller software for this vehicle is a knowledge-based artificial intelligence (AI) system requiring thorough analysis and testing before the AUV is operational. The manner in which rapid prototyping of this software has been demonstrated by developing a controller code on a LISP machine and using an Ethernet link with a graphics workstation to simulate the controller's environment is discussed. The development of a testing simulator using a knowledge engineering environment (KEE) expert system shell that examines AUV controller subsystems and vehicle models before integrating them with the full AUV for its test environment missions is discussed. This AUV simulator utilizes an interactive mission planning control console and is fully autonomous once initial parameters are selected     

Impact of natural (waves and currents) and anthropogenic (trawl) resuspension on the export of particulate matter to the open ocean: Application to the Gulf of Lion (NW Mediterranean)

Modern sediment deposits on continental margins form a vast reservoir of particulate matter that is regularly affected by resuspension processes. Resuspension by bottom trawling on shelves with strong fishing activity can modify the scale of natural disturbance by waves and currents. Recent field data show that the impact of bottom trawls on fine sediment resuspension per unit surface is comparable with that of the largest storms.