Postglacial isostatic adjustment in a self-gravitating spherical earth with power-law rheology

Since microphysics cannot say definitively whether the rheology of the mantle is linear or non-linear, the aim of this paper is to constrain mantle rheology from observations related to the glacial isostatic adjustment (GIA) process—namely relative sea-levels (RSLs), land uplift rate from GPS and gravity-rate-of-change from GRACE. We consider three earth model types that can have power-law rheology (n = 3 or 4) in the upper mantle, the lower mantle or throughout the mantle. For each model type, a range of A parameter in the creep law will be explored and the predicted GIA responses will be compared to the observations to see which value of A has the potential to explain all the data simultaneously. The coupled Laplace finite-element (CLFE) method is used to calculate the response of a 3D spherical self-gravitating viscoelastic Earth to forcing by the ICE-4G ice history model with ocean loads in self-gravitating oceans. Results show that ice thickness in Laurentide needs to increase significantly or delayed by 2 ka, otherwise the predicted uplift rate, gravity rate-of-change and the amplitude of the RSL for sites inside the ice margin of Laurentide are too low to be able to explain the observations. However, the ice thickness elsewhere outside Laurentide needs to be slightly modified in order to explain the global RSL data outside Laurentide. If the ice model is modified in this way, then the results of this paper indicate that models with power-law rheology in the lower mantle (with A ∼ 10⁻³⁵ Pa⁻³ s⁻¹ for n = 3) have the highest potential to simultaneously explain all the observed RSL, uplift rate and gravity rate-of-change data than the other model types.     

Using postglacial sea level, crustal velocities and gravity-rate-of-change to constrain the influence of thermal effects on mantle lateral heterogeneities

Lateral heterogeneities in the mantle can be caused by thermal, chemical and non-isotropic pre-stress effects. Here, we investigate the possibility of using observations of the glacial isostatic adjustment (GIA) process to constrain the thermal contribution to lateral variations in mantle viscosity. In particular, global historic relative sea level, GPS in Laurentide and Fennoscandia, altimetry together with tide-gauge data in the Great Lakes area, and GRACE data in Laurentide are used. The lateral viscosity perturbations are inferred from the seismic tomography model S20A by inserting the scaling factor β to determine the contribution of thermal effects versus compositional heterogeneity and non-isotropic pre-stress effects on lateral heterogeneity in mantle viscosity. When β = 1, lateral velocity variations are caused by thermal effects alone. With β < 1, the contribution of thermal effect decreases, so that for β = 0, there is no lateral viscosity variation and the Earth is laterally homogeneous. These lateral viscosity variations are superposed on four different reference models which differ significantly in the lower mantle viscosity. The Coupled Laplace Finite Element method is used to predict the GIA response on a spherical, self-gravitating, compressible, viscoelastic Earth with self-gravitating oceans, induced by the ICE-4G deglaciation model.Results show that the effect of β on uplift rates and gravity rate-of-change is not simple and involves the trade-off between the contribution of lateral viscosity variations in the transition zone and in the lower mantle. Models with small viscosity contrast in the lower mantle cannot explain the observed uplift rates in Laurentide and Fennoscandia. However, the RF3S20 model with a reference viscosity profile simplified from Peltier's VM2 with the value of β around 0.2–0.4 is found to explain most of the global RSL data, the uplift rates in Laurentide and Fennoscandia and the BIFROST horizontal velocity data. In addition, the changes in GIA signals caused by changes in the value of β are large enough to be detected by the data, although uncertainty in other parameters in the GIA models still exists. This may encourage us to further utilize GIA observations to constrain the thermal effect on mantle lateral heterogeneity as geodetic and satellite gravity measurements are improved.     

Using postglacial sea level,crustal velocities and gravity-rate-of-change to constrain the influence of thermal effects on mantle lateral heterogeneities

Lateral heterogeneities in the mantle can be caused by thermal, chemical and non-isotropic pre-stress effects. Here, we investigate the possibility of using observations of the glacial isostatic adjustment (GIA) process to constrain the thermal contribution to lateral variations in mantle viscosity. In particular, global historic relative sea level, GPS in Laurentide and Fennoscandia, altimetry together with tide-gauge data in the Great Lakes area, and GRACE data in Laurentide are used. The lateral viscosity perturbations are inferred from the seismic tomography model S20A by inserting the scaling factor β to determine the contribution of thermal effects versus compositional heterogeneity and non-isotropic pre-stress effects on lateral heterogeneity in mantle viscosity. When β = 1, lateral velocity variations are caused by thermal effects alone. With β < 1, the contribution of thermal effect decreases, so that for β = 0, there is no lateral viscosity variation and the Earth is laterally homogeneous. These lateral viscosity variations are superposed on four different reference models which differ significantly in the lower mantle viscosity. The Coupled Laplace Finite Element method is used to predict the GIA response on a spherical, self-gravitating, compressible, viscoelastic Earth with self-gravitating oceans, induced by the ICE-4G deglaciation model.Results show that the effect of β on uplift rates and gravity rate-of-change is not simple and involves the trade-off between the contribution of lateral viscosity variations in the transition zone and in the lower mantle. Models with small viscosity contrast in the lower mantle cannot explain the observed uplift rates in Laurentide and Fennoscandia. However, the RF3S20 model with a reference viscosity profile simplified from Peltier's VM2 with the value of β around 0.2–0.4 is found to explain most of the global RSL data, the uplift rates in Laurentide and Fennoscandia and the BIFROST horizontal velocity data. In addition, the changes in GIA signals caused by changes in the value of β are large enough to be detected by the data, although uncertainty in other parameters in the GIA models still exists. This may encourage us to further utilize GIA observations to constrain the thermal effect on mantle lateral heterogeneity as geodetic and satellite gravity measurements are improved.     

Glacial isostatic adjustment in 3-D earth models: Implications for the analysis of tide gauge records along the U.S. east coast

Tide gauge records of recent sea-level change along the U.S. east coast have received significant attention within the literature of glacial isostatic adjustment (GIA). Geographic trends in these tide gauge rates are not reduced by a GIA correction based on a commonly adopted radial viscosity profile (characterized, in particular, by a lower mantle viscosity 1−2×10²¹ Pa s), and this has led to speculation that the residual trends reflect contributions from neotectonics or oceanographic processes. While the trends can be significantly reduced by adopting an Earth model with a stiffer lower mantle, such a model appears to be incompatible with independent constraints from post-glacial decay times in Hudson Bay. We use a finite-element model of the GIA process to investigate whether 3-D viscosity variations superimposed onto the “common” radial viscosity profile may provide a route to reconciling the east coast sea-level trends. We find that the specific 3-D structure we impose has little impact on the geographic trends in the GIA-corrected rates. However, we do find that the imposed lateral variations in lower mantle viscosity introduce a nearly uniform upward shift of 0.5 mm/yr in GIA-induced sea-level rates along the U.S. east coast. Thus, inferences of regional (U.S. east coast) sea-level rise due to modern melting of ice reservoirs, based on tide gauge rates corrected using 1-D GIA models, may be significantly biased by this simplifying assumption.     

Glacial isostatic adjustment at the Laurentide ice sheet margin: Models and observations in the Great Lakes region

Glacial Isostatic Adjustment (GIA) modelling in North America relies on relative sea level information which is primarily obtained from areas far away from the uplift region. The lack of accurate geodetic observations in the Great Lakes region, which is located in the transition zone between uplift and subsidence due to the deglaciation of the Laurentide ice sheet, has prevented more detailed studies of this former margin of the ice sheet. Recently, observations of vertical crustal motion from improved GPS network solutions and combined tide gauge and satellite altimetry solutions have become available. This study compares these vertical motion observations with predictions obtained from 70 different GIA models. The ice sheet margin is distinct from the centre and far field of the uplift because the sensitivity of the GIA process towards Earth parameters such as mantle viscosity is very different. Specifically, the margin area is most sensitive to the uppermost mantle viscosity and allows for better constraints of this parameter. The 70 GIA models compared herein have different ice loading histories (ICE-3/4/5G) and Earth parameters including lateral heterogeneities. The root-mean-square differences between the 6 best models and the two sets of observations (tide gauge/altimetry and GPS) are 0.66 and 1.57 mm/yr, respectively. Both sets of independent observations are highly correlated and show a very similar fit to the models, which indicates their consistent quality. Therefore, both data sets can be considered as a means for constraining and assessing the quality of GIA models in the Great Lakes region and the former margin of the Laurentide ice sheet.     

Glacial isostatic adjustment at the Laurentide ice sheet margin: Models and observations in the Great Lakes region

Glacial Isostatic Adjustment (GIA) modelling in North America relies on relative sea level information which is primarily obtained from areas far away from the uplift region. The lack of accurate geodetic observations in the Great Lakes region, which is located in the transition zone between uplift and subsidence due to the deglaciation of the Laurentide ice sheet, has prevented more detailed studies of this former margin of the ice sheet. Recently, observations of vertical crustal motion from improved GPS network solutions and combined tide gauge and satellite altimetry solutions have become available. This study compares these vertical motion observations with predictions obtained from 70 different GIA models. The ice sheet margin is distinct from the centre and far field of the uplift because the sensitivity of the GIA process towards Earth parameters such as mantle viscosity is very different. Specifically, the margin area is most sensitive to the uppermost mantle viscosity and allows for better constraints of this parameter. The 70 GIA models compared herein have different ice loading histories (ICE-3/4/5G) and Earth parameters including lateral heterogeneities. The root-mean-square differences between the 6 best models and the two sets of observations (tide gauge/altimetry and GPS) are 0.66 and 1.57 mm/yr, respectively. Both sets of independent observations are highly correlated and show a very similar fit to the models, which indicates their consistent quality. Therefore, both data sets can be considered as a means for constraining and assessing the quality of GIA models in the Great Lakes region and the former margin of the Laurentide ice sheet.     

Laurentia crustal motion observed using TOPEX/POSEIDON radar altimetry over land

A new method to estimate the vertical crustal motion from satellite altimetry over land was developed. The method was tested around Hudson Bay, where the observed vertical motion is largely caused by the incomplete glacial isostatic adjustment (GIA) as a result of the Laurentide ice sheet deglaciation since the last glacial maximum (LGM). Decadal (1992–2003) TOPEX/POSEIDON radar altimetry data over land surfaces were used. The results presented here are improved compared to a previous study (Lee, H., Shum, C.K., Kuo, C.Y., Yi, Y., Braun, A., 2008. Application of TOPEX altimetry for solid Earth deformation studies. Terr. Atmos. Ocean. Sci. 19, 37–46. doi:10.3319/TAO.2008.19.1-2.37(SA).) which estimated vertical motion only over relatively flat land surfaces (standard deviation of the height variation <40 cm). In this study, we extended the concept of traditional 1-Hz (one-per-frame) radar altimeter ocean stackfiles to build 10-Hz (10-per-frame) land stackfiles over Hudson Bay land regions, and succeeded in obtaining vertical motion estimates over much rougher surfaces (standard deviation of the height variation <2 m). 90-m C-band Shuttle Radar Topography Mission (SRTM) Digital Elevation Model (DEM) is used as a reference surface to select an optimal waveform retracker, to correct surface gradient errors, and to calculate land surface anomalies. Here, we developed an alternative retracker, called the modified threshold retracker, resulting in decadal vertical motion time series over a 1500 km by 1000 km region covering northern Ontario, northeastern Manitoba, and the Great Lakes region which is at the margin of the former Laurentide ice sheet. The average of the estimated uncertainties for the vertical motion is 2.9 mm/year which is comparable to 2.1 mm/year of recent GPS solutions. The estimated vertical motion is compared with other geodetic observations from GPS, tide gauge/altimetry, GRACE, and several GIA models. The data agree best with the laterally varying 3D GIA model, RF3S20 (β = 0.4) whereas the combination of land altimetry solution with other measurements match best with the models RF3S20 (β = 0.0) or RF3S20 (β = 0.2) in terms of mean and standard deviation of the differences. It is anticipated that this innovative technique could potentially be used to provide additional constraints for GIA model improvement, and be applied to other geodynamics studies.     

Efficient trapping of organic carbon in sediments on the continental margin with high fluvial sediment input off southwestern Taiwan

Accumulation rates of marine and terrigenous organic carbon in the continental margin sediments off southwestern Taiwan were estimated from the measured concentrations and isotopic compositions of total organic carbon (TOC) and previously reported sedimentation rates. Surficial sediments were collected from the study area spanning from the narrow shelf near the Kaoping River mouth to the deep slope with depths reaching almost 3000 m. The average sediment loading of Kaoping River is 17 Mt/yr, which yields high sediment accumulation rates ranging from 0.08 to 1.44 g cm⁻² yr⁻¹ in the continental margin. About half of the discharged sediments were deposited on the margin within 120 km of the river mouth. Carbon isotopic compositions of terrestrial and marine end-members of organic matter were determined, respectively, based on suspended particulate matter (SPM) collected from three major rivers in the southwestern Taiwan and from an offshore station. All samples were analyzed for the TOC content and its isotopic composition (δ¹³C_org). The SPM samples were also analyzed for the total nitrogen (TN) content. TOC content in marine sediments ranges from 0.45% to 1.35% with the highest values on the upper slope near the Kaoping River mouth. The TOC/TN ratio of the SPM samples from the offshore station is 6.8±0.6, almost identical to the Redfield ratio, indicating their predominantly marine origin; their δ¹³C_org values are also typically marine with a mean of −21.5±0.3‰. The riverine SPM samples exhibit typical terrestrial δ¹³C_org values around −25‰. The δ¹³C_org values of surficial sediments range from −24.8‰ to −21.2‰, showing a distribution pattern influenced by inputs from the Kaoping River. The relative contributions from marine and terrestrial sources to sedimentary organic carbon were determined by the isotope mixing model with end-member compositions derived from the riverine and marine SPM. High fluvial sediment inputs lead to efficient trapping of organic carbon over a wide range of water depth in this continental margin. The marine organic accumulation rate ranges from 1.6 to 70 g C m⁻² yr⁻¹ with an area weighted mean of 4.2 g C m⁻² yr⁻¹, which is on a par with the mean terrestrial contribution and accounts for 2.3% of mean primary production. The depth-dependent accumulation rate of marine organic carbon can be simulated with a function involving primary productivity and mineral accumulation rate, which may be applicable to other continental margins with high sedimentation rates. Away from the nearshore area, the content of terrigenous organic carbon in surficial sediments decreases with distance from the river mouth, indicating its degradation in marine environments.     

Organic biomarkers for tracing carbon cycling in the Gulf of Papua (Papua New Guinea)

Sediment traps were deployed in the Gulf of Papua in June–July 1997, to determine fluxes of organic matter and inorganic elements from the photic zone to deeper waters at the base of the continental slope and in the northern Coral Sea. Three stations, ranging from 900 to 1500 m depth, had “shallow” traps at 300 m below the water surface and “deep” traps set 100 m above the bottom. Infiltrex II water samplers collected particulate and dissolved organic matter from the Fly, Purari and Kikori rivers, and near-surface water from the shelf of the Gulf of Papua. Samples were analysed for molecular organic biomarkers to estimate the sources of organic carbon and its cycling processes.Dry weight fluxes from the shallow traps ranged from 115 to 181 mg m⁻² day⁻¹ and particulate organic carbon (POC) fluxes ranged from 1.2 to 1.9 mM OC m⁻² d⁻¹ with molar organic carbon to particulate nitrogen ratios (C/N) ranging from 6.0 to 6.5. Fluxes in deep traps were likely influenced by both early diagenesis and entrapment of resuspended shelf sediments. Dry weight fluxes in deep traps ranged from 106 to 574 mg m⁻² day⁻¹ and POC fluxes ranged from 0.6 to 1.5 mM OC m⁻² d⁻¹, with C/N ratios ranging from 8.5 to 10.8. ¹³C/¹²C ratios were −20.2‰ to −21.7‰ in all trap samples, indicating that most of the settling POC was “marine-derived”. Shallow traps had δ¹⁵N values of 6.3‰ to 7.2‰ while the values in deep traps were 4.9–5.0‰, indicating the N-rich near-surface OC was less degraded than that in the deep traps. The biogenic lipids consisted of hydrocarbon, sterol and fatty acid biomarkers indicative of marine zooplankton, phytoplankton and bacteria. Sterol markers for diatoms and dinoflagellates were abundant in the water samples. Highly branched isoprenoid alkenes, usually attributable to diatoms, were also detected in both water and shallow traps. Traces of C₂₆–C₃₄ n-alcohols indicative of land–plant biomarkers, were found in river water samples and in the shallow sediment traps. A large unresolved complex mixture (UCM) of hydrocarbons, and a uniform distribution of n-alkanes, indicative of petroleum hydrocarbons, were also detected in the traps. Hopane and sterane biomarkers detected in the trap oil were characteristic of a marine carbonate source, and the aromatic hydrocarbon composition distinguished at least two different oil signatures.We concluded that mass and POC fluxes were similar to those reported for other continental shelves and marginal oceans in tropical and subtropical regions. There was a dramatic decrease in POC as particles sank, due to zooplankton repackaging and photochemical and bacterial decomposition. Carbon isotopic and biomarker patterns showed most of the POC in the sediment traps was marine-sourced with only traces of terrestrial input. There was a significant flux of petroleum, which may signal the existence of natural petroleum seeps in this region.     

Transport of oceanic nitrate from the continental shelf to the coastal basin in relation to the path of the Kuroshio

Hydrographic and biogeochemical observations were conducted along the longitudinal section from Ise Bay to the continental margin (southern coast of Japan) to investigate changes according to the Kuroshio path variations during the summer. The strength of the uplift of the cold deep water was influenced by the surface intrusion of the Kuroshio water to the shelf region. When the intrusion of the Kuroshio surface water to the shelf region was weak in 2006, the cold and NO₃⁻-rich shelf water intruded into the bottom layer in the bay from the shelf. This bottom intrusion was intensified by the large river discharge. The nitrogen isotope ratio (δ¹⁵N) of NO₃⁻ (4–5‰) in the bottom bay water was same as that in the deeper NO₃⁻ over the shelf, indicating the supply of new nitrogen to the bay. The warm and NO₃⁻-poor shelf water intruded into the middle layer via the mixing region at the bay mouth when the Kuroshio water distributed in the coastal areas off Ise Bay in 2005. The regenerated NO₃⁻ with isotopically light nitrogen (δ¹⁵N=−1‰) was supplied from the shelf to the bay. This NO₃⁻ is regenerated by the nitrification in the upper layer over the shelf. The contribution rate of regenerated NO₃⁻ over the shelf to the total NO₃⁻ in the subsurface chlorophyll maximum layer in the bay was estimated at 56% by a two-source mixing model coupled with the Rayleigh equation.     

Mantle rheology: radial and lateral viscosity variations inferred from microphysical creep laws

Inferences on the rheology of the mantle based on theoretical and experimental rate equations for steady state creep are discussed and compared with results from geophysical models. The radial increase of viscosity by one to three orders of magnitude across the mantle, required by inversion of postglacial rebound and geodynamic data, is confirmed by microphysical models based on the estimation of continuous and discontinuous changes of creep parameters with depth. The upper mantle (viscosity 10²⁰–10²¹ Pa s) is likely to show non-Newtonian rheology (power-law creep) for average grain sizes larger than 0.1 mm as an order of magnitude. Given the variability of both grain size and stress conditions, local regions of linear rheology can be present. The rheology of transition zone and lower mantle (viscosity 10²²–10²⁴ Pa s) cannot be definitely resolved at present. Estimation of creep parameters leads to possible nonlinear or mixed rheology, if grain sizes are not lower than 0.1 mm and flow conditions can be approximated by a constant strain rate of about 10⁻¹⁵ s⁻¹. This conclusion can be modified by different flow conditions (e.g. a decrease in strain rate or constant viscous dissipation). Furthermore, experiments on fine-grained garnetites and perovskite analogues have shown that diffusion creep is predominant at laboratory conditions. However, the pressure dependence of creep in these phases is unknown, and therefore direct extrapolation to lower mantle conditions is necessarily speculative. Lateral variations of viscosity, largest in the upper and lowermost mantle (up to 2–4 orders of magnitude) are predicted by models based on lateral temperature anomalies derived from seismic tomographic models.     

Hemolytic activity and fatty acids composition in the ichthyotoxic dinoflagellate Cochlodinium polykrikoides isolated from Bahía de La Paz, Gulf of California

The hemolytic activity of the dinoflagellate Cochlodinium polykrikoides from Bahía de La Paz, Gulf of California was investigated as part of the ichthyotoxic mechanism of this microalga. Two different kinds of erythrocytes, fish and human, were tested for the hemolytic assay. Since fatty acids have been associated with hemolytic activity in C. polykrikoides, the composition of fatty acids of this dinoflagellate was also analyzed. The concentration of C. polykrikoides causing 50% hemolysis (HE₅₀) was 4.88 and 5.27 × 10⁶ cells L⁻¹, for fish and human erythrocytes, respectively. According to the standard curve of saponin, an equivalence between the hemolytic activity of saponin and the dinoflagellate concentration was found with 1 μg saponin mL⁻¹ equivalent to 1 × 10⁶ cells L⁻¹ of C. polykrikoides. The polyunsaturated fatty acids: hexadecaenoic (16:0), docosahexaenoic (22:6 n3) and octadecapentaenoic (18:5 n3) were found in an abundance of 62% of total fatty acids.     

Thermal equation of state of magnesite to 32 GPa and 2073 K

Pressure–volume–temperature relations have been measured to 32 GPa and 2073 K for natural magnesite (Mg_0.975Fe_0.015Mn_0.006Ca_0.004CO₃) using synchrotron X-ray diffraction with a multianvil apparatus at the SPring-8 facility. A least-squares fit of the room-temperature compression data to a third-order Birch–Murnaghan equation of state (EOS) yielded K₀ = 97.1 ± 0.5 GPa and K′ = 5.44 ± 0.07, with fixed V₀ = 279.55 ± 0.02 Å³. Further analysis of the high-temperature compression data yielded the temperature derivative of the bulk modulus (∂K_T/∂T)_P = −0.013 ± 0.001 GPa/K and zero-pressure thermal expansion α = a₀ + a₁T with a₀ = 4.03 (7) × 10⁻⁵ K⁻¹ and a₁ = 0.49 (10) × 10⁻⁸ K⁻². The Anderson–Grüneisen parameter is estimated to be δ_T = 3.3. The analysis of axial compressibility and thermal expansivity indicates that the c-axis is over three times more compressible (K_Tc = 47 ± 1 GPa) than the a-axis (K_Tc = 157 ± 1 GPa), whereas the thermal expansion of the c-axis (a₀ = 6.8 (2) × 10⁻⁵ K⁻¹ and a₁ = 2.2 (4) × 10⁻⁸ K⁻²) is greater than that of the a-axis (a₀ = 2.7 (4) × 10⁻⁵ K⁻¹ and a₁ = −0.2 (2) × 10⁻⁸ K⁻²). The present thermal EOS enables us to accurately calculate the density of magnesite to the deep mantle conditions. Decarbonation of a subducting oceanic crust containing 2 wt.% magnesite would result in a 0.6% density reduction at 30 GPa and 1273 K. Using the new EOS parameters we performed thermodynamic calculations for magnesite decarbonation reactions at pressures to 20 GPa. We also estimated stability of magnesite-bearing assemblages in the lower mantle.     

Progressive cooling of the hyaloclastite ridge at Gjálp, Iceland, 1996–2005

In the subglacial eruption at Gjálp in October 1996 a 6 km long and 500 m high subglacial hyaloclastite ridge was formed while large volumes of ice were melted by extremely fast heat transfer from magma to ice. Repeated surveying of ice surface geometry, measurement of inflow of ice, and a full Stokes 2-D ice flow model have been combined to estimate the heat output from Gjálp for the period 1996–2005. The very high heat output of order 10⁶ MW during the eruption was followed by rapid decline, dropping to  2500 MW by mid 1997. It remained similar until mid 1999 but declined to 700 MW in 1999–2001. Since 2001 heat output has been insignificant, probably of order 10 MW. The total heat carried with the 1.2 × 10¹² kg of basaltic andesite erupted (0.45 km³ DRE) is estimated to have been 1.5 × 10¹⁸ J. About two thirds of the thermal energy released from the 0.7 km³ edifice in Gjálp occurred during the 13-day long eruption, 20% was released from end of eruption until mid 1997, a further 10% in 1997–2001, and from mid 2001 to present, only a small fraction remained. The post-eruption heat output history can be reconciled with the gradual release of 5 × 10¹⁷ J thermal energy remaining in the Gjálp ridge after the eruption, assuming single phase liquid convection in the cooling edifice. The average temperature of the edifice is found to have been approximately 240 °C at the end of the eruption, dropping to  110 °C after 9 months and reaching  40 °C in 2001. Although an initial period of several months of very high permeability is possible, the most probable value of the permeability from 1997 onwards is of order 10^− 12 m². This is consistent with consolidated/palagonitized hyaloclastite but incompatible with unconsolidated tephra. This may indicate that palagonitization had advanced sufficiently in the first 1–2 years to form a consolidated hyaloclastite ridge, resistant to erosion. No ice flow traversing the Gjálp ridge has been observed, suggesting that it has effectively been shielded from glacial erosion in its first 10 years of existence.     

Studies on turbulent diffusion processes and evaluation of diffusivity values from hydrodynamic observations in Corpus Christi Bay

The physical process of dispersion which can be attributed to turbulence (turbulent diffusion) or shear (shear-augmented diffusion) within the flow field is very important as it ultimately governs the distribution of constituents of interest within the environment. A series of diffusion experiments were conducted in Corpus Christi Bay, TX with the purpose of characterizing turbulent diffusion through dispersion coefficients or turbulent diffusivity, K_i (i=x, y, z) dependent on the degree of randomness or turbulence intensity, I.Measured with a boat-mounted acoustic doppler current profiler (ADCP), the Eulerian velocity time-series of fluid particles in random motion, u_i was used in the evaluation of the Eulerian time-scale of turbulence, T^E based on the velocity correlation function, R^E with T^E being related to the Lagrangian time-scale T^L through a scaling parameter, β(=T^L/T^E). Surface currents were obtained with high frequency (HF) Radar equipment deployed over the study area from which the horizontal velocity gradients were determined.Within the spatial scale of the experiment (1000 m), the observed low horizontal gradients (10⁻⁴ s⁻¹) allowed for the generation of velocity time-series from an ADCP mounted on a moving platform. A numerical scheme for evaluating turbulent diffusivity values was developed on the basis of Eulerian current measurements and calibrated through the statistics of an evolving dye patch for the scaling parameter β which in this scheme was found to be in the range 1–3.     

The Polar Cap (PC) indices: Relations to solar wind parameters and global magnetic activity

The polar geomagnetic activity resulting from solar wind–magnetosphere interactions can be characterized the Polar Cap (PC) indices, PCN and PCS. PC index values are derived from polar magnetic variations calibrated on a statistical basis such that the index approximate values in units of mV/m of the interplanetary “geo-effective” (or “merging”) electric field (E_M) conveyed by the solar wind. The timing and amplitude relations of the PC index to solar wind plasma and magnetic field parameters are reported. The solar wind effects are parameterized in terms of the geo-effective electric field (E_M) and the dynamical pressure (P_DYN). The PC index has a delayed and damped response to E_M variations and display saturation-like effects for E_M values exceeding 10 mV/m. Steady or slowly varying levels of solar wind dynamical pressure have little or no impact on the PC index above the effects related to E_M for which the solar wind velocity is also a factor. Sharp increases in the dynamical pressure generate impulsive variations in the PC index comprising a initial negative impulse of 5–10 min duration followed by a positive impulse lasting 10–20 min. Typical amplitudes of both the negative and the positive impulses are 0.2–0.5 units. A sharp decrease in the pressure produces the inverse sequence of pulses in the PC index. Auroral substorm activity represented by the AL index level has a marked influence on the average PC/E_M level at the transition from very quiet (AL0 nT) to disturbed conditions while more or less disturbed conditions (AL<−100 nT) have no systematic effect on the average PC/E_M values. At distinct substorm events the PC/E_M ratio has a minimum (0.8) in the pre-onset phase at around 20 min before substorm onset. The average ratio gradually increases in the expansion phase to reach a maximum value (1.1) at around 40 min after substorm onset (or 20 min after the largest (negative) peak in AL). At substorm recovery during the next 2 h the PC/E_M ratio decreases. Finally, we report on the application of polar magnetic variations to model the disturbance storm time (Dst) index development during magnetic storms by using the PC index as a source function to quantify the energy input to the ring current representing accumulated storm energy and characterized by the Dst index.     

The emplacement of an obsidian dyke through thin ice: Hrafntinnuhryggur, Krafla Iceland

An eruption along a 2.5 km-long rhyolitic dyke at Krafla volcano, northern Iceland during the last glacial period formed a ridge of obsidian (Hrafntinnuhryggur). The ridge rises up to 80 m above the surrounding land and is composed of a number of small-volume lava bodies with minor fragmental material. The total volume is < 0.05 km³. The lava bodies are flow- or dome-like in morphology and many display columnar-jointed sides typical of magma–ice interaction, quench-fragmented lower margins indicative of interaction with meltwater and pumiceous upper surfaces typical of subaerial obsidian flows. The fragmental material compromises poorly-sorted perlitic quench hyaloclastites and poorly-exposed pumiceous tuffs. Lava bodies on the western ridge flanks are columnar jointed and extensively hydrothermally altered. At the southern end of the ridge the feeder dyke is exposed at an elevation  95 m beneath the ridge crest and flares upwards into a lava body.Using the distribution of lithofacies, we interpret that the eruption melted through ice only 35–55 m thick, which is likely to have been dominated by firn. Hrafntinnuhryggur is therefore the first documented example of a rhyolitic fissure eruption beneath thin ice/firn. The eruption breached the ice, leading to subaerial but ice/firn-contact lava effusion, and only minor explosive activity occurred. The ridge appears to have been well-drained during the eruption, aided by the high permeability of the thin ice/firn, which appears not to have greatly affected the eruption mechanisms. We estimate that the eruption lasted between 2 and 20 months and would not have generated a significant jökulhlaup (< 70 m³ s^− 1).     

Hydrothermal activity on the southern Mid-Atlantic Ridge: Tectonically- and volcanically-controlled venting at 4–5°S

We report results from an investigation of the geologic processes controlling hydrothermal activity along the previously-unstudied southern Mid-Atlantic Ridge (3–7°S). Our study employed the NOC (UK) deep-tow sidescan sonar instrument, TOBI, in concert with the WHOI (USA) autonomous underwater vehicle, ABE, to collect information concerning hydrothermal plume distributions in the water column co-registered with geologic investigations of the underlying seafloor. Two areas of high-temperature hydrothermal venting were identified. The first was situated in a non-transform discontinuity (NTD) between two adjacent second-order ridge-segments near 4°02′S, distant from any neovolcanic activity. This geologic setting is very similar to that of the ultramafic-hosted and tectonically-controlled Rainbow vent-site on the northern Mid-Atlantic Ridge. The second site was located at 4°48′S at the axial-summit centre of a second-order ridge-segment. There, high-temperature venting is hosted in an  18 km² area of young lava flows which in some cases are observed to have flowed over and engulfed pre-existing chemosynthetic vent-fauna. In both appearance and extent, these lava flows are directly reminiscent of those emplaced in Winter 2005−06 at the East Pacific Rise, 9°50′N and reference to global seismic catalogues reveals that a swarm of large (M 4.6−5.6) seismic events was centred on the 5°S segment over a  24 h period in late June 2002, perhaps indicating the precise timing of this volcanic eruptive episode. Temperature measurements at one of the vents found directly adjacent to the fresh lava flows at 5°S MAR (Turtle Pits) have subsequently revealed vent-fluids that are actively phase separating under conditions very close to the Critical Point for seawater, at  3000 m depth and 407 °C: the hottest vent-fluids yet reported from anywhere along the global ridge crest.     

Nature of sediment dispersal off the Sepik River, Papua New Guinea: preliminary sediment budget and implications for margin processes

A sediment budget is constructed for the slope and narrow continental shelf off the Sepik River in order to estimate the relative importance of turbid plumes versus bottom gravity transport through a near-shore submarine canyon in the dispersal of sediment across this collision margin. ²¹⁰Pb geochronology and inventories of Kasten cores are consistent with the northwestward dispersal of sediment from the river mouth via hypopycnal and possible isopycnal plumes. Sediment accumulation rates are 5 cm yr⁻¹ on the upper slope just off of the Sepik mouth, decreasing gradually to 1 cm yr⁻¹ toward the northwest, and decreasing abruptly offshore (<0.2 cm yr⁻¹ at 1200 m water depth). A sediment budget indicates that only about 7–15% of the Sepik River sediment discharge accumulates on the adjacent open shelf and slope. The remainder presumably escapes offshore via gravity flows through a submarine canyon, the head of which extends into the river mouth. The divergent sediment pathways observed off the Sepik River (i.e., surface and subsurface plumes versus sediment gravity flows through a canyon) may be common along high-yield collision margins of the Indo–Pacific archipelago, and perhaps are analogous to most margins during Late Quaternary low sea-level conditions.     

The Deformational Response of a Viscoelastic Solid Earth Model Coupled to a Thermomechanical Ice Sheet Model

We apply a coupled thermomechanical ice sheet—self-gravitating viscoelastic solid Earth model (SGVEM), allowing for the dynamic exchange of ice thickness and bedrock deformation, in order to investigate the effect of viscoelastic deformation on ice dynamics and vice versa. In a synthetic glaciation scenario, we investigate the interaction between the ice sheet and the solid Earth deformation, the glacial-isostatic adjustment (GIA), accounting for an atmospheric forcing depending on the ice sheet surface altitude. We compare the results from the coupled model to runs with the common elastic lithosphere/relaxing asthenosphere (ELRA) model, where the lithosphere is represented by a thin plate and the mantle relaxes with one characteristic relaxation time, as well as to a rigid Earth without any deformation. We find that the deformational behaviour of the SGVEM on ice dynamics (i.e. stored ice volume, ice thickness and velocity field) is comparable to the ELRA for an optimal choice of the parameters in steady state, but exhibits differences in the transient behaviour. Beyond the ice sheet, in the region of peripheral forebulge, the differences in the transient surface deformation between ELRA and SGVEM are substantial, demonstrating the inadequacy of the ELRA model for interpreting constraints on GIA in the periphery of the ice sheet, such as sea-level indicators and GPS uplift rates.