Turbulence measurements in a submarine canyon

Profiles of velocity turbulence in Monterey Canyon, made with a recently developed expendable probe, show the existence of a very turbulent bottom boundary layer. The turbulent flow is up to 170 m thick and has peak microscale shears of 1 m s⁻¹ per meter. The rate of dissipation of kinetic energy, based on the observed shear variance, averaged over the depth of the turbulent boundary layer ranged from 70 to 500 × 10⁻⁶W m⁻³. Temperature measurements indicate that the flow was up canyon at a time of low tide. The upper bound for the vertical eddy viscosity is estimated to be17 × 10⁻⁴m²s⁻¹ and for the vertical eddy diffusivity is estimated to be 15 × 10⁻⁴m²s⁻¹. The large vertical scale and the intensity of the observed boundary layer suggest that the flow in Monterey Canyon may be important for the renewal and circulation of water over the continental shelf in the bay area.     

Application of a spectral model to the calculation of wind drift currents in an idealized stratified sea

The solution of the hydrodynamic equations describing wind drift currents in an idealized stratified sea is developed using a modal expansion in the vertical. The influence of the magnitude of wind induced surface eddy viscosity, within the surface mixed layer, and reduced eddy viscosity at depth associated with a stable pycnocline, upon wind induced current structure is examined.Changes in modal structure with eddy viscosity, mixed layer depth, pycnocline thickness and total water depth are considered. It is evident from the modal equations that the time variation of the contribution of each mode to current structure is determined by the modal eigenvalue, mean eddy viscosity, water depth and the vertical integral of the mode squared.Calculations show that the magnitude and direction of steady state surface currents under strongly stratified conditions, are significantly influenced by mixed layer depth and pycnocline thickness. Surface eddy viscosity is also important, and in weakly stratified conditions eddy viscosity within and below the pycnocline influences surface current.The presence of a thin near surface ‘wall layer’ is shown to influence significantly the direction of surface current, but to have little effect upon currents below this surface layer.     

Importance of subgrid-scale parameterization in numerical simulations of lake circulation

Two subgrid-scale modeling techniques––Smagorinsky’s postulation for the horizontal eddy viscosity and the Mellor–Yamada level-2 model for the vertical eddy viscosity––are applied as turbulence closure conditions to numerical simulations of resolved-scale baroclinic lake circulations. The use of the total variation diminishing (TVD) technique in the numerical treatment of the advection terms in the governing equations depresses numerical diffusion to an acceptably low level and makes stable numerical performances possible with small eddy viscosities resulting from the turbulence closure parameterizations. The results show that, with regard to the effect of an external wind stress, the vertical turbulent mixing is mainly restricted to the topmost epilimnion with the order of magnitude for the vertical eddy viscosity of 10⁻³ m² s⁻¹, whilst the horizontal turbulent mixing may reach a somewhat deeper zone with an order of magnitude for the horizontal eddy viscosity of 0.1–1 m² s⁻¹. Their spatial and temporal variations and influences on numerical results are significant. A comparison with prescribed constant eddy viscosities clearly shows the importance of subgrid-scale closures on resolved-scale flows in the lake circulation simulation. A predetermination of the eddy viscosities is inappropriate and should be abandoned. Their values must be determined by suitable subgrid-scale closure techniques.     

Sea breezes drive currents on the inner continental shelf off southwest Western Australia

In southwest Western Australia, strong and persistent sea breezes are common between September and February. We hypothesized that on the inner continental shelf, in the absence of tidal forcing, the depth, magnitude, and lag times of the current speed and direction responses to sea breezes would vary though the water column as a function of the sea breeze intensity. To test this hypothesis, field data were used from four sites were that were in water depths of up to 13 m. Sites were located on the inner continental shelf and were on the open coast and in a semi-enclosed coastal embayment. The dominant spectral peak in currents at all sites indicated that the majority of the spectral energy contained in the currents was due to forcing by sea breezes. Currents were aligned with the local orientation of the shoreline. On a daily basis, the sea breezes resulted in increased current speeds and also changed the current directions through the water column. The correlation between wind–current speeds and directions with depth, and the lag time between the onset of the sea breeze and the response of currents, were dependent on the intensity of the sea breezes. A higher correlation between wind and current speeds occurred during strong sea breezes and was associated with shorter lag times for the response of the bottom currents. The lag times were validated with estimates of the vertical eddy viscosity. Solar heating caused the water column to stratify in summer and the sea breezes overcame this stratification. Sea breezes caused the mixed layer to deepen and the intensity of the stratification was correlated to the strength of the sea breezes. Weak sea breezes of <5 m s⁻¹ were associated with the strongest thermal stratification of the water column, up to 1°C between the surface and bottom layers (6 and 10 m below the surface). In comparison, strong sea breezes of >14 m s⁻¹ caused only slight thermal stratification up to 0.5°C. Apart from these effects on the vertical structure of water column, the sea breezes also influenced transport and mixing in the horizontal dimension. The sea breezes in southwest Western Australia rotated in an anticlockwise direction each day and this rotation was translated into the currents. This current rotation was more prominent in surface currents and in the coastal embayment compared to the open coast.     

Application of a sigma coordinate sea model to the calculation of wind-induced currents

A three-dimensional numerical sea model is formulated in terms of sigma coordinates in the vertical. The vertical grid spacing in the model is arbitrary and can be refined to give enhanced resolution in high shear regions (e.g., close to the sea surface in wind-driven flows, and/or across the thermocline in stratified flows). A method of accurately determining surface currents and indicating how fine a grid is required in the surface layer is described.The problem of determining a suitable formulation of vertical eddy viscosity to use in a model of wind-induced flow in a tidal sea is considered in detail. A formulation in which surface eddy viscosity depends upon the roughness of the sea surface and the transfer of momentum to depth by surface waves appears reasonable. Below the surface layer turbulence is related to the current at depth.Idealized calculations are performed to demonstrate the accuracy and stability of the sigma coordinate model. Results of these calculations indicate that the formulation of eddy viscosity developed in this paper can explain the high surface shears reported in lake measurements of wind-induced surface currents, and the lack of shear under strong wind conditions in the open sea (GORDON, 1982, Journal of Geophysical Research, 87, 1939–1951).Surface current to surface wind ratio are also computed.     

Loading effects in Metsähovi from the atmosphere and the Baltic Sea

Loading by atmosphere and by the Baltic Sea cause gravity change at Metsähovi, located 15 km from the open sea. Gravity is changed by both the Newtonian attraction of the loading mass and by the crustal deformation. We have performed loading calculations using appropriate Green's function for both gravity and deformation, for both atmospheric and Baltic loading. The loading by atmosphere has been computed using a detailed surface pressure field from high resolution limited area model (HIRLAM) for north Europe up to 10° distances. Baltic Sea level is modelled using tide gauge records. Calculations show that 1 m of uniform layer of water corresponds to 31 nm s⁻² in gravity and −11 mm in height. Modelled loading is compared with observations of the superconducting gravimeter T020 for years 1994–2002. The combination of HIRLAM and a tide gauge record decreases RMS of gravity residuals by 14% compared to single admittance in air pressure corrections without sea level data. Regression of gravity residuals on the tide gauge record at Helsinki (at 30 km distance) gives a gravity effect of 26 nm s⁻² m⁻¹ for Baltic loading.The gravity station is co-located with a permanent GPS station. We have also associated the loading effects of the atmosphere and of the Baltic Sea with temporal height variations. The range of modelled vertical motion due to air pressure was 46 mm and that due to sea level 18 mm. The total range was 38 mm. The effects of the Baltic Sea and of the atmosphere partly cancel each other, since at longer periods the inverse barometer assumption is valid. Regression of the modelled height on local air pressure gives −0.37 mm hPa⁻¹, corresponding approximately to width 6° for pressure system.We have tested the models using one year of daily GPS data. Multilinear regression on local air pressure and sea level in Helsinki gives the coefficient −0.34 mm hPa⁻¹ for pressure, and −11 mm m⁻¹ for sea level. These match model values. Loading by air pressure and Baltic Sea explains nearly 40% of the variance of daily GPS height solutions.     

The effect of the Baltic Sea level on gravity at the Metsähovi station

The loading effect of the Baltic Sea is immediately recognizable in the gravity record of the superconducting gravimeter T020 in Metsähovi, Finland, by simply inspecting residual gravity together with the tide gauge record at Helsinki 30 km away. The station is 10 km from the nearest bay of the Baltic Sea and 15 km from the open sea. Sea level variations in the Baltic are non-tidal and driven at short periods primarily by wind stress, at longer periods by water exchange through the Danish straits. Locally they can have a range of 2–3 m. Loading calculations show that a uniform layer of water covering the complete Baltic Sea increases the gravity in Metsähovi by 31 nm/s² per 1 m of water, and the vertical deformation is −11 mm. The observed gravity response to the local sea level is generally less, since the variations at short periods are far from uniform areally, the same water volume just being redistributed to different places. Regression of the whole gravity record (1994-2001) on local sea level gives 50–70% of the uniform layer response, as do loading calculations using actual water distributions derived from 11 tide gauges. However, both fits are dominated by some extreme values of short duration, and parts of the gravity record with long-period variations in sea level are close to the uniform layer response. The gravity observations can be used to test corrections for other co-located geodetic observations (GPS, satellite laser ranging) which are influenced by the load effect but not sensitive enough to discriminate between models.     

Observations of Reynolds Stresses in wind waves

Summary Data are presented concerning Reynolds Stresses in wind waves obtained from time series records of horizontal and vertical velocity components of motion beneath the ocean free surface. The stresses, of the order of 25 dyne cm^–2, are generally positive indicating horizontal momentum transfer downward through the dynamic wind wave regime. The magnitude of the observed stress increases with wind speed and sea state. The co-spectra show strong negative peaks which appear at the ambient wave frequencies and indicate that the correlations or eddy stresses of the gross wave motions are responsible for the momentum flux. This is a corroboration of results reported previously by the writer in this journal.     

The polar spirals of Mars may be due to glacier surges deflected by Coriolis forces

All previous accounts of the spiral patterns at the Martian poles emphasize that the north polar spiral is centered about the geographic pole, whereas that of the south polar region is off-set by about 4°. This paper demonstrates that the patterns near both poles are centered on topographic highs rather than the spin poles themselves. This is circumstantial evidence in favour of the relatively unexplored mechanism of radial outflow of viscous rock by gravity spreading.The hypothesis developed here is that the spiral patterns are essentially due to crevasse patterns formed perpendicular to flow lines which are perturbed by Coriolis forces. In order to account for a crevasse pattern that has a form concave to the east the angular deflection of an hypothetical ice flow emanating from the topographic high centered about the geographical north pole, must be about 40° or 0.7 radians in a westward direction at 85°N latitude.The polar cap rock has previously been assumed to consist mainly of either frozen carbondioxide or water ice. Corresponding viscosities (at 190 K) allow for the occurrence of radial outflow or gravity driven tectonics at a maximum rate of 1 cm a⁻¹, but the flow pattern remains unaffected by Coriolis forces.The spiral patterns of the Martian poles can be explained if the flowing mass has an occasional effective kinematic viscosity as low as about 7 × 10⁶ m² s⁻¹, because gravity tectonics will then be deflected by Coriolis forces resulting in appropriately curved flowlines. A tensile fracture pattern, resembling an anticlockwise spiral pattern perpendicular to the clockwise deflected flowlines may subsequently form by local brittle failure.The occasional kinematic viscosity 7 × 10⁶ m² s⁻¹ would cause flow rates of 0.2 m s⁻¹ along the slopes of the topographic highs. This velocity and the corresponding viscosity is tentatively thought to be possible when thermal and pressure runaway occurs in the polar layered deposits. This would mean glacier surges on the Martian poles are two orders of magnitude faster than those hitherto observed on Earth.     

First observation of mesospheric wind shear as high as 330 m s−1 km−1

Mesospheric wind profiles with an altitude resolution of 25 m have been obtained by means of radar tracking of foil chaff clouds. Such experiments were performed during winter 1990 at Biscarrosse, France (44^°N, 1^°W). On one flight, a wind shear as high as 330 m s⁻¹ km⁻¹ at 87.4 km and a region of dynamical instability between 86 and 88 km was measured. This wind shear is believed to be the largest value ever measured in the mesosphere. The region of dynamical instability results from a superposition of two wave motions, and is found to link well with enhanced turbulence and small-scale wave activity.     

Numerical investigation of an oceanic resonant regime induced by hurricane winds

The oceanic mixed layer (OML) response to an idealized hurricane with different propagation speeds is investigated using a two-layer reduced gravity ocean model. First, the model performances are examined with respect to available observations relative to Hurricane Frances (2004). Then, 11 idealized simulations are performed with a Holland (Mon Weather Rev 108(8):1212–1218, 1980) symmetric wind profile as surface forcing with storm propagation speeds ranging from 2 to 12 m s⁻¹. By varying this parameter, the phasing between atmospheric and oceanic scales is modified. Consequently, it leads to different momentum exchanges between the hurricane and the OML and to various oceanic responses. The present study determines how OML momentum and heat budgets depend on this parameter. The kinetic energy flux due to surface wind stress is found to strongly depend on the propagation speed and on the cross-track distance from the hurricane center. A resonant regime between surface winds and near-inertial currents is clearly identified. This regime maximizes locally the energy flux into the OML. For fast-moving hurricanes (>6 m s⁻¹), the ratio of kinetic energy converted into turbulence depends only on the wind stress energy input. For slow-moving hurricanes (<6 m s⁻¹), the upwelling induced by current divergence enhances this conversion by shallowing the OML depth. Regarding the thermodynamic response, two regimes are identified with respect to the propagation speed. For slow-moving hurricanes, the upwelling combined with a sharp temperature gradient at the OML base formed in the leading part of the storm maximizes the oceanic heat loss. For fast propagation speeds, the resonance mechanism sets up the cold wake on the right side of the hurricane track. These results suggest that the propagation speed is a parameter as important as the surface wind speed to accurately describe the oceanic response to a moving hurricane.     

海面动量、感热和水汽输送的动力学分析

基于大气边界层动量、感热和水汽通量的基本方程,定量地计算了波动海面的动量、感热和水汽通量。首先,应用Prandtl的混合长概念,推导出贴海面大气层中风速、位温和比湿的涡动交换率及其贴海面层厚度,并且证明了波面上位温或比湿贴海面层厚度与速度贴海面层厚度的比值,和平面上的相应比值完全相等。随后,利用空气动力学理论讨论了贴海面动量、感热和水汽输送的参数化问题。最后,对现有五种理论模式进行了比较说明。     

Effects of shear in the convective boundary layer: analysis of the turbulent kinetic energy budget

Effects of convective and mechanical turbulence at the entrainment zone are studied through the use of systematic Large-Eddy Simulation (LES) experiments. Five LES experiments with different shear characteristics in the quasi-steady barotropic boundary layer were conducted by increasing the value of the constant geostrophic wind by 5 m s^-1 until the geostrophic wind was equal to 20 m s^-1. The main result of this sensitivity analysis is that the convective boundary layer deepens with increasing wind speed due to the enhancement of the entrainment heat flux by the presence of shear. Regarding the evolution of the turbulence kinetic energy (TKE) budget for the studied cases, the following conclusions are drawn: (i) dissipation increases with shear, (ii) the transport and pressure terms decrease with increasing shear and can become a destruction term at the entrainment zone, and (iii) the time tendency of TKE remains small in all analyzed cases. Convective and local scaling arguments are applied to parameterize the TKE budget terms. Depending on the physical properties of each TKE budget contribution, two types of scaling parameters have been identified. For the processes influenced by mixed-layer properties, boundary layer depth and convective velocity have been used as scaling variables. On the contrary, if the physical processes are restricted to the entrainment zone, the inversion layer depth, the modulus of the horizontal velocity jump and the momentum fluxes at the inversion appear to be the natural choices for scaling these processes. A good fit of the TKE budget terms is obtained with the scaling, especially for shear contribution.     

The aeolian dust accumulation curve

This article presents a simple physical concept of aeolian dust accumulation, based on the behaviour of the subprocesses of dust deposition and dust erosion. The concept is tested in an aeolian dust wind tunnel. The agreement between the accumulation curve predicted by the model and the accumulation curve obtained in the experiments is close to perfect and shows that it is necessary to discriminate between the processes of aeolian dust deposition and aeolian dust accumulation. Two important thresholds determine the accumulation process. For wind speeds below the deflation threshold, the aeolian accumulation of dust increases linearly with the wind speed. For wind velocities between the deflation threshold and the accumulation limit, the sedimentation balance is above unity and there is still accumulation, though it rapidly drops once the deflation threshold has been exceeded. At wind speeds beyond the accumulation limit, the sedimentation balance is below unity and there will no longer be an accumulation of dust. The thresholds have been determined in a wind tunnel test at friction velocity u_* = 0·34 m s^?1 (deflation threshold) and u_* = 0·43 m s^?1 (accumulation limit), but these values are only indicative since they depend heavily on the characteristics of the accumulation surface and of the airborne grains. Copyright © 2001 John Wiley & Sons, Ltd.     

Study on the rheology of subducting slabs

We calculate thermal and phase structures of subducting slabs for different subducting velocities by a modified coupling code of the kinetic phase-transformation equations and the heat-diffusion equation with latent-heat release. Whereafter, we estimate their rheology structures based on the thermal and phase structures from the mineral physical point of view. At shallow depth, the upper layer has a high effective viscosity greater than 10³⁴Pa · s; while the lower layer has a relatively low effective viscosity, which is greater than 10²⁶Pa · s nevertheless. The effective viscosities below the kinetic phase boundary of olivine to wadsleyite decrease obviously, and reach a minimum of 10²²Pa · s. Small areas with higher effective viscosities exist above the depth of about 700 km in subducting slabs, which are produced by lower temperatures that are related with endothermic phase transformation of spinel to perovskite and magnesiowustite. The 1% and 99% isograds of spinel proportion delineate tortuous belts with low effective viscosities, which would affect the geodynamic behavior of subducting slabs.     

On the attenuation of sea waves by rain

Abstract

Some conflicting evidence on Reynolds' (1900) hypothesis that rain should attenuate any wave motion on the sea surface is discussed. It is concluded that rain drops ought to produce vortex rings in the sea which mix the water to a sufficient depth to affect most waves, as asserted by Reynolds. By introducing vertical and horizontal eddy viscosities to model the mixing process, an estimate is found for the rate of attenuation of wave energy by the rain. Consequently, the net effect on the wave field of attenuation by rain and generation by wind is calculated.     

Convective instabilities in a variable viscosity fluid cooled from above

Whether in the mantle or in magma chambers, convective flows are characterized by large variations of viscosity. We study the influence of the viscosity structure on the development of convective instabilities in a viscous fluid which is cooled from above. The upper and lower boundaries of the fluid are stress-free. A viscosity dependence with depth of the form ν₀ + ν₁ exp(?γ.z) is assumed. After the temperature of the top boundary is lowered, velocity and temperature perturbations are followed numerically until convective breakdown occurs. Viscosity contrasts of up to 10⁷ and Rayleigh numbers of up to 10⁸ are studied.For intermediate viscosity contrasts (around 10³), convective breakdown is characterized by the almost simultaneous appearance of two modes of instability. One involves the whole fluid layer, has a large horizontal wavelength (several times the layer depth) and exhibits plate-like behaviour. The other mode has a much smaller wavelength and develops below a rigid lid. The “whole layer” mode dominates for small viscosity contrasts but is suppressed by viscous dissipation at large viscosity contrasts.For the “rigid lid” mode, we emphasize that it is the form of the viscosity variation which determines the instability. For steep viscosity profiles, convective flow does not penetrate deeply in the viscous region and only weak convection develops. We propose a simple method to define the rigid lid thickness. We are thus able to compute the true depth extent and the effective driving temperature difference of convective flow. Because viscosity contrasts in the convecting region do not exceed 100, simple scaling arguments are sufficient to describe the instability. The critical wavelength is proportional to the thickness of the thermal boundary layer below the rigid lid. Convection occurs when a Rayleigh number defined locally exceeds a critical value of 160–200. Finally, we show that a local Rayleigh number can be computed at any depth in the fluid and that convection develops below depth z_r (the rigid lid thickness) such that this number is maximum.The simple similarity laws are applied to the upper mantle beneath oceans and yield estimates of 5 × 10¹⁵?5 × 10¹⁶ m² s^?1 for viscosity in the thermal boundary layer below the plate.     

Solar wind velocity distribution on the heliospheric current sheet during Carrington rotations 1787-1795

The solar wind velocity distribution in the heliosphere is best represented using a v-map, where velocity contours are plotted in heliographic latitude-longitude coordinates. It has already been established that low-speed regions of the solar wind on the source surface correspond to the maximum bright regions of the K-corona and the neutral line of the coronal magnetic field. In this analysis, v-maps on the source surface for Carrington rotations (CRs) 1787-1795, during 1987, have been prepared using the interplanetary scintillation measurements at Research Institute of Atmospherics (RIA), Nagoya Univ., Japan. These v-maps were then used to study the time evolution of the low-speed (\leq450 km s⁻¹) belt of the solar wind and to deduce the distribution of solar wind velocity on the heliospheric current sheet. The low-speed belt of the solar wind on the source surface was found to change from one CR to the next, implying a time evolution. Instead of a slow and systematic evolution, the pattern of distribution of solar wind changed dramatically at one particular solar rotation (CR 1792) and the distributions for the succeeding rotations were similar to this pattern. The low-speed region, in most cases, was found to be close to the solar equator and almost parallel to it. However, during some solar rotations, they were found to be organised in certain longitudes, leaving regions with longitudinal width greater than 30^○ free of low-speed solar wind, i.e. these regions were occupied by solar wind with velocities greater than 450 km s⁻¹. It is also noted from this study that the low-speed belt, in general, followed the neutral line of the coronal magnetic field, except in certain cases. The solar wind velocity on the heliospheric current sheet (HCS) varied in the range 300–585 km s⁻¹ during the period of study, and the pattern of velocity distribution varied from rotation to rotation.     

Three-dimensional structure of mesoscale eddies in the western tropical Pacific as revealed by a high-resolution ocean simulation

The three-dimensional structure of mesoscale eddies in the western tropical Pacific(6°S–20°N, 120°E–150°E)is investigated using a high-resolution ocean model simulation. Eddy detection and eddy tracking algorithms are applied to simulated horizontal velocity vectors, and the anticyclonic and cyclonic eddies identified are composited to obtain their three-dimensional structures. The mean lifetime of all long-lived eddies is about 52 days, and their mean diameter is 147 km. Two typical characteristics of mesoscale eddies are revealed and possible dynamic explanations are analyzed. One typical characteristic is that surface eddies are generally separated from subthermocline eddies along the bifurcation latitude(~13°N) of the North Equatorial Current in the western tropical Pacific, which may be associated with different eddy energy sources and vertical eddy energy fluxes in subtropical and tropical gyres. Surface eddies have maximum swirl velocities of 8–9 cm s~(-1) and can extend to about 1500 m depth. Subthermocline eddies occur below 200 m, with their cores at about 400–600 m depth, and their maximum swirl velocities can reach 10 cm s~(-1). The other typical characteristic is that the meridional velocity component of the eddy is much larger than the zonal component. This characteristic might be due to more zonal eddy pairs(two eddies at the same latitude),which is also supported by the zonal wavelength(about 200 km) in the high-frequency meridional velocity component of the horizontal velocity.     

Inertial oscillations in the Irish Sea

An inertial current with a maximum speed of about 0.20 m s¹ was observed in June 1981, by a current meter 17 m below the surface of the thermally stratified waters of the western Irish Sea. The dynamics of the surface waters have been modelled using contemporary wind data recorded at the site to drive a horizontally unbounded surface slab 20 m deep with a velocity decay time of about three inertial periods (1.86 days). From thermistor chain data, inertial frequency potential energy was observed to radiate beneath the thermocline, descending at a rate of about 0.06 cm s¹. The energy available for internal mixing from the observed oscillation may have been as much as 80 J m², compared with 150 J m² from wind-induced surface mixing. The mean surface mixing power in summer is about7.1 × 10⁴W m², while the mean shear-induced mixing power could be as much as2.0 × 10⁴W m². The likely occurrence of strong inertial events has been estimated using the wind record fcr 10 successive summers at Dublin airport. Large oscillations with speeds 0.20m s¹, or mixed layer kinetic energy(MLKE) > 400 J m², have an average duration of about 28 h. They occur about once every 2 months between mid-May and mid-August but seem to be more frequent from then until mid-September. Moderate oscillations (MLKE > 200 J m²) are more evenly distributed, occurring about once every 2 weeks. Small oscillations (MLKE > 40 J m²) occur nearly 50% of the time.