Accuracy assessment of ocean tide loading computations for precise geodetic observations

The increasing resolution of ground based gravity measurements (e.g. by superconducting gravimeters) as well as satellite based gravity field studies allows to study very small signals, globally as well as local. On the other hand, this requires the correction of such signals to uncover others. To study the Earth’s deep interior and the on-going dynamic processes requires the correction of disturbing signals, and one of these signals is related to ocean tidal loading. Although new ocean tide models are being derived from current satellite missions, there are still uncertainties.In this paper we present an intercomparison ocean tide models to test their fit to world-wide observations. Therefore, three TOPEX/POSEIDON (T/P) satellite derived models (CSR3.0, FES95.2 and TPXO.2) beside the classical SCHW80 model were selected for an accuracy assessment study. The selected models have been subjected to an intercomparison test, tide gauge validation test and comparison to 59 tidal gravity stations.The intercomparison test shows a good agreement between the T/P-based models for the open ocean and remarkable disagreement between the selected models in the coastal regions indicating that such models are still problematic in these regions. The tide gauge validation shows that the T/P derived models fit tide gauges better than SCHW80, with a better fit for the semidiurnal constituents than for the diurnal constituents. Comparing the gravimetric ocean-tide loading computed from the selected models with the residuals from a set of 59 tidal gravity stations shows that there is an improvement of the T/P derived models with respect to the Schwiderski model, especially in M2. However, this improvement is not as significant as the result of the comparison with the pelagic data. The procedure developed for the comparison of T/P derived models with SCHW80 is presented. The results provide not only information and improvement with regard to SCHW80, but also information about the properties of the new models. It is intended to continue this work applying the very recent models to see how they perform compared to this study.With this study we provide boundary conditions for the improvement of new ocean-tide models in order to benefit from the gravity measurements now possible regarding the evaluation of Earth structures and dynamic processes.     

Accuracy assessment of ocean tide loading computations for precise geodetic observations

The increasing resolution of ground based gravity measurements (e.g. by superconducting gravimeters) as well as satellite based gravity field studies allows to study very small signals, globally as well as local. On the other hand, this requires the correction of such signals to uncover others. To study the Earth’s deep interior and the on-going dynamic processes requires the correction of disturbing signals, and one of these signals is related to ocean tidal loading. Although new ocean tide models are being derived from current satellite missions, there are still uncertainties.In this paper we present an intercomparison ocean tide models to test their fit to world-wide observations. Therefore, three TOPEX/POSEIDON (T/P) satellite derived models (CSR3.0, FES95.2 and TPXO.2) beside the classical SCHW80 model were selected for an accuracy assessment study. The selected models have been subjected to an intercomparison test, tide gauge validation test and comparison to 59 tidal gravity stations.The intercomparison test shows a good agreement between the T/P-based models for the open ocean and remarkable disagreement between the selected models in the coastal regions indicating that such models are still problematic in these regions. The tide gauge validation shows that the T/P derived models fit tide gauges better than SCHW80, with a better fit for the semidiurnal constituents than for the diurnal constituents. Comparing the gravimetric ocean-tide loading computed from the selected models with the residuals from a set of 59 tidal gravity stations shows that there is an improvement of the T/P derived models with respect to the Schwiderski model, especially in M2. However, this improvement is not as significant as the result of the comparison with the pelagic data. The procedure developed for the comparison of T/P derived models with SCHW80 is presented. The results provide not only information and improvement with regard to SCHW80, but also information about the properties of the new models. It is intended to continue this work applying the very recent models to see how they perform compared to this study.With this study we provide boundary conditions for the improvement of new ocean-tide models in order to benefit from the gravity measurements now possible regarding the evaluation of Earth structures and dynamic processes.     

Morphotectonics of the Iberian Peninsula

A morphotectonic interpretation of the relief of the Iberian Peninsula (IP) is presented based on an analysis of its gross forms. Several geological–geophysical and geomorphological methods were used in order to build up a morphostructural sketch. Three main categories were established using the approach of Rantsman (1979): Territorial Units (1-megablock, 9-macroblocks, 34-mesoblocks, 1.374-blocks and 2.523-microblocks); Morphostructural Alignments (2-first, 8-second, 20-third and 1-fourth rank; and 43 knots between Morphostructural Alignments (second-fourth rank). The main seismic activity is concentrated on the first- and second-rank lineaments, and some important epicenters are located near the lineament intersections. The origin of the earthquakes in the vicinity of such knots can be explained by the forcing/pushing of macroblocks westwards. The existence of earthquakes along the lineaments may be explained by tension. From the present study it appears that earthquake occurrence in the IP is due principally to stress concentrations around morphotectonic zones. A seismotectonic interpretation of the IP is also presented. On this map three zones are distinguished (A, B, C), each of them with a different active level and dimensions.Received: 25 November 2000     

印尼地区地震断层面上的海潮负荷应力特征

对印尼地区3574次M≥5.0地震震源断层面上的海潮负荷应力进行计算,研究了地震断层上的负荷应力随深度变化以及随断层倾角和走向的变化特征.结果表明,地震断层上的负荷正应力大于负荷切应力.逆断层上的负荷应力最大,其它依次为斜滑断层、正断层和走滑断层.倾角范围在0°—50°的断层上的负荷正应力随倾角增大有减小趋势,负荷切应力在30°—50°倾角范围的断层上较大.负荷应力随断层走向的趋势性变化不明显.不同类型断层上的负荷正应力和切应力在震源深度15km处均达到最大,之后随深度增加而减小.当震源深度小于40km时,断层面上的海潮负荷应力的量级与固体潮应力相当.因此指出,在潮汐应力对地震影响的研究中,对震源深度小于40km的海洋和沿海地震,除固体潮汐应力外还需要考虑海潮负荷应力的影响.     

Assessing regional geodiversity: the Iberian Peninsula

Geodiversity is a landscape characteristic related to the heterogeneity of the physical properties of the earth surface. In this work, we quantify and compare geodiversity in several geodynamic zones of the Iberian Peninsula. For this purpose, we have developed a geographic information system (GIS) procedure to carry out a regional terrain classification based on geodiversity factors. A classification process helped to produce a morphometric map (10 classes), a morphoclimatic map (five classes) and a geological map (15 classes). These three maps were combined using an overlay operation (union) to obtain the final terrain classification (419 classes), which was then applied to calculate diversity landscape indices. The latter were calculated using common landscape diversity indices (Patch Richness Density, Shannon's Diversity Index, Shannon's Evenness Index, Simpson's Diversity Index and Simpson's Evenness Index), provided by FRAGSTATS free software. These indices were calculated for the whole landscape of the main Iberian geological regions, thus revealing a close relationship between some index values and the geological and geomorphological characteristics. The highest diversity values are associated with Alpine collisional orogens and reactivated chains of the Precambrian‐Palaeozoic massif. Intraplate orogen with sedimentary cover, characterized by extensive planation surfaces, have lower values. Mesozoic areas with no significant tectonic deformation and Cenozoic basins are characterized by the lowest diversity values. Amongst the latter, the major diversity is associated with the most dissected basins, which also present higher morphoclimatic variety. Though depending on the chosen scale and the landscape classification criteria, these indices provide an objective assessment of the regional geodiversity of Iberia. Copyright © 2009 John Wiley & Sons, Ltd.     

Seismic Hazard Estimate at the Iberian Peninsula

—?Seismic hazard at the Iberian Peninsula has been evaluated by using a methodology which combines both zonified and non-zonified probabilistic methods. Seismic sources are used when considering zones where certain calculation parameters may be considered homogeneous, as in zonified methods, while, on the other hand, earthquakes are considered wherever it has taken place, as in non-zonified methods. The methodology which is applied in this paper has been originally used to calculate the seismic hazard maps in the United States. In our case, it has been necessary to adapt the method to the specific features of the seismicity in the Iberian Peninsula and its geographical surroundings, not only with respect to its distribution and characteristics, but also with respect to the properties of the seismic catalog used.¶Geographically, the main feature of the result is the fact that it reflects both historical seismicity and current seismic clusters of the region. Despite the smoothing, maps show marked differences between several seismic zones; these differences becoming more noticeable as exposure time increases. Maximum seismic hazard is found to be in the southwestern region of the Peninsula, especially in the area of the Cape St. Vicent, and around Lisbon. The uncertainty of the results, without considering that due to the attenuation laws, as deduced from the other evaluation parameters, is quite stable, being more sensitive to the parameters b and m _max of the Gutenberg-Richter relation.     

Determination of vertical displacements over the coastal area of Korea due to the ocean tide loading using GPS observations

This paper describes the GPS applicability for detecting the vertical displacements of ground stations caused by ocean tide loading effects. An experiment was carried out using 12 permanent GPS stations located in the coastal area of Korea using data in the period 1 July until 26 August 2003. The relative height differences were calculated from hourly DGPS data processing based on the carrier-phase observation. The power spectra of the M₂ and N₂ constituents of ocean tide loading were derived using the CLEAN algorithm. The differential vertical displacements generated by the ocean tide loading effect are typically 3–25 mm in coastal area of the Korea. We compared the results from GPS with those of the ocean tide models, NAO.99Jb regional model and GOT00.2, FES99 global models. The M₂ (N₂) amplitude differences of vertical displacements between GPS and GOT00.2 is 1.22 ± 3.61 mm (1.01 ± 1.48 mm), and that of the M₂ (N₂) amplitude difference between GPS and FES99 is 0.04 ± 4.64 mm (0.64 ± 1.75 mm), whereas the M₂ (N₂) amplitude difference between GPS and NAO.99Jb are 0.05 ± 1.03 mm (0.86 ± 1.18 mm). The highest vertical displacements at the PALM station are found for 24.5 ± 0.7 mm from GPS observation, and 22.9 mm from the regional model NAO.99Jb and 13.17 and 10.00 mm from the global models GOT00.2 and FES99, respectively. These values show that the vertical displacements derived from GPS are in good agreement with those of the regional model NAO.99Jb around Korea, more than with the global models. This result indicated that GPS is an effective tool to measure the vertical displacement caused by the ocean tide loading effect in the coastal area, and we need to use the NAO.99Jb ocean tide model rather than the global ocean tide models in and around the Korean peninsula for position determination with permanent GPS installations. This work demonstrates that vertical displacement caused by the M₂ and N₂ constituents of ocean tide loading can be measured by carrier-phase DGPS.     

Elastic-anelastic Regional Structures for the Iberian Peninsula Obtained from a Rayleigh Wave Tomography and a Causal Uncoupled Inversion

The elastic and anelastic structure of the lithosphere and asthenosphere of the Iberian Peninsula is derived by means of tomographic techniques applied to local phase and group velocities and local attenuation coefficients of Rayleigh wave fundamental mode. The database consists of surface wavetrains recorded at the broadband stations located in the Iberian Peninsula on the occasion of the ILIHA project. Path-averaged phase and group velocities and attenuation coefficients were previously obtained by standard filtering techniques of surface wavetrains and, subsequently, local dispersion curves were computed according to the Yanovskaya-Ditmar formulation. First, a principal component analysis (PCA) and the average linkage (AL) clustering algorithm are applied to these local values in order to classify the Iberian Peninsula in several rather homogeneous domains from the viewpoint of the similarity of the corresponding local dispersion curves, without previous seismotectonic constraints. Second, averaged phase and group velocities and attenuation coefficients representing each homogeneous region are used to derive the respective elastic and anelastic models of the lithosphere and asthenosphere. This purpose is achieved by using the uncoupled causal inversion of phase and group velocities and attenuation coefficients. The main features of the homogeneous regions are discussed by taking as reference the Hercynic, Alpine and Neogene domains of the Iberian Peninsula, and two questions affecting the reliability of the elastic-anelastic models are revised. First, the coherence of the shear-velocity and Q_β⁻¹ models obtained by causal uncoupled inversion for each region is analysed. Second, the influence of the causal phase and group velocities on the shear-velocity models is evaluated by comparing elastic and anelastic models derived from causal uncoupled inversion with those deduced from non-causal inversion.     

Recent trends of extreme temperature indices for the Iberian Peninsula

Climate change and extreme climate events have a significant impact on societies and ecosystems. As a result, climate change projections, especially related with extreme temperature events, have gained increasing importance due to their impacts on the well-being of the population and ecosystems. However, most studies in the field are based on coarse global climate models (GCMs). In this study, we perform a high resolution downscaling simulation to evaluate recent trends of extreme temperature indices. The model used was Weather Research and Forecast (WRF) forced by MPI-ESM-LR, which has been shown to be one of the more robust models to simulate European climate. The domain used in the simulations includes the Iberian Peninsula and the simulation covers the 1986–2005 period (i.e. recent past). In order to study extreme temperature events, trends were computed using the Theil-Sen method for a set of temperature indexes defined by the Expert Team on Climate Change Detection and Indices (ETCCDI). For this, daily values of minimum and maximum temperatures were used. The trends of the indexes were computed for annual and seasonal values and the Mann-Kendall Trend test was used to evaluate their statistical significance. In order to validate the results, a second simulation, in which WRF was forced by ERA-Interim, was performed. The results suggest an increase in the number of warm days and warm nights, especially during summer and negative trends for cold nights and cold days for the summer and spring. For the winter, contrary to the expected, the results suggest an increase in cold days and cold nights (warming hiatus). This behavior is supported by the WRF simulation forced by ERA-Interim for the autumn days, pointing to an extension of the warming hiatus phenomenon to the remaining seasons. These results should be used with caution since the period used to calculate the trends may not be long enough for this purpose. However, the general sign of trends are similar for both simulations despite some differences in their magnitudes.     

Coastal and oceanic SST variability along the western Iberian Peninsula

The inter-annual variability of the sea surface temperature (SST) was analyzed along the western Iberian Peninsula in the region ranging from 9.5 °W to 21.5 °W and from 37.5 °N to 42.5 °N with a spatial resolution of 1°×1° from 1900 to 2008. Both coastal and oceanic SST showed an overall increase with warming and cooling cycles similar to those observed in the North Atlantic region and in previous regional studies. In addition, the evolution of coastal and ocean water has been observed to be different. In general, ocean water is more affected by the different warming–cooling cycles than coastal water. In spite of coast and ocean are highly influenced by global changes affecting the whole North Atlantic region, near shore SST has been observed to be correlated with local wind regime, which is itself a manifestation of the Eastern Atlantic (EA) teleconnection pattern.     

Regionalisation of precipitation for the Iberian Peninsula and climate change

Temporal variability of precipitation over the Iberian Peninsula (IP) has high spatial gradients. Therefore, statistics of the temporal behaviour of precipitation and derived quantities over the IP must be estimated taking into account these spatial gradients. Some statistics can be displayed over a map. However there are statistics, such as Probability Density Functions at each location of the IP, that are impossible to display in a map. Because of this, it is mandatory to reduce the number of degrees of freedom which, in this case, consists of a reduction of the time series representative of the IP domain. In this work, we present a spatial partition of the IP region into areas of similar precipitation. For that, an observed dataset of daily-total precipitation for the years between 1951 and 2003 was used. The land-only high resolution data was obtained on a regular grid with 0.2° resolution in the IP domain. This data was subjected to a k-means Cluster Analysis in order to divide the IP into K regions. The clustering was performed using the squared Euclidean distance. Four clusters of IP grid points, defining 4 IP regions, were identified. The grid points in each region share the same time-varying behaviour which is different from region to region. The annual precipitation discriminates the following regions: (1) north Iberia, (2) a large region extending from the centre to the Mediterranean shores of the IP, (3) a large region ranging from the centre to the western and southwestern shores of the Iberia, and (4) northwest Iberia. The regions obtained for the four seasons of the year are similar. These results are consistent with the thermodynamic characteristics described in the available literature. These Iberian regions were used to assess climate change of seasonal precipitation from the multi-model ensemble of the fifteen simulations provided by the European project ENSEMBLES. Probability Density Functions of annual- and seasonal-total precipitation, consecutive dry days, and total precipitation above the 95th percentile, averaged in each region were estimated for a reference climate (1961–1960), a near-future climate (2021–2050), and a distant-future climate (2069–2098). Climate change projections are based on comparisons of these functions between each future climate and the reference climate.Finally, we emphasize that: (i) the methodology used here, based on Cluster Analysis, can be used to regionalise other areas of the world, and (ii) the identified regions of the IP can be used to represent the Iberian precipitation by four time series that can be subjected to further analysis, whose results can be presented in a concise manner.     

Ten-year study of cloud-to-ground lightning activity in the Iberian Peninsula

The characteristics (annual and diurnal cycle, polarity, multiplicity and first stroke peak current) of ∼4.3×10⁶ cloud-to-ground flashes recorded in the Iberian Peninsula during the first decade of measurements of the lightning detection network installed in Spain are analyzed. The mean monthly variation shows maximum lightning activity between May and September, while minimum values are observed in January and February. The mean diurnal cycle shows maximum values at 1700 LT and minimum values at around 1000 LT. The average maximum flash density (not corrected for detection efficiency) is 2.1 flashes km⁻² year⁻¹. Maximum lightning activity is associated with mountainous areas. The effect of the Mediterranean Sea is also seen. The percentage of positive flashes is 9%, although this changes over the year from 6.5% in June to 22.6% in January. The average multiplicity is found to be 2.0 for the negative flashes and 1.1 for the positive flashes, and the percentages of single-stroke flashes are 53.6% and 89%, respectively. The monthly distribution of multiplicity for negative flashes peaks in the summer and minimum values are found in the winter. The multiplicity of the positive flashes does not seem to be function of the month. The median (mean) first stroke peak currents are found to be 23.5 kA (27.3 kA) for the negative flashes and 35.3 kA (47.1 kA) for the positive flashes. For both polarities, the peak current is higher in the summer than in the winter. The percentage of positive flashes and the mean peak currents for negative flashes are higher over the sea areas than over land.     

南极中山与昭和站重力海潮负荷效应及背景噪声研究

利用南极中山站LCR-ET21重力仪器与昭和站GWR058仪器获得的重力潮汐观测资料,采用最新的三个全球海潮模型（Dtu10,Eot11A和HAM11A）研究了南极地区的海潮负荷效应和背景噪声.结果表明,由三个海潮模型计算的重力负荷均值改正后,中山站O1和M2振幅观测残差分别由13.83%和20.55%下降到5.32%和5.95%,昭和站O1和M2振幅观测残差分别由10.84%和21.52%下降到1.91%和3.40%,说明海潮负荷改正的有效性.利用加汉宁窗的FFT变换,获得了地震频段的地震噪声等级（SNM）,其值分别为1.574（中山站）和1.289（昭和站）.而在潮汐频段,中山站的背景噪声比昭和站高一个数量级,主要由不同观测仪器和台站局部环境所致.本文结果可为进一步利用南极重力资料研究局部环境和全球动力学问题提供有效参考.     

Experiments in ocean circulation modelling

Abstract

In a laboratory model ocean, fluid in a rotating tank of varying depth is subjected to “wind-stress”, For a certain range of the parameters, Ekman number E and Rossby number R, a homogeneous fluid displays steady, westward intensified flow. For the same range of E and R, a two-layer fluid can have baroclinic instabilities. The parameter range for the various kinds of instabilities is mapped in a regime diagram. The northward transport in the western boundary current is measured as it varies with Rossby number for both homogeneous and two-layer fluid.     

Sky luminance distribution in the North of Iberian Peninsula during winter

This work is based on the study of outdoor daylight conditions during the winter period in Pamplona (the Iberian Peninsula). The 15 standard sky types proposed by the Commission Internationale de l’Eclairage (CIE) are used. The selected sky type in each moment is the one exhibiting the lowest root mean standard deviation (RMSD) when the theoretical and experimental luminance distributions in the sky hemisphere are compared. The selection was derived from luminance distribution data from 145 patches of the sky hemisphere and registered every 10 min from October to March from 2007 and 2008.The most frequent sky type observed in winter in Pamplona is V.5 (cloudless polluted with a broad solar corona), with a frequency of occurrence of 21%. Notwithstanding, the group of overcast skies exhibits a slightly higher frequency (40%) than the clear skies (37%). Seven out of the 15 standard sky types, i.e.V.5, I.1, III.4, II.1, IV.4, II.2 and III.3, are practically 80% of those under study. In order to provide a possible application in daylight climates studies, the frequency of occurrence of the 15 sky types is included in eleven intervals of solar elevation on the threshold corresponding to the time period and the location considered.     

Internal tide modelling and the influence of wind effects

Initially the development of shallow sea three-dimensional barotropic tidal models is briefly reviewed with a view to determining what were the key measurements that allowed progress in this field and rigorous model validation. Subsequently this is extended to a brief review of baroclinic tidal models to try to determine a “way forward” for baroclinic model development. The difficulty of high spatial variability, and wind influence are identified as possibly important issues that must be considered in validating baroclinic tidal models. These are examined using a three-dimensional unstructured grid model of the M₂ internal tide on the shelf edge region off the west coast of Scotland. The model is used to investigate the spatial variability of the M₂ internal tide, and associated turbulence energy and mixing in the region. Initial calculations are performed with tidal forcing only, with subsequent calculations briefly examining how the tidal distribution is modified by down-welling and up-welling favourable winds. Calculations with tidal forcing only, show that there is significant spatial variability in the internal tide and associated mixing in the region. In addition, these are influenced by wind effects which may have to be taken into account in any model validation exercise. The paper ends by discussing the comprehensive nature of data sets that need to be collected to validate internal tidal models to the same level currently attained with three dimensional barotropic tidal models.     

Solid earth tide and Arctic oceanic loading tide at Longyearbyen (Spitsbergen)

The three components of the indirect oceanic effect are calculated for the M ₂, K ₁ and O ₁ waves and compared with the observed earth tide. The vertical component of the near Arctic oceanic load explains fairly well the large 45° observed phase lag of M ₂. The results for the horizontal components are satisfactory; the discrepancies between different tiltmeters are not due to the oceanic perturbations (unless some local or cavity effects are supposed).