高效优化非分裂PML边界二阶标量波方程数值模拟方法

卷积完全匹配层（convolution perfectly matched layer,CPML）吸收边界是一种高效处理波动方程数值模拟中人工边界反射波的方法。本文基于传统的一阶系统CPML吸收边界条件推广并推导了新的二阶系统CPML边界条件（NCPML）。与常规二阶系统CPML边界条件不同,新边界条件推导的核心思想是在复数-频率域中忽略部分衰减因子空间变化特性,避免其在时间域产生复杂卷积算子,然后反变换至时间域得到基于CPML吸收条件的二阶标量波方程,并应用于二阶标量波方程数值模拟。均匀介质模型测试验证了NCPML吸收条件在内存使用上相对于常规二阶系统CPML与SPML（split PML）吸收条件更少。在对人工边界反射的吸收效果上,NCPML稍逊色于常规二阶系统CPML,但二者均相对于SPML优势明显。最后通过层状模型和Marmousi模型测试验证了NCPML的稳定性及其在效率上的优势。     

基于NAD算法的声波方程时间四阶差分解法

波动方程数值模拟是研究地震波传播机理的重要工具,有限差分求解波动方程是当前地震波数值模拟的主要方法之一。当地下介质中的地震波速度较低或地震波高频成分丰富时,常规有限差分技术常常产生严重的数值频散误差,这种误差会降低数值模拟的精度,影响对地震波传播机理的分析。为压制地震波数值模拟时产生的数值频散误差,提高波场模拟精度,提出了基于NAD算子的时间四阶精度波动方程差分格式。根据对应的差分格式,分析了该差分格式的数值频散关系。与常规四阶精度差分算法的频散曲线相比,基于NAD时间四阶精度差分方法不但能够实现时间频散的有效压制,同时其基于更多网格点的位移分量和位移梯度分量空间微分求解方法还能够实现空间频散的有效压制。另外在相同模型条件下,基于NAD算法的声波方程时间四阶差分解法可采用大网格对模拟空间进行差分离散,减少网格数,提高计算效率。     

Numerical study for vegetation effects on coastal wave propagation by using nonlinear Boussinesq model

The vegetation has important impacts on coastal wave propagation. In the paper, the sensitivities of coastal wave attenuation due to vegetation to incident wave height, wave period and water depth, as well as vegetation configurations are numerically studied by using the fully nonlinear Boussinesq model. The model is based on the implementation of drag resistances due to vegetation in the fully nonlinear Boussinesq equation where the drag resistance is provided by the Morison’s formulation for rigid structure induced drag stresses. The model is firstly validated by comparing with the experimental results for wave propagation in vegetation zones. Subsequently, the model is used to simulate waves with different height, period propagating on vegetation zones with different water depth and vegetation configurations. The sensitivities of wave attenuation to incident wave height, wave period, water depth, as well as vegetation configurations are investigated based on the numerical results. The numerical results indicate that wave height attenuation due to vegetation is sensitive to incident wave height, wave period, water depth, as well as vegetation configurations, and attenuation ratio of wave height is increased monotonically with increases of incident wave height and decreases of water depth, while it is complex for wave period. Moreover, more vegetation segments can strengthen the interaction of vegetation and wave in a certain range.     

Australian ‘gayborhoods’ and ‘lesborhoods’: a new method for estimating the number and prevalence of adult gay men and lesbian women living in each Australian postcode

ABSTRACT

Gay men and lesbian women often demonstrate unique settlement patterns, forming what have been referred to as ‘gayborhoods’. This study sought to provide the first postcode-level estimates of population size and prevalence of gay and lesbian people in Australia. Data on same-gender-partnered households from the Australian Census were combined with information from six different surveys conducted from 2011 to 2017. We estimated that in 2016 there were 132,203 gay men (1.5% of adult males; 95% CI: 1.4–1.6) and 79,931 lesbian women (0.9% of adult females; 95% CI: 0.8–1.0) in Australia. While many postcodes were sparsely populated by gay and lesbian people (40.1% had prevalences of <0.1%), 24.6% were moderately populated (prevalences in the 50-95th percentile) and 2.7% were highly populated (95th percentile). By jurisdiction, the Australian Capital Territory had the highest prevalences of gay men (2.1%; 95% CI: 2.0–2.2) and lesbian women (1.5%; 95% CI: 1.4–1.6). Although the majority of highly populated postcodes were found in major cities (83.7%), some were also found in regional and remote area (16.3%). This method can be applied in other countries to enhance populate estimates. The accompanying dataset can be used to guide service delivery, conduct geographically contextualised research and develop policies relevant to gay men and lesbian women in Australia.     

传统聚落景观基因的地理信息特征及其理解

传统聚落是民族传统文化遗产的重要组成部分,对社会经济和文化建设具有重要的价值。针对现有研究缺乏探讨传统聚落景观基因（简称景观基因）蕴藏的地理信息特征及理解方法,论文从以下方面开展了深入分析：① 地理信息包含了语义描述、几何形态、属性特征、维度、时空框架、尺度、要素相互关系（空间关联）、演化过程（存在状态）等属性;② 景观基因是一种特殊的文化因子,蕴含着丰富的哲理,是认识传统聚落特征的分析方法,也是文化符号的集合;③ 景观基因包含了空间定位、载体特征、历史与文化特征等丰富的地理信息。根据前述分析结果,论文从符号化、数据挖掘、谱系分析与空间格局制图等途径系统地探讨了景观基因地理信息特征的理解方法。论文认为：在地理时空大数据与地理服务日益深化发展的背景下,结合GIS原理探讨景观基因的地理信息特征对于延伸地理信息科学的内涵,深化传统聚落的地学认识,促进人文GIS的发展具有积极的意义,今后应该继续加强相关方法探索。     

Advances in pan‐European flood hazard mapping

Flood hazard maps at trans‐national scale have potential for a large number of applications ranging from climate change studies, reinsurance products, aid to emergency operations for major flood crisis, among others. However, at continental scales, only few products are available, due to the difficulty of retrieving large consistent data sets. Moreover, these are produced at relatively coarse grid resolution, which limits their applications to qualitative assessments. At finer resolution, maps are often limited to country boundaries, due to limited data sharing at trans‐national level. The creation of a European flood hazard map would currently imply a collection of scattered regional maps, often lacking mutual consistency due to the variety of adopted approaches and quality of the underlying input data. In this work, we derive a pan‐European flood hazard map at 100 m resolution. The proposed approach is based on expanding a literature cascade model through a physically based approach. A combination of distributed hydrological and hydraulic models was set up for the European domain. Then, an observed meteorological data set is used to derive a long‐term streamflow simulation and subsequently coherent design flood hydrographs for a return period of 100 years along the pan‐European river network. Flood hydrographs are used to simulate areas at risk of flooding and output maps are merged into a pan‐European flood hazard map. The quality of this map is evaluated for selected areas in Germany and United Kingdom against national/regional hazard maps. Despite inherent limitations and model resolution issues, simulated maps are in good agreement with reference maps (hit rate between 59% and 78%, critical success index between 43% and 65%), suggesting strong potential for a number of applications at the European scale. Copyright © 2013 John Wiley & Sons, Ltd.     

An ensemble approach for the analysis of extreme rainfall under climate change in Naples (Italy)

In the present paper, an ensemble approach is proposed to estimate possible modifications caused by climate changes in the extreme precipitation regime, with the rain gauge Napoli Servizio Idrografico (Naples, Italy) chosen as test case. The proposed research, focused on the analysis of extremes on the basis of climate model simulations and rainfall observations, is structured in several consecutive steps. In the first step, all the dynamically downscaled EURO‐CORDEX simulations at about 12 km horizontal resolution are collected for the current period 1971–2000 and the future period 2071–2100, for the RCP4.5 and the RCP8.5 concentration scenarios. In the second step, the significance of climate change effects on extreme precipitation is statistically tested by comparing current and future simulated data and bias‐correction is performed by means of a novel approach based on a combination of simple delta change and quantile delta mapping, in compliance with the storm index method. In the third step, two different ensemble models are proposed, accounting for the variabilities given by the use of different climate models and for their hindcast performances. Finally, the ensemble models are used to build novel intensity–duration–frequency curves, and their effects on the early warning system thresholds for the area of interest are evaluated.     

Regional‐scale mapping of groundwater discharge zones using thermal satellite imagery

Mapping groundwater discharge zones at broad spatial scales remains a challenge, particularly in data sparse regions. We applied a regional scale mapping approach based on thermal remote sensing to map discharge zones in a complex watershed with a broad diversity of geological materials, land cover and topographic variation situated within the Prairie Parkland of northern Alberta, Canada. We acquired winter thermal imagery from the USGS Landsat archive to demonstrate the utility of this data source for applications that can complement both scientific and management programs. We showed that the thermally determined potential discharge areas were corroborated with hydrological (spring locations) and chemical (conservative tracers of groundwater) data. This study demonstrates how thermal remote sensing can form part of a comprehensive mapping framework to investigate groundwater resources over broad spatial scales. Copyright © 2013 John Wiley & Sons, Ltd.     

Influence of dissolution/reprecipitation reactions on metamorphic greenschist to amphibolite facies mica 40Ar/39Ar ages in the Longmen Shan (eastern Tibet)

Linking ages to metamorphic stages in rocks that have experienced low‐ to medium‐grade metamorphism can be particularly tricky due to the rarity of index minerals and the preservation of mineral or compositional relicts. The timing of metamorphism and the Mesozoic exhumation of the metasedimentary units and crystalline basement that form the internal part of the Longmen Shan (eastern Tibet, Sichuan, China), are, for these reasons, still largely unconstrained, but crucial for understanding the regional tectonic evolution of eastern Tibet. In situ core‐rim ⁴⁰Ar/³⁹Ar biotite and U–Th/Pb allanite data show that amphibolite facies conditions (~10–11 kbar, 530°C to 6–7 kbar, 580°C) were reached at 210–180 Ma and that biotite records crystallization, rather than cooling, ages. These conditions are mainly recorded in the metasedimentary cover. The ⁴⁰Ar/³⁹Ar ages obtained from matrix muscovite that partially re‐equilibrated during the post peak‐P metamorphic history comprise a mixture of ages between that of early prograde muscovite relicts and the timing of late muscovite recrystallization at c. 140–120 Ma. This event marks a previously poorly documented greenschist facies metamorphic overprint. This latest stage is also recorded in the crystalline basement, and defines the timing of the greenschist overprint (7 ± 1 kbar, 370 ± 35°C). Numerical models of Ar diffusion show that the difference between ⁴⁰Ar/³⁹Ar biotite and muscovite ages cannot be explained by a slow and protracted cooling in an open system. The model and petrological results rather suggest that biotite and muscovite experienced different Ar retention and resetting histories. The Ar record in mica of the studied low‐ to medium‐grade rocks seems to be mainly controlled by dissolution–reprecipitation processes rather than by diffusive loss, and by different microstructural positions in the sample. Together, our data show that the metasedimentary cover was thickened and cooled independently from the basement prior to c. 140 Ma (with a relatively fast cooling at 4.5 ± 0.5°C/Ma between 185 and 140 Ma). Since the Lower Cretaceous, the metasedimentary cover and the crystalline basement experienced a coherent history during which both were partially exhumed. The Mesozoic history of the Eastern border of the Tibetan plateau is therefore complex and polyphase, and the basement was actively involved at least since the Early Cretaceous, changing our perspective on the contribution of the Cenozoic geology.     

Urban forest loss using a GIS-based approach and instruments for integrated urban planning:A case study of Jo(a)o Pessoa,Brazil

Urban forests play an important role in the thermal comfort and overall life of local populations in large-and medium-sized cities.This study analyzes urban forest loss and maps land use and land cover (LULC) changes between 1991 and 2018 by evaluating the use of urban planning instruments for the mitigation of urban forest loss in Jo～o Pessoa,Brazil.For this purpose,satellite-derived LULC images from 1991,2006,2010 and 2018 and data on urban forest loss areas obtained using the Google Earth Engine were used.In addi-tion,this paper also discusses the instruments used for integrated urban planning,which are(a) the legal sector,responsibility and nature;(b) the urban expansion process;and (c) the elements of urban infrastructure.The results show a clear shift in land use in the study area.The major changes in LULC classes occurred in urban areas and herbaceous vegetation,while the greatest loss was in arboreal/shrub vegetation.Thus,an increase in the pressure to occupy zones intended for environmental preservation could be estimated.Our results showed similar accuracies with other studies and more spatial details.The characteristics of the patterns,traces,and hotspots of urban expansion and forest cover loss were explored.We highlighted the potential use of this proposed framework to be applied and validated in other parts of the world to help better understand and quantify various aspects of ur-ban-related problems such as urban forest loss mapping using instruments for integrated urban planning and low-cost approaches.