On the rifting dynamics of plate divergence and magma accumulation at oceanic ridge axes

Rifting dynamics at spreading axes is governed by two processes: the large-scale plate divergence and the local magma accumulation in the crust-mantle transition layer. Both evolve simultaneously. A model is developed particularly for the situation in Iceland where a well studied rifting episode occurred in the Krafla volcanic system 1975–1984. Both the divergence and the buoyant rise of magma create tensile deviatoric stress in the axial region, but while divergence generates an altogether extensional stress field, uprising of buoyant melt produces tension only near the axis but compression of the sides. The buoyant rise is driven by the differential pressure gradient in rock and melt. The processes are studied with a two-dimensional finite-element routine. Presently thermal effects are neglected. The model parameters are density difference, size of the buoyant body, externally applied stress field, mechanical properties of rock and melt. Relatively small amounts of divergence and small increases of buoyancy are shown to generate axial tension which can overcome the tensile strength. Axial tension produced by buoyant bodies can even overcome lateral compression. Observed long intervals of quiescence require either large rock strength, quasi-continuous stress relaxation, small buoyant bodies, and/or a compressive deviatoric stress normal to the axis during much of the time between rifting episodes. Buoyant rise and injection of melt must be important in generating compression.     

Micro-scale modelling in the study of plume evolution in heterogeneous media

The migration of contaminants in heterogeneous aquifers involves dispersive processes that act at different scales. The interaction of these processes as a plume evolves can be studied by micro-scale modelling whereby two scales, a local- or micro-scale and an aquifer- or macro-scale, are covered simultaneously. Local-scale dispersive processes are represented through the local dispersion coefficient in the transport equation, while large-scale dispersion due to heterogeneities is represented through the resolution of the flow field and the diffusive exchange between streamtubes. The micro-scale model provides both the high degree of resolution compatible with local-scale processes, and the extent required for the approach to asymptotic conditions, using grids of up to a million nodal points. The model is based on the dual potential-streamfunction formulation for flow, and the transport problem is formulated in a natural coordinate system provided by the flownet. Simulations can be used to verify stochastic theories of dispersion, without the restrictive assumptions inherent in the theory. For the two-dimensional case, results indicate convergence of the effective dispersivity to the theoretical macrodispersivity value. Convergence takes place within a travel distance of about 50 correlation lengths of the hydraulic conductivity field. However, the approach taken to asymptotic conditions, as well as the macrodispersivity value, may differ for different realizations of the same medium. The influence of early-time events such as plume splitting on the asymptotic convergence remains to be investigated.     

Geochemistry and spatial distribution of OIB and MORB in A’nyemaqen ophiolite zone: Evidence of Majixueshan ancient ridge-centered hotspot

The mafic volcanic association is made up of OIB, E-MORB and N-MORB in the A'nyemaqen Paleozoic ophiolites. Compared with the same type rocks in the world, the mafic rocks generally display lower Nb/U and Ce/Pb ratios and some have Nb depletion and Pb enrichment. The OIB are LREE-enriched with (La/Yb)N =5―20, N-MORB are LREE-depleted with (La/Yb)N = 0.41―0.5. The OIB are featured by incompatible element enrichment and the N-MORB are obviously depleted with some metasomatic effect, and E-MORB are geochemically intermediated. These rocks are distributed around the Majixueshan OIB and gabbros in a thickness greater than a thousand meters and transitionally change along the ophiolite extension in a west-east direction, showing a symmetric distribution pattern as centered by the Majixueshan OIB, that is, from N-MORB, OIB and E-MORB association in the Dur'ngoi area to OIB in the Majixueshan area and then to N-MORB, OIB and E-MORB assemblage again in the Buqingshan area. By consideration of the rock association, the rock spatial distribution and the thickness of the mafic rocks in the Majixueshan, coupled with the metasomatic relationship between the OIB and MORB sources, it can be argued that the Majixueshan probably corresponds to an ancient hotspot or an ocean island formed by mantle plume on the A'nyemaqeh ocean ridge, that is the ridge-centered hotspot, tectonically similar to the present-day Iceland hotspot.     

Surface current variability in the Keum River Estuary (South Korea) during summer 2002 as observed by high-frequency radar and coastal monitoring buoy

High-frequency (HF) radar observations of surface currents were conducted for 3 months during summer 2002 in the Keum River estuary. A comparison between HF radar-derived currents and directly measured ones form a buoy showed that the regression slope is close to 1 and the correlation coefficient greater than 0.86, with an RMS difference less than 13 cm/s which is less than 17% of the tidal current. This fairly good agreement allows us to use HF radar observation in investigating the surface flow and circulation in this tidal-current-dominant coastal-plume area. To examine the spatial variation in tidal current characteristics, as well as currents associated with non-tidal forcing, the HF radar-derived currents were separated into tidal and sub-tidal frequency currents. The overall pattern of M₂-current ellipse distribution in the study area showed a counterclockwise rotation, with the offshore maximum current direction to the northeast. Eccentricity, the direction of maximum current, and the phase of net motion of the ellipse changed near the estuary mouth and near the gap of the Saemangeum reclamation tide dyke due to the complex coastal geometry and the out-flowing jet during the ebb period.     

Tidal effects on estuarine circulation and outflow plume in the Chesapeake Bay

The response of the Chesapeake Bay to river discharge under the influence and absence of tide is simulated with a numerical model. Four numerical experiments are examined: (1) response to river discharge only; (2) response to river discharge plus an ambient coastal current along the shelf outside the bay; (3) response to river discharge and tidal forcing; and (4) response to river discharge, tidal forcing, and ambient coastal current. The general salinity distribution in the four cases is similar to observations inside the bay. Observed features, such as low salinity in the western side of the bay, are consistent in model results. Also, a typical estuarine circulation with seaward current in the upper layer and landward current in the lower layer is obtained in the four cases. The two cases without tide produce stronger subtidal currents than the cases with tide owing to greater frictional effects in the cases with tide. Differences in salinity distributions among the four cases appear mostly outside the bay in terms of the outflow plume structure. The two cases without tide produce an upstream (as in a Kelvin wave sense) or northward branch of the outflow plume, while the cases with tide produce an expected downstream or southward plume. Increased friction in the cases with tide changes the vertical structure of outflow at the entrance to the bay and induces large horizontal variations in the exchange flow. Consequently, the outflow from the bay is more influenced by the bottom than in the cases without tide. Therefore, a tendency for a bottom-advected plume appears in the cases with tide, rather than a surface-advected plume, which develops in the cases without tide. Further analysis shows that the tidal current favors a salt balance between the horizontal and vertical advection of salinity around the plume and hinders the upstream expansion of the plume outside the bay.     

The gas content and buoyancy of strombolian ash plumes

Plinian plumes erupt with a bulk density greater than that of air, and depend upon air entrainment during their gas-thrust phase to become buoyant; if entrainment is insufficient, the column collapses into a potentially deadly pyroclastic flow. This study shows that strombolian ash plumes can be erupted in an initially buoyant state due to their extremely high initial gas content, and in such cases are thus impervious to column collapse. The high gas content is a consequence of decoupled gas rise in the conduit, in which particles are ultimately incidental. The relations between conduit gas flow, eruption style and plume density are explored here for strombolian scenarios and contrasted with conventional wisdom derived from plinian eruptions. Considering the inherent relation between gas content and initial plume density together with detailed measurements of plume velocities can help unravel ambiguities surrounding conduit processes, eruption styles and hazards at poorly understood volcanoes. Analysis of plume dynamics at Santiaguito volcano, Guatemala adds further support for a model involving decoupled gas rise in the conduit.     

Azimuthal variation of the P phase in Icelandic receiver functions

A curious observation has been made on radial receiver functions calculated from teleseisms recorded by 29 broad-band seismometers distributed over Iceland. The arrival time of the direct P phase of the radial receiver functions depends critically upon the azimuth of the teleseismic source. For a seismic station in West Iceland, the direct P  phase of the radial receiver function arrives consistently later for easterly source azimuths than for westerly source azimuths. The reverse applies for stations in East Iceland. In the original seismograms, the delayed P phase of the receiver function appears up to 450 ms later on the radial than on the vertical component. The seismometer locations in East and West Iceland are separated by the Neovolcanic Zone, a constructive plate boundary. The delayed P phases occur for seismic rays travelling across this zone. However, it is not obvious how wave propagation across the plate boundary zone could cause the observed delays. The tentative explanation proposed here involves the regional dip of the Icelandic lava sequences towards the Neovolcanic Zone. A dipping interface at shallow depth results in a P–S converted phase arriving shortly after the P phase. These phases cannot be separated in the radial receiver functions, given the bandwidth of the observed signals. However, a calculation of receiver functions from estimates of the P , SV and SH wavefields clearly reveals a P–S converted phase at about 500 ms for easterly source azimuths in West Iceland and for westerly source azimuths in East Iceland. The amplitudes of the direct P phase and the P–S phase converted at a dipping interface would be expected to vary strongly with azimuth in accordance with the observed behaviour.     

Geochemistry and magnetostratigraphy of deccan flows at Anjar, Kutch

Chemical analysis of nine Deccan flow basalts at Anjar, Kutch, western India, indicates that all, except the uppermost flow F-9, are alkaline. In their major and trace element composition, the alkali basalts resemble Ocean island basalts (OIB). Similarities of many diagnostic trace element ratios (e.g. Sm/Nd, Ba/Nb,Y/Nb and Zr/Nb) are similar to those found in the Réunion Island basalts. The uppermost basalt is tholeiitic and chemically resembles the least contaminated Deccan basalt (Ambenali type). The Anjar basalts have iridium concentration ranging between 2 and 178 pg/g. Some of these values are higher by about an order of magnitude compared to the Ir concentration in other basalts of the Deccan. A synthesis of chemical, palaeomagnetic and geochronologic data enables us to construct a chemical and magnetic stratigraphy for these flows. The three flows below the iridium enriched intertrappean bed (IT III) show normal magnetic polarity whereas all except one of the upper basalts show reversed magnetic polarity. The sequence seems to have started in polarity zones 31N and probably continued up to 28R or 27R. The results presented here support the view that Deccan volcanism in Kutch occurred on a time span of a few million years.