The effect of a shallow low viscosity zone on the apparent compensation of mid-plate swells

A simple model for mid-plate swells is that of convection in a fluid which has a low viscosity layer lying between a rigid bed and a constant viscosity region. Finite element calculations have been used to determine the effects of the viscosity contrast, the layer thickness and the Rayleigh number on the flow and on the perceived compensation mechanism for the resulting topographic swell. As the viscosity decreases in the low viscosity zone, the effective local Rayleigh number for the top boundary layer of the convecting cell increases. Also, because the lower viscosity facilitates greater velocities in the low viscosity zone, the low viscosity layer produces proportionally greater horizontal flow near the conducting lid, causing the base of the conducting lid to appear like a free boundary. The change in the local Rayleigh number and in the effective boundary condition both cause the top boundary layer to thin. Through a Green's function analysis, we have found that the low viscosity zone damps the response of the surface topography to the temperature anomalies at depth, whereas it causes the gravity and geoid response functions to change sign at depth counteracting the positive contributions from the shallower temperature variations. By increasing the viscosity contrast, the conbined effects of the thinning of the boundary layer and the behaviour of the response functions allow the apparent depth of compensation to become arbitrarily small. Therefore, shallow depths of compensation cannot be used to argue against dynamic support of mid-plate swells. Furthermore, we compared the distribution of the effective compensating densities, which is used to obtain the geoid, to that of Pratt compensation, which is often used to calculate the depth of compensation from geoid and topography data for mid-plate swells. For all of our calculations including those with no low viscosity layer, the effective gravitational mass distribution is more complex than assumed in simple Pratt models, so that the Pratt models are not an appropriate gauge of the compensation mechanism.     

利用多种地球物理观测资料直接反演地幔对流模型

假定地幔为一个均匀的、粘滞系数为常数、同时均匀分布放射性热源的流体球层,其内部存在的对流则由流体力学3个基本方程:运动方程、能量方程和连续性方程确定.如果假定地幔处于低瑞利数的状态（临界瑞利数1.5倍左右）,那么上述方程中的非线性项可以忽略不计.作为一类可能的模型,本文计算一组用6个边界条件确定6个未知数的线性方程组.这些条件包括板块绝对运动极型场、地球大地水准面异常和地震层析结果提供的地幔密度分布横向不均匀相应的“刚性地球”水准面异常等.模型计算表明:1.地幔中流体运动格局不仅受地幔热动力学参数（瑞利数）控制,而且强烈地受边界条件的影响.2.若不限定下边界为等温边界,则上、下地幔之间并不呈现出活动性明显差异;但是在模型瑞利数加大到一定值时,核-幔边界附近将出现一些局部的小尺度对流环.3.当模型瑞利数从很小增加时,对流格局将发生变化,这些格局可能反应由地幔热动力学参数决定的地幔固有特性.4.当瑞利数为50000和80000时,核-幔边界形变与PcP波得到的结果吻合较好.     

Dynamos with weakly convecting outer layers: implications for core-mantle boundary interaction

Convection in the Earth's core is driven much harder at the bottom than the top. This is partly because the adiabatic gradient steepens towards the top, partly because the spherical geometry means the area involved increases towards the top, and partly because compositional convection is driven by light material released at the lower boundary and remixed uniformly throughout the outer core, providing a volumetric sink of buoyancy. We have therefore investigated dynamo action of thermal convection in a Boussinesq fluid contained within a rotating spherical shell driven by a combination of bottom and internal heating or cooling. We first apply a homogeneous temperature on the outer boundary in order to explore the effects of heat sinks on dynamo action; we then impose an inhomogeneous temperature proportional to a single spherical harmonic Y ₂² in order to explore core-mantle interactions. With homogeneous boundary conditions and moderate Rayleigh numbers, a heat sink reduces the generated magnetic field appreciably; the magnetic Reynolds number remains high because the dominant toroidal component of flow is not reduced significantly. The dipolar structure of the field becomes more pronounced as found by other authors. Increasing the Rayleigh number yields a regime in which convection inside the tangent cylinder is strongly affected by the magnetic field. With inhomogeneous boundary conditions, a heat sink promotes boundary effects and locking of the magnetic field to boundary anomalies. We show that boundary locking is inhibited by advection of heat in the outer regions. With uniform heating, the boundary effects are only significant at low Rayleigh numbers, when dynamo action is only possible for artificially low magnetic diffusivity. With heat sinks, the boundary effects remain significant at higher Rayleigh numbers provided the convection remains weak or the fluid is stably stratified at the top. Dynamo action is driven by vigorous convection at depth while boundary thermal anomalies dominate in the upper regions. This is a likely regime for the Earth's core.     

Heat flow and depth versus age for the Mesozoic northwest Atlantic Ocean: results from the Sohm abyssal plain and implications for the Bermuda Rise

Heat flow in the Sohm abyssal plain is measured to be 53 mW/m² at an age of 163 Ma. This is 25% higher than predicted by conductive cooling models, even though the sediment-corrected basement depth of 6.5 km at this location is normal for its age. An analysis of existing heat flow, depth and geoid anomalies in the northwest Atlantic shows that there is little correlation between heat flow and depth throughout the entire region. Depth and geoid are clearly related to the Bermuda swell while the associated heat flow anomaly, once adjusted for variations with age, is limited to 5 mW/m² and only decays to the south. This means that the Bermuda swell is probably not caused by extensive thermal reheating within the lithosphere, but instead by dynamic uplift at its lower boundary due to the convective upwelling of a mantle plume. The regionally high heat flow in the northwest Atlantic may be a thermal remanent of previous plumes which passed beneath this region early in its history. Therefore, depth and heat flow anomalies from this region cannot be used to provide constraints on steady-state parameters of the lithosphere, such as the presence or absence of a long-term boundary layer at its base.     

Convective instabilities in a variable viscosity fluid cooled from above

Whether in the mantle or in magma chambers, convective flows are characterized by large variations of viscosity. We study the influence of the viscosity structure on the development of convective instabilities in a viscous fluid which is cooled from above. The upper and lower boundaries of the fluid are stress-free. A viscosity dependence with depth of the form ν₀ + ν₁ exp(?γ.z) is assumed. After the temperature of the top boundary is lowered, velocity and temperature perturbations are followed numerically until convective breakdown occurs. Viscosity contrasts of up to 10⁷ and Rayleigh numbers of up to 10⁸ are studied.For intermediate viscosity contrasts (around 10³), convective breakdown is characterized by the almost simultaneous appearance of two modes of instability. One involves the whole fluid layer, has a large horizontal wavelength (several times the layer depth) and exhibits plate-like behaviour. The other mode has a much smaller wavelength and develops below a rigid lid. The “whole layer” mode dominates for small viscosity contrasts but is suppressed by viscous dissipation at large viscosity contrasts.For the “rigid lid” mode, we emphasize that it is the form of the viscosity variation which determines the instability. For steep viscosity profiles, convective flow does not penetrate deeply in the viscous region and only weak convection develops. We propose a simple method to define the rigid lid thickness. We are thus able to compute the true depth extent and the effective driving temperature difference of convective flow. Because viscosity contrasts in the convecting region do not exceed 100, simple scaling arguments are sufficient to describe the instability. The critical wavelength is proportional to the thickness of the thermal boundary layer below the rigid lid. Convection occurs when a Rayleigh number defined locally exceeds a critical value of 160–200. Finally, we show that a local Rayleigh number can be computed at any depth in the fluid and that convection develops below depth z_r (the rigid lid thickness) such that this number is maximum.The simple similarity laws are applied to the upper mantle beneath oceans and yield estimates of 5 × 10¹⁵?5 × 10¹⁶ m² s^?1 for viscosity in the thermal boundary layer below the plate.     

基于背景噪声和地震面波联合反演华北克拉通中部岩石圈结构

基于ChinArray三期项目布设于华北克拉通中部的流动台阵观测数据,利用背景噪声互相关和地震面波层析成像获取了研究区内6—140 s周期的瑞雷面波频散,使用蒙特卡罗非线性反演方法获得了华北克拉通中部岩石圈的高分辨率三维S波速度结构。结果显示华北克拉通不同地块的岩石圈速度结构存在显著的横向差异：其中鄂尔多斯盆地腹地整体表现为高速特征,延伸至200 km以下,但其东南缘存在小范围的低速异常;东部的华北盆地整体表现为低速特征,具有较薄的地壳和岩石圈厚度;中部造山带南北两端以及南北重力梯度线下方存在相连接的低速区域,在深处延伸至华北盆地下方;在下地壳和上地幔顶部,大同火山群区域的低速体逐渐向西偏移至鄂尔多斯盆地东北角下方;而在上地幔中,该区域的低速异常随深度增加而逐渐减弱,低速体延伸至东南方向的华北盆地下方。基于本研究获得的S波速度模型,我们认为：鄂尔多斯盆地腹地保持了克拉通特性,但其东南缘存在局部的岩石圈改造作用;华北盆地发生了强烈的岩石圈破坏减薄和地壳伸展变形;中部造山带南北端以及南北重力梯度线下方的岩石圈发生了局部的改造减薄,其机制可能都来源于华北盆地下方地幔热物质的上涌;大同火山群下...     

蒙古及周边地区重力异常和地壳不均匀体分布

基于全球EGM2008自由空气重力异常模型,本文计算了蒙古及周边地区的布格重力异常和AiryHeiskanen均衡重力异常.在此基础上,本文采用Crust 1.0地壳模型为参考,通过重力正演方法,对蒙古及周边地区不同深度地壳密度结构模型的重力异常进行了计算,并对得到的正演布格重力异常与实际重力异常进行了对比和分析.研究结果表明:蒙古西部杭爱山地区与阿尔泰山地区的构造变形差异性明显,现今均衡重力异常中杭爱山周边没有明显的均衡异常高值区,而阿尔泰山地区西南方向存在均衡重力异常高值分布,分析与新构造运动密切相关;Crust 1.0模型给出的壳幔横向密度不均匀体分布对于计算Moho面起伏引起的重力异常作用明显;Crust 1.0给出的地壳内界面变形可以反映深大活动断裂的深部构造变形.研究结果对于认识蒙古东西部构造特征差异,以及现今西部活动断裂的地球物理场特征具有参考意义,也可以为进一步应用Crust 1.0模型为参考开展三维密度结构反演提供一定帮助.     

21世纪初青藏高原感热年代际增强对中国东部季风雨带关键区夏季降水年代际转折的影响

基于中国气象局提供的气象站点月值资料,NOAA、CMAP降水格点月值资料,NDVI卫星资料及再分析资料,利用统计方法分析了1961-2014年青藏高原感热与中国东部季风雨带关键区夏季降水的年代际变化,并根据热动力平衡方程结合CESM模式试验解释了21世纪初高原感热异常对关键区夏季降水的影响机理.结果表明:21世纪初,黄淮、江淮地区降水增加,而长江以南地区降水减少.同时,高原感热也发生年代际增强,当春季感热增强后,大气热能上传导致夏季高原近地面产生气旋性环流异常,大气辐合;高层产生反气旋性环流异常,大气辐散.黄淮、江淮地区在对流层中低层受异常偏南风控制,高层受高原上空的大尺度反气旋环流影响产生异常偏北风.此外,高原感热增强通过影响黄淮、江淮地区产生暖平流输送和非绝热加热正异常,该区域产生异常的上升运动,降水量增加.长江以南地区在对流层中低层存在一个异常的反气性环流,有来自海洋的冷平流输送,同时大气非绝热加热在该地区为负异常,产生异常的下沉运动,降水量减少.模式敏感性试验的结果证实了当高原感热发生年代际增强,黄淮、江淮地区水平温度平流及非绝热加热为正异常,而在华南地区为负异常,从而导致黄淮、江淮地区大气上升运动增强,降水增加;而华南地区下沉运动增强,降水减少.     

Comparison of global geopotential models from the champ and grace missions for regional geoid modelling in Turkey

The continuous efforts on establishment and modernization of the geodetic control in Turkey include a number of regional geoid models that have been determined since 1976. The recently released gravimetric Geoid of Turkey, TG03, is used in geodetic applications where GPS-heights need to be converted to the local vertical datum. To reach a regional geoid model with improved accuracy, the selection of the appropriate global geopotential model is of primary importance. This study assesses the performance of a number of recent satellite-only and combined global geopotential models (GGMs) derived from CHAMP and GRACE missions’ data in comparison to the older EGM96 model, which is the underlying reference model for TG03. In this respect, gravity anomalies and geoid heights from the global geopotential models were compared with terrestrial gravity data and low-pass filtered GPS/levelling data, respectively. Also, five new gravimetric geoid models, computed by the Fast Fourier Transform technique using terrestrial gravity data and the geopotential models, were validated at the GPS/levelling benchmarks. The findings were also compared with the validation results of the TG03 model. The tests showed that as it was expected any of the high-degree combined models (EIGEN-CG03C, EIGEN-GL04C, EGM96) can be employed for determining the gravity anomalies over Turkey. In the west of Turkey, EGM96 and EIGEN-CHAMP03S fit the GPS/levelling surface better. However, all the tested GGMs revealed equal performance when they were employed in gravimetric geoid modelling after de-trending the gravimetric geoid model with corrector surface fitting. The new geoid models have improved accuracy (after fit) compared to TG03.     

Subducted lithosphere,hotspots, and the geoid

The earth's largest positive geoid height anomalies are associated with subduction zones and hotspots. Although the correlation with subduction has been noted for many years, the correlation with hotspots is fully evident only when the subduction-related geoid highs are removed from the observed field. Using the assumption that subducted lithospheric slabs are uncompensated and are thermally re-equilibrated with the asthenosphere at the maximum depth of earthquakes, the expected geoid anomaly over subduction zones is calculated. This field provides a satis-factory fit to the observed circum-Pacific high. Subtraction of this predicted anomaly leaves a residual field which is correlated, at greater than the 99% confidence level, with the distribution of hotspots. Broad residual geoid highs occur over the central Pacific and the Africa/eastern Atlantic regions, the same areas where the hotspots are concentrated. The mass anomalies associated with hotspots and subducted slabs apparently control the location of the earth's spin axis.     

Numerical simulations of thermal convection in a rapidly rotating spherical shell cooled inhomogeneously from above

Abstract

Numerical simulations of thermal convection in a rapidly rotating spherical fluid shell with and without inhomogeneous temperature anomalies on the top boundary have been carried out using a three-dimensional, time-dependent, spectral-transform code. The spherical shell of Boussinesq fluid has inner and outer radii the same as those of the Earth's liquid outer core. The Taylor number is 10⁷, the Prandtl number is 1, and the Rayleigh number R is 5R_c (R_c is the critical value of R for the onset of convection when the top boundary is isothermal and R is based on the spherically averaged temperature difference across the shell). The shell is heated from below and cooled from above; there is no internal heating. The lower boundary of the shell is isothermal and both boundaries are rigid and impermeable. Three cases are considered. In one, the upper boundary is isothermal while in the others, temperature anomalies with (l,m) = (3,2) and (6,4) are imposed on the top boundary. The spherically averaged temperature difference across the shell is the same in all three cases. The amplitudes of the imposed temperature anomalies are equal to one-half of the spherically averaged temperature difference across the shell. Convective structures are strongly controlled by both rotation and the imposed temperature anomalies suggesting that thermal inhomogeneities imposed by the mantle on the core have a significant influence on the motions inside the core. The imposed temperature anomaly locks the thermal perturbation structure in the outer part of the spherical shell onto the upper boundary and significantly modifies the velocity structure in the same region. However, the radial velocity structure in the outer part of the shell is different from the temperature perturbation structure. The influence of the imposed temperature anomaly decreases with depth in the shell. Thermal structure and velocity structure are similar and convective rolls are more columnar in the inner part of the shell where the effects of rotation are most dominant.     

Coastal upwelling along the north coast of Papua New Guinea and El Niño events during 1981–2005

We investigate the relationship between sea surface temperature (SST) cooling and upwelling along Papua New Guinea’s (PNG) north coast before the onset of El Niño events using a hindcast experiment with a high-resolution ocean general circulation model. Coastal upwelling and related SST cooling appear along PNG north coast during the boreal winter before the onsets of six El Niño events occurring during 1981–2005. Relatively cool SSTs appear along PNG north coast during that time, when anomalous northwesterly surface wind stress, which can cause coastal upwelling by offshore Ekman transport appearing over the region. In addition, anomalous cooling tendencies of SST are observed, accompanying anomalous upward velocities at the base of the mixed layer and shallow anomalies of 27°C isotherm depth. It is also shown that entrainment cooling plays an important role in the cooling of the mixed layer temperature in this region.     

On the effect of a low viscosity asthenosphere on the temporal change of the geoid—A challenge for future gravity missions

New satellite technology to measure changes in the Earth’s gravity field gives new possibilities to detect layers of low viscosity inside the Earth. We used density models for the Earth mantle based on slab history as well as on tomography and fitted the viscosity by comparison of predicted gravity to the new CHAMP gravity model. We first confirm that the fit to the observed geoid is insensitive to the presence of a low viscosity anomaly in the upper mantle as long as the layer is thin ( 200 km) and the viscosity reduction is less than two orders of magnitude. Then we investigated the temporal change in geoid by comparing two stages of slablet sinking based on subduction history or by advection of tomography derived densities and compared the spectra of the geoid change for cases with and without a low viscosity layer, but about equal fit to the observed geoid. The presence of a low viscosity layer causes relaxation at smaller wavelength and thus leads to a spectrum with relatively stronger power in higher modes and a peak around degrees 5 and 6. Comparing the spectra to the expected degree resolution for GRACE data for a 5 years mission duration shows a weak possibility to detect changes in the Earth’s gravity field due to large scale mantle circulation, provided that other causes of geoid changes can be taken into account with sufficient accuracy. A discrimination between the two viscosity cases, however, demands a new generation of gravity field observing satellites.     

Geoid anomalies over Gorringe Ridge,North Atlantic Ocean

Gorringe Ridge is a strong uplifted block of oceanic crust and upper mantle lying at the eastern end of the Azores-Gibraltar plate boundary. The geoid over this structure derived from Seasat altimeter data exhibits a 9-m height anomaly with a north-south lateral extension smaller than 200 km. An attempt is made to interpret this geoid together with the gravity anomalies and with the seismicity, which has been compiled as a function of depth.It is first shown that the flexure of the oceanic lithosphere due to the ridge loading does not provide a good fit of the geoid anomalies and probably should be discarded, as it assumes a continuous unfractured elastic plate.Models involving local heterogeneities are then tested. The comparison of the observed geoid anomalies with the anomalies due to the uncompensated relief indicates that the topographic high has no shallow compensation.Uncompensated models, previously proposed to explain the gravity anomalies, are tested using the geoid. One model (Purdy and Bonnin, in Bonnin [11]), which involves an uplift of upper mantle material at depth, generates too strong geoid anomalies and must be discarded. Another model, which represents a nascent subduction zone (Le Pichon et al. [25]), fits both the gravity and geoid anomalies, but leads to difficulties in explaining the deep seismicity north of Gorringe Ridge.A model in isostatic equilibrium is also able to fit both gravity and geoid anomalies. This model involves a deep root of density 3.0 g cm^?3, as has been previously proposed for many oceanic ridges and plateaus. This model is compatible with the deep seismicity, but the origin of this low-density material at great depth is up to now an unresolved question.More likely, dynamical models taking into account the forces induced by the convection related to the slow plate convergence in this area will have to be considered.     

A comparison of GOCE gravitational models with EGM2008

Four new gravity field models from GOCE, two of them combined with GRACE, are compared here with EGM2008. The objectives are to look into the differences in consecutive ranges of the spherical harmonic expansion globally as well as in selected geographical regions and in the regions of the various data sources used for EGM2008. In general, GOCE is able to contribute to improved global gravity models in the spherical harmonic range between 120 and 200 (and above). The agreement between EGM2008 and the GOCE models is very good in well-surveyed regions such as North America, Europe and Australia, with geoid RMS-differences on the order of 4–6 cm. In other regions, where the surface gravity data available for the development of EGM2008 were poor, such as South America, Africa, South-East Asia or China the RMS-differences are on a level of 30 cm. Here GOCE leads to a significant improvement. These findings are confirmed by the analysis of the areas of the various EGM2008 data sources. In the regions of the so-called “fill-in” data of EGM2008 RMS-geoid height differences are high. In Antarctica GOCE also gives important improvements in terms of spatial resolution and accuracy. In general, the agreement between EGM2008 and the GOCE-models up to degree and order (d/o) 200 is good, with a global (excluding the polar gaps of GOCE orbits, throughout) geoid difference RMS of 11 cm, in the ocean areas 8 cm and 20 cm in the continental areas. GOCE models are better suited for ocean circulation studies because no prior ocean information enters into the data reduction process, as it is the case when deducing gravity anomalies from an altimetric mean sea surface. On the other hand, the good consistency between GOCE-models and EGM2008 in ocean areas very likely indicates that the influence of ocean circulation information on EGM2008 is rather small. The four tested GOCE models behave similarly except at the highest latitudes where GOCE lacks data due to its orbit inclination of 96.5° and some form of regularization which has to be applied.     

Detection of Local Geoid Deformations by Gravity Disturbances from GPS/Gravimetric Observations

This paper deals with a method for detection of local geoid deformations; as a consequence, the methods main application concerns geoid adjustment to GPS/levelling points. This is based on the fact that these points should present no local geoid deformation to avoid errors in the adjustments. These type of miscalculations would lead to an incorrect adjustment and result in further errors in subsequent studies with GPS in the proximity at the point with local deformation.The method proposed is based on predictions of gravity disturbance from geoid undulations using Poisson integral with modified kernel, and its comparison with the gravity disturbance from GPS and gravimetric observations.The use of gravity disturbance instead of gravity anomalies has been chosen since gravity disturbance is a quantity derived from GPS and not from levelling. The loss of accuracy arising with a local height reference system is therefore theoretically avoided as far as the differences in geodetic reference systems regarding positions of gravity measurements and coefficients of the global models are accounted for.Extended numerical tests using computed geoidal undulations and the corresponding gravity disturbances obtained from the geopotential model GPM98cr computed up to degree 720 illustrate the validity of the proposed method and its usefulness as local geoid deformations detection tool.Finally, the method is tested using real GPS/Gravimetric data and geoid models IBERGEO95 and EGG97 with good results.     

On Gravity Inversion for Point Mass Anomalies by Means of the Truncated Geoid

The physical meaning of the truncated geoid, which is defined by the convolution of gravity anomalies with the Stokes function on a spherical cap of specified radius, has been studied by the authors. They investigated its relation to the density distribution, generating the surface gravity, and its potential use in inversion. Some progress results for simulated studies on point mass anomalies are presented. The behavior of the truncated geoid is controlled by the radius of the integration domain, hereinafter referred to as the truncation parameter, which is treated as a free parameter. The change of the truncated geoid in response to the change of the truncation parameter was studied in the context of the simulated mass distributions. By means of such computer simulations we have managed to demonstrate the clear sensitivity of the truncated geoid to the depths, in addition to the horizontal positions, of point mass anomalies generating the synthetic surface gravity. The objective of this paper is to illustrate, with the help of computer simulation as the method of our study, the contribution of the truncated geoid to the solution of the gravimetric inverse problem. Further work towards employing the truncated geoid in gravity exploration is being conducted.     

星载SAR遥感图像反演海洋大气边界层高度

本研究将边界层相似理论与对流理论应用到具有海洋大气边界层(Marine Atmospheric Boundary Layer, MABL)对流特征的星载合成孔径雷达(Synthetic Aperture Radar, SAR)遥感图像,探讨了星载SAR遥感图像描述海气应力作用下水平扰动尺度变化的潜在可能性.针对具有三维对流涡旋Cell和二维水平滚轴涡旋Roll特征的星载SAR遥感图像,反演了中国海海域MABL高度,并与同步实验获取的MABL高度结果进行对比.结果表明,利用具有对流特征的星载SAR遥感图像反演MABL高度是可行的,展示了以高分辨率、大面积观测为特点的星载SAR遥感图像探测MABL的广阔前景.     

The effect of coastal upwelling on the sea-breeze circulation at Cabo Frio, Brazil: a numerical experiment

The effect of coastal upwelling on sea-breeze circulation in Cabo Frio (Brazil) and the feedback of sea-breeze on the upwelling signal in this region are investigated. In order to study the effect of coastal upwelling on sea-breeze a non-linear, three-dimensional, primitive equation atmospheric model is employed. The model considers only dry air and employs boundary layer formulation. The surface temperature is determined by a forcing function applied to the Earths surface. In order to investigate the seasonal variations of the circulation, numerical experiments considering three-month means are conducted: January-February-March (JFM), April-May-June (AMJ), July-August-September (JAS) and October-November-December (OND). The model results show that the sea-breeze is most intense near the coast at all the seasons. The sea-breeze is stronger in OND and JFM, when the upwelling occurs, and weaker in AMJ and JAS, when there is no upwelling. Numerical simulations also show that when the upwelling occurs the sea-breeze develops and attains maximum intensity earlier than when it does not occur. Observations show a similar behavior. In order to verify the effect of the sea-breeze surface wind on the upwelling, a two-layer finite element ocean model is also implemented. The results of simulations using this model, forced by the wind generated in the sea-breeze model, show that the sea-breeze effectively enhances the upwelling signal.     

联合多代卫星测高和多源重力数据的局部大地水准面精化方法

本文研究了基于泊松小波径向基函数融合多代卫星测高及多源重力数据精化大地水准面模型的方法.分别以沿轨垂线偏差和大地水准面高高差作为卫星测高观测量,研究了使用不同类型测高数据对于大地水准面建模精度的影响.针对全球潮汐模型在浅水区域及部分开阔海域精度较低的问题,引入局部潮汐模型研究了不同潮汐模型对于大地水准面的影响.数值分析表明:相比于使用沿轨垂线偏差作为测高观测量,基于沿轨大地水准面高高差解算得到的大地水准面模型的精度更高,特别是在海域区域,其精度提高了2.3cm.由于使用沿轨大地水准面高高差作为测高观测量削弱了潮汐模型长波误差的影响,采用不同潮汐模型对大地水准面解算的影响较小.总体而言,船载重力及测高观测数据在海洋重力场的确定中呈现互补性关系,联合两类重力场观测量可以提高局部重力场的建模精度.