地球科学大数据的管理与共享：以英国地质调查局为例

大数据驱动的研究范式正在引起地学领域的革命,而海量大数据的有效管理和共享是数据高效利用的前提。英国地质调查局作为最早成立的国家地质调查局,拥有海量的地学数据资源,通过近年来对数字化工作的全面推进,在数据的开放共享方面走在了世界各国的最前沿。文章对英国地质调查局的数据资源管理和数据共享方式进行了分析调研,详细介绍了他们的一站式管理平台——开放地学的主体组成,以及他们与同行合作建设的数据库。开放地学全面汇总了地调局内的数据资源,并通过一系列数据共享将所有数据集有机链接,通过数据和模型的巧妙结合,在满足用户数据需求的同时,对数据的应用进行了全方位的拓展,是地球系统科学研究框架下地球科学数字化工作的良好典范。     

Integration of zircon and apatite U–Pb geochronology and geochemical mapping of the Wude basalts(Emeishan large igneous province):A tool for a better understanding of the tectonothermal and geodynamic evolution of the Emeishan LIP

Radiogenic isotopic dating and Lu–Hf isotopic composition using laser ablation-inductively coupled plasma-mass spectrometry(LA-ICP-MS)of the Wude basalt in Yunnan province from the Emeishan large igneous province(ELIP)yielded timing of formation and post-eruption tectonothermal event.Holistic lithogeochemistry and elements mapping of basaltic rocks were further reevaluated to provide insights into crustal contamination and formation of the ELIP.A zircon U–Pb age of 251.3±2.0 Ma of the Wude basalt recorded the youngest volcanic eruption event and was consistent with the age span of 251-263 Ma for the emplacement of the ELIP.Such zircons hadε_Hf(t)values ranging from7.3 to+2.2,identical to those of magmatic zircons from the intrusive rocks of the ELIP,suggesting that crust-mantle interaction occurred during magmatic emplacement,or crust-mantle mixing existed in the deep source region prior to deep melting.The apatite U–Pb age at 53.6±3.4 Ma recorded an early Eocene magmatic superimposition of a regional tectonothermal event,corresponding to the Indian–Eurasian plate collision.Negative Nb,Ta,Ti and P anomalies of the Emeishan basalt may reflect crustal contamination.The uneven Nb/La and Th/Ta values distribution throughout the ELIP supported a mantle plume model origin.Therefore,the ELIP was formed as a result of a mantle plume which was later superimposed by a regional tectonothermal event attributed to the Indian–Eurasian plate collision during early Eocene.     

Geochronological and geochemical evidence for a Late Ordovician to Silurian arc-back-arc system in the northern Great Xing'an Range,NE China

The early Paleozoic tectonic evolution of the Xing'an-Mongolian Orogenic Belt is dominated by two oceanic basins on the northwestern and southeastern sides of the Xing'an Block,i.e.,the Xinlin-Xiguitu Ocean and the Nenjiang Ocean.However,the early development of the Nenjiang Ocean remains unclear.Here,we present zircon U-Pb geochronology and whole-rock elemental and Sr-Nd isotopic data on the gabbros in the Xinglong area together with andesitic tuffs and basalts in the Duobaoshan area.LA-ICP-MS zircon U-Pb dating of gabbros and andesitic tuffs yielded crystallization ages of 443-436 Ma and 452-451 Ma,respectively.The Early Silurian Xinglong gabbros show calc-alkaline and E-MORB affinities but they are enriched in LILEs,and depleted in HFSEs,with relatively low U/Th ratios of 0.18-0.36 andεNd(t)values of-1.6 to+0.5.These geochemical features suggest that the gabbros might originate from a mantle wedge modified by pelagic sediment-derived melts,consistent with a back-arc basin setting.By contrast,the andesitic tuffs are characterized by high MgO(>5 wt.%),Cr(138-200 ppm),and Ni(65-110 ppm)contents,and can be termed as high-Mg andesites.Their low Sr/Y ratios of 15.98-17.15 and U/Th values of 0.24-0.25 and moderate(La/Sm)_n values of 3.07-3.26 are similar to those from the Setouchi Volcanic Belt(SW Japan),and are thought to be derived from partial melting of subducted sediments,and subsequent melt-mantle interaction.The Duobaoshan basalts have high Nb(8.44-10.30 ppm)and TiO2 contents(1.17-1.60 wt.%),typical of Nb-enriched basalts.They are slightly younger than regional adakitic rocks and have positiveεNd(t)values of+5.2 to+5.7 and are interpreted to be generated by partial melting of a depleted mantle source metasomatized by earlier adakitic melts.Synthesized with coeval arc-related igneous rocks from the southeastern Xing'an Block,we propose that the Duobaoshan high-Mg andesitic tuffs and Nbenriched basalts are parts of the Late Ordovician and Silurian Sonid Zuoqi-Duobaoshan arc belt,and they were formed by the northwestern subduction of the Nenjiang Ocean.Such a subduction beneath the integrated Xing'an-Erguna Block also gave rise to the East Ujimqin-Xinglong igneous belt in a continental back-arc basin setting.Our new data support an early Paleozoic arc-back-arc model in the northern Great Xing'an Range.     

Hydrodynamic performance of a T-shaped floating breakwater

The comprehensive utilization of floating breakwaters, specially acting as a supporting structure for offshore marine renewable energy explorations, has received more and more attention recently. Based on linear water-wave theory, the hydrodynamic performance of a T-shaped floating breakwater is semi-analytically investigated through the matched eigenfunction expansion method (MEEM). Auxiliary functions, to speed up the convergence and improve the accuracy in the numerical computations, are introduced to represent the singular behavior of fluid field near the lower salient corners of the structure. The effects of the height and installation position of the vertical screen on the reflection and transmission coefficients, dynamic response and wave forces are examined. It is found that the presence of the screen shifts the resonance frequency of RAO for both surge and pitch modes to the low-frequency area, while has no effect on heave mode. The identical added masses, damping and transmission coefficients can be obtained in the cases where the screen holds the same distance away from the longitudinal central axis of the upper box-type structure. Moreover, a relatively small pitch response can be achieved in a wide wave–frequency range, when the breakwater is Γ-shaped.     

An efficient and robust approach to predict ship self-propulsion coefficients

Achieving a reliable and accurate numerical prediction of the self-propulsion performance of a ship is still an open problem that poses some relevant issues. Several CFD methods, ranging from boundary element methods (BEM) to higher-fidelity viscous Reynolds averaged Navier–Stokes (RANS) based solvers, can be used to accurately analyze the separate problems, i.e. the open water propeller and the hull calm water resistance. However, when the fully-coupled self-propulsion problem is considered, i.e. the hull advancing at uniform speed propelled by its own propulsion system, several complexities rise up. Typical flow simplifications adopted to speed-up the simulations of the single analysis (hull and propeller separately) lose their validity requiring a more complex solver to tackle the fully-coupled problem. The complexity rises up further when considering a maneuver condition. This aspect increases the computational burden and, consequently, the required time which becomes prohibitive in a preliminary ship design stage.The majority of the simplified methods proposed in literature to include propeller effects, without directly solve the propeller flow, in a high-fidelity viscous solver are not able to provide all the commonly required self-propulsion coefficients. In this work, a new method to enrich the results from a body force based approach is proposed and investigated, with the aim to reduce as much as possible the computational burden without losing any useful result. This procedure is tested for validation on the KCS hull form in self-propulsion and maneuver conditions.     

CFD modelling of a small–scale fixed multi–chamber OWC device

Wave Energy Converters (WECs) have excellent potential as a source of renewable energy that is yet to be commercially realised. Recent attention has focused on the installation of Oscillating Water Column (OWC) devices as a part of harbor walls to provide advantages of cost–sharing structures and proximity of power generation facilities to existing infrastructure. In this paper, an incompressible three–dimensional CFD model is constructed to simulate a fixed Multi–Chamber OWC (MC–OWC) device. The CFD model is validated; the simulation results are found to be in good agreement with experimental results obtained from a scale physical model tested in a wave tank. The validated CFD model is then used for a benchmark study of 96 numerical tests. These investigate the effects of the PTO damping caused by the power take–off (PTO) system on device performance. The performance is assessed for a range of regular wave heights and periods. The results demonstrate that a PTO system with an intermediate damping can be used for all chambers in the MC–OWC device for most wave period ranges, except for the long wave periods. These require a higher PTO damping. An increased incident wave height reduces the device capture width ratio, but there is a noticeable improvement for long wave periods.     

Cartographic Design and the Space–Time Cube

This paper reports on the results of an empirical evaluation that aimed to define the effectiveness and efficiency of different visual variables in depicting the Space–Time Cube’s (STC) content. Existing STC applications demonstrate that the most used visual variables are size and colour hue. Less is known, however, about their usability metrics. The research sets design criteria for STC contents, such as space–time paths, based on the cartographic design theory. The visual variables colour hue, colour value, colour saturation, size and orientation have been applied in two different use case studies. Besides, to support the three-dimensional visual environment, depth cues such as shading and transparency were considered too. User tests have been executed based on real-world problems with particular attention for the visualization strategy and data complexity. The outcomes revealed the most efficient and effective visual variables to represent data of various complexities in the STC.     

Cartographic Symbol Design Considerations for the Space–Time Cube

The cartographic representation of geographic phenomena in the space–time cube comes with special challenges and opportunities when compared with two-dimensional maps. While the added dimension allows the display of attributes that vary with time, it is difficult to display rapidly varying temporal data given the limited display height. In this study, we adapt 2D cyclic point symbols to construct 3D surfaces designed along a helical path for the space–time cube. We demonstrate how these complex?3D helical surfaces can display detailed data, including data reported daily over 100 years and data reported in four-hour intervals over a year. To create the point symbols, each value is plotted along the curve of a helix, with each turn of the helix representing one year or week, respectively. The model is modified by varying the radii from the time axis to all points using the attribute value, in these cases maximum daily temperature and four-hourly ridership, and then creating a triangulated surface from the resulting points. Using techniques common to terrain representation, we apply hue and saturation to the surface based on attribute values, and lightness based on relief shading. Multiple surfaces can be displayed in a space–time cube with a consistent time interval facing the viewer, and the surfaces or viewer perspective can be rotated to display synchronized variations. We see this method as one example of how cartographic design can refine or enhance operations in the space–time cube.     

Re-examining urban region and inferring regional function based on spatial–temporal interaction

Urban system is shaped by the interactions between different regions and regions planned by the government, then reshaped by human activities and residents’ needs. Understanding the changes of regional structure and dynamics of city function based on the residents’ movement demand are important to evaluate and adjust the planning and management of urban services and internal structures. This paper constructed a probabilistic factor model on the basis of probabilistic latent semantic analysis and tensor decomposition, for purpose of understanding the higher order interactive population mobility and its impact on urban structure changes. First, a four-dimensional tensor of time (T)?×?week (W)?×?origin (O)?×?destination (D) was constructed to identify the day-to-day activities in three time modes and weekly regularity of weekday/weekend pattern. Then we reclassified the urban regions based on the space clustering formed by the space factor matrix and core tensor. Finally, we further analysed the space–time interaction on different time scales to deduce the actual function and connection strength of each region. Our research shows that the application of individual-based spatial–temporal data in human mobility and space–time interaction study can help to analyse urban spatial structure and understand the actual regional function from a new perspective.     

Future Changes in Extreme High Temperature over China at 1.5℃-5℃ Global Warming Based on CMIP6 Simulations

Extreme high temperature(EHT)events are among the most impact-related consequences related to climate change,especially for China,a nation with a large population that is vulnerable to the climate warming.Based on the latest Coupled Model Intercomparison Project Phase 6(CMIP6),this study assesses future EHT changes across China at five specific global warming thresholds(1.5℃-5℃).The results indicate that global mean temperature will increase by 1.5℃/2℃ before 2030/2050 relative to pre-industrial levels(1861-1900)under three future scenarios(SSP1-2.6,SSP2-4.5,and SSP5-8.5),and warming will occur faster under SSP5-8.5 compared to SSP1-2.6 and SSP2-4.5.Under SSP5-8.5,global warming will eventually exceed 5℃ by 2100,while under SSP1-2.6,it will stabilize around 2℃ after 2050.In China,most of the areas where warming exceeds global average levels will be located in Tibet and northern China(Northwest China,North China and Northeast China),covering 50%-70%of the country.Furthermore,about 0.19-0.44 billion people(accounting for 16%-41%of the national population)will experience warming above the global average.Compared to present-day(1995-2014),the warmest day(TXx)will increase most notably in northern China,while the number of warm days(TX90p)and warm spell duration indicator(WSDI)will increase most profoundly in southern China.For example,relative to the present-day,TXx will increase by 1℃-5℃ in northern China,and TX90p(WSDI)will increase by 25-150(10-80)days in southern China at 1.5℃-5℃ global warming.Compared to 2℃-5℃,limiting global warming to 1.5℃ will help avoid about 36%-87%of the EHT increases in China.