南极纳尔逊冰帽的雪层及其成冰作用研究

纳尔逊冰帽雪—粒雪的演化依赖于融水渗浸冻结作用下的暖变质过程。积雪的密实化过程的快慢,取决于温度条件和融水的参与程度,以及自身的物理状况。粒雪的密实化过程表现为均匀且变幅小。纳尔逊冰帽成冰深度在23～25m,成冰历时17～19年。成冰带分为暖渗浸—重结晶带,渗浸—冻结带,消融带。     

南极洛多姆冰帽BJ钻孔粒雪层的热学性质及其温度分布

根据详细的密度资料,对洛多姆冰帽BJ钻孔粒雪层的热学性质进行了讨论。应用非均匀介质热传导方程描述活动层温度分布。通过两种简化的层状结构和竖向运动模型给出了解析解。计算和实测温度差异的分析表明边界条件的选取非常重要。     

Solar Wind Effects Associated with a Recurrent,High-Latitude Coronal Brightness Enhancement

During the descent of Ulysses following the 2001 solar north pole passage, the SOHO LASCO C2 telescope recorded a particularly strong sequence of recurrent polarization brightness (pB) features at latitudes of around 55°. As Ulysses passed overhead, solar rotation swept the interplanetary extensions of these persistent coronal structures over the spacecraft. Comparison of solar remote sensing and Ulysses in situ observations through 2002 reveals the solar wind effects of very bright and recurrent K-coronal structures at high solar latitudes and of a steeply inclined heliospheric neutral sheet (HNS). Despite the high level of solar activity, the HNS at high latitude still organizes solar wind stream structure much as it did near the previous solar minimum. The recurrent coronal streamers originate slow solar wind and mark the northern extremity of a very tilted HNS whose passage at Ulysses is accompanied by slow, dense solar wind, enhanced temperature, depressed α abundance, enhanced magnetic fields, and magnetic field directional changes that evolve with spacecraft latitude.     

The Number of Magnetic Null Points in the Quiet Sun Corona

The coronal magnetic field above a particular photospheric region will vanish at a certain number of points, called null points. These points can be found directly in a potential field extrapolation or their density can be estimated from the Fourier spectrum of the magnetogram. The spectral estimate, in which the extrapolated field is assumed to be random and homogeneous with Gaussian statistics, is found here to be relatively accurate for quiet Sun magnetograms from SOHO’s MDI. The majority of null points occur at low altitudes, and their distribution is dictated by high wavenumbers in the Fourier spectrum. This portion of the spectrum is affected by Poisson noise, and as many as five-sixths of null points identified from a direct extrapolation can be attributed to noise. The null distribution above 1500 km is found to depend on wavelengths that are reliably measured by MDI in either its low-resolution or high-resolution mode. After correcting the spectrum to remove white noise and compensate for the modulation transfer function we find that a potential field extrapolation contains, on average, one magnetic null point, with altitude greater than 1.5 Mm, above every 322 Mm² patch of quiet Sun. Analysis of 562 quiet Sun magnetograms spanning the two latest solar minima shows that the null point density is relatively constant with roughly 10% day-to-day variation. At heights above 1.5 Mm, the null point density decreases approximately as the inverse cube of height. The photospheric field in the quiet Sun is well approximated as that from discrete elements with mean flux 〈|φ|〉=1.0×10¹⁹ Mx distributed randomly with density n=0.007 Mm⁻².     

Comparison of Five Numerical Codes for Automated Tracing of Coronal Loops

The three-dimensional (3D) modeling of coronal loops and filaments requires algorithms that automatically trace curvilinear features in solar EUV or soft X-ray images. We compare five existing algorithms that have been developed and customized to trace curvilinear features in solar images: i) the oriented-connectivity method (OCM), which is an extension of the Strous pixel-labeling algorithm (developed by Lee, Newman, and Gary); ii) the dynamic aperture-based loop-segmentation method (developed by Lee, Newman, and Gary); iii) unbiased detection of curvilinear structures (developed by Steger, Raghupathy, and Smith); iv) the oriented-direction method (developed by Aschwanden); and v) ridge detection by automated scaling (developed by Inhester). We test the five existing numerical codes with a TRACE image that shows a bipolar active region and contains over 100 discernable loops. We evaluate the performance of the five codes by comparing the cumulative distribution of loop lengths, the median and maximum loop length, the completeness or detection efficiency, the accuracy, and flux sensitivity. These algorithms are useful for the reconstruction of the 3D geometry of coronal loops from stereoscopic observations with the STEREO spacecraft, or for quantitative comparisons of observed EUV loop geometries with (nonlinear force-free) magnetic field extrapolation models.     

Observations of Low-Latitude Coronal Plumes

Using Fe ix/x 17.1 nm observations from the Extreme-Ultraviolet Imaging Telescope (EIT) on the Solar and Heliospheric Observatory (SOHO), we have identified many coronal plumes inside low-latitude coronal holes as they transited the solar limb during the late declining phase of cycle 23. These diffuse, linear features appear to be completely analogous to the familiar polar plumes. By tracking them as they rotate from the limb onto the disk (or vice versa), we confirm that EUV plumes seen against the disk appear as faint, diffuse blobs of emission surrounding a brighter core. When the EIT images are compared with near-simultaneous magnetograms from the SOHO Michelson Doppler Imager (MDI), the low-latitude, on-disk plumes are found to overlie regions of mixed polarity, where small bipoles are in contact with unipolar flux concentrations inside the coronal hole. The birth and decay of the plumes are shown to be closely related to the emergence of ephemeral regions, their dispersal in the supergranular flow field, and the cancellation of the minority-polarity flux against the dominant-polarity network elements. In addition to the faint polar and nonpolar plumes associated with ephemeral regions, we note the existence of two topologically similar coronal structures: the giant plume-like features that occur above active regions inside coronal holes, and the even larger scale “pseudostreamers” that separate coronal holes of the same polarity. In all three cases, the basic structure consists of open field lines of a given polarity overlying a photospheric region of the opposite polarity; ongoing interchange reconnection at the X-point separating the open field domains from the underlying double-arcade system appears to result in the steady evaporation of material from the closed into the open region.     

新冠肺炎疫情扩散与人口流动的空间关系及对中国城市公共卫生分类治理启示

人口流动影响新冠肺炎疫情传播与风险扩散。基于百度迁徙大数据和各省市区卫生健康委员会数据，结合地理信息技术，研究了2020年1月1日至3月5日136个城市新冠肺炎疫情扩散与武汉市人口流出的空间关系及其对我国城市公共卫生治理启示。研究表明：1）中国新冠肺炎疫情扩散过程具有阶段性特征，经历了疫情发生与隐性扩散、快速扩散与暴发、扩散遏制和扩散衰减4个阶段。2）研究时间段武汉市人口主要流向湖北省境内以及周边省市和北京、上海、广州、深圳等一线城市，具有地理邻近性和倾向区域中心城市的人口流入特征。受地理距离、时间成本、社会经济联系、境外输入等因素的影响，新冠肺炎疫情空间分布的不平衡性明显，长江中游城市群、京津冀城市群、长三角城市群、粤港澳大湾区和成渝城市群成为新冠肺炎疫情集中分布的重点区域，一些重点出入境口岸城市的新冠疫情扩散风险较大。3）新冠肺炎疫情扩散与人口流出之间具有较强的正向等级相关性。两者之间的空间关系可以分为8种调控类型，近90%的城市具有人口流入多、确诊病例数高或人口流入少、确诊病例数低的特征。其中，人口流入多、确诊病例数高的城市主要集中分布在湖北省境内以及中国重点城市群的中心城市，其防控压力来自人口流入多、确诊病例数高带来的疫情扩散风险；而人口流入少、确诊病例数低的城市分布较为分散，其防控的难点在于提高防控对策的精准性。我国疫情防控取得了显著成效，但随着时间的演进全球疫情形势反弹的不确定性仍然存在，外防输入、内防反弹的压力依然很重，现阶段乃至未来一段时期，疫情防控将伴随中国经济社会发展成为新常态。建议针对8种调控类型，从人口流动、交通和资源等引导与管控方面分类提出精细化的疫情防控策略，提升城市公共卫生治理能力。     

稻城冰帽库照日冰碛垄的宇生核素10Be年代测定

沙鲁里山南端的稻城冰帽在第四纪曾发生过大规模的冰川作用,然而,关于冰川作用的年代数据相对较少。选择稻城冰帽南缘库照日附近的3套冰碛垄进行宇生核素10Be的暴露年代测定,研究结果表明:(1)库照日最老冰碛垄(E)的暴露年代为(714±8.0)～(529±8.0)kaBP,D垄的暴露年代为(121±2.9)～(114±2.9)kaBP,DC垄的暴露年代为(47±1.3)～(43±1.3)kaBP。由于宇生核素暴露年代所反映的是冰川消退时的年代,因此,我们结合深海氧同位素的气候特征以及其他测年数据推断,E垄、D垄、DC垄可能是MIS16、MIS6以及MIS3b阶段冰川作用的产物;(2)冰川退却是一种逐步退却的模式。样品E的两次测验结果比较吻合,表明宇生核素10 Be暴露测年方法对冰碛垄形成年代的测定是可靠的,该方法在地貌学和第四纪年代学研究中有着巨大的潜力。     

Hydrogeochemistry and groundwater salinization in an ephemeral coastal flood plain: Cap Bon,Tunisia

Abstract

The Wadi Al Ayn plain is a coastal system on the eastern coast of Cap Bon in northeastern Tunisia. The area is known for its intensive agriculture, which is based mainly on groundwater exploitation. The aim of this study is to identify the sources of groundwater salinization in the Wadi Al Ayn aquifer system and deduce the processes that drive the mineralization. Surface water and groundwater samples were taken and analysed for major ions and stable isotopes. The geochemical data were used to characterize and classify the water samples based on a variety of ion plots and diagrams. Stable isotopes are useful tools to help us understand recharge processes and to differentiate between salinity origins. The oilfield brines infiltrated from the sandy bed of Wadi Al Ayn comprise the main source of groundwater salinization in the central part of the plain, while seawater intrusion is mainly responsible for the increased salinity in the groundwater of the coastal part of the plain (at Daroufa).

Citation Chekirbane, A., Tsujimura, M., Kawachi, A., Isoda, H., Tarhouni, J., and Benalaya, A., 2013. Hydrogeochemistry and groundwater salinization in an ephemeral coastal flood plain: Cap Bon, Tunisia. Hydrological Sciences Journal, 58 (5), 1097–1110.     

Small-scale ductile shear zones as transposed rheologically weak domains: A numerical modeling investigation and practical application

We test the hypothesis that small ductile shear zones are developed from initial rheologically weak domains. We regard weak domains as ellipsoidal inhomogeneities and apply Eshelby's formalism extended for power-law viscous materials to investigate the kinematics and finite strain evolution of the partitioned flow fields in weak domains. We show that, under an imposed bulk flow field, weak domains, regardless of their initial shapes and orientations, deform into zone-like features at relatively small bulk strains and the deformation paths inside weak domains have all characteristics expected in ductile shear zones. We apply our model to the Cap de Creus area Spain, where abundant small ductile shear zones exist. To relate the fabrics inside these shear zones with the regional deformation, we take a multi-scale approach. We assume that the area contained many weak domains which were randomly shaped and oriented initially and were deformed into shear zones eventually in response to the regional flow. We constrain the regional flow field by the fabric patterns across the area and compute numerically the partitioned flow fields in individual weak domains. The latter are related to fabrics inside shear zones. Our model reproduces first-order features of field-observed fabrics. Although the deformation path of each shear zone was close to simple shearing, the deformation of the whole belt was not. Our approach also resolves the strain compatibility problem for a finite-sized shear zone embedded in a far less deformed country rock.