Coastal and global sea level change

Both coastal and global mean sea level rise by about 3.0 ± 0.5 mm/year from January 1993 to December 2004. Over shorter intervals the coastal sea level rises faster and over longer intervals slowly than the global mean, which trend is almost constant for each interval and is equal to 2.9 ± 0.5 mm/year in 1993–2008. The different trends are due to the higher interannual variability of coastal sea level, caused by the sea level regional variability, that is further averaged out when computing the global mean.Coastal sea level rise is well represented by a selected set of 267 stations of the Permanent Service for Mean Sea Level and by the corresponding co-located altimeter points. Its departure from coastal sea level computed from satellite altimetry in a 150 km distance from coast, dominated by a large rise in the Eastern Pacific, is due to the regional interannual variability.Regionally the trends of the coastal and open-ocean sea level variability are in good agreement and the main world basins have a positive averaged trend. The interannual variability is highly correlated with the El Nino Southern Oscillation (SO) and the North Atlantic Oscillation (NAO) climatic indices over both the altimeter period and the interval 1950–2001. Being the signal of large scale a small number of stations with good spatial coverage is needed. The reconstruction of the interannual variability using the spatial pattern from altimetry and the temporal patterns from tide gauges correlated to NAO and SOI restitutes about 50% of the observed interannual variability over 1993–2001.     

Sea level time series and ocean tide analysis from multipath signals at five GPS sites in different parts of the world

We present sea level observations derived from the analysis of signal-to-noise ratio (SNR) data recorded by five coastal GPS stations. These stations are located in different regions around the world, both in the northern and in the southern hemisphere, in different multipath environments, from rural coastal areas to busy harbors, and experience different tidal ranges.The recorded SNR data show periodic variations that originate from multipath, i.e. the interference of direct and reflected signals. The general assumption is that for satellite arcs facing the open sea, the rapid SNR variations are due to reflections off the sea surface. The SNR data recorded from these azimuth intervals were analyzed by spectral analysis with two methods: a standard analysis method assuming a static sea level during a satellite arc and an extended analysis method assuming a time dependent sea level during a satellite arc.The GPS-derived sea level results are compared to sea level records from co-located traditional tide gauges, both in the time and in the frequency domain. The sea level time series are highly correlated with correlation coefficients to the order of 0.89–0.99. The root-mean-square (RMS) difference is 6.2 cm for the station with the lowest tidal range of 165 cm and 43 cm for the station with the highest tidal range of 772 cm. The relative accuracy, defined as the ratio of RMS and tidal range, is between 2.4% and 10.0% for all stations.Comparing the standard analysis method and the extended analysis method, the results based on the extended analysis method agree better with the independent tide gauge records for the stations with a high tidal range. For the station with the highest tidal range (772 cm), the RMS is reduced by 47% when using the extended analysis method. Furthermore, the results also indicate that the standard analysis method, assuming a static sea level, can be used for stations with a tidal range of up to about 270 cm, without performing significantly worse than the extended analysis method.Tidal amplitudes and phases are derived by harmonic analysis of the sea level records. Again, a high level of agreement is observed between the tide gauge and the GPS-derived results. Comparing the GPS-derived results, the results based on the extended analysis method show a higher degree of agreement with the traditional tide gauge results for stations with larger tidal ranges. Spectral analysis of the residuals after the harmonic analysis reveals remaining signal power at multiples of the draconitic day. This indicates that the observed SNR data are to some level disturbed by additional multipath signals, in particular for GPS stations that are located in harbors.     

Geodetic observation of sea-level change and crustal deformation in the Baltic Sea region

Based on tide gauge observations spanning almost 200 years, homogeneous time series of the mean relative sea level were derived for nine sites at the southern coast of the Baltic Sea. Our regionally concentrated data were complemented by long-term relative sea-level records retrieved from the data base of the Permanent Service for Mean Sea Level (PSMSL). From these records relative sea-level change rates were derived at 51 tide gauge stations for the period between 1908 and 2007. A minimum observation time of 60 years is required for the determination of reliable sea-level rates. At present, no anthropogenic acceleration in sea-level rise is detected in the tide gauge observations in the southern Baltic. The spatial variation of the relative sea-level rates reflects the fingerprint of GIA-induced crustal uplift. Time series of extreme sea levels were also inferred from the tide gauge records. They were complemented by water level information from historic storm surge marks preserved along the German Baltic coast. Based on this combined dataset the incidence and spatial variation of extreme sea levels induced by storm surges were analysed yielding important information for hazard assessments. Permanent GPS observations were used to determine recent crustal deformation rates for 44 stations in the Baltic Sea region. The GPS derived height change rates were applied to reduce the relative sea-level changes observed by tide gauges yielding an estimate for the eustatic sea-level change. For 13 tide gauge-GPS colocation sites a mean eustatic sea-level trend of 1.3 mm/a was derived for the last 100 years.     

Water mass variation in the Mediterranean and Black Seas

The mass-induced sea level variability and the net mass transport between Mediterranean Sea and Black Sea are derived for the interval between August 2002 and July 2008 from satellite-based observations and from model data. We construct in each basin two time series representing the basin mean mass signal in terms of equivalent water height. The first series is obtained from steric-corrected altimetry while the other is deduced from GRACE data corrected for the contamination by continental hydrology. The series show a good agreement in terms of annual and inter-annual signals, which is in line with earlier works, although different model corrections influence the consistency in terms of seasonal signal and trend.In the Mediterranean Sea, we obtain the best agreement using a steric correction from the regional oceanographic model MFSTEP and a continental hydrological leakage correction derived from the global continental hydrological model WaterGAP2. The inter-annual time series show a correlation of 0.85 and a root mean square (RMS) difference of 15 mm. The two estimates have similar accuracy and their annual amplitude and phase agree within 3 mm and 23 days respectively. The GRACE-derived mass-induced sea level variability yields an annual amplitude of 27 ± 5 mm peaking in December and a trend of 5.3 ± 1.9 mm/yr, which deviates within 3 mm/yr from the altimetry-derived estimate.In the Black Sea, the series are less consistent, with lower accuracy of the GRACE-derived estimate, but still show a promising agreement considering the smaller size of the basin. The best agreement is realized choosing the corrections from WaterGAP2 and from the regional oceanographic model NEMO. The inter-annual time series have a correlation and RMS differences of 0.68 and 55 mm, their annual amplitude and phase agree within 4 mm and 6 days respectively. The GRACE-derived seawater mass signal has an annual amplitude of 32 ± 4 mm peaking in April. On inter-annual time scales, the mass-induced sea level variability is stronger than in the Mediterranean Sea, with an increase from 2003 to 2005 followed by a decrease from 2006 to 2008.Based on mass conservation, the mass-induced sea level variations, river runoff and precipitation minus evaporation are combined to derive the strait flows between the basins and with the Atlantic Ocean. At the Gibraltar strait, the net inflow varies annually with an amplitude of 52 ± 10 × 10⁻³ Sv peaking end of September (1 Sv = 10⁶ m³ s⁻¹). The inflow through the Bosphorus strait displays an annual amplitude of 13 ± 3 ×10⁻³ Sv peaking in the middle of March. Additionally, an increase of the Gibraltar net inflow (3.4 ± 0.8 × 10⁻³ Sv/yr) is detected.     

Suspended microplastics in the surface water of the Yangtze Estuary System,China: First observations on occurrence,distribution

Levels of microplastics (MPs) in China are completely unknown. This study characterizes suspended MPs quantitatively and qualitatively for the Yangtze Estuary and East China Sea. MPs were extracted via a floatation method. MPs were counted and categorized according to shape and size under a stereomicroscope. The MP densities were 4137.3 ± 2461.5 and 0.167 ± 0.138 n/m³, respectively, in the estuarine and the sea samples. Plastic abundances varied significantly in the estuary. Higher densities in three sea trawls confirmed that rivers were the important sources of MP to the marine environment. Plastic particles (>5 mm) were observed with a maximum size of 12.46 mm, but MPs (0.5–5 mm) constituted more than 90% by number of items. The most frequent geometries were fibres, followed by granules and films. Plastic spherules occurred sparsely. Transparent and coloured plastics comprised the majority of the particles. This study provides clues in understanding the fate and potential sources of MPs.     

In situ cosmogenic 10Be production-rate calibration from the Southern Alps,New Zealand

We present a ¹⁰Be production-rate calibration derived from an early Holocene debris-flow deposit at about 1000 m above sea level in the central Southern Alps, New Zealand, in the mid-latitude Southern Hemisphere. Ten radiocarbon ages on macrofossils from a soil horizon buried by the deposit date the deposit to 9690 ± 50 calendar years before AD2008. Surface ¹⁰Be concentrations of seven large boulders partially embedded in the stable surface of the deposit are tightly distributed, yielding a standard deviation of ～2%. Conversion of the ¹⁰Be measurements to sea level/high-latitude values using each of five standard scaling methods indicates ¹⁰Be production rates of 3.84 ± 0.08, 3.87 ± 0.08, 3.83 ± 0.08, 4.15 ± 0.09, and 3.74 ± 0.08 atoms g^?1 a^?1, relative to the ‘07KNSTD’ ¹⁰Be AMS standard, and including only the local time-integrated production-rate uncertainties. When including a sea level high-latitude scaling uncertainty the overall error is ～2.5% (1σ) for each rate. To test the regional applicability of this production-rate calibration, we measured ¹⁰Be concentrations in a set of nearby moraines deposited before 18 060 ± 200 years before AD2008. The ¹⁰Be ages are only consistent with minimum-limiting ¹⁴C age data when calculated using the new production rates. This also suggests that terrestrial in situ cosmogenic-nuclide production did not change significantly from Last Glacial Maximum to Holocene time in New Zealand. Our production rates agree well with those of a recent calibration study from northeastern North America, but are 12–14% lower than other commonly adopted values. The production-rate values presented here can be used elsewhere in New Zealand for rock surfaces exposed during or since the last glacial period.     

GPS constraints on the deformation of Azerbaijan and surrounding regions

In this study, we present new GPS observations in Azerbaijan to provide an improved basis for determining the distribution of crustal deformation throughout the country and surrounding areas. The deformation field in the region has been analyzed with a dense GPS network configuration and a reliable quantification of the ongoing deformation was achieved. Results show that while contraction is dominant over the whole region, it is mostly concentrated on the middle and eastern parts of Caucasus Thrust Fault reaching up to 6.4 ± 0.2 mm/yr and Lesser Caucasus Fault does no accommodate more than 1–2 mm/yr of contraction. New network also clearly substantiates that the West Caspian Fault, which is a continuation of Caucasus Thrust Fault in the south, accommodates right-lateral slip rates of 7.1 ± 0.3 mm/yr in addition to 5.5 ± 0.3 mm/yr contraction rates.     

OSL chronology of Quaternary terraced deposits outcropping between Mt. Etna volcano and the Catania Plain (Sicily,southern Italy)

In this paper we tested the applicability of the Optically Stimulated Luminescence (OSL) technique through Single-Aliquot Regenerative-dose (SAR) protocol, on single grain quartz extracted from alluvial–coastal sediments. Five samples were collected from deposits belonging to a flight of seven orders of coastal–alluvial terraces outcropping in the area between Mt. Etna volcano and the Catania Plain (Sicily, southern Italy), at the front of the Sicilian fold and thrust system. After various performance tests, we obtained OSL ages ranging between 240 ± 12 and 80 ± 4 ka, consistent with the normal evolutionary model of a terraced sequence, moving from the highest to the lowest elevation. Obtained data allowed us to determine a mean uplift rate of 1.2 mm/year during the last 240 ka, mostly related to regional uplift processes coupled with sea-level changes. Moreover, terraces belonging to the two highest orders are folded, forming a large anticline. According to our results, the frontal thrust of the Sicilian chain was active between 236 and 197 ka ago, even though seismological and geodetic data suggest current activity to the back.     

Determining sea level change in the German Bight

Regional mean sea level changes in the German Bight are considered. Index time series derived from 15 tide gauge records are analysed. Two different methods for constructing the index time series are used. The first method uses arithmetic means based on all available data for each time step. The second method uses empirical orthogonal functions. Both methods produce rather similar results for the time period 1924–2008. For this period, we estimate that regional mean sea level increased at rates between 1.64 and 1.74 mm/year with a 90% confidence range of 0.28 mm/year in each case. Before 1924, only data from a few tide gauges are available with the longest record in Cuxhaven ranging back till 1843. Data from these tide gauges, in particular from Cuxhaven, thus receive increasingly more weight when earlier years are considered. It is therefore analysed to what extent data from Cuxhaven are representative for the regional sea level changes in the German Bight. While this cannot be clarified before 1924, it is found that this is not the case from 1924 onwards when changes in Cuxhaven can be compared to that derived from a larger data set. Furthermore, decadal variability was found to be substantial with relatively high values towards the end of the analysis period. However, these values are not unusual when compared to earlier periods.     

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.     

Temporal changes in electrical resistivity at Sakurajima volcano from continuous magnetotelluric observations

Continuous magnetotelluric (MT) measurements were conducted from May 2008 to July 2009 at Sakurajima, one of the most active volcanoes in Japan. Two observation sites were established at locations 3.3 km east and 3 km west–northwest of the summit crater. At both observation sites, the high-quality component of the impedance tensor (Zyx) showed variations in apparent resistivity of approximately ± 20% and phase change of ± 2°, which continued for 20–180 days in the frequency range between 320 and 4 Hz. The start of the period of changes in apparent resistivity approximately coincided with the start of uplift in the direction of the summit crater, as observed by a tiltmeter, which is one of the most reliable pieces of equipment with which to detect magma intrusion beneath a volcano. A 2D inversion of MT impedance suggests that the resistivity change occurred at a depth around sea level. One of the possible implications of the present finding is that the degassed volatiles migrated not only vertically through the conduit but also laterally through a fracture network, mixing with shallow groundwater beneath sea level and thereby causing the observed resistivity change.     

Atmospheric polychlorinated biphenyl (pcb) inputs to a coastal city near the marmara sea

Recent studies have shown that polychlorinated biphenyl (PCB) levels are high in coastal sites. Various types of deposition events, including dry deposition, air–water exchange and wet deposition, were analyzed as part of the study to illustrate the pollution level. The atmospheric levels and deposition mechanisms of 82 PCBs were investigated in Mudanya, a coastal city in Turkey. A total of 175 samples and 112 blanks were collected. Air–water exchange and particle phase dry deposition fluxes were 120 ± 90 ng/m²-d and 125 ± 160 ng/m²-d, respectively. The wet-dry deposition sampler (WDDS) consisted of two reservoirs, wet and dry. Collected with the WDDS, the dry deposition flux was 21 ± 20 ng/m²-d, and the wet deposition fluxes (during the rainy period) were 145 ± 130 ng/m²-d and 195 ± 270 ng/m²-d for the dissolved and particulate phases, respectively. Mass transfer coefficients (MTCs), dry deposition velocities and washout ratios were calculated and compared with the literature values.     

Constraining the long-term evolution of the slip rate for a major extensional fault system in the central Aegean,Greece, using thermochronology

The brittle/ductile transition is a major rheologic boundary in the crust yet little is known about how or if rates of tectonic processes are influenced by this boundary. In this study we examine the slip history of the large-scale Naxos/Paros extensional fault system (NPEFS), Cyclades, Greece, by comparing published slip rates for the ductile crust with new thermochronological constraints on slip rates in the brittle regime. Based on apatite and zircon fission-track (AFT and ZFT) and (U–Th)/He dating we observe variable slip rates across the brittle/ductile transition on Naxos. ZFT and AFT ages range from 11.8 ± 0.8 to 9.7 ± 0.8 Ma and 11.2 ± 1.6 to 8.2 ± 1.2 Ma and (U–Th)/He zircon and apatite ages are between 10.4 ± 0.4 to 9.2 ± 0.3 Ma and 10.7 ± 1.0 to 8.9 ± 0.6 Ma, respectively. On Paros, ZFT and AFT ages range from 13.1 ± 1.4 Ma to 11.1 ± 1.0 Ma and 12.7 ± 2.8 Ma to 10.5 ± 2.0 Ma while the (U–Th)/He zircon ages are slightly younger between 8.3 ± 0.4 Ma and 9.8 ± 0.3 Ma. All ages consistently decrease northwards in the direction of hanging wall transport. Most of our new thermochronological results and associated thermal modeling more strongly support the scenario of an identical fault dip and a constant or slightly accelerating slip rate of ∼ 6–8 km Myr^− 1 on the NPEFS across the brittle/ductile transition. Even the intrusion of a large granodiorite body into the narrowing fault zone at ∼ 12 Ma on Naxos does not seem to have affected the thermal structure of the area in a way that would significantly disturb the slip rate. The data also show that the NPEFS accomplished a minimum total offset of ∼ 50 km between ∼ 16 and 8 Ma.     

Maternal-to-fetal transfer and concentration profiles of PCB congeners for Steller sea lions (Eumetopias jubatus) from Hokkaido,Japan

The concentrations of PCB congeners in the blubber and liver of mother Steller sea lions (Eumetopias jubatus; SSLs) and their fetuses from the coast of Hokkaido, Japan in 2008, 2010 and 2012 were analyzed by HRGC–HRMS, in order to elucidate PCB congener profiles and maternal-to-fetal transfer of PCBs in SSLs. ΣPCBs in the fetuses were 1400 ± 660 (the mean ± SD) ng/g-fat in the blubber and 570 ± 320 ng/g-fat in the liver, respectively. There was a concern that SSLs had been contaminated by PCBs during the fetal period. The concentrations of the major congeners in the blubber and liver were a correlation between the fetus and mother (blubber: r = 0.9934, liver: r = 0.9160; P ≦ 0.05). The results indicated that PCBs in the fetuses came from the mothers. #177 and #199 showed no correlations between the fetus and the mother in the blubber and liver. This indicated a selective capture by some natural protector such as the placenta.     

OSL-dating the Quaternary landscape evolution in the Vértes Hills forelands (Hungary)

Fluvial, colluvial, and aeolian sediments were dated by optically stimulated luminescence (OSL) on quartz to improve the chronological framework for Quaternary sedimentation and landscape evolution in the forelands of the Vértes Hills (central Hungary). The separated quartz was suitable for age determination based on an OSL SAR protocol. Most samples have asymmetric equivalent dose distributions and OSL ages were calculated by the mean, central, and minimum D_e values. Considering geomorphology and earlier age data from the area, the central D_e values seem most appropriate for age calculation. A fan on the geomorphological level QV in the western foreland of the Vértes Hills was deposited 79–75 (±8) ka ago. In the south-eastern foreland an alluvial fan on level QIIb is at most 42 ± 4 ka old. Fluvial incision and aggradation occurred 16–10 (±1) ka ago on the geomorphic surface QIIa. Loess is 14 ± 1 ka old, and slope sedimentation was active 11–9 (±1) ka ago. Our OSL data demonstrate that in the north-western foreland of the Vértes Hills wind remained an important agent after the last glacial times, into the early Holocene (9–8 ± 1 ka) and was able to accumulate large aeolian dunes.     

High resolution mapping of Earth tide response based on GPS data in Japan

We observe the Earth tidal fields at diurnal and semi-diurnal periods using Kinematic Precise Point Positioning (KPPP) GPS analysis. Our KPPP GPS solutions compare well with super-conducting gravimeter (SG) observations and a theoretical Earth tidal model, that includes both ocean tide loading model and body tides. We make a high resolution map of the observed Earth tidal response fields using the Japanese GEONET GPS network which consists of 1200 sites. We find that: (1) the average phase of GPS data lags 0.11 ± 0.04° from our theoretical Earth tidal model, (2) the average amplitude ratio between GPS and the theoretical Earth tidal model is 1.007 ± 0.003, (3) the amplitude in the Kyushu district is about 1.0–1.5 ± 0.3% larger than in the Hokkaido district, and (4) the amplitude at the Japan Sea side is about 0.5 ± 0.2% larger than that at the Pacific Ocean side. These results suggest that we may be able to place constraints on Earth structure using GPS-derived tidal information.     

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.     

Luminescence dating of baked earth and sediments from the Qujialing archaeological site,China

The Qujialing site is a representative Neolithic archaeological site in the middle reaches of the Yangtze River, China. Absence of suitable material for radiocarbon dating in this region makes the timing of the similar sites difficult. Here we applied optically stimulated luminescence (OSL–SAR) and thermoluminescence (TL–SAR) techniques to date the archaeological and natural deposits from the Qujialing site with known age, testing the techniques on samples at archaeological sites in this region. The results showed that the luminescence properties of quartz from sediment and baked earth samples are very similar. The quartz OSL ages obtained for a sediment sample and a baked earth sample from the cultural layer are 5.4 ± 0.3 and 5.1 ± 0.3 ka, respectively. The quartz TL age of the baked earth sample is 5.6 ± 0.5 ka. These dates are consistent with the calibrated radiocarbon ages (4.9 ± 0.1 and 5.1 ± 0.1 ka cal BP (±1σ)) of the two charcoal samples from the cultural layer at a nearby locality, and are also in agreement with the age of Qujialing culture period. The results indicate that the OSL dating techniques can be applied to date similar archaeological sites in the middle reaches of the Yangtze River, China.     

Crustal structure of the Reykjanes Ridge near 62°N,on the basis of seismic refraction and gravity data

Explosion deep seismic sounding data sections of high quality had been obtained with RV Meteor in the Reykjanes Iceland Seismic Project (RRISP77 [Angenheister, G., Gebrande, H., Miller, H., Goldflam, P., Weigel, W., Jacoby, W.R., Pálmason, G., Björnsson, S., Einarsson, P., Pavlenkova, N.I., Zverev, S., Litvinenko, I.V., Loncarecic, B., Solomon, S., 1980. Reykjanes Ridge Iceland Seismic Experiment (RRISP 77). J. Geophys. 47, 228–238]) which close an information gap near 62°N. Preliminary results were presented by Weigel [Weigel, W., 1980. Aufbau des Reykjanes Rückens nach refraktionsseismischen Messungen. In: Weigel, W. (Ed.), Reykjanes Rücken, Island, Norwegischer Kontinentalrand. Abschlusskolloquium, Hamburg zur Meteor-Expedition, vol. 45. DFG, Bonn, pp. 53–61], and here we report on the data and results of interpretation. Clear refracted phases to 90 km distance permit crustal and uppermost mantle structure to be modelled by ray tracing. The apparent P-wave velocities are around 4.5, 6–6.5, 7–7.6 and 8.2–8.7 km/s, but no wide-angle reflections have been clearly seen. Accompanying sparker reflection data reveal thin sediment ponds in the axial zone and up to 400 m thick sediments at 10 Ma crustal age. Ray tracing reveals the following model below the sediments: (1) a distinct, 1–2 km thick upper crust (layer 2A) with Vp increasing with age (to 10 Ma) from <3.4 to 4.9 km/s and with a vertical gradient of 0.1–0.2 km/s/km, (2) a lower crust or layer 3 beginning at depths of 2 (axis) to 4 km (10 Ma age) below sea level with 6.1–6.8 km/s and similar vertical gradients as above, (3) the lower crust bottoms at 5.2–9.5 km depth below sea level (0–10 Ma) with a marked discontinuity, underneath which (4) Vp rises from about 7.5–7.8 km/s (0–10 Ma) with a positive vertical gradient of, again, 0.1–0.2 km/s/km such that 8 km/s would be reached at 12 km and deeper near the axis. Our preferred interpretation is that the mantle begins at the distinct discontinuity (“Moho”), but a deeper “Moho” of Vp ≈ 8 km/s cannot be excluded. From Iceland southward to 60°N several experiments show a decrease of crustal thickness from 14 to 8 km. Velocity trends with age across the ridge reflect cooling and filling of cracks, and thickness trends probably suggest volcanic productivity variations as previously suggested.Gravity inversion concentrates on a profile across the ridge with the above seismic a priori information; with 0.2–0.5 km depth uncertainty it leads to a good fit (±2.5 mGal where seismic data exist). Best fitting densities are (in kg/m³) for sediments, 2180; upper crust, 2450–2570; lower crust, 2850–2940; mantle lithosphere, 3215–3240 with a deficit for an asthenospheric wedge of no more than −100 kg/m³. The morphological ridges and troughs superimposed on the SE ridge flank are partly correlated, partly anti-correlated with the Bouguer anomaly and suggest that variable crustal density variations accompany the morphology variations.