The 1957 great Aleutian earthquake

The 9 March 1957 Aleutian earthquake has been estimated as the third largest earthquake this century and has the longest aftershock zone of any earthquake ever recorded—1200 km. However, due to a lack of high-quality seismic data, the actual source parameters for this earthquake have been poorly determined. We have examined all the available waveform data to determine the seismic moment, rupture area, and slip distribution. These data include body, surface and tsunami waves. Using body waves, we have estimated the duration of significant moment release as 4 min. From surface wave analysis, we have determined that significant moment release occurred only in the western half of the aftershock zone and that the best estimate for the seismic moment is 50–100×10²⁰ Nm. Using the tsunami waveforms, we estimated the source area of the 1957 tsunami by backward propagation. The tsunami source area is smaller than the aftershock zone and is about 850 km long. This does not include the Unalaska Island area in the eastern end of the aftershock zone, making this area a possible seismic gap and a possible site of a future large or great earthquake. We also inverted the tsunami waveforms for the slip distribution. Slip on the 1957 rupture zone was highest in the western half near the epicenter. Little slip occurred in the eastern half. The moment is estimated as 88×10²⁰ Nm, orM _w =8.6, making it the seventh largest earthquake during the period 1900 to 1993. We also compare the 1957 earthquake to the 1986 Andreanof Islands earthquake, which occurred within a segment of the 1957 rupture area. The 1986 earthquake represents a rerupturing of the major 1957 asperity.     

Convective and absolute instability of baroclinic eddies

Abstract

The long time response to an arbitrary unstable disturbance in a two layer model on an f-plane is sought. It has been found that depending on the ratio of the shear to the average speed of the mean flow two types of baroclinic instabilities exist: convective and absolute. When the system supports convective instabilities the long time response to an initial pulse excitation decreases with time at a fixed point in space. When such a system is excited by a wave maker the steady state frequency of response of the system corresponds to a spatially amplifying wave oscillating with the frequency of the wave maker. If the dispersion relation yields a saddle point of the frequency in the wave number complex plane with positive imaginary part of the frequency the system supports absolute instabilities. The response of the system at any point in space excited by an arbitrary signal grows exponentially with time at a rate determined by the properties of the system at the saddle point. This response is different from that of unstable normal modes. It is conjectured that absolute instability may be responsible for local cyclogenesis activity in certain geographical regions.     

Paramater regimes for studying isolated eddies

Abstract

The parameter ranges relevant to propagating oceanic eddies are explored with a view to simplifying the equation sets previously used to simulate such eddies. The reduced gravity shallow-water model is used as a starting point. We consider eddies whose width is of the same order as, or somewhat larger than, the Rossby deformation radius. Two model equations emerge from these scalings, both permitting modon solutions and easy numerical integration. Both, too, give fair simulations of the shallow-water model, but disagree in details.

Using another approach, an ad hoc vortex model is used to predict the path of an eddy. This also gives a fair but not perfect simulation of the one-layer results. It appears that the only ideal way to reproduce primitive equation results accurately is to use the primitive equations themselves.     

阿留申岛弧的非火山颤动

在阿留申岛弧观测到了一些非火山颤动。在数座火山上的大跨区的地震台网也监测到颤动状的事件,因此不可能与任何一座火山体有关。虽然目前还不能给这些颤动定位,但我们根据多个台网的反射时差和振幅给出了大概的位置。多数非火山颤动是在所推断的太平洋消减板块闭锁区或从蠕动到闭锁转换区记录到的。有一次非火山颤动事件是在拉特群岛和安德烈亚诺夫群岛附近地区地震活动和火山活动增强时期记录到的。类似的同步增加以往在阿留申群岛记录到过,并被假设是由慢滑动事件引起的大范围应力变化造成的。2008年7月12日奥克莫克火山喷发约1小时后在科罗文和大锡特金地震台网记录到显著的非火山颤动。非火山颤动并不是分散地出现在整个阿留申岛弧,而似乎是集中在几个地区。用2500km的消减带,我们有机会检查这些地区并将其与不产生非火山颤动地区加以比较,以便确定产生该颤动的作用因素。     

Lens eddies in rotating two-fluid systems

Abstract

Laboratory experiments are conducted on the instantaneous release of a constant volume of intermediate density fluid along the interface of a two-layer fluid system in rigid body rotation about a vertical axis. Such a system leads to the development of a thin lens of fluid with anticyclonic motion which grows in radial extent to reach an equilibrium radius as the radial motion is captured by Coriolis effects, becomes unstable to non-symmetric dusturbances, and finally decays due to the action of fluid viscosity. Scaling and dimensional analysis arguments are advanced for the growth, equilibrium and decay phases of the motion. The scaling analysis for the initial growth and equilibrium phases are shown to be in good agreement with the experimental observations. The lens decay data do not collapse under the assumption of a simple Ekman spindown model assuming immiscible fluids. An empirical fit of the data for the decay phase is presented.     

Observation of baroclinic eddies southeast of Okinawa Island

In the region southeast of Okinawa, during May to July 2001, a cyclonic and an anticyclonic eddy were observed from combined measurements of hydrocasts, an upward-looking moored acoustic Doppler current profiler (MADCP), pressure-recording inverted echo sounders (PIESs), satellite altimetry, and a coastal tide gauge. The hydrographic data showed that the lowest/highest temperature (T) and salinity (S) anomalies from a 13-year mean for the same season were respectively -3.0/ 2.5℃ and -0.20/ 0.15 psu at 380/500 dbar for the cyclonic/anticyclonic eddies. From the PIES data, using a gravest empirical mode method, we estimated time-varying surface dynamic height (D) anomaly referred to 2000 dbar changing from -20 to 30 cm, and time-varying T and S anomalies at 500 dbar ranging through about ±2 ℃ and ±0.2 psu, respectively. The passage of the eddies caused variations of both satellite-measured sea surface height anomaly (SSHA) and tide-gauge-measured sea level anomaly to change from about –20 to 30 cm, consistent with the D anomaly from the PIESs. Bottom pressure sensors measured no variation related to these eddy activities, which indicated that the two eddies were dominated by baro-clinicity. Time series of SSHA map confirmed that the two eddies, originating from the North Pacific Subtropical Countercurrent region near 20°―30°N and 150°―160°E, traveled about 3000 km for about 18 months with mean westward propagation speed of about 6 cm/s, before arriving at the region southeast of Okinawa Island.     

Topographic eddies in temporally varying oceanic flows

Abstract

In this paper we examine the behaviour of oceanic unsteady flow impinging on isolated topography by means of numerical simulation. The ocean model is quasigeostrophic and forced by an oscillatory mean flow. The fluid domain is of the channel type and open-boundary numerical conditions are used to represent downstream and upstream flow.

In certain cases, vortex shedding, either cyclonic or anticyclonic, is observed in the lee of obstacles. Such shedding can be explained as the consequence of both an enhanced process of vorticity dissipation over the topography which locally affects the balance of potential vorticity on the advective timescale, and a periodic dominance of advective effects which sweep the fluid particles trapped on the seamount. For refined resolution and smallest viscosity the model will predict flows in which the shed eddies are coherent structures with closed streamlines.

The model suggests a mechanism by which topographically generated eddies may be swept away from a seamount in the ocean.     

Existence of eddies at crossroad of the Indonesian seas

An eddy-resolving Indo-Pacific ocean circulation model was applied to highlight the behavior of eddies throughout the Indonesian seas. The complexity of the topography and coastline at the entrance of the Makassar Strait induce an eddy-type throughflow, instead of a straightforward flow. A sill and a narrow passage in the Makassar strait creates a barrier and impedes the continuation of eddies from the Pacific ocean, but the existence of a steep deep basin (>500 m depth) between the Java and Flores seas indicates a possible area for eddies. Based on our numerical results, we described the presence of a unique eddy structure north of Lombok Island, which we designated the “Lombok Eddy” and verified it by performing a drifter release field experiment and reviewing monthly mean climatology data from the World Ocean Atlas 2001 and XBT PX2 track data. NCEP/NCAR reanalysis, satellite observation data, and mixed layer depth analysis were also used to confirm these processes. By analyzing numerical simulation results and available temperature datasets, two additional eddies were found. All eddies form primarily due to eastward local winds correlated with seasonal monsoon winds during the austral summer. These eddies vary synchronously at an interannual time scale. Together, they are referred to as the Lombok Eddy Train (LET), which affects the surface layer down to a depth of 60 m, and the intensity of the eddy system is strongly affected by mixed layer depth variability from December to February.     

An analytical study of general hyper-diffusivity and barotropic eddies

An exact solution to the barotropic potential vorticity equation is used to examine the properties of barotropic vortices under arbitrary nth-order hyper-diffusivity. Analytical expressions are derived for an eddy's lifetime, meridional drift, decay in size, and energy, as functions of the Coriolis parameter, order and magnitude of diffusivity, and the eddy's size, shape and strength. These expressions provide a simple explanation for many observed features of oceanic and atmospheric vortices. For example, the competition between the Coriolis effect and eddy strength in giving permitted eddy geometries; the bias towards a zonal anisotropy for large vortices but not for small ones; energetic preference for axisymmetry; poleward meridional drift of cyclonic vortices; and meridional speed variation depending on eddy geometry and strength.     

海洋中尺度涡与内波的地震图像

海洋反射地震通常用于调查、研究海底地质构造,勘探油气与天然气水合物资源.近期研究表明多道反射地震方法也可以对水柱的热盐细结构成像.中尺度涡与内波是重要的物理海洋现象,但是常规的物理海洋调查是在间隔若干公里的离散测站上进行的,水平分辨率较低,因此对中尺度涡的结构与内波的横向分布了解较差.本文利用在大西洋东部、南海采集的地震数据给出了低频反射地震可以对中尺度涡与内波清晰成像的新的证据.反射地震方法较传统海洋观测手段,具有明显的优势,主要体现在高的水平分辨率和短时间内对整个海水剖面进行成像方面．从地震剖面上,能够清楚地观测到中尺度涡、内波造成的反射特征变化,从而有助于改进对能量在不同尺度的海水运动之间传递过程的认识.     

Hydrodynamic modelling of mesoscale eddies in the Black Sea

A three dimensional structure of mesoscale circulation in the Black Sea is simulated using the Proudman Oceanographic Laboratory Coastal Ocean Modelling System. A number of sensitivity tests reveal the response of the model to changes in the horizontal resolution, time steps, and diffusion coefficients. Three numerical grids are examined with x-fine (3.2 km), fine (6.7 km) and coarse (25 km) resolution. It is found that the coarse grid significantly overestimates the energy of the currents and is not adequate even for the study of basin-scale circulation. The x-fine grid, on the other hand, does not give significant advantages compared to the fine grid, and the latter is used for the bulk of simulations. The most adequate parameters are chosen from the sensitivity study and used to model both the basin-scale circulation and day-to-day variability of mesoscale currents for the months of May and June of 2000. The model is forced with actual wind data every 6 h and monthly climatic data for evaporation, precipitation, heat fluxes and river run-off. The results of the fine grid model are compared favourably against the satellite imagery. The model adequately reproduces the general circulation and many mesoscale features including cyclonic and anticyclonic eddies, jets and filaments in different parts of the Black Sea. The model gives a realistic geographical distribution and parameters of mesoscale currents, such as size, shape and evolution of the eddies.     

A vortex-tube model of eddies in the inertial range

Abstract

A vortex-tube geometry of the cascade of energy to small-scale eddies, in the inertial range of fully-developed turbulence, is proposed. The model is a special case of the beta model of Frisch, Sulem and Nelkin (1978). We require that the cascade conserve the principal invariants of inviscid, incompressible flow, namely volume, topological knottedness, circulation, and, at discrete times marking the termination of steps in the cascade, energy. The process terminates in a finite time, as in any beta model, leaving behind a self-similar network of “inactive” tubes. We associate a self-similar scaling dimension D with the structure, equal to the Hausdorff dimension of the set of “active” tubes at the termination of the cascade. Because circulation Λ plays a key role in the analysis of the cascade, we refer to these vortex-tube geometries as “gamma models”. The viewpoint throughout is entirely deterministic.

We describe two examples of gamma models. In the ring geometry, an eddy is a vortex ring, and the cascade produces “rings upon rings”, so we allow cutting and fusing of tubes while conserving total helicity. In the preferred helical model, no cutting is needed, and the cascade produces an infinite progression of braided “coils upon coils”. We suggest that latter geometry as a candidate for the topology of a singularity of the inviscid limit of a Navier-Stokes flow, when modeled by discrete vortex tubes.

A crucial ingredient of a gamma model, not explicitly present in a beta model, is the possibility of “splitting” a vortex tube into sub-tubes carrying smaller circulation. We suggest a dynamical basis for this process, as an instability of tubes whose cores violate the Rayleigh criterion.

The parameters describing a gamma model are not uniquely determined by our study, but there is a “simplest” helical gamma model, involving minimal splitting and distortion of tubes. The dimension D of the structure is 13/5, with a scale factor Λ = 2^?5/4. This value of D agrees with that suggested by Hentschel and Procaccia (1982), by analogy with established results for certain branched polymers.     

Application of altimetry data assimilation on mesoscale eddies simulation

Mesoscale eddy plays an important role in the ocean circulation. In order to improve the simulation accuracy of the mesoscale eddies, a three-dimensional variation (3DVAR) data assimilation system called Ocean Variational Analysis System (OVALS) is coupled with a POM model to simulate the mesoscale eddies in the Northwest Pacific Ocean. In this system, the sea surface height anomaly (SSHA) data by satellite altimeters are assimilated and translated into pseudo temperature and salinity (T-S) profile data. Then, these profile data are taken as observation data to be assimilated again and produce the three-dimensional analysis T-S field. According to the characteristics of mesoscale eddy, the most appropriate assimilation parameters are set up and testified in this system. A ten years mesoscale eddies simulation and comparison experiment is made, which includes two schemes: assimilation and non-assimilation. The results of comparison between two schemes and the observation show that the simulation accuracy of the assimilation scheme is much better than that of non-assimilation, which verified that the altimetry data assimilation method can improve the simulation accuracy of the mesoscale dramatically and indicates that it is possible to use this system on the forecast of mesoscale eddies in the future.     

Seismicity and tectonics of the far Western Aleutian Islands

Study of focal mechanisms of earthquakes in the Near and Komandorsky Islands indicate that there are several distinct zones of tectonic activity. South of the Near Islands, normal faulting occurs in the trench east of 172°E and low-angle thrusting dominates the Aleutian ridge. Mechanisms indicate underthrusting as far west as Mednyy Island with strike-slip faulting restricted to the south and west of Beringa Island. A zone of northeast striking left-lateral faulting near 1645.°E is proposed to separate the Aleutian Ridge from Kamchatka Peninsula. This motion, as well as faulting north of the Komandorsky Islands, may be related to the existance of a buffer plate comprising the Aleutian Ridge in the Komandorsky Islands. Active subduction terminates near 173°E and the faulting north of the Komandorsky Islands may, in part, be due to the bouyancy of a remnant slab. Depth phase modelling indicates bulletin-reported depths are overestimated due to a misidentification of depth phases.     

Eddy energy sources and mesoscale eddies in the Sea of Okhotsk

Based on eddy-permitting ocean circulation model outputs, the mesoscale variability is studied in the Sea of Okhotsk. We confirmed that the simulated circulation reproduces the main features of the general circulation in the Sea of Okhotsk. In particular, it reproduced a complex structure of the East-Sakhalin current and the pronounced seasonal variability of this current. We established that the maximum of mean kinetic energy was associated with the East-Sakhalin Current. In order to uncover causes and mechanisms of the mesoscale variability, we studied the budget of eddy kinetic energy (EKE) in the Sea of Okhotsk. Spatial distribution of the EKE showed that intensive mesoscale variability occurs along the western boundary of the Sea of Okhotsk, where the East-Sakhalin Current extends. We revealed a pronounced seasonal variability of EKE with its maximum intensity in winter and its minimum intensity in summer. Analysis of EKE sources and rates of energy conversion revealed a leading role of time-varying (turbulent) wind stress in the generation of mesoscale variability along the western boundary of the Sea of Okhotsk in winter and spring. We established that a contribution of baroclinic instability predominates over that of barotropic instability in the generation of mesoscale variability along the western boundary of the Sea of Okhotsk. To demonstrate the mechanism of baroclinic instability, the simulated circulation was considered along the western boundary of the Sea of Okhotsk from January to April 2005. In April, the mesoscale anticyclonic eddies are observed along the western boundary of the Sea of Okhotsk. The role of the sea ice cover in the intensification of the mesoscale variability in the Sea of Okhotsk was discussed.     

Impact of eddies on ocean diapycnal mixing in Gulf Stream region

Based on high-resolution,Array for Real-time Geostrophic Oceanography(Argo)profiles and Sea Level Anomaly(SLA)data,this study statistically analyzes and compares turbulent diapycnal mixing profiles inside and outside mesoscale eddies in the Gulf Stream region.The result indicates that average diapycnal diffusivity at 300–540 m depths in anticyclonic eddies reaches4.0×10-5 m2 s-1.This is significantly higher than the 1.6×10-5 m2 s-1 outside eddies and 0.8×10-5 m2 s-1 in cyclonic eddies.Probabilities of diapycnal diffusivity greater than 10-4 m2 s-1 within anticyclonic and cyclonic eddies and outside eddies are29%,5%and 12%,respectively.However,magnitudes of average diapycnal diffusivity at 540–900 m depths in these three cases are of the same order,10-5 m2 s-1.Twenty-four of a total 38 anticyclonic eddies had enhanced mixing in the ocean interior,and 22 were observed during or shortly after strong winds.The coincidence between enhanced mixing and strong wind stress indicates that more wind-induced,near-inertial wave energy propagates downward in anticyclonic eddies.The deeper part of 12 profiles(below 540 m)in anticyclonic eddies had vertical overturns with Thorpe scale exceeding 5 m,among which three profiles had overturns reaching 20 m.Enhanced mixing may have occurred in deep layers of some profiles,although it was not evident in average conditions.     

Decadal variability of the Kuroshio Extension: mesoscale eddies and recirculations

An eddy-resolving multidecadal ocean model hindcast simulation is analyzed to investigate time-varying signals of the two recirculation gyres present respectively to the north and south of the Kuroshio Extension (KE) jet. The northern recirculation gyre (NRG), which has been detected at middepth recently by profiling float and moored current meter observations, is a major focus of the present study. Low-frequency variations in the intensity of the recirculation gyres are overall highly correlated with decadal variations of the KE jet induced by the basin-wide wind change. Modulation of the simulated mesoscale eddies and its relationship with the time-varying recirculation gyres are also evaluated. The simulated eddy kinetic energy in the upstream KE region is inversely correlated with the intensity of the NRG, consistent with previous observational studies. Eddy influence on the low-frequency modulation of the NRG intensity at middepth is further examined by a composite analysis of turbulent Sverdrup balance, assuming a potential vorticity balance between the mean advection and the convergent eddy fluxes during the different states of the recirculation gyre. The change in the NRG intensity is adequately explained by that inferred by the turbulent Sverdrup balance, suggesting that the eddy feedback triggers the low-frequency modulation of the NRG intensity at middepth.     

High-resolution ensemble forecasting for the Gulf of Mexico eddies and fronts

High-resolution models and realistic boundary conditions are necessary to reproduce the mesoscale dynamics of the Gulf of Mexico (GOM). In order to achieve this, we use a nested configuration of the Hybrid Coordinate Ocean Model (HYCOM), where the Atlantic TOPAZ system provides lateral boundary conditions to a high-resolution (5 km) model of the GOM . However, such models cannot provide accurate forecasts of mesoscale variability, such as eddy shedding event, without data assimilation. Eddy shedding events involve the rapid growth of nonlinear instabilities that are difficult to forecast. The known sources of error are the initial state, the atmospheric condition, and the lateral boundary condition. We present here the benefit of using a small ensemble forecast (10 members) for providing confidence indices for the prediction, while using a data assimilation scheme based on optimal interpolation. Our set of initial states is provided by using different values of a data assimilation parameter, while the atmospheric and lateral boundary conditions are perturbed randomly. Changes in the data assimilation parameter appear to control the main position of the large features of the GOM in the initial state, whereas changes in the boundary conditions (lateral and atmospheric) appears to control the propagation of cyclonic eddies at their boundary. The ensemble forecast is tested for the shedding of Eddy Yankee (2006). The Loop Current and eddy fronts observed from ocean color and altimetry are almost always within the estimated positions from the ensemble forecast. The ensemble spread is correlated both in space and time to the forecast error, which implies that confidence indices can be provided in addition to the forecast. Finally, the ensemble forecast permits the optimization of a data assimilation parameter for best performance at a given forecast horizon.