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 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.     

Rapid crustal uplift in Patagonia due to enhanced ice loss

A vertical crustal uplift rate of 39 mm yr^? 1 is measured between 2003 and 2006 using Global Positioning System (GPS) measurements at the northeastern edge of the Southern Patagonia Icefield (SPI). This is the largest present-day glacial isostatic rate ever recorded. The combination of SPI's rapid melting and the unique regional slab-window tectonics that promotes a relatively low viscosity, is central to our interpretation of the observations. The two effects lead to a strong interaction of short relaxation times with ice loads that change on a comparable time scale. The profile of GPS observations link ice loss to the soft viscoelastic isostatic flow response over the time scale of the Little Ice Age (LIA), including ice loss in the period of observation. The agreement of the results with our model predictions strongly suggests the large crustal uplift in Patagonia is due an accelerated glacier wasting since the termination of the LIA and that the effective regional mantle viscosity is near 4.0–8.0 × 10¹⁸ Pa s. A century-long diminution of the icefields, at rates that are about 1/4 – 1/2 the contemporary loss rates, is consistent with multidecadal-scale temperature trends estimated for the past 50–100 years and is, in fact, a key feature in any model capable of explaining the uplift data.     

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.     

Changes in ice phenology on polish lakes from 1961 to 2010 related to location and morphometry

Observations of lake ice at the shore, complete ice cover, ice duration, ice thickness and other measures for 18 Polish lakes were collected for the 50 year period (1961–2010). Average ice dates in early winter became later: first appearance of ice along shore 2.3 days decade⁻¹ and complete ice cover 1.2 days decade⁻¹ while complete ice cover disappeared earlier (5.6 days decade⁻¹) as did last ice at the shore (4.3 days decade⁻¹). The duration of ice cover decreased by 5.6 days decade⁻¹ and average ice thickness declined by 6.1 cm decade⁻¹. The magnitude of these values for individual lakes decreased from eastern to western Poland. This geographic gradient is likely related to regional atmospheric circulation because in winter this part of Europe is strongly affected by continental air, an influence that is greater in the east. A multivariate redundancy analysis (RDA), used in order to examine the dependence of ice measures on lake physical properties and location, indicated longitude and altitude as key factors explaining lake ice dynamics such as the disappearance of ice and ice cover, ice cover duration and thickness. Lake volume and average depth influenced mostly the appearance of ice and ice cover.     

Scattering coefficients in the mantle revealed from the seismogram envelope analysis based on the multiple isotropic scattering model

At longer periods, scattered ScS waves sometimes dominate over coda waves at large lapse times. Examining recordings of seismic envelopes at 9 IRIS seismic stations of regional earthquakes with focal depths deeper than 150 km in periods from 1 to 20 s for a wide lapse time range up to 2000 s, we found significant frequency dependence. The coda decay gradient at short periods is steeper than that at longer periods; however, the change of coda gradient associated with the ScS arrival becomes distinct as the period becomes longer. In particular, a clear offset of coda amplitude appears in central Asia for 10 and 15 s period bands. The multiple isotropic scattering process of S-waves in the heterogeneous mantle can be simply simulated by using the Monte Carlo simulation method based on the radiative transfer theory in scattering media. Assuming a two-plane-layer attenuation structure and smoothed velocity model of the PREM, we estimated the average total scattering coefficients of S-waves such as 7.52 × 10^− 4∼1.32 × 10^− 3 km^− 1 and 2.08 × 10^− 4∼6.23 × 10^− 4 km^− 1 at 4 s, and 4.51 × 10^− 4∼7.37 × 10^− 4 km^− 1 and 2.80 × 10^− 5∼2.71 × 10^− 4 km^− 1 at 10 s, for the lithosphere and the upper mantle and for the lower mantle, respectively. Our results indicate that scattering occurs mostly in the lithosphere and the upper mantle and support that medium heterogeneity spreads over the whole mantle though its scattering power is small. Strong scattering occurs beneath central Asia and Papua New Guinea, whereas the scattering beneath Italy and regions of east Russia is much weaker. The numerical calculation enables us to confirm that much stronger scattering than intrinsic attenuation causes the offset behavior with coda decay gradient change after the ScS arrival for 4 and 10 s period bands in some regions.     

Evapotranspiration in the Pampean Region using field measurements and satellite data

Evapotranspiration (LE) is an important factor for monitoring crops, water requirements, and water consumption at local and regional scale. In this paper, we applied the semi-empirical model to estimate the daily latent heat flux (LE_d = Rn_d + A − B(Ts − Ta)). LE_d has been estimated using satellite images (Thematic Mapper sensor) and a local dataset (incoming and outgoing short- and long-wave radiation) measured during three years. We first estimated the daily net Radiation (Rn_d) from a linear equation derived from the instantaneous net Radiation (Rn_d = CRn_i + D). Subsequently, coefficients A and B have been estimated for two different cover vegetations (pasture and soybean). For each vegetation cover, an error analysis combining Rn_d, A, B, and surface and air temperatures has been calculated. Results showed that Rn_d had good performance (nonbias and low RMSE). LE_d errors for pasture and soybean were ±28 W m⁻² and ±40 W m⁻² respectively.     

Volume and mass changes of the Greenland ice sheet inferred from ICESat and GRACE

This study examines the recent evolution of the Greenland ice sheet and its six major drainage basins. Based on laser altimetry data acquired by the Ice, Cloud and Land Elevation Satellite (ICESat), covering the period September–November 2003 to February–March 2008, ice surface height changes and their temporal variations were inferred. Our refined repeat track analysis is solely based on ICESat data and is independent of external elevation models, since it accounts for both ice height changes and the local topography. From the high resolution ice height change pattern we infer an overall mean surface height trend of −0.12 ± 0.006 m yr⁻¹. Furthermore, the largest changes could be identified at coastal margins of the ice sheet, exhibiting rates of more than −2 m yr⁻¹. The total ice volume change of the entire ice sheet amounts to −205.4 ± 10.6 km³ yr⁻¹. In addition, we assessed mass changes from 78 monthly Gravity Recovery and Climate Experiment (GRACE) solutions. The Release-04 gravity field solutions of GeoForschungsZentrum Potsdam cover the period between August 2002 and June 2009. We applied an adjusted regional integration approach in order to minimize the leakage effects. Attention was paid to an optimized filtering which reduces error effects from different sources. The overall error assessment accounts for GRACE errors as well as for errors due to imperfect model reductions. In particular, errors caused by uncertainties in the glacial isostatic adjustment models could be identified as the largest source of errors. Finally, we determined both seasonal and long-term mass change rates. The latter amounts to an overall ice mass change of −191.2 ± 20.9 Gt yr⁻¹ corresponding to 0.53 ± 0.06 mm yr⁻¹ equivalent eustatic sea level rise. From the combination of the volume and mass change estimates we determined a mean density of the lost mass to be 930 ± 11 kg m⁻³. This value supports our applied density assumption 900 ± 30 kg m⁻³ which was used to perform the volume–mass-conversion of our ICESat results. Hence, mass change estimates from two independent observation techniques were inferred and are generally in good agreement.     

The thermal state and strength of the lithosphere in the Spanish Central System and Tajo Basin from crustal heat production and thermal isostasy

In this work we have modeled the thermal structure of the lithosphere of the Spanish Central System and the Tajo Basin, and their implications for lithospheric strength. For this, we have used refined heat-producing elements (HPE) values to obtain new estimates of heat production rates in the Spanish Central System and Tajo Basin areas, which have been used joined to the relation between topography and thermal structure of the lithosphere to calculate the best-fit surface heat flows in the study area. Moreover, we have implemented a temperature-dependent thermal conductivity (appropriate for olivine) for the lithospheric mantle to improve the calculations of temperature profiles in the mantle. The geotherms so obtained, together with the implementation of a new rheological law for the upper lithospheric mantle, have been used to calculate refined estimations of the strength and effective elastic thickness of the lithosphere. We have obtained surface heat flow values of 84 mW m⁻² and ∼82 mW m⁻² for the Spanish Central System and the Tajo Basin, respectively. The thermal state of the lithosphere affects mantle temperatures, and hence may be playing an important role in the uplift and maintenance of the Spanish Central System.     

The replacement of native plants by exotic species may affect the colonization and reproduction of aquatic hyphomycetes

The tropical riparian zone has a high diversity of plant species that produce a wide variety of chemical compounds, which may be released into streams. However, in recent decades there has been an extensive replacement of tropical native vegetation by Eucalyptus monocultures. Our objective was to compare fungal colonization of Eucalyptus camaldulensis leaves with fungal colonization of native plant species from riparian zones in Brazilian Cerrado (savannah) streams. The fungal colonization and enzymatic activity significantly influenced leaf litter decomposition. Fungal sporulation rates from leaf litter varied significantly with leaf species, with E. camaldulensis showing the highest sporulation rate (1226 conidia mg⁻¹AFDM day⁻¹) and leaf mass loss (23.2 ± 0.9%). This species has the lowest lignin content and highest N concentration among the studied species. Among the studied native species, we observed the highest sporulation rate for Protium spruceanum (271 conidia mg⁻¹AFDM day⁻¹), Maprounea guianensis (268 conidia mg⁻¹AFDM day⁻¹) and Copaifera langsdorffii (196 conidia mg⁻¹AFDM day⁻¹). Overall, native plant species of the Brazilian Cerrado exhibited recalcitrant characteristics and a higher lignin:N ratio. Therefore, variations in the physical and chemical characteristics of the leaf litter could explain the higher decay rate and reproductive activity observed for E. camaldulensis. However, the detritus of this species were colonized almost exclusively by Anguillospora filiformis (99.6 ± 0.4%) and exhibited a reduction in aquatic hyphomycetes species diversity. Our results suggest that the disturbance in the composition of riparian vegetation and consequently, in the diversity of leaf litter input into streams, could change the patterns and rates of leaf litter utilization by microbial decomposers. These changes may have important consequences in the processing of organic matter and, consequently, in the functioning of freshwater ecosystems.     

Resolving sea level contributions by identifying fingerprints in time-variable gravity and altimetry

We have studied the ability of the GRACE gravimetry mission and Jason-1 altimetry to resolve ice and glacier induced contributions to sea level rise, by means of a fingerprint method. Here, the signals from ice sheet and land glacier changes, steric changes, glacial isostatic adjustment and terrestrial hydrology are assumed to have fixed spatial patterns. In a joint inversion using GRACE and Jason-1 data the unknown temporal components can then be estimated by least-squares. In total, we estimate temporal components for up to ∼ 80 individual patterns. From a propagation of the full error-covariance from GRACE and a diagonal error-covariance from Jason-1 altimetry we find that: (1) GRACE almost entirely explains the mass related parameters in the joint inversion, (2) an inversion using only Jason-1 data has a marginal ability to estimate the mass related parameters, while the steric parameters have much better formal accuracy. In terms of mean sea level rise the steric patterns have a maximum formal accuracy of 0.01 mm for an 11 week running mean. In general, strong negative error correlations (ρ < − 0.9) exists between the high and low elevation parts of the ice sheet drainage basins, when those are estimated independently. The largest formal errors found are in the order of 40 Gton for small high elevation subbasins in the southern Greenland ice sheet, which are difficult to separate. In a simplified joint inversion, merging high and low elevation basins, we have investigated the ability of the GRACE and Jason-1 data to separate the geocenter motion into a present-day contribution and a contribution from glacial isostatic adjustment (GIA). We find that the GIA related signal is larger than the present-day component with a maximum of −0.71 mm/year in the Z direction. Total geocenter motion rates are found to be −0.28, 0.43, −1.08 mm/year for the X, Y and Z components, respectively. The inversion results have been propagated to the Jason-1 along-track measurements. Over the time period considered, we see that a large part of the variability in the Pacific, Atlantic and Indian ocean can be explained by our inversion results. The applied inversion method therefore seems a feasible way to separate steric from mass induced sea level changes. At the same time, the joint inversion would benefit from more advanced parameterizations, which may aid in fitting remaining signal from altimetry.     

Stability of massive ground ice bodies in University Valley,McMurdo Dry Valleys of Antarctica: Using stable O–H isotope as tracers of sublimation in hyper-arid regions

To date, studies of the stability of subsurface ice in the McMurdo Dry Valleys of Antarctica have been mainly based on climate-based vapor diffusion models. In University Valley (1800 m), a small glacier is found at the base of the head of the valley, and adjacent to the glacier, a buried body of massive ice was uncovered beneath 20–40 cm of loose cryotic sediments and sandstone boulders. This study assesses the origin and stability of the buried body of massive ice by measuring the geochemistry and stable O–H isotope composition of the ice and applies a sublimation and molecular diffusion model that accounts for the observed trends. The results indicate that the buried massive ice body represents an extension of the adjacent glacier that was buried by a rock avalanche during a cold climate period. The contrasting δ¹⁸O profiles and regression slope values between the uppermost 6 cm of the buried massive ice (upward convex δ¹⁸O profile and S_D-18O = 5.1) and that below it (progressive increase in δ¹⁸O and S_D-18O = 6.4) suggest independent post-depositional processes affected the isotope composition of the ice. The upward convex δ¹⁸O profile in the uppermost 6 cm is consistent with the ice undergoing sublimation. Using a sublimation and molecular diffusion model, and assuming that diffusion occurred through solid ice, the sublimation rate needed to fit the measured δ¹⁸O profile is 0.2 ? 10^? 3 mm yr^? 1, a value that is more similar to net ice removal rates derived from ³He data from cobbles in Beacon Valley till (7.0 ? 10^? 3 mm yr^? 1) than sublimation rates computed based on current climate (0.1–0.2 mm yr^?1). We suggest that the climate-based sublimation rates are offset due to potential ice recharge mechanisms or to missing parameters, particularly the nature and thermo-physical properties of the overlying sediments (i.e., temperature, humidity, pore structure and ice content, grain size).     

Expected temporal absolute gravity change across the Taiwanese Orogen,a modeling approach

The island of Taiwan is located on the convergent boundary between the Philippine Sea plate and the Chinese continental margin. It offers very active mountain building and collapsing processes well illustrated by the rugged topography, rapid uplift and denudation, young tectonic landforms, active faulting and numerous earthquakes. In this paper, using simple models, we have estimated vertical movements and associated absolute gravity variations which can be expected along a profile crossing the southern part of the island and probably suffering the highest rates of rising. The two different tectonic styles proposed for the island, thin-skinned and thick-skinned, were taken into account. Horizontal and vertical movements were modeled by an elastic deformation code. Gravity variations due to these deformations are then modeled at a second step. They are dominated by plate and free-air effects, i.e. elevation of the topography, with several μGal yr⁻¹. By comparison, gravity changes generated by mass transfers are weak: maximum 0.1 μGal yr⁻¹ with the thin-skinned tectonic and 0.3 μGal yr⁻¹ with the thick-skinned tectonic. Though elastic rheology has limitations, this modeling offers interesting results on what gravity signal can be expected from the AGTO project (Absolute Gravity in the Taiwanese Orogen), which proposes to study the dynamic of these mountain ranges using absolute gravimetry (AG) and also including relative gravimetry (RG) and GPS measurements.     

Mercury levels in selected bycatch fish species from industrial shrimp-trawl fishery in the SE Gulf of California

Baseline Hg concentration in bycatch fish from the SE Gulf of California were determined in muscle and liver of 19 species. Levels of Hg in muscle were compared with legal limits of this element in national and international legislation. Considering all fish species, mean concentrations in liver (2.458 ± 1.997 μg g⁻¹) were significantly higher (p < 0.05) than in muscle (0.993 ± 0.670 μg g⁻¹). The sequence of averaged Hg concentrations in most ichthyofauna was liver > muscle. Highest level of Hg in muscle (2.556 μg g⁻¹) and liver (7.515 μg g⁻¹) corresponded to Diapterus peruvianus and Ophioscion strabo, respectively. Considering muscle samples, none of the species had levels of Hg above the limit (1.0 μg g⁻¹ wet weight) in the Mexican legislation; with respect to the Japanese (0.4 μg g⁻¹ wet weight) and British (0.3 μg g⁻¹ wet weight) legislations, 26.3% and 31.6% of the species respectively, were above the corresponding limits.     

Natural antifouling compounds produced by a novel fungus Aureobasidium pullulans HN isolated from marine biofilm

A fungus, Aureobasidium pullulans, was isolated from marine biofilm and identified. A bioassay-guided fractionation procedure was developed to isolate and purify antifouling compounds from A. pullulans HN. The procedure was: fermentation broth—aeration and addition of sodium thiosulfate–graduated pH and liquid–liquid extraction—SPE purification—GC–MS analysis. Firstly, the fermentation broth was tested for its toxicity. Then it was treated with aeration and addition of sodium thiosulfate, and its toxicity was almost not changed. Lastly, antifouling compounds were extracted at different pH, the extract had high toxicity at pH 2 but almost no toxicity at pH 10, which suggested the toxicants should be fatty acids. The EC₅₀ of the extract against Skeletonema costatum was 90.9 μg ml⁻¹, and its LC₅₀ against Balanus amphitrete larvae was 22.2 μg ml⁻¹. After purified by HLB SPE column, the EC₅₀ of the extract against S. costatum was 49.4 μg ml⁻¹. The myristic and palmitic acids were found as the main toxicants by GC–MS.     

Comparative study of top soil magnetic susceptibility variation based on some human activities

An investigation of the effect of some human activities on the magnetic susceptibility and frequency dependent susceptibility was conducted on top soil samples from, a commercial area, a motor park and a school environment in Jalingo, Taraba State, N-E Nigeria. The purpose was to assess the variation of magnetic susceptibility with different land use, detect pollution hotspots using magnetic proxy parameters and evaluate the contribution of superparamagnetic (SP) grain size contribution to the magnetic susceptibility from calculation of the frequency dependence of magnetic susceptibility (MS). The results of the mass specific low frequency magnetic susceptibility measurements showed significant enhancement with values ranging from 67.8 − 495.3 × 10⁻⁸ m³kg⁻¹ with a mean value of 191.61 × 10⁻⁸ × m³kg⁻¹ for the Jalingo College of Education (JCOE) data; 520.1 − 1612.8 × 10⁻⁸ m³kg⁻¹ with a mean value of 901.34 × 10⁻⁸ m³kg⁻¹ for the Jalingo main Market (JMM) and 188.5 − 1203.6 × 10⁻⁸m³kg⁻¹ with an average value of 574 92 × 10⁻⁶ m³kg⁻¹ for the Jalingo Motor Park (JMP). The significant magnetic enhancement indicates high concentration of ferrimagnetic minerals in the soil and hence increased pollution. The magnetic susceptibility of the different land use studied decreased in the order commercial area (market) > motor park > school premises. The results of the percentage frequency dependence susceptibility showed that most of the samples had a mixture of SP and coarse multi domain grains or SP grains < 0.05 µm. The value of χ_fd% range from 2.68 to 13.80% with an average value of 8.67% in the JCOE samples, 0.49 to 10.04% with an a-verage of 5.05% in the JMM samples and 0.56 to 13.04% with an average value of 5.86% in the JMP samples.     

Helium isotopes in the Izu Peninsula,Japan: Relation of magma and crustal activity

The relation of magma and crustal activity has been studied from spatial distribution of ³He/⁴He ratios of gas and/or water samples over the Izu Peninsula, where significant crustal deformation associated with seismic swarm activities has been observed since 1970s. The air-corrected values of ³He/⁴He ratios ranged from 3.5 to 8.2 R_A, where R_A is the atmospheric ³He/⁴He ratio = 1.4 × 10^? 6, indicating that helium is mostly of magmatic origin. Among the three pressure sources proposed to explain the crustal deformation, two inflation sources beneath the inland of northeast and the mid east coast of the Izu Peninsula locate in the broad distribution of high ³He/⁴He ratios, which supports relation of magma to the crustal uplift. In contrast, the distribution of ³He/⁴He ratios around the tensile fault assumed in the area of seismic swarms appears not to indicate existence of significant amount of magma below the tensile fault. Alternatively, the results suggest magma below a point several kilometers south of the tensile fault. The seismic swarms are explained either by fluid pressurization of thermal water heated by this magma or by intrusion of magma to the tensile fault moved obliquely from the deep magma reservoir.     

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.