Ambipolar diffusion in the solar atmosphere

The process of non-linear ambipolar diffusion in the region overlying the solar surface can be an effective mechanism for producing sharp magnetic structures and current sheets. These may be the sites responsible for the occurrence of connectivity of magnetic field lines, and the subsequent explosive input of energy for heating of some of the features in the atmosphere of the Sun..     

Asymmetrically zoned reaction rims: assessment of grain boundary diffusivities and growth rates related to natural diffusion-controlled mineral reactions

This study explores garnet coronas around hedenbergite, which were formed by the reaction plagioclase + hedenbergite→garnet + quartz, to derive information about diffusion paths that allowed for material redistribution during reaction progress. Whereas quartz forms disconnected single grains along the garnet/hedenbergite boundaries, garnet forms ～20‐μm‐wide continuous polycrystalline rims along former plagioclase/hedenbergite phase boundaries. Individual garnet crystals are separated by low‐angle grain boundaries, which commonly form a direct link between the reaction interfaces of the plagioclase|garnet|hedenbergite succession. Compositional variations in garnet involve: (i) an overall asymmetric compositional zoning in Ca, Fe²⁺, Fe³⁺ and Al across the garnet layer; and (ii) micron‐scale compositional variations in the near‐grain boundary regions and along plagioclase/garnet phase boundaries. These compositional variations formed during garnet rim growth. Thereby, transfer of the chemical components occurred by a combination of fast‐path diffusion along grain boundaries within the garnet rim, slow diffusion through the interior of the garnet grains, and by fast diffusion along the garnet/plagioclase and the garnet/hedenbergite phase boundaries. Numerical simulation indicates that diffusion of Ca, Al and Fe²⁺ occurred about three to four, four and six to seven orders of magnitude faster along the grain boundaries than through the interior of the garnet grains. Fast‐path diffusion along grain boundaries contributed substantially to the bulk material transfer across the growing garnet rim. Despite the contribution of fast‐path diffusion, bulk diffusion through the garnet rim was too slow to allow for chemical equilibration of the phases involved in garnet rim formation even on a micrometre scale. Based on published garnet volume diffusion data the growth interval of a 20‐μm‐wide garnet rim is estimated at ～10³–10⁴ years at the inferred reaction conditions of 760 ± 50 °C at 7.6 kbar. Using the same parameterization of the growth law, 100‐μm‐ and 1‐mm‐thick garnet rims would grow within 10⁵–10⁶ and 10⁶–10⁷ years respectively.     

THERIA_G: a software program to numerically model prograde garnet growth

We present the software program THERIA_G, which allows for numerical simulation of garnet growth in a given volume of rock along any pressure–temperature–time (P–T–t) path. THERIA_G assumes thermodynamic equilibrium between the garnet rim and the rock matrix during growth and accounts for component fractionation associated with garnet formation as well as for intracrystalline diffusion within garnet. In addition, THERIA_G keeps track of changes in the equilibrium phase relations, which occur during garnet growth along the specified P–T–t trajectory. This is accomplished by the combination of two major modules: a Gibbs free energy minimization routine is used to calculate equilibrium phase relations including the volume and composition of successive garnet growth increments as P and T and the effective bulk rock composition change. With the second module intragranular multi-component diffusion is modelled for spherical garnet geometry. THERIA_G allows to simulate the formation of an entire garnet population, the nucleation and growth history of which is specified via the garnet crystal size frequency distribution. Garnet growth simulations with THERIA_G produce compositional profiles for the garnet porphyroblasts of each size class of a population and full information on equilibrium phase assemblages for any point along the specified P–T–t trajectory. The results of garnet growth simulation can be used to infer the P–T–t path of metamorphism from the chemical zoning of garnet porphyroblasts. With a hypothetical example of garnet growth in a pelitic rock we demonstrate that it is essential for the interpretation of the chemical zoning of garnet to account for the combined effects of the thermodynamic conditions of garnet growth, the nucleation history and intracrystalline diffusion. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.

Circulation in the Gulf of Aqaba (Red Sea) during winter—spring

This study is devoted to oceanographic features of the semi-enclosed Gulf of Aqaba, Red Sea. The data were recorded in winter—spring 1999 on the R/V Meteor cruise leg 44/2. Temperature and salinity profiles were measured at six positions (I—VI). The shipboard NarrowBand Acoustic Doppler Current Profiler (NB ADCP) 150 kHz continuously recorded current profiles down to 350 m en route. The research revealed that the current near the Strait of Tiran front (position VI) represents a semidiurnal signal of an internal tide wave (∼12 h period; 0.2 ms⁻¹ amplitude) that might be generated by the barotropic tide at the sill of the Strait. A sequence of cyclonic and anticyclonic eddy pairs is found along the axis of the Gulf of Aqaba during winter-spring seasons. These sub-mesoscale signals are dominant above the main thermocline and might be caused by wind forcing and the narrowness of the Gulf; it might remain in other seasons with different dimensions in relation to the depth of thermocline. The total diameter of each pair was twice the baroclinic Rossby radius (R ≈ 10 km). A single anti-cyclonic eddy was observed in the upper 300 m in the northern tip of the Gulf with a diameter of about 5–8 km.     

CO<Subscript>2</Subscript> Flux Measurements over the Sea Surface by Eddy Correlation and Aerodynamic Techniques

The measurements of the vertical transport of CO₂ were carried out over the Sea of Japan using the specially designed pier of Kyoto University on September 20 to 22, 2000. CO₂ fluxes were measured by the eddy correlation and aerodynamic techniques. Both techniques showed comparable CO₂ fluxes during sea breeze conditions: −0.001 to −0.08 mg m⁻²s⁻¹ with the mean of −0.05 mg m⁻²s⁻¹. This means that the measuring site satisfies the fetch requirement for meteorological observations under sea breeze conditions. Moreover, the eddy diffusivity coefficient used in the aerodynamic technique is found to be consistent with the coefficient used in the eddy correlation technique. The present result leads us to conclude that the aerodynamic technique may be applicable to underway CO₂ flux measurements over the ocean and may be used in place of the bulk technique. The important point is the need to maintain a measuring accuracy of CO₂ concentration difference of the order of 0.1 ppmv on the research vessels or the buoys.     

浅海流体动力学一种含变湍粘性系数的流速分解模型(Ⅱ)——模型在超浅海风暴潮零阶问题中的应用

本文将流速分解模型应用于作为超浅海风暴潮的渤海风潮,并讨论了变湍粘性系数的确定。作为一个初步的,但较为成功的数值试验例子,描述了实际风场作用下的渤海风潮,比较了变湍粘性系数模型与常湍粘性系数模型的计算结果间的差异。     

Three-dimensional Layer-integrated Modelling of Estuarine Flows with Flooding and Drying

Details are given of the refinement and application of a thee-dimensional (3-D) layer-integrated numerical model of tidal circulation, with the aim of simulating severe tidal conditions for practical engineering applications. The mode splitting strategy has been used in the model. A set of depth-integrated 2-D equations are first solved to give the pressure gradient, and the layer-integrated 3-D equations are then solved to obtain the vertical distributions of the flow velocities. Attention has been given to maintaining consistency of the physical quantities derived from the 2-D and 3-D equations. A TWO=layer mixing length turbulence model for the vertical shear stress distribution has been included in the model. Emphasis has been focused on applying the model to a real estuary, which is geometrically complicated and has large tidal ranges giving rise to extensive flooding and drying. The model has been applied to three examples, including: wind-driven flow in a rectangular lake, tidal circulation in a model rectangular harbour, and tidal circulation in a large estuary. Favourable results have been obtained for both the simple and complex flow beds.     

A model of the formation and development of ocean shear wave

Ａｍｏｄｅｌｏｆｔｈｅｆｏｒｍａｔｉｏｎａｎｄｄｅｖｅｌｏｐｍｅｎｔｏｆｏｃｅａｎｓｈｅａｒｗａｖｅ￥ＱｉａｏＦａｎｇｌｉａｎｄＺｈａｎｇＱｉｎｇｈｕａ（ＦｉｒｓｔＩｎｓｔｉｔｕｔｅｏｆＯｃｅａｎｏｇｒａｐｈｙ，ｓｔａｔｅｏｃｅａｎｉｃＡｄｍｉｎｉｓ...     

Monitoring the Kuroshio Extension with Dynamically Constrained Synthesis of the Acoustic Tomography, Satellite Altimeter and in situ Data

A finite-difference quasigeostrophic (QG) model of an open ocean region has been employed to produce a dynamically constrained synthesis of acoustic tomography and satellite altimetry data with in situ observations. The assimilation algorithm is based upon the 4D variational data interpolation scheme controlled by the model's initial and boundary conditions. The data sets analyzed include direct and differential travel times measured at the array of five acoustic transceivers deployed by JAMSTEC in the region of the Kuroshio Extension in 1997, Topex/Poseidon altimetry, CTD soundings, and ADCP velocity profiles. The region monitored is located within the area 27.5°–36.5°N, 143°–155°. The results of assimilation show that mesoscale variability can be effectively reconstructed by five transceivers measuring direct and reciprocal travel times supported by relatively sparse in situ measurements. The misfits between model and data lie within the observational error bars for all the data types used in assimilation. We have compared the results of assimilation with the statistical inversion of travel time data and analyzed energy balances of the optimized model solution. Energy exchange between the depth-averaged and shear components of the observed currents reveals a weak decay of the barotropic mode at the rate of 0.2 ± 0.7⋅10⁻⁵ cm²/s³ due to topographic interaction. Mean currents in the region are unstable with an estimate of the available potential energy flux from the mean current to the eddies of 4.7 ± 2.3⋅10⁻⁵ cm²/s³. Kinetic energy transition has the same sign and is estimated as 2.8 ± 2.5⋅10⁻⁵ cm²/s³. Potential enstrophy is transferred to the mesoscale at a rate of 5.5 ± 2.7⋅10⁻¹⁸ s⁻³. These figures provide observational evidence of the properties of free geostrophic turbulence which were predicted by theory and observed in numerical experiments. This revised version was published online in July 2006 with corrections to the Cover Date.     

Eddy Shedding from the Kuroshio Bend at Luzon Strait

TOPEX/POSEDIENT-ERS satellite altimeter data along with the mean state from the Parallel Ocean Climate Model result have been used to investigate the variation of Kuroshio intrusion and eddy shedding at Luzon Strait during 1992–2001. The Kuroshio penetrates into the South China Sea and forms a bend. The Kuroshio bend varies with time, periodically shedding anticyclonic eddies. Criteria of eddy shedding are identified: 1) When the shedding event occurs, there are usually two centers of high Sea Surface Height (SSH) together with negative geostrophic vorticity in the Kuroshio Bend (KB) area. 2) Between the two centers of high SSH there usually exists positive geostrophic vorticity. These criteria have been used to determine the eddy shedding times and locations. The most frequent eddy shedding intervals are 70, 80 and 90 days. In both the winter and summer monsoon period, the most frequent locations are 119.5°E and 120°E, which means that the seasonal variation of eddy shedding location is unclear.