磁宁静期磁尾爆发性整体流持续时间多点卫星研究

本文利用星簇CLUSTER的三颗卫星数据分析了磁宁静期磁尾爆发性整体流（BBFs, Bursty Bulk Flows）的时间尺度, 并与单个卫星的结果做了比较. 事例研究表明, 利用三颗卫星观测数据判断的BBFs的时间尺度比单个卫星的大一倍左右. 对于三颗卫星观测到的同一个BBFs, BBFs在晨昏方向上的摆动决定了CLUSTER的三个卫星观测到BBFs的先后次序. 三颗卫星的观测也显示了BBFs的高度局域化特征. 磁宁静期磁尾BBFs寿命的增大, 使得BBFs携带的质量和能量的地向输运增加. 这种地向输运增加的结果是: 磁尾储存的能量得到较为平稳的释放, 改变了亚暴起始产生的时间, 为解决磁层压力平衡矛盾（PBI, Pressure Balance Inconsistency）问题提供了新的思路.     

地球磁层里磁泡的跨尾尺度研究

爆发流(Busty Bulk Flows)事件是发生在地球磁层里的等离子体输运现象,磁泡模型能很好地解释这一过程.现有的理论和观测事实已给出了对磁泡在跨尾方向上的尺度以及其他重要的物理参数的估计,但由于观测手段的限制这些参数并不十分精确,而只是个数值域.本文从最新的Tsyganenko磁层模型出发,利用磁泡在极光区根部的位置参数和自编的磁力线跟踪程序,通过映射给出了对应的磁泡在磁赤道面跨尾方向上的尺度.计算结果与理论预计和观测事实相符.     

利用GPS网观测反射海啸波引发的电离层扰动

震中产生的海啸波传播到海岸或者遇到水下障碍时会发生反射,从而形成反射海啸波.本文利用稠密的日本GEONET网,首次在电离层扰动图中观测到2011年3月11日Tohoku地震引发的反射海啸波信号.观测到的电离层扰动与海平面的反射海啸波具有相似的波形、水平速度、方向、周期以及到达时间等传播特征,表明观测的电离层扰动为反射海啸波所引起,本文的观测结果表明反射海啸产生的大气内重力波也能向上传播到电离层与等离子体作用.     

利用低阶卫星重力位系数研究下地幔横向密度异常分布

本文试图采用卫星重力资料和一种新的反演方法来研究地幔的横向密度异常分布.先将密度异常△（r,,（?））在一个三维正交函数系下进行展开,其展开系数待定.然后,根据密度异常与重力扰动位之间的关系建立观测方程组,其中未知向量由密度异常展开系数组成,重力扰动观测向量由 GEM10B 重力模型中的位系数计算而得,并通过适当选取重力位系数的阶数,对观测向量进行滤波.最后,就下地幔（670km——CM 界面）作了实际计算.计算中,重力扰动位阶数取为2——11阶,密度异常展开式的截断阶数取为 K=4和 L=6,求解观测方程组时采用阻尼最小二乘法.结果表明:密度扰动值在670km 不连续面及核幔界面处达到极大值,且在环太平洋地区存在一高密度带,太平洋中部对应于一低密度区,这些特征与 Dziewonski 得到的下地幔三维波速异常分布特征相一致.但是,在南极地区、大西洋及印度洋部分地区,所得的密度异常分布与三维波速异常分布呈负相关,文章就其原因作了初步分析.      

不同水动力扰动下太湖沉积物的悬浮特征

利用2002年7月及2003年4月分别在太湖乌龟山和梅梁湾的水动力观测资料,计算了波浪和湖流产生的湖底切应力,分析了二者对太湖沉积物悬浮的不同贡献.结果发现,在小扰动情况下,流切应力大于波切应力,但不能引起沉积物的悬浮;在强扰动情况下,底流对沉积物悬浮的影响可以忽略,且扰动越强,其作用越小,而波浪的影响却愈加突出;在中等扰动下,二者的相互作用产生沉积物的悬浮,但波浪的作用较为显著.该结论能为计算太湖水体悬浮物浓度和动态内源释放提供重要的参考依据.     

磁尾爆发性整体流与亚暴的关系

利用星簇Cluster的三颗卫星(C1, C3和C4)在2001年和2002年的数据,研究磁尾爆发性整体流（BBF, Bursty Bulk Flow）事件及其与亚暴的关系. 三颗卫星可以同时观测到同一次BBF事件, 有时只有一颗或两颗观测到BBF,其原因有:（1）等离子体整体流的速度峰值满足BBF选取原则中的峰值要求, 但卫星的运行轨道不满足;（2）卫星处于所要求的区域内,等离子体整体流的速度峰值不满足;（3）中性片的复杂结构及其运动使得选取条件不能同时满足;（4）BBF空间分布高度局域化. 统计研究结果表明：单颗卫星观测存在局限性, C1、C3和C4卫星独立观测到BBF的持续总时间分别占它们联合观测到的5507%、7748%和5552%; 大部分亚暴爆发期间都能观测到BBF, 甚至在一次亚暴爆发期间观测到多次BBF; 少数亚暴爆发期间没有观测到BBF.     

极隙区电子束流模电磁不稳定性

解析地研究了中高度（离他心3－4个地球半径）极隙区极低密度上行电子束流引起的沿磁力线传播的电磁个稳定性，上行电子束流和背景等离子体都考虑成冷等离子体．结果表明，上行电子引起的左旋和右旋圆偏振电磁来流模是不稳定的，当它与离子回旋模耦合时增长率达到最大值，频率色散关系仍为电子束流模特征这此结果对解释权隙区纬度地面站低频电磁波观测资料和理解极隙区动力学过程是很有益的．     

内磁层电流体系的地磁效应

对Ⅱ区场向电流及其伴随的部分环电流和电离层电流组成的内磁层三维电流体系（PRFI电流系）的磁场效应进行了数值计算.这一三维电流体系在中低纬度地面产生的磁场呈现出特殊的纬度分布:X分量几乎不随纬度变化,Y分量随纬度增高近似呈线性变化.这些特征明显不同于对称环电流的磁场分布特征（X∝ cosφ,φ是纬度,Y=0）,也不同于DP2、S_q、L等电流体系的磁场分布特征.利用这一特征我们可以从地磁台子午链观测到的磁场扰动中分离出PRFI电流系的贡献.用1989年3月磁暴的实例检验了上述模型,观测结果与理论结果符合得很好.分析结果还表明,最大的D_st指数并不一定对应着最强的对称环电流.     

亚暴膨胀相近磁尾位形不稳定性模型Ⅰ．近磁尾位形不稳定性

用单流体和双流体ＭＨＤ近似，研究了近磁尾位形不稳定性（ＮＥＴＣ）．分析表明，ＮＥＴＣ可能存在两种漂移不稳定情况Ｃ１和Ｃ２与卫星观测资料对比显示，Ｃ２较容易在亚暴膨胀相前夕出现，它可以解释亚暴膨胀相期间磁场和等离子体扰动的特征周期、尾向传播速度、磁场扰动和等离子体压强扰动之间的位相关系，场向电流的周期性结构，西向涌浪头部的电子沉降和极光隆起等观测特性和现象．薄电流片的极端情况（Ｒｃ≈ｒｉ）不在本文的讨论范围之内．     

亚暴膨胀相近磁尾位形不稳定性模型Ⅰ．近磁尾位形不稳定性

用单流体和双流体ＭＨＤ近似，研究了近磁尾位形不稳定性（ＮＥＴＣ）．分析表明，ＮＥＴＣ可能存在两种漂移不稳定情况Ｃ１和Ｃ２与卫星观测资料对比显示，Ｃ２较容易在亚暴膨胀相前夕出现，它可以解释亚暴膨胀相期间磁场和等离子体扰动的特征周期、尾向传播速度、磁场扰动和等离子体压强扰动之间的位相关系，场向电流的周期性结构，西向涌浪头部的电子沉降和极光隆起等观测特性和现象．薄电流片的极端情况（Ｒｃ≈ｒｉ）不在本文的讨论范围之内．     

Jet drop enrichment of bacteria,virus, and dissolved organic material

Four parameters that control the enrichment of bacteria in jet drops are bubble scavenging, drop size, drop position in the jet set, and the type of bacteria. Without the scavenging of bacteria as a bubble rises through the water it is doubtful that the observed enrichment factors, EF, greater than 1000 could be obtained. There is a maximum in EF as a function of top jet drop size, and the EF decreases from the top to the bottom drop of the jet set. The efficiency by which bubbles scavenge bacteria varies with species. Presumably these parameters apply in some degree to the EF of virus in jet drops. Dissolved organic material in natural waters can adsorb to bubbles and contribute to a large EF in jet drops, but there is a feedback mechanism whereby changes in bubble surface free energy modify the jet drop-size distribution. However, there is reason to believe this will not significantly influence the jet drop-size distribution produced by bubbles in the sea.     

基于含气泡液体声波方程的海底冷泉数值模拟

海底冷泉羽状流与海底天然气水合物的分布密切相关,对水合物稳定带的边界具有指示作用,是未来能源勘探的重要领域.研究海底冷泉羽状流的地震响应特征,对确定天然气水合物的储集区域及成藏环境等均有重要意义.当前获得海底冷泉羽状流的地震响应主要通过数值模拟进行,然而该过程所依据的含气泡介质声速模型及随机介质理论不能完整描述海底冷泉的物理性质,采用的声波方程也不适用于高频地震波数值模拟.为了准确地实现海底冷泉羽状流地震波数值模拟,精确分析其地震响应特征,本文提出利用Keller-Miksis气泡振动模型来描述气泡在声波作用下的运动状态,同时考虑气泡间的相互作用,建立海底冷泉气泡模型.在此基础上,本文创新性地采用含气泡液体声波方程进行海底冷泉高频地震波数值模拟.数值模拟结果表明,本文提出的方法能够实现海底冷泉羽状流地震响应的高精度数值模拟.     

Ionospheric Signatures Of Bursty Bulk Flows

While the plasma convection in the Earth's magnetosphere was for a long timeconsidered to consist mostly of laminar flows with wide spatial extents, about adecade ago the phenomenon of bursty bulk flows (BBFs), which now could beunderstood as long, but narrow channels of fast earthward plasma flow in thecentral plasma sheet of the magnetospheric tail, was discovered. Soon after thisit became clear that such events are not exceptional, but a large portion of theearthward plasma flow in the inner part of the magnetospheric tail is organisedin this bursty, intermittent mode. Since the Earth's magnetosphere is connectedby highly conducting magnetic field lines with its ionosphere, the next logicalstep was the search for the ionospheric signatures of BBFs. We review the resultsobtained so far in this young field of Space Physics, in terms of theauroral and ground magnetic signatures caused by BBFs, the particle precipitationto the ionosphere, as well as of the ionospheric electrodynamics of the processesassociated with the BBFs in the magnetosphere. Finally, we briefly review somemodels of ionosphere-magnetosphere coupling with respect to their ability to explainthe ionospheric signatures observed.     

Comparison of methods for modelling the behaviour of bubbles produced by marine seismic airguns

Three models for the dynamics of seismic airgun‐generated bubbles and their associated far‐field signals are developed and compared with geophysical data. The first model of an airgun‐generated bubble uses a spherical approximation, the second is an approximate Lagrangian model which allows for small deformations from a spherical shape, whilst the final model is an axisymmetric boundary‐integral method which permits the bubble to evolve into highly non‐spherical geometries. The boundary‐integral method also allows both geometric interference and strong dynamic interactions in multi‐bubble studies. When comparing the spherical model to experimental data there are three apparent, significant differences: the magnitude of the primary pressure peak, which is greater in the model; the subsequent decay of the pressure peaks and motion – the experimental data demonstrating greater decay and a slower rise rate; and the frequency of oscillation, which is slower in the experimental data. It is believed that the first discrepancy is due to the initial stages of expansion where the compressed air is forced to sparge through the airgun ports. The other differences indicate that there is some other energy‐loss mechanism which is not accounted for in the spherical bubble model. Non‐spherical bubble behaviour is investigated through the use of two different deformable many‐bubble codes and their predictions are compared with the spherical model and experimental data. The Lagrangian model predicts the formation of a buoyancy‐driven liquid jet on the first collapse of a typical airgun bubble; however, the model breaks down when the bubble becomes significantly deformed, due to a low‐order spherical‐harmonic approximation for the potential. The axisymmetric boundary‐integral code models the jet shape accurately and it is found that these bubbles evolve to toroidal geometries when the jet impacts on the opposite surface of the bubble. This highly non‐spherical behaviour is readily observed on high‐speed films of airgun bubbles, and is one key source of energy loss; it damps the pulsations of the bubble and slows its rise speed. Inter‐bubble interactions are investigated using the two deformable bubble models, and the predictions are compared to field data. It was found that as the bubbles approach each other, their periods of oscillation increase in accordance with observations, and jets are formed in the direction of motion upon collapse.     

The influence of gas bubbles on sediment acoustic properties: in situ, laboratory, and theoretical results from Eckernförde Bay, Baltic sea

Acoustic turbidity caused by the presence of gas bubbles in seafloor sediments is a common occurrence worldwide,but is as yet poorly understood. The Coastal Benthic Boundary Layer experiment in the Baltic off northern Germany was planned to better characterize the acoustic response of a bubbly sediment horizon. In this context, in situ measurements of compressional wave speed and attenuation were made over the frequency range of 5–400 kHz in gassy sediments of Eckernförde Bay. Dispersion of compressional speed data was used to determine the upper limit of the frequency of methane bubble resonance at between 20 and 25 kHz. These data, combined with bubble size distributions determined from CT scans of sediments in cores retained at ambient pressure, yield estimates of effective bubble sizes of 0.3–5.0 mm equivalent radius. The highly variable spatial distribution of bubble volume and bubble size distribution is used to reconcile the otherwise contradictory frequency-dependent speed and attenuation data with theory. At acoustic frequencies above resonance (>25 kHz) compressional speed is unaffected by bubbles and scattering from bubbles dominates attenuation. At frequencies below resonance (<1 kHz) ‘compressibility effects’ dominate, speed is much lower (250 m s^-1) than bubble-free sediments, and attenuation is dominated by scattering from impedance contrasts. Between 1.5 and 25 kHz bubble resonance greatly affects speed and attenuation. Compressional speed in gassy sediments (1100–1200 m s^-1) determined at 5–15 kHz is variable and higher than predicted by theory (<250 m s^-1). These higher measured speeds result from two factors: speeds are an average of lower speeds in gassy sediments and higher speeds in bubble-free sediments; and the volume of smaller-sized bubbles which contribute to the lower observed speeds is much lower than total gas volume. The frequency-dependent acoustic propagation is further complicated as the mixture of bubble sizes selectively strips energy near bubble resonance frequencies (very high attenuation) allowing lower and higher frequency energy to propagate. It was also demonstrated that acoustic characterization of gassy sediments can be used to define bubble size distribution and fractional volume.     

Bubble coalescence in basalt flows: comparison of a numerical model with natural examples

Gas accumulation in magma may be aided by coalescence of bubbles because large coalesced bubbles rise faster than small bubbles. The observed size distribution of gas bubbles (vesicles) in lava flows supports the concept of post-eruptive coalescence. A numerical model predicts the effects of rise and coalescence consistent with observed features. The model uses given values for flow thickness, viscosity, volume percentage of gas bubbles, and an initial size distribution of bubbles together with a gravitational collection kernel to numerically integrate the stochastic collection equation and thereby compute a new size spectrum of bubbles after each time increment of conductive cooling of the flow. Bubbles rise and coalesce within a fluid interior sandwiched between fronts of solidification that advance inward with time from top and bottom. Bubbles that are overtaken by the solidification fronts cease to migrate. The model predicts the formation of upper and lower vesicle-rich zones separated by a vesicle-poor interior. The upper zone is broader, more vesicular, and has larger bubbles than the lower zone. Basaltic lava flows in northern California exhibit the predicted zonation of vesicularity and size distribution of vesicles as determined by an impregnation technique. In particular, the size distribution at the tops and bottoms of flows is essentially the same as the initial distribution, reflecting the rapid initial solidification at the bases and tops of the flows. Many large vesicles are present in the upper vesicular zones, consistent with expected formation as a result of bubble coalescence during solidification of the lava flows. Both the rocks and model show a bimodal or trimodal size distribution for the upper vesicular zone. This polymodality is explained by preferential coalescence of larger bubbles with subequal sizes. Vesicularity and vesicle size distribution are sensitive to atmospheric pressure because bubbles expand as they decompress during rise through the flow. The ratio of vesicularity in the upper to that in the lower part of a flow therefore depends not only on bubble rise and coalescence, but also on flow thickness and atmospheric pressure. Application of simple theory to the natural basalts suggests solidification of the basalts at 1.0±0.2 atm, consistent with the present atmospheric pressure. Paleobathymetry and paleoaltimetry are possible in view of the sensitivity of vesicle size distributions to atmospheric pressure. Thus, vesicular lava flows can be used to crudely estimate ancient elevations and/or sea level air pressure.     

Modeling discrete gas bubble formation and mobilization during subsurface heating of contaminated zones

During thermal remediation the increase in subsurface temperature can lead to bubble formation and mobilization. In order to investigate the effect of gas formation on resulting aqueous concentrations, a 2D finite difference flow and mass transport model was developed which incorporates a macroscopic invasion percolation (MIP) model to simulate bubble expansion and movement. The model was used to simulate three soil scenarios with different permeabilities and entry pressures at various operating temperatures and groundwater velocities. It was observed that discrete bubble formation occurred in all three soils, upward mobility being limited by lower temperatures and higher entry pressures. Bubble mobilization resulted in a different aqueous mass distribution than if no discrete gas formation was modeled, especially at higher temperatures. This was a result of bubbles moving upwards to cooler areas, then collapsing, and contaminating previously clean zones. The cooling effect also led to possible non-aqueous phase liquid (NAPL) formation which was not predicted using a model without discrete bubble formation.     

Dynamics of the gas flux from shallow gas hydrate deposits: interaction between oily hydrate bubbles and the oceanic environment

Decomposition of methane hydrates on the continental margins is a potentially significant source of atmospheric methane, but the input depends upon the poorly understood fate of the hydrocarbon bubbles rising from the sea floor. During a field trip to the Gulf of Mexico, three different seepages were imaged and analyzed. Three different imaging techniques were tried (side, front, and back illumination), of which back illumination produced the best results. The images were analyzed and the size-dependent bubble distribution, mass flux, and rise speeds determined. The total observed gas flux was 62.3×10⁻³ mol s⁻¹, primarily methane, of which a single vent produced seven times the next largest vent. Of this major vent, 50% of the bubble mass was contained in the largest bubbles, r>5500 μm. The vertical velocities demonstrated that these bubbles were heavily contaminated with oil, which was also corroborated by bubble shape and oscillation observations.