GPS/LEO掩星观测的变分同化技术   总被引：6，自引：0，他引：6  

刘敏  郭鹏 《天文学进展》2006,24(1):27-42

在简单介绍GPS/LEO掩星探测大气的发展历史和科学意义之后,详细阐述了反演的基本原理;分析了标准反演中存在的问题,并说明一维变分同化(1DVAR)在反演方法中的重要性;给出了一维变分同化中价值函数的求解,以及各种同化因子;简单介绍了对当前气象学中普遍使用的四维变分同化(4DVAR);重点讨论了各种同化方法,以及使用各种同化因子的优缺点。最后,通过CHAMP卫星的观测实例分析,验证了GPS数据在数值天气预报(NWP)中的作用,以及相对于标准反演法一维变分对气象要素的改进。  相似文献   

Strain localization in sand: an overview of the experimental results obtained in Grenoble using stereophotogrammetry     

Jacques Desrues  Gioacchino Viggiani 《国际地质力学数值与分析法杂志》2004,28(4):279-321

Experimental results are presented from the extensive program of drained plane strain compression tests on sand carried out in Grenoble over the last two decades. Systematic analysis of photographs of the deforming specimen allowed for measuring deformations and determining strain fields throughout the test, that is: prior to, at, and after the onset of strain localization. The principles, details and accuracy of the procedure are described, as well as its suitability to properly depict the patterns of deformation. Findings concerning the occurrence and progression of strain localization are discussed. The issues of shear band orientation and thickness are addressed, as well as temporary and persistent complex localization patterns, and the volumetric behaviour inside a band after its formation. The influence of such variables as initial state of the sand (effective stress and relative density), specimen size and slenderness, as well as grain size, is discussed. Copyright © 2004 John Wiley & Sons, Ltd  相似文献   

Bond rolling resistance and its effect on yielding of bonded granulates by DEM analyses     

M. J. Jiang  H. S. Yu  D. Harris 《国际地质力学数值与分析法杂志》2006,30(8):723-761

A discrete element modelling of bonded granulates and investigation on the bond effect on their behaviour are very important to geomechanics. This paper presents a two‐dimensional (2‐D) discrete element theory for bonded granulates with bond rolling resistance and provides a numerical investigation into the effect of bond rolling resistance on the yielding of bonded granulates. The model consists of mechanical contact models and equations governing the motion of bonded particles. The key point of the theory is that the assumption in the original bond contact model previously proposed by the authors (55th CSCE‐ASCE Conference, Hamilton, Ont., Canada, 2002; 313–320; J. Eng. Mech. (ASCE) 2005; 131 (11):1209–1213) that bonded particles are in contact at discrete points, is here replaced by a more reliable assumption that bonded particles are in contact over a width. By making the idealization that the bond contact width is continuously distributed with the normal/tangential basic elements (BE) (each BE is composed of spring, dashpot, bond, slider or divider), we establish a bond rolling contact model together with bond normal/tangential contact models, and also relate the governing equations to local equilibrium. Only one physical parameter β needs to be introduced in the theory in comparison to the original bond discrete element model. The model has been implemented into a 2‐D distinct element method code, NS2D. Using the NS2D, a total of 86 1‐D, constant stress ratio, and biaxial compressions tests have been carried out on the bonded granular samples of different densities, bonding strengths and rolling resistances. The numerical results show that: (i) the new theory predicts a larger internal friction angle, a larger yielding stress, more brittle behaviour and larger final broken contact ratio than the original bond model; (ii) the yielding stress increases nonlinearly with the increasing value of β, and (iii) the first‐yield curve (initiation of bond breakage), which define a zone of none bond breakage and which shape and size are affected by the material density, is amplified by the bond rolling resistance in analogous to that predicted by the original bond model. Copyright © 2006 John Wiley & Sons, Ltd.  相似文献   

Studying Hydrodynamic Instability Using Shock-Tube Experiments     

O. Sadot  K. Levy  A. Yosef-Hai  D. Cartoon  Y. Elbaz  Y. Srebro  G. Ben-Dor  D. Shvarts 《Astrophysics and Space Science》2005,298(1-2):305-312

The hydrodynamic instability, which develops on the contact surface between two fluids, has great importance in astrophysical phenomena such as the inhomogeneous density distribution following a supernova event. In this event acceleration waves pass across a material interface and initiate and enhance unstable conditions in which small perturbations grow dramatically. In the present study, an experimental technique aimed at investigating the above-mentioned hydrodynamic instability is presented. The experimental investigation is based on a shock-tube apparatus by which a shock wave is generated and initiates the instability that develops on the contact surface between two gases. The flexibility of the system enables one to vary the initial shape of the contact surface, the shock-wave Mach number, and the density ratio across the contact surface. Three selected sets of shock-tube experiments are presented in order to demonstrate the system capabilities: (1) large-initial amplitudes with low-Mach-number incident shock waves; (2) small-initial amplitudes with moderate-Mach-number incident shock waves; and (3) shock bubble interaction. In the large-amplitude experiments a reduction of the initial velocity with respect to the linear growth prediction was measured. The results were compared to those predicted by a vorticity-deposition model and to previous experiments with moderate- and high-Mach number incident shock waves that were conducted by others. In this case, a reduction of the initial velocity was noted. However, at late times the growth rate had a 1/t behavior as in the small-amplitude low-Mach number case. In the small-amplitude moderate-Mach number shock experiments a reduction from the impulsive theory was noted at the late stages. The passage of a shock wave through a spherical bubble results in the formation of a vortex ring. Simple dimensional analysis shows that the circulation depends linearly on the speed of sound of the surrounding material and on the initial bubble radius.  相似文献   

Experimental solidification of anhydrous latitic and trachytic melts at different cooling rates: The role of nucleation kinetics   总被引：1，自引：0，他引：1  

Gianluca Iezzi  Silvio Mollo  Guido Ventura  Andrea Cavallo  Claudia Romano 《Chemical Geology》2008,253(3-4):91-101

Two sets of cooling experiments were run at atmospheric conditions for two anhydrous starting latitic and trachytic melts: 1) five cooling rates (25, 12.5, 3, 0.5, and 0.125 °C/min) between 1300° and 800 °C, and 2) a 0.5 °C/min cooling rate from 1300 °C with quench temperatures at 1200°, 1100°, 1000° and 900 °C. Trachytic run-products are invariably glassy. Nucleation is also suppressed in the latitic run-products at the three highest cooling rates. Conversely, in the 0.5 and 0.125 °C/min runs, latites have a crystal content of  90 vol.%. The phases are: plagioclase, clinopyroxene, glass and iron-bearing oxide (in order of abundance). The variable quench temperatures, investigated by coupling experiments with Pt wire and Pt capsule sample containers in set 2, again did not produce crystallization of trachyte, whereas latitic samples are characterized by  10 vol.% of oxides, pyroxenes and plagioclase (in order of appearance), at temperature < 1000 °C. Effects of (preferential) heterogeneous nucleation on sample holders, of superheating degree, and chemical species loss during cooling are absent for both melt compositions. The difference of solidification paths between these two silicate melts can be ascribed only to their small chemical differences. In comparison with calculated equilibrium conditions all the experimental latitic and trachytic run-products revealed strong kinetic effects, interpretable in the light of the nucleation theory. The glass-forming ability (GFA) of trachyte is higher, whereas their critical cooling rate (Rc) is lower (< 0.125 °C/min), in comparison to latitic melts (Rc > 0.5 °C/min). The experimental results carried out in this study can be applied to lava flows and domes; trachytic lavas are able to flow for longer period with respect to latitic ones in a metastable condition. Glass-rich terrestrial lavas, i.e. obsidians, can be the result of sluggish nucleation kinetics due to the relative high polymerisation of evolved silicate melts.  相似文献   

The Ca-Eskola component in eclogitic clinopyroxene as a function of pressure, temperature and bulk composition: an experimental study to 15 GPa with possible implications for the formation of oriented SiO2-inclusions in omphacite     

Jürgen Konzett  Daniel J. Frost  Alexander Proyer  Peter Ulmer 《Contributions to Mineralogy and Petrology》2008,155(2):215-228

Experiments have been conducted in the P-T range 2.5–15 GPa and 850–1,500°C using bulk compositions in the systems SiO₂–TiO₂–Al₂O₃–Fe₂O₃–FeO–MnO–MgO–CaO–Na₂O–K₂O–P₂O₅ and SiO₂–TiO₂–Al₂O₃–MgO–CaO–Na₂O to investigate the Ca-Eskola (CaEs Ca_0.5□_0.5AlSi₂O₆) content of clinopyroxene in eclogitic assemblages containing garnet + clinopyroxene + SiO₂ ± TiO₂ ± kyanite as a function of P, T, and bulk composition. The results show that CaEs_ss in clinopyroxene increases with increasing T and is strongly bulk composition dependent whereby high CaEs-contents are favoured by bulk compositions with high normative anorthite and low diopside contents. In this study, a maximum of 18 mol% CaEs_ss was found at 6 GPa and 1,350°C in a kyanite-eclogite assemblage garnet + clinopyroxene + kyanite + rutile + coesite. By comparison, no significant increase in CaEs_ss with increasing P could be observed. If the formation of oriented SiO₂-rods frequently observed in eclogititc clinopyroxenes is due to the retrogressive breakdown of a CaEs-component then these textures are a cooling rather than a decompression phenomenon and are most likely to be found in kyanite-bearing eclogites cooled from temperatures ≥750°C. The presence of clinopyroxene with approx. 4 mol% CaEs_ss in an experiment conducted at 2.5 GPa/850°C confirms earlier suggestions based on field data that vacancy-rich clinopyroxenes are not necessarily restricted to ultrahigh pressure metamorphic conditions. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.  相似文献   

The impact of the M₂ internal tide on data-assimilative model estimates of the surface tide     

E.D. Zaron  G.D. Egbert   《Ocean Modelling》2007,18(3-4):210-216

We use a synthetic data experiment to assess the accuracy of ocean tides estimated from satellite altimetry data, with emphasis on the impact of the phase-locked internal tide, which has a surface expression of several centimeters near its sites of genesis. Previous tidal estimates have regarded this signal as a random measurement error; however, it is deterministic and not scale-separated from the barotropic (surface) tide around complex bathymetric features. The synthetic data experiments show that the internal tide has a negligible impact on the barotropic tidal fields inferred under these circumstances, and the barotropic dissipation (a quadratic functional of the tidal fields) is in good agreement with the energetics of the three-dimensional regional primitive equations model which is the source of the synthetic data.  相似文献   

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

Konstantin V. Lebedev  Max Yaremchuk  Humio Mitsudera  Iwao Nakano  Gang Yuan 《Journal of Oceanography》2003,59(6):751-763

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.  相似文献   

Data assimilation for a coastal area morphodynamic model: Morecambe Bay     

T.R. Scott  D.C. Mason 《Coastal Engineering》2007

This paper investigates the use of data assimilation in coastal area morphodynamic modelling using Morecambe Bay as a study site. A simple model of the bay has been enhanced with a data assimilation scheme to better predict large-scale changes in bathymetry observed in the bay over a 3-year period. The 2DH decoupled morphodynamic model developed for the work is described, as is the optimal interpolation scheme used to assimilate waterline observations into the model run. Each waterline was acquired from a SAR satellite image and is essentially a contour of the bathymetry at some level within the inter-tidal zone of the bay. For model parameters calibrated against validation observations, model performance is good, even without data assimilation. However the use of data assimilation successfully compensates for a particular failing of the model, and helps to keep the model bathymetry on track. It also improves the ability of the model to predict future bathymetry. Although the benefits of data assimilation are demonstrated using waterline observations, any observations of morphology could potentially be used. These results suggest that data assimilation should be considered for use in future coastal area morphodynamic models.  相似文献   

Kuroshio extension variability explored through assimilation of TOPEX/POSEIDON altimeter data into a quasi-geostrophic model     

Yousuke Nishihama  Motoyoshi Ikeda 《Journal of Oceanography》2007,63(6):879-895

Mesoscale features in the eastward extension of the Kuroshio were investigated using assimilation of TOPEX/POSEIDON (T/P) data into a three-layer quasi-geostrophic model. The T/P data exhibited an elongated state of the southern recirculation gyre in 1993–95 and 1997, between whose two periods the gyre had a contracted state in 1995–96. A few stationary eddies were located in the southern gyre during the contracted state. The baroclinic instability, which was indicated by the phase shift from the uppermost-to the lowest-layer anomalies toward the downstream side, was evident near the Kuroshio Extension (KE) path. Since the instability never appeared in the artificial model without bottom topography, the topographic barrier for the eastward flow in the lowest layer was a necessary condition for the instability. The instability synchronized with the transition in the western region of the KE axis from the elongated to the contracted states. This evolution was interpreted as if the baroclinic instability played some part in the KE states and was a trigger for the transition from the elongated to the contracted states.  相似文献