19802000北京市农业土地生产性的变动分析

根据1980、1990、2000年的统计资料和实地走访观察,本文分析了北京市乡镇一级的粮食作物、蔬菜作物和果品作物的分布格局及农业土地生产性近20年的变化。研究认为:1980年北京农业生产分布格局与土地生产性为比较典型的杜能模式;随着农业经营的多样化,1990年和2000年北京农业土地生产性呈现出多样的分布格局;北京粮食作物、蔬菜作物、果品作物的分布格局及农业土地生产性近20年的变化与此阶段的农业生产政策、农业以外的社会经济活动变化有较强的关联性;对北京以上问题的研究为验证杜能、辛克利亚、布莱昂特城市周边农业地域研究经典理论的演变提供了实证案例。     

A multidimensional magnetohydrodynamics code based on semidiscrete central WENO and PPM schemes

We present a multidimensional magnetohydrodynamics code based on semidiscrete central Godunov-type schemes and high order weighted essentially non-oscillatory data reconstruction and incorporate elements of the piecewise parabolic method into the reconstruction schemes. This algorithm avoids spectral decomposition into characteristic waves and is easily implemented. We have implemented this algorithm and carried out one and two dimensional tests, which verifies that this code can achieve the accuracy and efficiency to solve the astrophysical problems.     

频率域全波形反演方法研究进展

全波形反演是一种利用地震波传播的动力学特征来获取地下介质物性参数的反演方法,可为揭示地下精细结构提供重要依据。本文以弹性波方程作为数学模型来模拟地震波传播规律并进行相应的反演方法研究。为提高计算效率与反演结果的准确性,可将近似解析离散化（NAD）算子用于频率域弹性波方程的正演模拟。本文在频率域NAD离散的基础上推导阻抗矩阵的稀疏分块结构与反演目标函数对模型参数的梯度计算公式,由此建立基于NAD算子的频率域弹性波全波形反演方法。为验证该方法的有效性,文中通过数值实验对多种典型介质模型进行反演计算,均得到了理想的反演结果。

     

福建中南部沿岸海域潮汐特征分析

利用2004、2007、2010年在福建中南部沿海长期验潮站和临时验潮站获取的实测逐时潮位资料,采用基于最小二乘原理的潮汐调和分析方法,分析了福建中南部沿岸海域18个验潮站的潮汐调和常数,并根据调和常数分析了沿岸海域潮汐特征.结果表明:福建中南部沿岸海域以向南传播的半日分潮为主,传播至浮头湾以南海域后衰减明显,对应平潭至浮头湾海域的潮汐性质为半日潮,古雷半岛以南海域为不规则半日潮混合潮;福建中南部沿岸海域自北部平潭向南至诏安,平均海平面升高约25 cm,潮差减小超过3 m,平均大潮差与平均小潮差之差也减小2 m左右,理论高潮位降低约1.5m,而理论低潮位升高约2 m,潮时推迟约2h.     

厦蓉高速公路某隧道洞口段开挖对仰坡变形影响的数值分析

针对一座浅埋偏压隧道,采用FLAC3D对该隧道进口段进洞开挖进行动态施工三维数值模拟。基于围岩应力分布特征,仰坡坡面轴向和横向位移的变化特征,分析了偏压浅埋隧道洞口段开挖引起的仰坡变形规律。计算结果揭示：仰坡后缘下沉,前缘向洞心外有移动趋势;隧道开挖引起隧洞洞身附近岩体出现较大应力集中和变形现象,洞口段洞身以上仰坡坡面主要以竖向沉降为主,洞身两侧向洞内挤压。在此基础上提出了保证仰坡稳定和安全进洞的一些建议。     

腾格里沙漠南部格状沙丘的形态演变及移动特征

沙丘的形态变化与移动蕴含区域风沙环境和地貌演化的关键信息,是风沙地貌研究的重要内容。以腾格里沙漠南缘的长格状和方格状沙丘为研究对象,利用风况资料和Google Earth卫星影像,监测2009—2020年两种格状沙丘的形态变化并分析其移动特征。结果表明：(1)腾格里沙漠西南缘和东南缘的主风向均为西北风,长格状沙丘分布区的次风向为东南风,方格状沙丘分布区次风向为东风和东南风,都属于低风能环境、中变率风况。西南部风能环境大于东南部,研究区近10年风动力呈衰减趋势。(2)格状沙丘的主副梁长度和间距在增加,其中,长格状沙丘高度增加,方格状沙丘高度在降低。沙丘主梁向东偏移,副梁向南偏移,形态整体保持稳定。(3)长格状沙丘平均移动速率为1.57～1.71 m·a^-1,方格状沙丘平均移动速率为1.63～2.01 m·a^-1,沙丘平均移动方向与合成输沙方向基本一致,沙丘的体积是造成移动速率差异的主因。     

中国海高分辨率业务化风暴潮模式的业务化预报检验

国家海洋环境预报中心建立了中国海高分辨率风暴潮数值预报模式,模式在水平分辨率和网格的嵌套方面都较以往的业务化模式有改进和提高.自2003年起将模式投入业务化运行以来,连续三年共对11个台风风暴潮过程进行了跟踪预报,并将数值预报结果与实测资料相对比.本文将对预报模式三年来的预报结果进行检验.     

Sedimentary evolution of the Holocene subaqueous clinoform off the southern Shandong Peninsula in the Western South Yellow Sea

Based on the stratigraphic sequence formed since the last glaciation and revealed by 3000 km long high-resolution shallow seismic profiles and the core QDZ03 acquired recently off the southern Shandong Peninsula, we addressed the sedimentary characteristics of a Holocene subaqueous clinoform in this paper. Integrated analyses were made on the core QDZ03, including sedimentary facies, sediment grain sizes, clay minerals, geochemistry, micro paleontology, and AMS 14 C dating. The result indicates that there exists a Holocene subaqueous clinoform, whose bottom boundary generally lies at 15–40 m below the present sea level with its depth contours roughly parallel to the coast and getting deeper seawards. The maximum thickness of the clinoform is up to 22.5 m on the coast side, and the thickness contours generally spread in a banded way along the coastline and becomes thinner towards the sea. At the mouths of some bays along the coast, the clinoform stretches in the shape of a fan and its thickness is evidently larger than that of the surrounding sediments. This clinoform came into being in the early Holocene(about 11.2 cal kyr BP) and can be divided into the lower and upper depositional units(DU 2 and DU 1, respectively). The unit DU 2, being usually less than 3 m in thickness and formed under a low sedimentation rate, is located between the bottom boundary and the Holocene maximum flooding surface(MFS), and represents the sediment of a post-glacial transgressive systems tract; whereas the unit DU 1, the main body of the clinoform, sits on the MFS, belonging to the sediment of a highstand systems tract from middle Holocene(about 7–6 cal kyr BP) to the present. The provenance of the clinoform differs from that of the typical sediments of the Yellow River and can be considered as the results of the joint contribution from both the Yellow River and the proximal coastal sediments of the Shandong Peninsula, as evidenced by the sediment geochemistry of the core. As is controlled mainly by coactions of multiple factors such as the Holocene sea-level changes, sediment supplies and coastal dynamic conditions, the development of the clinoform is genetically related with the synchronous clinoform or subaqueous deltas around the northeastern Shandong Peninsula and in the northern South Yellow Sea in the spatial distribution and sediment provenance, as previously reported, with all of them being formed from the initial stage of the Holocene up to the present.     

Vortex-Induced Vibrations of A Long Flexible Cylinder in Linear and Exponential Shear Flows

A numerical study based on a wake oscillator model was conducted to determine the response performance of vortex-induced vibration(VIV) on a long flexible cylinder with pinned-pinned boundary conditions subjected to linear and exponential shear flows. The coupling equations of a structural vibration model and wake oscillator model were solved using a standard central finite difference method of the second order. The VIV response characteristics including the structural displacement, structural frequency, structural wavenumber, standing wave behavior,travelling wave behavior, structural velocity, lift force coefficient and transferred energy from the fluid to the structure with different flow profiles were compared. The numerical results show that the VIV displacement is a combination of standing waves and travelling waves. For linear shear flow, standing waves and travelling waves dominate the VIV response within the low-velocity and high-velocity zones, respectively. The negative values of the transferred energy only occur within the low-velocity zone. However, for exponential shear flow, travelling waves dominate the VIV response and the negative energy occurs along the entire length of the cylinder.     

The Magnetic Topological Structure and Energy of the 2B/X2 Flare in NOAA 8100

In this paper, we reconstruct the finite energy force-free magnetic field of the active region NOAA 8100 on 4 November 1997 above the photosphere. In particular, the 3-D magnetic field structures before and after a 2B/X2 flare at 05:58 UT in this region are analyzed. The magnetic field lines were extrapolated in close coincidence with the Yohkoh soft X-ray (SXR) loops accordingly. It is found that the active region is composed of an emerging flux loop, a complex loop system with differential magnetic field shear, and large-scale, or open field lines. Similar magnetic connectivity has been obtained for both instants but apparent changes of the twisting situations of the calculated magnetic field lines can be observed that properly align with the corresponding SXR coronal loops. We conclude that this flare was triggered by the interaction of an emerging flux loop and a large loop system with differential magnetic field shear, as well as large-scale, or open field lines. The onset of the flare was at the common footpoints of several interacting magnetic loops and confined near the footpoints of the emerging flux loop. The sheared configuration remained even after the energetic flare, as demonstrated by calculated values of the twist for the loop system, which means that the active region was relaxed to a lower energy state but not completely to the minimum energy state (two days later another X-class flare occurred in this region).