Present day geodynamics in Iceland monitored by a permanent network of continuous GPS stations

Iceland is located on the Mid-Atlantic Ridge and thereby offers a rare opportunity to study crustal movements at a divergent plate boundary. Iceland is not only characterized by the divergence of the Eurasian and North American Plates, as several active volcanoes are located on the island. Moderate size earthquakes occur in the transform zones, causing measurable crustal deformation. In 1999 the installation of a permanent network of continuous GPS stations (ISGPS) was initiated in order to observe deformation due to unrest in the Hengill volcanic system and at the Katla volcano. The ISGPS network has been enlarged over the years and consists today of more than 25 CGPS stations. Most of the stations are located along the plate boundary, where most of the active deformation takes place. Uplift due to post-glacial rebound due to the melting of the largest glacier in Europe, Vatnajökull, is also detected by the ISGPS network. This study presents results from analysis of 9 years of data from the ISGPS network, in the global reference frame PDR05, which has been evaluated by the Potsdam-Dresden-Reprocessing group with reprocessed GPS data only. We thus determine subsidence or land uplift in a global frame. The horizontal station velocities clearly show spreading across the plate boundary of about 20 mm/a. Stations in the vicinity of the glacier Vatnajökull indicate uplift in the range of 12 mm/a, while a station in the central part of Iceland shows uplift rates of about 25 mm/a. Tide gauge readings in Reykjavik and current subsidence rates observed with CGPS agree also quite well.     

The GGOS as the backbone for global observing and local monitoring: A user driven perspective

A key geodetic contribution to both the three Global Observing Systems and initiatives like the European Global Monitoring for Environment and Security is an accurate, long-term stable, and easily accessible reference frame as the backbone. Many emerging scientific as well as non-scientific high-accuracy applications require access to an unique, technique-independent reference frame decontaminated for short-term fluctuations due to global Earth system processes. Such a reference frame can only be maintained and made available through an observing system such as the Global Geodetic Observing System (GGOS), which is currently implemented and expected to provide sufficient information on changes in the Earth figure, its rotation and its gravity field. Based on a number of examples from monitoring of infrastructure, point positioning, maintenance of national references frames to global changes studies, likely future accuracy requirements for a global terrestrial reference frame are set up as function of time scales. Expected accuracy requirements for a large range of high-accuracy applications are less than 5 mm for diurnal and sub-diurnal time scales, 2–3 mm on monthly to seasonal time scales, better than 1 mm/year on decadal to 50 years time scales. Based on these requirements, specifications for a geodetic observing system meeting the accuracy requirements can be derived.     

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 rise in the north-western part of the Arabian Gulf

Relative sea level variations in the north-western part of the Arabian Gulf have been estimated in the past using no more than 10 to 15 years of observations. In this study, we have almost doubled the period to 28.7 years by examining all available tide gauge data in the area and constructing a mean gauge time-series from seven coastal tide gauges. We found for the period 1979–2007 a relative sea level rise of 2.2 ± 0.5 mm/year. Using the subsidence observed at 6 GPS stations within a radius of 100 km of the tide gauges as an indication of the vertical land motion, the corresponding absolute sea level rise is 1.5 ± 0.8 mm/year that is in agreement with the global estimate of 1.9 ± 0.1 mm/year (Church and White, 2011) for the same studied period. By taking into account the temporal correlations we conclude that previous published results underestimate the true sea level rate uncertainty in this area by a factor of 5–10.     

Geological,tomographic, kinematic and geodynamic constraints on the dynamics of sinking slabs

We use geodynamic models with imposed plate velocities to test the forward-modeled history of subduction based on a particular plate motion model against alternative seismic tomography models. We utilize three alternative published reference frames: a hybrid moving hotspot-palaeomagnetic, a hybrid moving hotspot-true polar wander corrected-palaeomagnetic, and a Subduction Reference Frame, a plate model including longitudinal shifts of subduction zones by matching subduction volumes imaged by P-wave tomography, to assess which model best predicts present day mantle structure compared with seismic tomography and volumetrically derived subduction history. Geodynamic modeling suggests paleo-longitudinal corrections applied to the Subduction Reference Frame result in lower mantle slab material beneath North America and East Asia accumulating up to 10–15° westward of that imaged by tomography, whereas the hybrid models develop material offset by 2–9°. However, the Subduction Reference Frame geodynamic model produces slab material beneath the Tethyan Domain coinciding with slab volumes imaged by tomography, whereas the hybrid reference frame models do not, suggesting regional paleo-longitudinal corrections are required to constrain slab locations. We use our models to test inferred slab sinking rates in the mantle focusing on well-constrained regions. We derive a globally averaged slab-sinking rate of 13 ± 3 mm/yr by combining the ages of onset and cessation of subduction from geological data and kinematic reconstructions with images of subducted slabs in the mantle. Our global average slab-sinking rate overlaps with the 15–20 mm/yr rate implied by mantle convection models using a lower mantle viscosity 100 times higher than the upper mantle.     

Kinematics of the eastern part of the North Anatolian Fault Zone

The North Anatolian Fault Zone (NAFZ), which marks the boundary between Anatolia and the Eurasian plate, is one of the world's most seismically active structures. Although the eastern part of NAFZ has high seismic hazard, there is a lack of geodetic information about the present tectonics of this region. Even though many scientists would like to study this area, geographical and logistical problems make performing scientific research difficult. In order to investigate contemporary neotectonic deformation on the eastern NAFZ and in its neighborhood, a relatively dense Global Positioning System (GPS) monitoring network was established in 2003. Geodetic observations were performed in three GPS campaigns in an area of 350 km × 200 km with 12-month intervals. In addition, 14 new GPS stations were measured far from the deforming area. Since this region includes the intersection of the NAFZ and the East Anatolian Fault Zone (EAFZ), deformation is complex and estimating seismic hazard is difficult. One important segment is the Yedisu segment and it has not broken since the 1784 earthquake. After the 1992 Erzincan and 2003 Pulumur earthquakes, the Coulomb stress loading on the Yedisu segment of the NAFZ has increased significantly, emphasizing the need to monitor this region. We computed the horizontal velocity field with respect to Eurasia and strain rates field as well. GPS-derived velocities relative to Eurasia are in the range of 16–24 mm/year, which are consistent with the regional tectonics. The principal strain rates were derived from the velocity field. Results show that strain is accumulating between the NAFZ and EAFZ along small secondary fault branches such as the Ovacik Fault (OF).     

Intra-plate deformation in west-central Europe

The central part of Europe north of the alpine orogenic belt is generally seen as a relatively stable area of the western tip of the Eurasian plate. Indeed, up to now, no geodetically significant motions have been detected although an active rift system running roughly in SSE–NNW direction along the Rhine valley could have some effect on the stability of this region. Presently, the increasing accuracy of geodetic point motions should allow the study of small motions at levels down to nearly 0.1 mm/yr. We start our investigation with a closer look at the ‘true’ accuracy and significance of GPS derived point velocities of permanent stations. We compare and discuss the different levels of formal errors obtained by the three analysis centers considered in this study (EPN, JPL and SOPAC) and present additional ways of assessing the accuracy using the redundancy offered by different independent analyses and multiple systems operating at one site. On the average, all results indicate that a one-sigma level of ±0.3 mm/yr can be seen as a conservative estimate for the horizontal accuracy of point motions in central Europe. On the basis of this assumption we find that at present, the actual velocity field as determined by different analysis groups and centers does not show any significant east–west extensional deformation. We do however see a prominent north–south compressional velocity gradient of about 1 mm/yr/1000 km (1 nanostrain/yr) which could be associated with the Alpine thrust in conjunction with a south-directed horizontal component of the Fennoscandian Glacial Isostatic Adjustment.     

New insights in the geodynamics of the Lipari–Vulcano area (Aeolian Archipelago,southern Italy) from geological,geodetic and seismological data

Geological, geodetic and seismological data have been analyzed in order to frame the Lipari–Vulcano complex (Aeolian archipelago, southern Italy) into the geodynamic context of the southeastern Tyrrhenian Sea. It is located at the northern end of a major NNW–SSE trending right-lateral strike-slip fault system named “Aeolian–Tindari–Letojanni” which has been interpreted as a lithospheric discontinuity extending from the Aeolian Islands to the Ionian coast of Sicily and separating two different tectonic domains: a contractional one to the west and an extensional one to the north-east. Structural field data consist of structural measurements performed on well-exposed fault planes and fractures. The mesostructures are mostly represented by NW–SE striking normal faults with a dextral-oblique component of motion. Minor structures are represented by N–S oriented joints and tension gashes widespread over the whole analyzed area and particularly along fumarolized sectors. The analyzed seismological dataset (from 1994 to 2013) is based on earthquakes with magnitude ranging between 1.0 and 4.8. The hypocenter distribution depicts two major alignments corresponding to the NNW–SSE trending Aeolian–Tindari–Letojanni fault system and to the WNW–ESE oriented Sisifo–Alicudi fault system. GPS data analysis displays ∼3.0 mm/yr of active shortening between the two islands, with a maximum shortening rate of about 1.0 × 10⁻¹³ s⁻¹, between La Fossa Caldera and south of Vulcanello. This region is bounded to the north by an area where the maximum values of shear strain rates, of about 0.7 × 10⁻¹³ s⁻¹ are observed. This major change occurs in the area south of Vulcanello that is also characterized by a transition in the way of the vertical axis rotation. Moreover, both the islands show a clear subsidence process, as suggested by negative vertical velocities of all GPS stations which exhibit a decrease from about −15 to −7 mm/yr from north to south. New data suggest that the current kinematics of the Lipari–Vulcano complex can be framed in the tectonic context of the eastward migrating Sisifo–Alicudi fault system. This is dominated by transpressive tectonics in which contractional and minor extensional structures can coexist with strike-slip motion.     

Sea level: A review of present-day and recent-past changes and variability

In this review article, we summarize observations of sea level variations, globally and regionally, during the 20th century and the last 2 decades. Over these periods, the global mean sea level rose at rates of 1.7 mm/yr and 3.2 mm/yr respectively, as a result of both increase of ocean thermal expansion and land ice loss. The regional sea level variations, however, have been dominated by the thermal expansion factor over the last decades even though other factors like ocean salinity or the solid Earth's response to the last deglaciation can have played a role. We also present examples of total local sea level variations that include the global mean rise, the regional variability and vertical crustal motions, focusing on the tropical Pacific islands. Finally we address the future evolution of the global mean sea level under on-going warming climate and the associated regional variability. Expected impacts of future sea level rise are briefly presented.     

Global-scale validation of model-based load deformation of the Earth's crust from continental watermass and atmospheric pressure variations using GPS

Temporal mass variations in the continental hydrosphere and in the atmosphere lead to changes in the gravitational potential field that are associated with load-induced deformation of the Earth’s crust. Therefore, models that compute continental water storage and atmospheric pressure can be validated by measured load deformation time series. In this study, water mass variations as computed by the WaterGAP Global Hydrology Model (WGHM) and surface pressure as provided by the reanalysis product of the National Centers for Environmental Prediction describe the hydrological and atmospheric pressure loading, respectively. GPS observations from 14 years at 208 stations world-wide were reprocessed to estimate admittance factors for the associated load deformation time series in order to determine how well the model-based deformation fits to real data. We found that such site-specific scaling factors can be identified separately for water mass and air pressure loading. Regarding water storage variation as computed by WGHM, weighted global mean admittances are 0.74 ± 0.09, 0.66 ± 0.10, 0.90 ± 0.06 for the north, east and vertical component, respectively. For the dominant vertical component, there is a rather good fit to the observed displacements, and, averaged over all sites, WGHM is found to slightly overestimate temporal variations of water storage. For Europe and North America, with a dense GPS network, site-specific admittances show a good spatial coherence. Regarding regional over- or underestimation of WGHM water storage variations, they agree well with GRACE gravity field data. Globally averaged admittance estimates of pre-computed atmospheric loading displacements provided by the Goddard Geodetic VLBI Group were determined to be 0.88 ± 0.04, 0.97 ± 0.08, 1.13 ± 0.01 for the north, east and vertical, respectively. Here, a relatively large discrepancy for the dominant vertical component indicates an underestimation of corresponding loading predictions.     

Velocity field and tectonic strain in Southern Spain and surrounding areas derived from GPS episodic measurements

The goal of this paper is to study the velocity field and deformation parameters in Southern Spain and surrounding areas (Ibero-Maghrebian region) using GPS episodic measurements. Results are compared to those previously published as well as deformation parameters derived from seismic data. For this purpose, a geodetic GPS network of 12 stations was observed during eight field campaigns from 1998 to 2005 by the San Fernando Naval Observatory (ROA), Spain. Relative GPS velocities in the Gulf of Cadiz with respect to the stable part of Eurasia are ～4.1 mm/yr in a NW–SE to NNW–SSE direction. In the Betics, Alboran Sea and North of Morocco, velocities are ～4.4 mm/yr in a NW–SE direction, and they are ～2.3 mm/yr in a N–S direction in the eastern part of the Iberian Peninsula. These results are in agreement with the anticlockwise rotation of the African plate. GPS strain tensors are determined from the velocity model, to obtain a more realistic crustal deformation model. The Gulf of Cadiz is subjected to uniform horizontal compression in a NNW–SSE direction, with a rotation to N–S in the Alboran Sea and Northern Morocco. An extensional regime in a NW–SE direction, which rotates to W–E, is present in the Internal Betics area. In the Betic, Alboran Sea and North of Morocco regions we compare seismic deformation rates from shallow earthquakes with the determined GPS deformation rates. The comparison indicates a seismic coupling of 27%, while the remaining 73% might be generated in aseismic processes. Deformations measured in the Ibero-Maghrebian region with GPS could be interpreted in terms of either elastic loading or ductile deformation.     

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.     

Imaging Pn velocities beneath the Pannonian basin

The Pannonian depression is an extensional back-arc basin in central Europe and is an integral part of the Alpine–Carpathian orogenic mountain belts. It can be characterized by thinned lower crust, shallow Moho discontinuity, high surface heat flow and Moho temperature, implying recent active tectonic processes. Imaging the velocity structure of the upper mantle may help us to better understand the structure and formation of the Pannonian region.In this paper, P_n traveltimes from regional earthquakes are used to tomographically image the lateral velocity variations in the uppermost mantle beneath the Pannonian basin. The set of linear tomographic equations, built up of the time term equation for each source–receiver pair, is solved by a truncated singular value decomposition algorithm. The explicit computation of the generalized inverse of the tomographic equations makes it possible to deduce both the resolution matrix and the model covariance matrix, allowing us to estimate the resolution and reliability of the solution.The mean compressional wave velocity in the uppermost mantle beneath the Pannonian basin is 7.9 km/s, substantially lower than the average continental P_n velocity of 8.1 km/s. It is mostly due to the high Moho temperature having values on average 400–500 °C more than those in the surrounding areas. The velocity anomalies range from −0.3 to 0.3 km/s relative to the mean velocity of 7.9 km/s. Due to high Moho temperature, below the North Hungarian range low (7.6–7.7 km/s) velocities can be found. High-velocity anomalies of around 8.1 km/s can be detected along the W-SW boundaries of Hungary and at the junction of the Pannonian basin and the Southern Carpathians. The Great Hungarian Plain shows average (7.9 km/s) P_n velocities.     

Crustal deformation before the 2008 Wenchuan MS8.0 earthquake studied using GPS data

We used GPS velocities from approximately 700 stations in western China to study the crustal deformation before the Wenchuan M_S8.0 earthquake. The processing methods included analyses of the strain rate field, inversion of fault locking and the GPS velocity profiles. The GPS strain rate in the E-W direction in the Qinghai-Tibet block shows that extensional deformation was dominant in the western region of the block (west of 92.5° E), while compressive deformation predominated in the eastern region of the block (from 92.5° E to 100° E). On a regional scale, the hypocentral region of the Wenchuan earthquake was located at the edge of an intense compression deformation zone of about 1.9 × 10⁻⁸/a in an east-west direction. The characteristic deformation in the seismogenic fault was compressive with a dextral component. The compression deformation rate was greater in the fault's western region than in its eastern region, and the strain accumulation was very slow on the fault scale. The results of a fault locking inversion show that the locking fraction and slip deficit was greater in the middle-northern section of the seismogenic fault than in the southern section. The GPS velocity profile before the Wenchuan earthquake shows that the compression deformation was smaller than the dextral deformation, which is asymmetrical with respect to the distribution of co-seismic displacement. These deformation characteristics should provide some clues to the Wenchuan earthquake which occurred in the later period of the earthquake cycle.     

Estimation of Shallow S-Wave Velocity Structure in Damascus City,Syria, Using Microtremor Exploration

Array measurements of microtremors were carried out at thirty sites in Damascus city, Syria to estimate S-wave velocity structures of shallow soil formations for site effect analysis. The microtremor data were recorded by 6 vertical-component seismometers distributed along the circumferences of two circles as well as a 3-component seismometer deployed in the center. The phase velocities were estimated at each site from the vertical components of recorded microtremor data by using the Spatial Autocorrelation method. Then, Genetic Simulated Annealing Algorithm technique was applied for inversion of the phase velocities to estimate 1-D S-wave velocity structures beneath the sites. The inverted Vs profiles are not uniform in Damascus city and the results show that a shallow soft layer (∼200 m/s) appears in the eastern part of the city as well as the central part along Barada River. This layer controls the amplification distribution in the city with a high amplification mainly observed at the locations having this layer. The inversion results also show that the depth to the engineering bedrock (∼750 m/s) is very shallow along the foothills of Mt. Qasyoun in the north-west. Then the depth increases towards the east and the south. The maximum depth to the engineering bedrock (∼80 m) was observed in the southern part of Damascus. To validate the results of the inversions, the spectral ratios between the horizontal and vertical components (H/V) of the recorded microtremor data at the central seismometer were compared with the computed ellipticities of the fundamental-mode Rayleigh-waves based on the respective Vs structure. The results show a good agreement in a period range of 0.05 s to 0.5 s. In this period range, the dominant peaks of the H/V ratios are due to the overall effect of the velocity contrasts between the shallow layers representing the subsurface S-wave velocity structure. Moreover, the average S-wave velocity for the top 10 m of soils (V_S10) shows a better correlation with the averaged site amplification in a period range of 0.05 s to 0.5 s than V_S30 which indicates that V_S10 can be a better proxy for high-frequency site amplification in the case of Damascus city.     

Depth constraints on azimuthal anisotropy in the Great Basin from Rayleigh-wave phase velocity maps

We present fundamental-mode Rayleigh-wave azimuthally anisotropic phase velocity maps obtained for the Great Basin region at periods between 16 s and 102 s. These maps offer the first depth constraints on the origin of the semi-circular shear-wave splitting pattern observed in central Nevada, around a weak azimuthal anisotropy zone. A variety of explanations have been proposed to explain this signal, including an upwelling, toroidal mantle flow around a slab, lithospheric drip, and a megadetachment, but no consensus has been reached. Our phase velocity study helps constrain the three-dimensional anisotropic structure of the upper mantle in this region and contributes to a better understanding of the deformation mechanisms taking place beneath the western United States. The dispersion measurements were made using data from the USArray Transportable Array. At periods of 16 s and 18 s, which mostly sample the crust, we find a region of low anisotropy in central Nevada coinciding with locally reduced phase velocities, and surrounded by a semi-circular pattern of fast seismic directions. Away from central Nevada the fast directions are ～ N–S in the eastern Great Basin, NW–SE in the Walker Lane region, and they transition from E–W to N–S in the northwestern Great Basin. Our short-period phase velocity maps, combined with recent crustal receiver function results, are consistent with the presence of a semi-circular anisotropy signal in the lithosphere in the vicinity of a locally thick crust. At longer periods (28–102 s), which sample the uppermost mantle, isotropic phase velocities are significantly reduced across the study region, and fast directions are more uniform with an ～ E–W fast axis. The transition in phase velocities and anisotropy can be attributed to the lithosphere–asthenosphere boundary at depths of ～ 60 km. We interpret the fast seismic directions observed at longer periods in terms of present-day asthenospheric flow-driven deformation, possibly related to a combination of Juan de Fuca slab rollback and eastward-driven mantle flow from the Pacific asthenosphere. Our results also provide context to regional SKS splitting observations. We find that our short-period phase velocity anisotropy can only explain ～ 30% of the SKS splitting times, despite similar patterns in fast directions. This implies that the origin of the regional shear-wave splitting signal is complex and must also have a significant sublithospheric component.     

Velocity dispersion in fractured rocks in a wide frequency range

Experimental measurements of fracture-induced seismic waves velocity variations at frequencies ~ 1 kHz, ~ 40 kHz and ~ 1 MHz were performed directly in the field at the rocky outcrop and in the laboratory on specific rock samples collected from the outcrops. The peridotite–lherzolite outcrop appeared macroscopically uniform and contained three systems of visible parallel sub-vertical fractures. This rock has substantial bulk density and higher than average value of seismic wave velocity. The presence of fracture systems gives rise to its velocity anisotropy. The seismic waves passing through the rock fractures are subject to velocity dispersion and frequency dependent attenuation. Our data, obtained from field and laboratory measurements, were compared with theoretical model predictions. In this model we successfully used displacement discontinuity approach. For the velocity dispersion evaluation we used multi-frequency measurements. The a priori observation of orientations and densities of fracture sets allowed evaluation of their stiffness. Our approach revealed that the first arrivals of seismic waves can be used for evaluation of P-wave group velocities, the specific case, in which we expect anomalous velocity dispersion. Our observations contribute to the issue of up-scaling of well-log derived velocities in fractured rock to the scale of standard seismic exploration frequencies.     

One-minute rain rate distribution in Nigeria derived from TRMM satellite data

Data from the Tropical Rainfall Measuring Mission (TRMM) satellite sensors, the Microwave Imager (TMI, 3A12 V6) and other satellite sources (3B43 V6) have been used to derive the thunderstorm ratio β, total rain accumulation M, and 1-min rainfall rates, R1min, for 37 stations in Nigeria, for 0.001–1% of an average year, for the period 1998–2006. Results of the rain accumulations from the TRMM satellite (1998–2006) were compared with the data collected from 14 ground stations in Nigeria for the period 1991–2000. The two data sets are reasonably positively correlated, with correlation coefficients varying from 0.64 to 0.99. Deduced 1-min rainfall rates compared fairly well with the previous ground data of Ajayi and Ezekpo (1988. Development of climatic maps of rainfall rate and attenuation for microwave applications in Nigeria. The Nigerian Engineering 23(4), 13–30) with correlation coefficients varying from 0.17 to 0.97 in all 37 stations. The agreement was much better when compared with the International Telecommunications Union Radio communication Study group 3 digital maps with correlation coefficients varying from 0.84 to 0.98 in 23 locations; however there were negative correlation coefficient (of 0.2 in 7 stations) in the Middle Belt and a weak positive coefficient (of 0.09 in 6 stations) in the South South. Regionally the inferred mean annual 1-min rainfall rates are the highest in the South-East region with values between 111 and 125 mm/h throughout the 9 years, followed by the South-South region (105–124 mm/h). The lowest rainfall rate and rainfall accumulation occur in the North-West region (60–86 mm/h) followed, in ascending order, by the North-East (66–95 mm/h) region, the Middle-Belt region (76–102 mm/h) and the South-West region (77–110 mm/h). The present results were also compared with 9 tropical stations around the world and there was positive correlation between the results. The present results will be very useful for satellite rain attenuation modeling in the tropics and subtropical stations around the world.It is useful to note that one country, particularly one as large as Nigeria, can have significant variations in its rainfall characteristics for a variety of reasons, and this is borne out by the results presented.     

Transition zone structure under a stationary hot spot: Cape Verde

We report on a two-year seismic deployment in the Cape Verde Islands, one goal of which was to study the upper mantle to determine its structure under a hot spot that is stationary in the hot spot reference frame. We find from analysis of P-to-S receiver functions estimated from broadband seismic recordings that, within uncertainty, the time separation between the 410 and 660 km discontinuities is normal compared to radial earth models. Thus, to exist, even stationary hot spots do not require vertical thermal anomalies from deep melting sources anchored in the lower mantle or at the core–mantle boundary or their anomalies are narrower than ～ 250 km in the upper mantle.     

Slope angle and aspect as influencing factors on the accuracy of the SRTM and the ASTER GDEM databases

Presently, the application of digital elevation or surface models have increasing relevance in all areas of scientific research and in practical engineering applications. The ASTER GDEM and SRTM databases are the most widely used digital surface models, due to their free accessibility and global coverage. The SRTM model was produced using a radar-based technique and the ASTER GDEM was developed using optical stereo image-pairs. Therefore, as all models contain errors (i.e. differences stemming from real surface or vertical biases), errors in these models will also differ. Our aim was to examine these vertical biases and to calculate the rate of error variance. A TIN (Triangulated Irregular Network) model was used as a reference surface, derived from the contour lines of a large scale topographic map. Errors were evaluated with statistical and geoinformation techniques. We discovered significant differences between the surfaces. The mean difference between topographic elevations minus the SRTM-V2 is +2.6 ± 4 m, while the mean difference between topographic elevations minus the SRTM-V3 is +2.7 ± 2.5 m. With the GDEM, the mean difference was 2.7 ± 9.1 m. Furthermore, we found that in the case of SRTM, the differences were significant considering the aspects and the steepness of the slopes: southern and eastern directions and larger slope angles showed greater differences compared to the reference data. The GDEM V2 DEM had a larger error variance, but the error did not vary significantly with slope angle.