Long-Term Stability of Sea Surface Topography from Topex/Poseidon Altimetry

The stability of the mean ocean level was investigated using the T/P altimeter data of 1993-1997 in 39 blocks of about 30° by 30°: 20 blocks forming the Pacific Ocean, 10 the Atlantic, and 9 blocks in the Indian Ocean. The 1993-1997 yearly means were found to be nearly constant, the computed linear terms came out as: (0.9±1.3) mm/year for the Pacific, (0.3±1.1) mm/year for the Atlantic, (–0.7 ± 1.4) mm/year for the Indian Ocean. No SST model was used in the solution.     

Validation of AVHRR and TMI-derived sea surface temperature in the northern South China Sea

Satellite-derived SSTs are validated in the northern South China Sea (NSCS) using in situ SSTs from the drifting buoys and well-calibrated sensors installed on Research/Vessel(R/V) Shiyan 3. The satellite SSTs are Advanced Very High Resolution Radiometer (AVHRR) daytime SST, AVHRR nighttime SST, Tropical rainfall Measuring Mission Microwave Imager (TMI) daytime SST and TMI nighttime SST. Availability of satellite SST, which is the ratio that the number of available satellite SST to the total ocean pixels in NSCS is calculated; annual average SST availabilities of AVHRR daytime SST, AVHRR nighttime SST, TMI daytime SST and TMI nighttime SST are 68.42%, 69.99%, 56.57% and 52.80%, respectively. Though the TMI SST availability is nearly constant throughout the year, the variations of the AVHRR SST availability are much larger because of seasonal variations of cloud cover in NSCS. Validation of the satellite-derived SSTs shows that bias±standard deviation (STD) of AVHRR SST is −0.43±0.76 and −0.33±0.79 °C for daytime and nighttime, respectively, and bias±STD of TMI SSTs is 0.07±1.11 and 0.00±0.97 °C for daytime and nighttime, respectively. It is clear that AVHRR SSTs have significant regional biases of about −0.4 °C against the drifting buoy SSTs. Differences between satellite-derived−in situ SSTs are investigated in terms of the diurnal SST cycle. When satellite-derived wind speeds decrease down below 6 m/s, the satellite SSTs become higher than the corresponding in situ SSTs, which means that the SST difference (satellite SST−Buoy SST) is positive. This wind-speed dependence of the SST difference is consistent with the previous results, which have mentioned that low wind speed coupled with clear sky conditions (high surface solar radiation) enhance the diurnal SST amplitude and the bulk-skin temperature difference.     

Computation of Ray Amplitudes in Inhomogeneous Media with Curved

The TOPEX/POSEIDON (T/P) satellite altimeter data from January 1, 1993 to January 3, 2001 (cycles 11–305) was used for investigating the long-term variations of the geoidal geopotential W ₀ and the geopotential scale factor R ₀ = GM÷W ₀ (GM is the adopted geocentric gravitational constant). The mean values over the whole period covered are W ₀ = (62 636 856.161 ± 0.002) m²s^-2, R ₀ = (6 363 672.5448 ± 0.0002) m. The actual accuracy is limited by the altimeter calibration error (2–3 cm) and it is conservatively estimated to be about ± 0.5 m²s^-2 (± 5 cm). The differences between the yearly mean sea surface (MSS) levels came out as follows: 1993–1994: –(1.2 ± 0.7) mm, 1994–1995: (0.5 ± 0.7) mm, 1995–1996: (0.5 ± 0.7) mm, 1996–1997: (0.1 ± 0.7) mm, 1997–1998: –(0.5 ± 0.7) mm, 1998–1999: (0.0 ± 0.7) mm and 1999–2000: (0.6 ± 0.7) mm. The corresponding rate of change in the MSS level (or R ₀) during the whole period of 1993–2000 is (0.02 ± 0.07) mm÷y. The value W ₀ was found to be quite stable, it depends only on the adopted GM, and the volume enclosed by surface W = W ₀. W ₀ can also uniquely define the reference (geoidal) surface that is required for a number of applications, including World Height System and General Relativity in precise time keeping and time definitions, that is why W ₀ is considered to be suitable for adoption as a primary astrogeodetic parameter. Furthermore, W ₀ provides a scale parameter for the Earth that is independent of the tidal reference system. After adopting a value for W ₀, the semi-major axis a of the Earth's general ellipsoid can easily be derived. However, an a priori condition should be posed first. Two conditions have been examined, namely an ellipsoid with the corresponding geopotential which fits best W ₀ in the least squares sense and an ellipsoid which has the global geopotential average equal to W ₀. It is demonstrated that both a-values are practically equal to the value obtained by the Pizzetti's theory of the level ellipsoid: a = (6 378 136.7 ± 0.05) m.     

Impact of Potential Climate Warming on Evaporation from Forest Drainage Basins in the Eastern Fennoscandia

The evaporation from forest area is shown to depend to a large degree on the growth conditions and productivity of forest stand; under the climatic conditions of the Eastern Fennoscandia, the variations in the productivity of forest stand of a particular forest type are mainly governed by the temperature. It is found that at the likely rise of the mean annual air temperature by 1 and 2°C, the average productivity of coniferous forest would increase by 0.3–0.7 (depending on the tree species and physiographic subzone) and 0.6–1.1 of the forest stand quality class, respectively; the increase in the annual evapotranspiration from a relatively large forest drainage basin could average 20–25 mm at a 2°C-temperature rise.     

Mean Earth'S Equipotential Surface From Topex/Poseidon Altimetry

The geopotential value of W ₀ = (62 636 855.611 ± 0.008) m ² s ^–2 which specifies the equipotential surface fitting the mean ocean surface best, was obtained from four years (1993 - 1996) of TOPEX/POSEIDON altimeter data (AVISO, 1995). The altimeter calibration error limits the actual accuracy of W ₀ to about (0.2 - 0.5) m ² s ^–2 (2 - 5) cm. The same accuracy limits also apply to the corresponding semimajor axis of the mean Earth's level ellipsoid a = 6 378 136.72 m (mean tide system), a = 6 378 136.62 m (zero tide system), a = 6 378 136.59 m (tide-free). The variations in the yearly mean values of the geopotential did not exceed ±0.025 m ² s ^–2 (±2.5 mm).     

On the Dynamical Effects of the Mantle in the Theory of Nutations

The influence of inerial forces in the mantle on its tidal deformations is considered. It is shown that the main role is played the corrections, described by spheroidal deformations, which result in corrections of nutational amplitudes of about –0.076 mas for the prograde semi-annual component, –0.056 mas for the retrograde annual component, 0.099 mas for the retrograde 19-year component, and –0.024 mas for the prograde 14-day component. These values exceed the errors of modem VLBI-measurements significantly (being of the order of 0.02-0.03 mas with the exception of the main nutational component); the effect of toroidal deformations is negligibly smalt in comparison with the accuracy of the VLBI-measurements, and for any practical purposes need not be considered.     

ANNUAL SEA LEVEL VARIABILITY INDUCED BY CHANGES IN SEA ICE EXTENT AND ACCUMULATION ON ICE SHEETS: AN ASSESSMENT BASED ON REMOTELY SENSED DATA

Changes of mean annual net accumulation at the surface on the grounded ice sheets of East Antarctica, West Antarctica and Greenland in response to variations in sea ice extent are estimated using grid-point values 100 km apart. The data bases are assembled principally by bilinear interpolation of remotely sensed brightness temperature (Nimbus-5 ESMR, Nimbus-7 SMMR), surface temperature (Nimbus-7 THIR), and surface elevation (ERS-1 radar altimeter). These data, complemented by field data where remotely sensed data are not available, are used in multivariate analyses in which mean annual accumulation (derived from firn emissivity) is the dependent variable; the independent variables are latitude, surface elevation, mean annual surface temperature, and mean annual distance to open ocean (as a source of energy and moisture). The last is the shortest distance measured between a grid point and the mean annual position of the 10% sea ice concentration boundary, and is used as an index of changes in sea ice extent as well as of mean concentration. Stepwise correlation analyses indicate that variations in sea ice extent of ± 50 km would lead to changes in accumulation inversely of ± 4% on East Antarctica, ± 10% on West Antarctica, and ±4% on Greenland. These results are compared with those obtained in a previous study using visually interpolated values from contoured compilations of field data; they substantiate the findings for the Antarctic ice sheets (±4% on East Antarctica, ±9% in West Antarctica), and suggest a reduction by one half of the probable change of accumulation on Greenland (from ±8%). The results also suggest a reduction of the combined contribution to sea level variability to ±0.19 mm a-1 (from ±0.22 mm a-1).     

Shallow seismicity, stress distribution and crustal deformation pattern in the Andaman-West Sunda arc and Andaman Sea, northeastern Indian Ocean

Shallow seismicity and available source mechanisms in the Andaman–westSunda arc and Andaman sea region suggest distinct variation in stressdistribution pattern both along and across the arc in the overriding plate.Seismotectonic regionalisation indicates that the region could be dividedinto eight broad seismogenic sources of relatively homogeneousdeformation. Crustal deformation rates have been determined for each oneof these sources based on the summation of moment tensors. The analysisshowed that the entire fore arc region is dominated by compressive stresseswith compression in a mean direction of N23^°, and the rates ofseismic deformation velocities in this belt decrease northward from 5.2± 0.65 mm/yr near Nias island off Sumatra and 1.12 ±0.13 mm/yr near Great Nicobar islands to as much as 0.4 ±0.04 mm/yr north of 8^°N along Andaman–Nicobar islandsregion. The deformation velocities indicate, extension of 0.83 ±0.05 mm/yr along N343^° and compression of 0.19 ±0.01 mm/yr along N73^° in the Andaman back arc spreadingregion, extension of 0.18 ± 0.01 mm/yr along N125^° andcompression of 0.16 ± 0.01 mm/yr along N35^° in NicobarDeep and west Andaman fault zone, compression of 0.84 ±0.12 mm/yr N341^° and extension of 0.77 ± 0.11 mm/yralong N72^° within the transverse tectonic zone in the Andamantrench, N-S compression of 3.19 ± 0.29 mm/yr and an E-Wextension of 1.24 ± 0.11 mm/yr in the Semangko fault zone ofnorth Sumatra. The vertical deformation suggests crustal thinning in theAndaman sea and crustal thickening in the fore arc and Semangko faultzones. The apparent stresses calculated for all major events range between0.1–10 bars and the values increase with increasing seismic moment.However, the apparent stress estimates neither indicate any significantvariation with faulting type nor display any variation across the arc, incontrast to the general observation that the fore arc thrust events showhigher stress levels in the shallow subduction zones. It is inferred that theoblique plate convergence, partial subduction of 90^°E Ridge innorth below the Andaman trench and the active back arc spreading are themain contributing factors for the observed stress field within the overridingplate in this region.     

The use of GPS arrays in detecting shock-acoustic waves generated during rocket launchings

This paper is concerned with the parameters of shock-acoustic waves (SAW) generated during rocket launchings. We have developed the interferometric method for determining SAW parameters (including angular characteristics of the wave vector, and the SAW phase velocity, as well as the direction towards the source) using GPS-arrays. Contrary to the conventional radio-probing techniques, the proposed method provides an estimate of SAW parameters without a priori information about the site and time of a rocket launching. The application of the method is illustrated by a case study of ionospheric effects from launchings of rockets PROTON, SOYUZ and SPACE SHUTTLE from Baikonur and Kennedy Space Center cosmodromes in 1998–2000. In spite of a difference of rocket characteristics, the ionospheric response for all launchings had the character of an N-wave corresponding to the form of a shock wave. The SAW period T is 270–360 s, and the amplitude exceeds the standard deviation of total electron content background fluctuations in this range of periods under quiet and moderate geomagnetic conditions by factors of 2–5 as a minimum. The angle of elevation of the SAW wave vector varies from 30° to 60°, and the SAW phase velocity (900–1200 m/s) approaches the sound velocity at heights of the ionospheric F-region maximum. The position of the SAW source, inferred by neglecting refraction corrections, corresponds to the segment of the rockets path at a distance no less than 200–900 km from the launch pad, and to the rocket flying altitude no less than 100 km. Our data are consistent with the existing view that SAW are generated during a nearly horizontal flight of the rocket with its engine in operation in the acceleration segment of the path at 100–130 km altitudes in the lower atmosphere.     

Quaternary vertical movements of the earth's surface in East-Central Europe

Fluvial terraces in mountain territories and granulometric cycles in fluvial sediment complexes show 300–400 m uplifts in mountain regions and about the same sized subsidences in major East-Central European basins during the Quaternary (−2.4 Ma). The vertical movements are not at all regular and equal in the mountains nor in the great basins. The number of the terraces is different in the valleys and their mutual proportions also differ. Similarly, some local Quaternary basins were developed deeper in the peripheries of the large basins while the central parts are shallow. There are also some stable blocks between the mountains and the lowlands, which did not move essentially during the last two and a half million years.The course and velocity of Pliocene and Quaternary subsidence were proved in a local sub-basin of the Carpathian mountain arch by paleomagnetic measurements of cores from two deep boreholes. In the Körös Basin, filled with fine-grained fluvial sediment, the sedimentation rate was 0.16–0.19 mm/y in the last 700.000 years; 0.15–0.16 mm/y during the proceding one million years; and 0.22–0.28 mm/y during the Pliocene (from −5.2 to −2.4 Ma).In sub-basins filled with coarse-grained sediments, the sedimentation rate was 0.3–0.4 mm/y during the Quaternary.     

Auroral pulsations and accompanying VLF emissions

Results of simultaneous TV observations of pulsating auroral patches and ELF-VLF-emissions in the morning sector carried out in Sodankylä (Finland) on February 15, 1991 are presented. Auroral pulsating activity was typical having pulsating patches with characteristic periods of about 7 s. Narrow-band hiss emissions and chorus elements at intervals of 0.3–0.4 s formed the main ELF-VLF activity in the frequency range 1.0–2.5 kHz at the same time. The analysis of auroral images with time resolution of 0.04 s allowed perfectly separate analysis of spatial and temporal variations in the auroral luminosity. Mutual correspondence between the behaviour of the luminous auroral patches and the appearance of ELF noise type hiss emissions and VLF chorus trains was found in two intervals chosen for analysis. While the hiss emissions were associated with the appearance of luminosity inside a limited area close to the zenith, the structured VLF emissions were accompanied by rapid motion of luminosity inside the area. The spatial dimension of the pulsating area was about 45–50 km and luminosity propagated inside it with velocity of about 10–12 kms. We discuss a new approach to explain the 5–15 s auroral pulsation based on the theory of flowing cyclotron maser and relaxation characteristics of ionosphere.     

Lunar and solar daily variations of the geomagnetic field at Toolangi

Summary Mean hourly values of magnetic declination, horizontal intensity and vertical intensity observed at Toolangi during two ten year periods (1924–1933 and 1949–1958) have been analysed to determine their solar and luni-solar diurnal components. The results, showing the variations of the first four harmonic components with season, degree of magnetic activity and annual sunspot number, are tabulated and discussed. It is shown that there are marked differences in the dependence ofS andL on the various parameters and a tentative explanation of this phenomenon is given.     

High-precision 40Ar/39Ar sanidine geochronology of ignimbrites in the Mogollon-Datil volcanic field,southwestern New Mexico

⁴⁰Ar/³⁹Ar age spectra have been obtained from 85 sanidine separates from 36 ignimbrites and one rhyolitic lava in the latest Eocene-Oligocene Mogollon-Datil volcanic field of southwestern New Mexico. Of the 97 measured age spectra, 94 yield weighted-mean plateau ages each giving single-spectrum 1 precision of±0.25%–0.4% (±0.07–0.14 Ma). Replicate plateau age determinations for eight different samples show within-sample 1 precisions averaging ±0.25%. Plateau ages from multiple (n=3–8) samples of individual ignimbrites show 1 within-unit precision of ±0.1%–0.4% (±0.04–0.13 Ma). This within-unit precision represents a several-fold improvement over published K-Ar data for the same ignimbrites, and is similar to the range of precisions reported from single-crystal laser fusion studies. A further indication of the high precision of unit-mean ⁴⁰Ar/³⁰Ar ages is their close agreement with independently established stratigraphic order. Two samples failed to meet plateau criteria, apparently due to geologic contamination by older feldspars. Effects of minor contamination are shown by six other samples, which yielded slightly anomalous plateau ages. ⁴⁰Ar/³⁹Ar plateau ages permit resolution of units differing in age by 0.5% (0.15 Ma) or less. This high resolution, combined with paleomagnetic studies, has helped to correlate ignimbrites among isolated ranges and has allowed development of an integrated timestratigraphic framework for the volcanic field. Mogollon-Datil ignimbrites range in age from 36.2 to 24.3 Ma. Ignimbrite activity was strongly episodic, being confined to four brief (<2.6 m.y.) eruptive episodes separated by 1–3 m.y. gaps. Ignimbrite activity generally tended to migrate from the southeast toward the north and west.     

On the 18-day quasi-periodic oscillation in the ionosphere

The presence and persistence of an 18-day quasi-periodic oscillation in the ionospheric electron density variations were studied. The data of lower ionosphere (radio-wave absorption at equivalent frequency near 1 MHz), middle and upper ionosphere (critical frequencies f₀E and f₀F2) for the period 1970–1990 have been used in the analysis. Also, solar and geomagnetic activity data (the sunspot numbers Rz and solar radio flux F10.7 cm, and a_N index respectively) were used to compare the time variations of the ionospheric with the solar and geomagnetic activity data. Periodogram, complex demodulation, auto- and cross-correlation analysis have been used. It was found that 18-day quasi-periodic oscillation exists and persists in the temporal variations of the ionospheric parameters under study with high level of correlation and mean period of 18–19 days. The time variation of the amplitude of the 18-day quasi-periodic oscillation in the ionosphere seems to be modulated by the long-term solar cycle variations. Such oscillations exist in some solar and geomagnetic parameters and in the planetary wave activity of the middle atmosphere. The high similarities in the amplitude modulation, long-term amplitude variation, period range between the oscillation of investigated parameters and the global activity of oscillation suggests a possible solar influence on the 18-day quasi-periodic oscillation in the ionosphere.     

A fracture mechanics study of subcritical tensile cracking of quartz in wet environments

Load relaxation and cross-head displacement rate-change experiments have been used to establish log₁₀ stress intensity factor (K) versus log₁₀ crack velocity (v) diagrams for double torsion specimens, of synthetic quartz cracked on thea plane in liquid water and moist air.For crack propagation normal toz and normal tor at 20°C,K _Ic (the critical stress intensity factor) was found to be 0.852±0.045 MN·m^–3/2 and 1.002±0.048 MN·m^–3/2, respectively.Subcritical crack growth at velocities from 10^–3 m·s^–1 to 10^–9 m·s^–1 at temperatures from 20°C to 80°C is believed to be facilitated by chemical reaction between the siloxane bonds of the quartz and the water or water vapour of the environment (stress corrosion). The slopes, of isotherms in theK-v diagrams are dependent upon crystallographic orientation. The isotherms have a slope of 12±0.6 for cracking normal tor and 19.9±1.7 for cracking normal toz. The activation enthalpy for crack propagation in the former orientation in liquid water at temperatures from 20°C to 80°C is 52.5±3.8 kJ·mole^–1.A discussion is presented of the characteristics of theK-v diagrams for quartz.     

The lunar tide in theE-layer above Istanbul

Summary The lunar tide in theE-layer of the ionosphere above Istanbul has been determined by the analysis of thefo E parameter from 1964–1967. Semi-diurnal variations were found to be significant. The seasonal and monthly variations of the amplitude and phase of the lunar tide show very little differences.     

Loading effect of a self-consistent equilibrium ocean pole tide on the gravimetric parameters of the gravity pole tides at superconducting gravimeter stations

The gravimetric parameters of the gravity pole tide are the amplitude factor δ, which is the ratio of gravity variations induced by polar motion for a real Earth to variations computed for a rigid one, and the phase difference κ between the observed and the rigid gravity pole tide. They can be estimated from the records of superconducting gravimeters (SGs). However, they are affected by the loading effect of the ocean pole tide. Recent results from TOPEX/Poseidon (TP) altimeter confirm that the ocean pole tide has a self-consistent equilibrium response. Accordingly, we calculate the gravity loading effects as well as their influence on the gravimetric parameters of gravity pole tide at all the 26 SG stations in the world on the assumption of a self-consistent equilibrium ocean pole tide model. The gravity loading effect is evaluated between 1 January 1997 and 31 December 2006. Numerical results show that the amplitude of the gravity loading effect reaches 10⁻⁹ m s⁻², which is larger than the accuracy (10⁻¹⁰ m s⁻²) of a SG. The gravimetric factor δ is 1% larger at all SG stations. Then, the contribution of a self-consistent ocean pole tide to the pole tide gravimetric parameters cannot be ignored as it exceeds the current accuracy of the estimation of the pole tide gravity factors. For the nine stations studied in Ducarme et al. [Ducarme, B., Venedikov, A.P., Arnoso, J., et al., 2006. Global analysis of the GGP superconducting gravimeters network for the estimation of the pole tide gravimetric amplitude factor. J. Geodyn. 41, 334–344.], the mean of the modeled tidal factors δ_m = 1.1813 agrees very well with the result of a global analysis δ_CH = 1.1816 ± 0.0047 in that paper. On the other hand, the modeled phase difference κ_m varies from −0.273° to 0.351°. Comparing to the two main periods of the gravity pole tide, annual period and Chandler period, κ_m is too small to be considered. Therefore, The computed time difference κ_L induced by a self-consistent ocean pole tide produces a negligible effect on κ_m. It confirms the results of Ducarme et al., 2006, where no convincing time difference was found in the SG records.     

Average crustal structure of Sweden

Summary Records obtained at the permanent stations of the Swedish seismograph network from explosions carried out in Scandinavian waters in June 1969 are evaluated. The study includes determination of velocities for all crustal phases observed, furthermore of layer thicknesses, Poisson ratios and amplitude ratios. The purpose of the study is partly to provide a first approximation to the crustal structure in Sweden, partly to provide regional data for location of earthquakes and explosions in the area in the future. Average velocities (km/sec) are forPn 7.88±0.05,Pg1 6.25±0.08,Pg2 5.70,Sn 4.58±0.04,S ^* 3.70±0.04,Sg1 (Lg1) 3.58±0.03,Sg2 (Sg) 3.40±0.03,Rg 3.02±0.07. The average thickness is 12 km for the granitic layer, and 23 km for the basaltic layer, thus making the average crustal thickness equal to 35 km. Relative amplitudes plotted versus distance complete the dynamical side of the study and they are useful for identification of waves. A regional travel-time table is presented for the distance range 0°–10° with entries for each 0.1° and including all crustal phases read.     

Evaluation of magnitudes using USSR and US (NOS) techniques

Summary WWSSN seismograms for the Caucasus earthquake of May 14, 1970, and the Alma-Ata earthquake of June 5, 1970, were scaled form _b using techniques employed by USSR (P amplitude during first 20 seconds) and USA (P amplitude during first few cycles). USSR seismograms will be used in a second phase of this study. Using both techniques, the scatter at each station form _b (SPZ) was 0.0–0.3 with half of the stations having no scatter. The scatter among all of the stations was over 1 unit. Using both techniques, the scatter at each station form _b (LPZ) was zero except for NUR on May 14. The scatter among all of the stations was 1 unit on May 14 and 0.4 on June 5. This preliminary study indicates the scaling techniques have no effect in computingm _b (SPZ) andm _b (LPZ) but the differences are more likely related to source functions, crust and upper mantle response functions, and instrumental band pass.NOS Seismology personnel were transferred to NOAA-ESL in 1971.     

Trends in athens temperature

Summary Temperature records of a long period (1858–1952) for Athens, in the form of overlapping 30-year intervals, are analyzed with a view to be investigated temperature trends.While the winter temperature has been rising almost continously during the period under consideration, the increase of summer temperature — to which the annual temperature and the annual temperature range as well are rather similar — has occured altogether during the second half of the record. During the first half the summer and annual temperatures and the annual range were all decreasing. The oceanicity increased during the first half and decreased during the second.The total change has been a rise of 0.90 deg. in winter temperature, a fall of 0.7 followed by a rise of 1.00 deg. in summer temperature, and a fall of 0.3 then a rise of 0.70 deg. in annual temperature.The rapid increase of population and also variations in solar activity seem to have partly influenced the Athens temperature trends.Presented at the tenth General Assembly of the International Union of Geodesy and Geophysics in Rome (Italy), September 1954.