Reflection of equatorial Kelvin waves at eastern ocean boundaries Part I: hypothetical boundaries

A baroclinic shallow-water model is developed to investigate the effect of the orientation of the eastern ocean boundary on the behavior of equatorial Kelvin waves. The model is formulated in a spherical polar coordinate system and includes dissipation and non-linear terms, effects which have not been previously included in analytical approaches to the problem. Both equatorial and middle latitude response are considered given the large latitudinal extent used in the model. Baroclinic equatorial Kelvin waves of intraseasonal, seasonal and annual periods are introduced into the domain as pulses of finite width. Their subsequent reflection, transmission and dissipation are investigated. It is found that dissipation is very important for the transmission of wave energy along the boundary and for reflections from the boundary. The dissipation was found to be dependent not only on the presence of the coastal Kelvin waves in the domain, but also on the period of these coastal waves. In particular the dissipation increases with wave period. It is also shown that the equatorial β-plane approximation can allow an anomalous generation of Rossby waves at higher latitudes. Nonlinearities generally have a small effect on the solutions, within the confines of this model.     

Analysis of Rayleigh waves in saturated porous elastic media by finite element method

A numerical procedure for the analysis of Rayleigh waves in saturated porous elastic media is proposed by use of the finite element method. The layer stiffness matrix, the layer mass matrix and the layer damping matrix in a layered system are presented for the discretized form of the solid-fluid equilibrium equation proposed by Biot. In order to consider the influence of the permeability coefficient on the behavior of Rayleigh waves, attention is focused on the following states: ‘drained’ state, ‘undrained’ state and the states between two extremes of ‘drained’ and ‘undrained’ states. It is found from computed results that the permeability coefficient exerts a significant effect on dispersion curves and displacement distributions of Rayleigh waves in saturated porous media.     

A three-dimensional model study of warm core ring interaction with continental shelf and slope

The interaction of warm core rings with a western wall and shelf/slope is examined with a three-dimensional primitive-equation model. The model ring is initialized with an axisymmetric Gaussian-type anticyclonic eddy placing far from the coastal boundary to allow the ring to freely propagate towards the wall and shelf/slope. The ring initially propagates steadily to the southwest at about 3 km/day under the combined planetary β and nonlinear effects. When colliding with a wall, the ring adjusts into a ‘D’ shape and moves poleward under primarily the image effect. When colliding with a shelf and slope, the ring however becomes stalled and bounces on and off the shelf/slope with little net movement. Small cyclones marked by strong upwelling are generated near the shelfbreak. Cyclones and anticyclones also are spawned at the periphery of the ring. Satellite SST images and concurrent ADCP transects are used to illustrate the strong interaction of a Gulf Stream warm ring (99B) with the Middle Atlantic Bight.     

Internal soil moisture and piezometric responses to rainfall-induced shallow slope failures

Better knowledge regarding internal soil moisture and piezometric responses in the process of rainfall-induced shallow slope failures is the key to an effective prediction of the landslide and/or debris flow initiation. To this end, internal soil moisture and piezometric response of 0.7-m-deep, 1.5-m-wide, 1.7-m-high, and 3.94-m-long semi-infinite sandy slopes rested on a bi-linear impermeable bedrock were explored using a chute test facility with artificial rainfall applications. The internal response time defined by the inflection point of the soil moisture and piezometric response curves obtained along the soil–bedrock interface were closely related to some critical failure states, such as the slope toe failure and extensive slope failures. It was also found that the response times obtained at the point of abrupt bedrock slope decrease can be used as indicators for the initiation of rainfall-induced shallow slope failures. An investigation of spatial distributions of soil water content, ω (or degrees of saturation, S_r), in the slope at critical failure states shows that the 0.2 m – below – surface zone remains unsaturated with S_r 40–60%, regardless of their distances from the toe and the rainfall intensity. Non-uniform distributions of ω (or S_r) along the soil–bedrock interface at critical failure states were always associated with near-saturation states (S_r 80–100%) around the point of bedrock slope change or around the transient ‘toe’ upstream of the slumped mass induced by the retrogressive failure of the slope. These observations suggest the important role of the interflow along the soil–bedrock interface and the high soil water content (or high porewater pressure) around the point of bedrock slope deflection in the rainfall-induced failure of sandy slopes consisting of shallow impermeable bedrocks. The present study proposes an ‘internal response time’ criterion to substantiate the prediction of rainfall-induced shallow slope failures. It is believed that the ‘internal response time’ reflects the overall characteristics of a slope under rainfall infiltration and can be as useful as the conventional meteorology-based threshold times. The ‘internal response time’ theory can be generalized via numerical modeling of slope hydrology, slope geology and slope stability in the future.     

The Modified Mercalli intensity and the geometry of the sedimentary basin as scaling parameters of the frequency dependent duration of strong ground motion

Several new empirical equations of the frequency dependent duration of strong earthquake ground motion are presented. The duration is considered as being composed of two parts: (1) the duration of stong motion as it is observed at recording stations located on basement rocks, and (2) the prolongation of this duration for stations located on sediments. The first part, called the ‘basic duration’, is modelled in terms of the Modified Mercalli intensity and (in some cases) the hypocentral distance. The depth of the sediments under the station, the distance from the station to the rocks surrounding it, and the angular measure of the size of those rocks (as seen from the station) are chosen as the parameters for modelling the prolongation of the duration. The new empirical equations are compared (a) with each other, (b) with our previous models which used similar ‘prolongation’ terms, but the ‘basic duration’ was expressed in terms of the magnitude of the earthquake and the source-to-station distance, and (c) with models with ‘intensity-type’ ‘basic duration’, but with a simplified ‘prolongation’ term (the geological conditions at the stations are modeled by lumping all the sites into three groups: basement rock, sediments and intermediate geology). This collection of models is found to have good internal consistency. The choice of the proper model depends on the availability of the earthquake and site parameters. The residuals of the empirical regression equations are found to have similar distribution functions for all the models. An explicit functional form for such distributions is proposed, and the frequency dependent coefficients are found for all the models of duration. This allows one to predict (for each set of earthquake and site parameters) the probability of exceedance of any given level of duration of strong ground motion at a given frequency.     

The evolution of an internal bore at the Malin shelf break

Observations of internal waves were made at the Malin shelf edge during SESAME (Shelf Edge Studies Acoustic Measurement Experiment), a part of the NERC LOIS-SES experiment, in August-September 1996. These measurements provide a high resolution dataset demonstrating internal wave generation and propagation. This note presents observations of the evolution of an internal bore. The process is shown clearly in a sequence of thermistor chain tows across the shelf break covering a complete tidal cycle, as the double-sided bore transforms into a group of undulations and eventually into more distinct solitary waveforms. Current structures associated with the bore and waves were also observed by ship-mounted ADCP. Analysis of the waveforms in terms of the linear modes and empirical orthogonal functions (EOFs) indicate the dominance of the first mode, which is typical of a shallow water seasonal thermocline environment. Determination of the phase speed of the waves from the consecutive ship surveys enabled the Doppler shift in the towed data to be removed, allowing analysis of the real length scales of the waves. The bore evolution has been modelled using a first order non-linear KdV model for the first mode, initialised with the waveform in the first survey. Comparison of the model and the observations show close agreement in the amplitudes, length scales, phase speeds and separations of the leading internal waves as they evolve. Finally, analysis of the observed internal wave shapes indicates that, within the uncertainties of measurement, the wave-lengths lie between those predicted by first and second order soliton theory.     

In-plane foundation-soil interaction for embedded circular foundations

Foundation soil interaction is studied using an analytical two-dimensional model, for circular foundations embedded in a homogeneous elastic half-space and for incident plane P- and SV- and for surface Rayleigh waves. The scattered waves are expanded in complete series of cyclindrical wave functions. A detailed analysis is presented of the foundation response to unit amplitude incident waves as a function of the type of incident waves and angle of incidence, the depth of the embedment and the foundation mass per unit length.It is shown that free-field translations and point rotation approximate well the foundation input motion only for very long incident waves. For shorter incident waves, those in general overestimate the foundation input motion. Neglecting the rotation of the foundation input motion (which is usually done in practice) may eliminate a major contribution to the base excitation of buildings and may cause nonconservative estimates of the forces in these buildings. Incident waves appear as ‘longer’ to a shallow foundation than to a deeper foundation. Therefore, deeper foundations are more effective in reflecting and scattering the short incident waves.     

Wave-induced setup of the mean surface over a sloping beach

A new theoretical approach for the wave-induced setup over a sloping beach is presented that takes into consideration the explicit variations of the surface waves due to bottom slope and viscosity. In this way, the wave forcing of the mean Lagrangian volume fluxes is calculated without assuming that the local depth is constant. The analysis is valid in the region outside the surf zone and is based on the shallow-water assumption. A novel approach for separating the viscous damping of the waves from the frictional damping of the mean flow is introduced, where the mean Eulerian velocity is applied in the bottom stress for the mean fluxes. In the case where the onshore Lagrangian mean transport is zero, a new formula is derived for the Eulerian mean free surface slope, in which the effects of bottom slope, viscous wave damping and frictional bottom drag on the mean flow are clearly identified. The analysis suggests that viscous damping of the waves and frictional dissipation of the Eulerian near-bed return flow could lead to setup outside the surf zone.     

Magnetosheath Interaction with the High Latitude Magnetopause

We present both statistical and case studies of magnetosheath interaction with the high-latitude magnetopause on the basis of Interball-1 and other ISTP spacecraft data. We discuss those data along with recently published results on the topology of cusp-magnetosheath transition and the roles of nonlinear disturbances in mass and energy transfer across the high-latitude magnetopause. For sunward dipole tilts, a cusp throat is magnetically open for direct interaction with the incident flow that results in the creation of a turbulent boundary layer (TBL) over an indented magnetopause and downstream of the cusp. For antisunward tilts, the cusp throat is closed by a smooth magnetopause; demagnetized ‘plasma balls’ (with scale ∼ few R_E, an occurrence rate of ∼25% and trapped energetic particles) present a major magnetosheath plasma channel just inside the cusp. The flow interacts with the ‘plasma balls’ via reflected waves, which trigger a chaotization of up to 40% of the upstream kinetic energy. These waves propagate upstream of the TBL and initiate amplification of the existing magnetosheath waves and their cascade-like decays during downstream passage throughout the TBL. The most striking feature of the nonlinear interaction is the appearance of magnetosonic jets, accelerated up to an Alfvenic Mach number of 3. The characteristic impulsive local momentum loss is followed by decelerated Alfvenic flows and modulated by the TBL waves; momentum balance is conserved only on time scales of the Alfvenic flows (1/f_A ∼12 min). Wave trains at f_A∼1.3 mHz are capable of synchronizing interactions throughout the outer and inner boundary layers. The sonic/Alfvenic flows, bounded by current sheets, control the TBL spectral shape and result in non-Gaussian statistical characteristics of the disturbances, indicating the fluctuation intermittency. We suggest that the multi-scale TBL processes play at least a comparable role to that of macro-reconnection (remote from or in the cusp) in solar wind energy transformation and population of the magnetosphere by the magnetosheath plasma. Secondary micro-reconnection constitutes a necessary chain at the small-scale (∼ion gyroradius) edge of the TBL cascades. The thick TBL transforms the flow energy, including deceleration and heating of the flow in the open throat, ‘plasma ball’ and the region downstream of the cusp.     

Drifter observations of the Hebrides slope current and nearby circulation patterns

The mean flow at and around the Hebrides and Shetland Shelf slope is measured with ARGOS tracked drifters. Forty-two drifters drogued at 50 m were deployed in three circles over the Hebrides slope at 56.15°N in two releases, one on 5th December, 1995 and the second on 5–9th May, 1996. The circles span a distance of some 20 km from water depths of 200 m to 1200 m. Drifters are initially advected poleward along-slope by the Hebrides slope current at between 0.05 and 0.70 m s^–1 in a laterally constrained (25–50 km wide) jet-like flow. Drifters released in winter remained in the slope current for over 2000 km whilst summer drifters were lost from the slope current beyond the Wyville-Thomson Ridge, a major topographic feature at 60°N. Dispersion from the slope region into deeper waters occurs at bathymetric irregularities, particularly at the Anton Dohrn Seamount close to which the slope current is found to bifurcate, both in summer and winter, and at the Wyville-Thomson Ridge where drifters move into the Faeroe Shetland Channel. Dispersion onto the continental shelf occurs sporadically along the Hebrides slope. The initial dispersion around the Hebrides slope is remarkably sensitive to initial position, most of the drifters released in shallower water moving onto the shelf, whilst those in 1000 m or more are mostly carried away from the slope into deeper water near the Anton Dohrn Seamount. The dispersion coefficients estimated in directions parallel and normal to the local direction of the 500 m contour, approximately the position of the slope current core, are approximately 8.8 × 10³ m² s^–1 and 0.36 × 10³ m² s^–1, respectively, during winter, and 11.4 × 10³ m² s^–1 and 0.36 x 10³ m² s^–1, respectively, during summer. At the slope there is a minimum in across-slope mean velocity, Reynolds stress, and across-slope eddy correlations. The mean across-slope velocity associated with mass flux is about 4 × 10^–3 m s^–1 shelfward across the shelf break during winter and 2 × 10^–3 m s^–1 during summer. The drifters also sampled local patterns of circulation, and indicate that the source of water for the seasonal Fair Isle and East Shetland currents are the same, and drawn from Atlantic overflows at the Hebrides shelf.     

Water mass distribution in Fram Strait and over the Yermak Plateau in summer 1997

The water mass distribution in northern Fram Strait and over the Yermak Plateau in summer 1997 is described using CTD data from two cruises in the area. The West Spitsbergen Current was found to split, one part recirculated towards the west, while the other part, on entering the Arctic Ocean separated into two branches. The main inflow of Atlantic Water followed the Svalbard continental slope eastward, while a second, narrower, branch stayed west and north of the Yermak Plateau. The water column above the southeastern flank of the Yermak Plateau was distinctly colder and less saline than the two inflow branches. Immediately west of the outer inflow branch comparatively high temperatures in the Atlantic Layer suggested that a part of the extraordinarily warm Atlantic Water, observed in the boundary current in the Eurasian Basin in the early 1990s, was now returning, within the Eurasian Basin, toward Fram Strait. The upper layer west of the Yermak Plateau was cold, deep and comparably saline, similar to what has recently been observed in the interior Eurasian Basin. Closer to the Greenland continental slope the salinity of the upper layer became much lower, and the temperature maximum of the Atlantic Layer was occasionally below 0.5 °C, indicating water masses mainly derived from the Canadian Basin. This implies that the warm pulse of Atlantic Water had not yet made a complete circuit around the Arctic Ocean. The Atlantic Water of the West Spitsbergen Current recirculating within the strait did not extend as far towards Greenland as in the 1980s, leaving a broader passage for waters from the Atlantic and intermediate layers, exiting the Arctic Ocean. A possible interpretation is that the circulation pattern alternates between a strong recirculation of the West Spitsbergen Current in the strait, and a larger exchange of Atlantic Water between the Nordic Seas and the inner parts of the Arctic Ocean.     

Experimental evidence for the exsolution of cratonic peridotite from high-temperature harzburgite

Phase equilibrium experiments were performed on typical ‘oceanic’ and ‘cratonic’ peridotite compositions and a Ca, Al-rich orthopyroxene composition, to test the proposal that garnet lherzolites exsolved from high-temperature harzburgites, and to further our understanding of the origin of ancient cratonic lithospheres. ‘Oceanic’ peridotites crystallize a garnet harzburgite assemblage at pressures above 5 GPa in the temperature range 1450–1600°C, but at 5 GPa and temperatures less than 1450°C, crystallize clinopyroxene to become true lherzolites. ‘Cratonic’ peridotites crystallize a garnet harzburgite assemblage at pressures above 5 GPa in the temperature range 1300–1600°C. Garnet-free harzburgite crystallizes from both ‘cratonic’ and ‘oceanic’ peridotite at temperatures above 1450°C and pressures below 4.5–5 GPa. Phase relations for the high Ca, Al-rich orthopyroxene composition essentially mirror those for ‘oceanic’ peridotite.The complete solution of garnet and clinopyroxene into orthopyroxene observed in all three starting compositions at temperatures near or above the mantle solidus at pressures less than 6 GPa supports the hypothesis that garnet lherzolite could have exsolved from harzburgite. The inferred cooling path for the original high-temperature harzburgite protoliths of garnet lherzolites differs depending on bulk composition. The precursor harzburgite protoliths of garnet lherzolites and harzburgites with ‘cratonic’ bulk compositions apparently experienced simple isobaric cooling from formation temperatures near the peridotite solidus to those at which most of these peridotites were sampled in the mantle (< 1200°C). The cooling histories for harzburgite protoliths of sheared garnet lherzolites with ‘oceanic’ compositional affinity are speculated to have involved convective circulation of mantle material to depths deeper than those at which it was originally formed.Phase equilibria and compositional relationships for orthopyroxenes produced in phase equilibrium experiments on peridotite and komatiite are consistent with an origin for ‘cratonic’ peridotite as a residue of Archean komatiite extraction, which has since cooled and exsolved clinopyroxene and garnet to become the common low-temperature, coarse-grained peridotite thought to comprise the bulk of the mantle lithosphere beneath the Archean Kaapvaal craton.     

A dynamically consistent analysis of circulation and transports in the southwestern Weddell Sea

An inverse model is applied for the analysis of hydrographic and current meter data collected on the repeat WOCE section SR4 in the Weddell Sea in 1989–1992. The section crosses the Weddell Sea cyclonic gyre from Kapp Norvegia to the northern end of the Antarctic Peninsula. The concepts of geostrophy, conservation of planetary vorticity and hydrostatics are combined with advective balances of active and passive properties to provide a dynamically consistent circulation pattern. Our variational assimilation scheme allows the calculation of three-dimensional velocities in the section plane. Current speeds are small except along the coasts where they reach up to 12 cm/s. We diagnose a gyre transport of 34 Sverdrup which is associated with a poleward heat transport of 28 × 10¹² W corresponding to an average heat flux of 15 Wm^–2 in the Weddell Sea south of the transect. This exceeds the estimated local flux on the transect of 2 Wm^–2. As the transect is located mostly in the open ocean, we conclude that the shelf areas contribute significantly to the ocean-atmosphere exchange and are consequently key areas for the contribution of the Weddell Sea to global ocean ventilation. Conversion of water masses occuring south of the section transform 6.6 ± 1.1 Sv of the inflowing warm deep water into approximately equal amounts of Weddell Sea deep water and Weddell Sea bottom water. The volume transport of surface water equals in the in-and outflow. This means that almost all newly formed surface water is involved in the deep and bottom water formation. Comparison with the results obtained by pure velocity interpolation combined with a hydrographic data subset indicates major differences in the derived salt transports and the water mass conversion of the surface water. The differences can be explained by deviations in the structure of the upper ocean currents to which shelf areas contribute significantly. Additionally a rigorous variance analysis is performed. When only hydrographic data are used for the inversion both the gyre transport and the poleward heat transport are substantially lower. They amount to less than 40% of our best estimate while the standard deviations of both quantities are 6.5 Sv and 37 × 10¹² W, respectively. With the help of long-term current meter measurements these errors can be reduced to 2 Sv and 8 × 10¹² W. Our result underlines the importance of velocity data or equivalent information that helps to estimate the absolute velocities.     

An Asymptotic Analysis of a Simple Model for the Structure and Dynamics of the Ramdas Layer

This paper presents a complete asymptotic analysis of a simple model for the evolution of the nocturnal temperature distribution on bare soil in calm clear conditions. The model is based on a simplified flux emissivity scheme that provides a nondiffusive local approximation for estimating longwave radiative cooling near ground. An examination of the various parameters involved shows that the ratio of the characteristic radiative to the diffusive timescale in the problem is of order 10⁻³, and can therefore be treated as a small parameter (μ). Certain other plausible approximations and linearization lead to a new equation whose asymptotic solution as μ → 0 can be written in closed form. Four regimes, consisting of a transient at nominal sunset, a radiative-diffusive boundary (‘Ramdas’) layer on ground, a boundary layer transient and a radiative outer solution, are identified. The asymptotic solution reproduces all the qualitative features of more exact numerical simulations, including the occurrence of a lifted temperature minimum and its evolution during night, ranging from continuing growth to relatively sudden collapse of the Ramdas layer.     

The characteristics of breaking waves, bubble clouds, and near-surface currents observed using side-scan sonar

An inverted 248 kHz two-component side-scan sonar mounted on the sea bed in a mean depth of 34 m has been used to detect the clouds of bubbles produced by breaking surface waves. The sonar has a range of about 150 m. The breaking waves appear on the sonograph records as short-lived intense echoes, and form patterns which can be explained in terms of the behaviour of groups of waves in which the highest are breaking.The bubble clouds are slightly elongated in the wave direction when they are first formed and, in winds of 5.1 m s⁻¹, have lifetimes of up to about 5 min. Soon after a wave breaks, the horizontal motion of the fluid in which the bubbles are formed becomes similar to that of the surroundings, and the bubbles continue to be advected by the near-surface currents. The rate of drift of bubble clouds along the directions of the sonar beams allows the components of the currents to be measured.The sonographs show that large changes in currents can occur over horizontal distances of as little as 5 to 10 m when fronts are passing. The motion of the fronts through the water can be measured. The sonographs have been used to measure surface currents induced by internal waves. Bands of bubbles associated with Langmuir circulation can be detected in strong winds or when moderate winds are accompanied by heavy rain.     

The errors in surface runoff prediction by neglecting the relationship between infiltration rate and overland flow depth

Many simplifications are used in modeling surface runoff over a uniform slope. A very common simplification is to determine the infiltration rate independent of the overland flow depth and to combine it afterward with the kinematic-wave equation to determine the overland flow depth. Another simplication is to replace the spatially variable infiltration rates along the slope i(x, t) due to the water depth variations h(x,t) with an infiltration rate that is determined at a certain location along the slope. The aim of this study is to evaluate the errors induced by these simplications on predicted infiltration rates, overland flow depths, and total runoff volume. The error analysis is accomplished by comparing a simplified model with a model where the interaction between the overland flow depth and infiltration rate is counted. In this model, the infiltration rate is assumed to vary along the slope with the overland flow depth, even for homogeneous soil profiles. The kinematic-wave equation with interactive infiltration rate, calculated along the slopy by Richard's equation, are then solved by a finite difference scheme for a 100-m-long uniform slope. In the first error analysis, we study the effect of combining an ‘exact’ and ‘approximate’ one-dimensional infiltration rate with the kinematic-wave equation for three different soil surface roughness coefficients. The terms ‘exact’ and ‘approximate’ stand for the solution of Richard's equation with and without using the overland flow depth in the boundary condition, respectively. The simulations showed that higher infiltration rates and lower overland flow depths are obtained during the rising stage of the hydrograph when overland flow depth is used in the upper boundary condition of the one-dimensional Richard's equation. During the recession period, the simplified model predicts lower infiltration rates and higher overland flow depths. The absolute relative errors between the ‘exact’ and ‘approximate’ solutions are positively correlated to the overland flow depths which increase with the soil surface roughness coefficient. For this error analysis, the relative errors in surface runoff volume per unit slope width throughout the storm are much smaller than the relative errors in momentary overland flow depths and discharges due to the alternate signs of the deviations along the rising and falling stages. In the second error analysis, when the spatially variable infiltration rate along the slope i(x, t) is replaced in the kinematic-wave equation by i(t), calculated at the slope outlet, the overland flow depth is underestimated during the rising stage of the hydrograph and overestimated during the falling stage. The deviations during the rising stage are much smaller than the deviations during the falling stage, but they are of a longer duration. This occurs because the solution with i(x, t) recognizes that part of the slope becomes dry after rainfall stops, while overland flow still exists with i(t) determined at the slope outlet. As obtained for the first error analysis, the relative errors in surface runoff volume per unit slope width are also much smaller than the relative errors in momentary overland flow depths and discharges. The relation between the errors in overland flow depth and discharge to different mathematical simplifications enables to evaluate whether certain simplifications are justified or more computational efforts should be used.     

Influence of the topography of the surface of the Earth on vertical sounding of the temperature profile

Atmospheric temperature and humidity fields as well as information on other meteorological parameters are nowadays retrieved from radiance measurements recorded by operational meteorological satellites. Up to now, the inversion procedures used only take into account crude information on the topography of the Earth’s surface. However, the applied radiative transfer codes have to consider the Earth’s surface as the lower boundary of the atmospheric model and, therefore, need a more precise mean elevation and a classification of the roughness of the Earth’s surface. The influence of the topography of the Earth surface on retrieved temperature profiles is studied by using a physico-statistical inversion method. An objective analysis is made of the more precise mean elevation and derivation of roughness parameters using a new high-resolution digital elevation model (DEM) with a resolution of 500 m×500 m. By means of a geomorphological process and a newly developed topography rejection test, areas with a high surface roughness are localized and singled out. The influence of topography on the retrieved temperature profiles is illustrated by case studies. Changes are found predominantly in areas with a high variation of topography. Using the new high-resolution DEM and the topography rejection test, the geographical position of the calculated temperature profiles tends to be shifted towards areas with a small vertical variation of topography. The mean elevation determined by the new elevation model better characterizes the area observed. Hence, the temperature profiles can be calculated down to lower atmospheric levels. Furthermore, a guess profile better describing the atmospheric situation is selected by the more precise elevation. In addition, the temperature profiles obtained near the coast are improved considerably by the more precise determination of the surface property ‘sea’ and ‘land,’ respectively. Integration of an independent physical information such as topography leads, on average, to a slight improvement of the results of the physico-statistical inversion procedure. In some cases, however, significant improvements have been achieved regarding the desired accuracy of temperature profiles of the order of 1 K. In future, the spatial resolution of new high-resolution sounding instrumentation on the next generation of operational meteorological satellites will be increased. To exploit the resolving power of this new instrumentation, the different variation of the topography of the Earth surface, especially in regions with a high variation of topography, can be taken into account more precisely by using a high-resolution DEM.     

Horizontal to vertical spectral ratio and geological conditions: an analysis of strong motion data from Greece and Taiwan (SMART-1)

Nowadays, several research efforts have been directed towards the evaluation of the horizontal to vertical, (h/v), spectral ratio technique, following the pioneering work of Nakamura originally proposed for microtremor measurements. In the present study, an extension of this technique based on Greek and Taiwanese (SMART-1) strong motion data is attempted. Whenever possible comparisons with results obtained by the standard spectral ratio technique for the same sites are made. Regarding the Greek data it is concluded that (h/v) spectral ratio clearly distinguishes the range of resonant peaks between ‘rock’ and ‘alluvium’ sites. On the other hand, nonlinear effects that take place beneath SMART-1 array most probably obscure the expected consistency between (h/v) and standard spectral ratio technique.     

Passive tracers in a general circulation model of the Southern Ocean

Passive tracers are used in an off-line version of the United Kingdom Fine Resolution Antarctic Model (FRAM) to highlight features of the circulation and provide information on the inter-ocean exchange of water masses. The use of passive tracers allows a picture to be built up of the deep circulation which is not readily apparent from examination of the veloCity or density fields. Comparison of observations with FRAM results gives good agreement for many features of the Southern Ocean circulation. Tracer distributions are consistent with the concept of a global “conveyor belt” with a return path via the Agulhas retroflection region for the replenishment of North Atlantic Deep Water.     

Effect of topography and subsurface inhomogeneity on seismic rayleigh waves

The effect of topography and subsurface inhomogeneity on surface motion is investigated in the case of Rayleigh waves. In the previous paper, the same effect was investigated in the case of SV waves. Several types of topography, such as cliffs both with and without a soft layer at the foot of the slope, are considered. Computations are made using a new hybrid method combining a particle model with a finite element method. In cases of harmonic Rayleigh waves, surface motions with amplitudes as large as 1.5 to 5 times the horizontal surface displacement of the incident Rayleigh waves are produced near the slope and the sloping interface. When a Rayleigh wave propagating through a hard single-layered ground encounters a sloping interface where hard ground and soft ground make contact with each other, Rayleigh waves having two different, phase velocities are produced and they correspond to the fundamental mode, and the first mode determined by Haskell's method. In addition, the transient response when Rayleigh waves propagate through the cliff is also simulated. Assuming the vertical component of the Tokachi-oki Earthquake (1968) measured on the surface to be a Rayleigh wave, the incident Rayleigh wave can be obtained by a Fourier synthesis of eigenfunctions of Rayleigh waves.