Reduced Hamiltonian formalism for general-relativistic adiabatic fluids and applications

For the dynamics of general-relativistic adiabatic fluids a reduced Hamiltonian formalism is presented. For slowly moving fluids an approximation scheme is defined. The first steps of the scheme are worked out explicitly. Radiation damping and emission, and backscattering are discussed. The usefulness of the formalism for numerical relativity is pointed out.     

西安东部夜间异常增温原因及其数值模式预报订正

利用陕西省气象监测站观测资料、NCEP/NCAR和ERA5全球再分析资料,对2019年12月10日凌晨西安东部一次夜间温度异常跃增过程进行深入分析.结果表明:①此次夜间增温过程发展迅速、局地性强、预报难度大.②500 hPa内蒙东部冷槽底后部干燥的西北气流和850 hPa干暖气团有利于产生大气下沉绝热增温.西安地区近地...     

湿位涡方程及其应用

本文在湿绝热、无摩擦条件下，推导了湿位涡方程及其近似表达式。然后通过尺度分析得到了大尺度条件下湿位涡守恒的性质，其形式与干位涡一样。由于考虑了水汽的作用，使用它较干位涡更为方便。最后，我们把湿位涡用于江南岭北3 月有无连阴雨的分析，其结果与预报员经验是一致的。     

Pre-eruptive and syn-eruptive conditions in the Black Butte, California dacite: Insight into crystallization kinetics in a silicic magma system

A series of experiments and petrographic analyses have been run to determine the pre-eruption phase equilibria and ascent dynamics of dacitic lavas composing Black Butte, a dome complex on the flank of Mount Shasta, California. Major and trace element analyses indicate that the Black Butte magma shared a common parent with contemporaneously erupted magmas at the Shasta summit. The Black Butte lava phenocryst phase assemblage (20 v.%) consists of amphibole, plagioclase (core An_77.5), and Fe–Ti oxides in a fine-grained (< 0.5 mm) groundmass of plagioclase, pyroxene, Fe–Ti oxides, amphibole, and cristobalite. The phenocryst assemblage and crystal compositions are reproduced experimentally between 890 °C and 910 °C, ≥ 300 MPa, X_H2O = 1, and oxygen fugacity = NNO + 1. This study has quantified the extent of three crystallization processes occurring in the Black Butte dacite that can be used to discern ascent processes. Magma ascent rate was quantified using the widths of amphibole breakdown rims in natural samples, using an experimental calibration of rim development in a similar magma at relevant conditions. The majority of rims are 34 ± 10 μm thick, suggesting a time-integrated magma ascent rate of 0.004–0.006 m/s among all four dome lobes. This is comparable to values for effusive samples from the 1980 Mount St. Helens eruption and slightly faster than those estimated at Montserrat. A gap between the compositions of plagioclase phenocryst cores and microlites suggests that while phenocryst growth was continuous throughout ascent, microlite formation did not occur until significantly into ascent. The duration of crystallization is estimated using the magma reservoir depth and ascent rate, as determined from phase equilibria and amphibole rim widths, respectively. Textural analysis of the natural plagioclase crystals yields maximum growth rates of plagioclase phenocryst rims and groundmass microlites of 8.7 × 10^− 8 and 2.5 × 10^− 8 mm/s, respectively. These rates are comparable to values determined from time-sequenced samples of dacite erupted effusively from Mount St. Helens during 1980 and indicate that syneruptive crystallization processes were important during the Black Butte eruptive cycle.     

Stepwise accumulation and ascent of magmas

One of the currently popular theories on magma ascent is that it mainly occurs by propagating hydrofractures (dykes) and that magma viscosity is the primary rate‐controlling factor. This theory is based on mathematical models for single hydrofractures under idealised conditions. We simulated magma ascent with air ascending through gelatine and observed that the air ascended in batches, following paths made by their predecessors. Multiple batches accumulate at obstacles along the path. Although magma viscosity may control ascent rate during movement, obstacles ultimately control the size and average ascent velocity of ascending batches. We propose that step‐wise movement of magma batches is the mechanism of primary accumulation and ascent from the partially molten source rock of a magma to its first emplacement site and therefore the main ascent mechanism for granitic magmas. ‘Classical’ dyking is the mechanism for secondary ascent from a magma chamber.     

Adiabatic invariants for the nonconservative Kepler's problem

This paper considers adiabatic invariants for the classical Kepler problem with resisting forces. The analysis is based on the theory of integrating factors and theory of adiabatic invariants in the Krylov-Bogoliubov-Mitropolski variables. The adiabatic invariants are series with respect to a small parameter. Also, for every particular case of nonconservative forces, it is shown that, with a complete set of adiabatic invariants, an approximate solution of the problem can be obtained. Four problems are analyzed in detail where approximate solutions are compared with numerical.     

Forced secondary resonances at the 3 : 1 commensurability

For the 3 : 1 Jovian resonance problem, the time scales of the two degrees of freedom of the resonant Hamiltonian are well-separated [5]. With the adiabatic approximation, the solution for the fast oscillations can be found in terms of the slowly varying variables. Thus the rapidly oscillating terms in the slow oscillation equations can be treated as forced terms. We refer to the resonance between the forcing and intrinsic frequencies as a forced secondary one in this paper. We discuss the forced secondary resonances in asteroidal motion at the 3 : 1 commensurability by using Wisdom's method. The results show that the orbits situated originally near the resonance will leave the neighbourhood of resonance and tend to the separatrices and critical points for different energies, respectively. We have not found any stochastic web as expected in this case. Moreover, we study the problem of validity on the approximation of a system.The Project Supported by the National Natural Science Foundation of China.     

IMPACT OF DIABATIC PROCESSES IN AGCM ON 4-DIMENSIONAL VARIATIONAL DATA ASSIMILATION*

The impact of diabatic processes on 4-dimensional variational data assimilation (4D-Var) was studied using the 1995 version of NCEP's global spectral model with and without full physics.The adjoint was coded manually.A cost function measuring spectral errors of 6-hour forecasts to "observation" (the NCEP reanalysis data) was minimized using the L-BFGS (the limited memory quasi-Newton algorithm developed by Broyden,Fletcher,Goldfard and Shanno) for optimizing parameters and initial conditions.Minimization of the cost function constrained by an adiabatic version of the NCEP global model converged to a minimum with a significant amount of decrease in the value of the cost function.Minimization of the cost function using the diabatic model, however,failed after a few iterations due to discontinuities introduced by physical parameterizations.Examination of the convergence of the cost function in different spectral domains reveals that the large-scale flow is adjusted during the first 10 iterations,in which discontinuous diabatic parameterizations play very little role.The adjustment produced by the minimization gradually moves to relatively smaller scales between 10-20th iterations.During this transition period,discontinuities in the cost function produced by "on-off" switches in the physical parameterizations caused the cost function to stay in a shallow local minimum instead of continuously decreasing toward a deeper minimum. Next,a mixed 4D-Var scheme is tested in which large-scale flows are first adiabatically adjusted to a sufficient level,followed by a diabatic adjustment introduced after 10 to 20 iterations. The mixed 4D-Var produced a closer fit of analysis to observations,with 38% and 41% more decrease in the values of the cost function and the norm of gradient,respectively,than the standard diabatic 4D-Var,while the CPU time is reduced by 21%.The resulting optimal initial conditions improve the short-range forecast skills of 48-hour statistics.The detrimental effect of parameterization discontinuities on minimization was also reduced.     

Adiabatic indices in a convecting anisotropic plasma

Adiabatic indices for a non-dissipative anisotropic convecting plasma are analyzed, and general expressions for the effective adiabatic index and the partial adiabatic indices parallel (γ_∥) and perpendicular (γ_⟂) to the magnetic field are obtained. It is shown that, in the general case, the value of the effective adiabatic index is not an universal constant and depends on the plasma temperature anisotropy and on the properties of the plasma motion. The values of γ_⟂ and γ_∥ are shown to be independent of the plasma parameters being completely determined by the characteristics of the plasma flow.     

黄东海大气边界层高度时空变化特征

根据2006-2007年大连、青岛和台州逐日高分辨率L波段二次测风雷达探空资料、地面观测资料以及再分析资料,采用干绝热法和位温梯度法分别计算了各个台站的边界层高度,继而分析了黄东海边界层高度多时间尺度变化规律和空间演变特征.结果表明:(1)用清晨探空资料计算每日边界层高度,干绝热法比位温梯度法更为合理可靠;(2)沿海台站边界层高度有显著季节变化,夏季最低,秋季最高,这有别于陆地边界层普遍存在的夏季高、秋冬季节低的季节变化;(3)通过小波分析,发现边界层高度有显著的2～4 d天气尺度振荡和10～20 d准双周振荡;(4)黄东海上空海洋大气边界层在时间尺度上同样存在着明显的春夏季节低、秋冬季节高的季节变化特征,大部分海区空间上呈现西北低、东南高的分布形态.