ALBERT

Description: Abstract The Miocene epoch, spanning 23.03–5.33 Ma, was a dynamic climate of sustained, polar amplified warmth. Miocene atmospheric CO2 concentrations are typically reconstructed between 300 and 600 ppm and were potentially higher during the Miocene Climatic Optimum (16.75–14.5 Ma). With surface temperature reconstructions pointing to substantial midlatitude and polar warmth, it is unclear what processes maintained the much weaker-than-modern equator-to-pole temperature difference. Here, we synthesize several Miocene climate modeling efforts together with available terrestrial and ocean surface temperature reconstructions. We evaluate the range of model-data agreement, highlight robust mechanisms operating across Miocene modeling efforts and regions where differences across experiments result in a large spread in warming responses. Prescribed CO2 is the primary factor controlling global warming across the ensemble. On average, elements other than CO2, such as Miocene paleogeography and ice sheets, raise global mean temperature by ∼2°C, with the spread in warming under a given CO2 concentration (due to a combination of the spread in imposed boundary conditions and climate feedback strengths) equivalent to ∼1.2 times a CO2 doubling. This study uses an ensemble of opportunity: models, boundary conditions, and reference data sets represent the state-of-art for the Miocene, but are inhomogeneous and not ideal for a formal intermodel comparison effort. Acknowledging this caveat, this study is nevertheless the first Miocene multi-model, multi-proxy comparison attempted so far. This study serves to take stock of the current progress toward simulating Miocene warmth while isolating remaining challenges that may be well served by community-led efforts to coordinate modeling and data activities within a common analytical framework.

Description: With the increasing pressure from population growth and economic development, northern China (NC) faces a grand challenge of water scarcity, which can be further exacerbated by climatic and societal changes. The South‐to‐North Water Diversion (SNWD) project is designed to mitigate the water scarcity in NC. However, few studies have quantified the impact of the SNWD on water scarcity within the context of climatic and societal changes and its potential effects on economic and agricultural food in the region. We used water supply stress index (WaSSI) to quantify water scarcity within the context of environmental change in NC, and developed a method to estimate the economic and agricultural impacts of the SNWD. Focuses were put on alleviating the water supply shortage and economic and agricultural benefits for the water‐receiving NC. We find that societal changes, especially economic growth, are the major contributors to water scarcity in NC during 2009–2099. To completely mitigate the water scarcity of NC, at least an additional water supply of 13 billion m3/year (comparable to the annual diversion water by SNWD Central Route) will be necessary. Although SNWD alone cannot provide the full solution to northern China's water shortage in next few decades, it can significantly alleviate the water supply stress in NC (particularly Beijing), considerably increasing the agricultural production (more than 115 Teracalories/year) and bringing economic benefits (more than 51 billion RMB/year) through supplying industrial and domestic water use. Additionally, the transfer project could have impacts on the ecological environment in the exporting regions.

Description: Eigen analysis has been a powerful tool to distinguish multiple processes into different simple principal modes in complex systems. For a non-equilibrium system, the principal modes corresponding to the non-equilibrium processes are usually evolving with time. Here, we apply the eigen analysis into the complex climate systems. In particular, based on the daily surface air temperature in the tropics (30° S–30° N, 0° E–360° E) between 1979-01-01 and 2016-12-31, we uncover that the strength of two dominated intra-annual principal modes represented by the eigenvalues significantly changes with the El Ni$\tilde {\mathrm{n}}$o/southern oscillation from year to year. Specifically, according to the 'regional correlation' introduced for the first intra-annual principal mode, we find that a sharp positive peak of the correlation between the El Ni$\tilde {\mathrm{n}}$o region and the northern (southern) hemisphere usually signals the beginning (end) of the El Ni$\tilde {\mathrm{n}}$o. We discuss the underlying physical mechanism and suppose that the evolution of the first intra-annual principal mode is related to the meridional circulations; the evolution of the second intra-annual principal mode responds positively to the Walker circulation. Our framework presented here not only facilitates the understanding of climate systems but also can potentially be used to study the dynamical evolution of other natural or engineering complex systems.

Description: Here we present a method to estimate the total number of nodes of a network using locally observed response dynamics. The algorithm has the following advantages: (a) it is data-driven. Therefore it does not require any prior knowledge about the model; (b) it does not need to collect measurements from multiple stimulus; and (c) it is distributed as it uses local information only, without any prior information about the global network. Even if only a single node is measured, the exact network size can be correctly estimated using a single trajectory. The proposed algorithm has been applied to both linear and nonlinear networks in simulation, illustrating the applicability to real-world physical networks.

Description: Distributed adaptive learning algorithms have played a critical role in signal processing and parameter estimation over networks. Most existing algorithms are based on the mean-square error (MSE) criterion, and they can achieve desirable performance when the noise is modeled as Gaussian. However, the performance of MSE-based algorithms may degrade dramatically with the impulsive noise. Therefore, the aim of this paper is to present a diffusion algorithm, named generalized correntropy-based logarithmic difference (d-GCLD) algorithm, for distributed estimation that incorporates robustness to wireless sensor networks (WSNs). By combining the logarithm operation and the correntropy criterion as the loss function, the proposed algorithm is robust to impulsive noise and achieves satisfactory performance in various situations. In addition, the stability problem is studied theoretically. Experimental results are given to demonstrate the validity of the new algorithm in different scenarios.

Description: Community detection is of great significance for understanding network functions and behaviors. With the successful application of deep learning in network-based analyses, recent studies have turned to utilizing graph convolutional networks (GCNs) to this problem due to their capability in capturing network attributes. Nevertheless, most existing GCN-based community detection approaches are semi-supervised and local structure-aware, even though community detection is an unsupervised learning problem essentially. In this paper, we develop a novel GCN method for unsupervised community detection under the framework of mutual information (MI) maximization, called UCDMI. Specifically, a novel MI maximization mechanism is developed to capture more fine-grained information of the global network structure in an unsupervised manner. Moreover, a new aggregation function is proposed for GCN to distinguish the importance between different neighboring nodes, which enables our method to identify more high-quality node representations and improve the community detection performance. Our extensive experiments demonstrate the effectiveness of our proposed UCDMI compared with several state-of-the-art community detection methods.

Description: Apatite is a ubiquitous accessory mineral in crustal rocks. The Sr-isotope record of apatite has a wide range of applications in earth science studies. However, apatite has been documented to be easily altered by certain fluids. Currently, the impact of metasomatic alteration on Sr-isotopic abundances in apatite is not well known. In order to better understand this issue, well-characterized fluorapatite grains have been metasomatized experimentally at a temperature of 800 or 600 °C and a pressure of 200 MPa. Fluids used included a CO2-H2O mixture and NaF-, CaCl2-, and HCl-bearing solutions, all of which were doped with a standard solution with a known 87Sr/86Sr ratio. In the fluorapatite + CO2-H2O experiments, the fluorapatite grains were not altered by the fluids; thus, their Sr isotope compositions were generally kept unchanged. However, the other fluids induced partial to complete alteration of fluorapatite. In experiments involving the NaF- or HCl-bearing solutions, the Sr content remained constant or was increased in altered areas of the fluorapatite, and the Sr isotopes underwent changes with partial isotopic signature of the reacting solutions. In experiments involving CaCl2, the Sr content was decreased in the altered fluorapatite because high Ca activity in the solution caused Ca to replace Sr on the Ca site. Notably, the Sr isotopic ratios are still changed, although to a relatively small extent. The efficient Sr isotopic exchange between fluid and apatite is attributed to the rapid transport rate of Sr between the reaction-interface fluid and bulk fluid surrounding the apatite. This experimental study demonstrates that the response of apatite Sr isotopes to metasomatic alteration is mainly controlled by the chemistry of fluids. Overall, Sr isotopes become susceptible to hydrothermal alteration once the apatite is chemically reactive with the fluids. Therefore, it is important to evaluate the fluid-rock history, especially the conjectured fluid composition before using Sr isotopes from apatite as a geochemical tracer.

Description: Dense hydrous magnesium silicates (DHMSs) with large water content and wide stability fields are a potential H2O reservoir in the deep Earth. Al-bearing superhydrous phase B (shy-B) with a wider stability field than the Al-free counterpart can play an important role in understanding H2O transport in the Earth’s transition zone and topmost lower mantle. In this study, a nominally Al-free and two different Al-bearing shy-B with 0.47(2) and 1.35(4) Al atoms per formula unit (pfu), were synthesized using a rotating multi-anvil press. The single-crystal structures were investigated by X-ray diffraction (XRD) complemented by Raman spectroscopy, and Fourier-transform infrared spectroscopy (FTIR). Single-crystal XRD shows that the cell parameters decrease with increasing Al-content. By combining X-ray diffraction and spectroscopy results, we conclude that the Al-poor shy-B crystallizes in the Pnn2 space group with hydrogen in two different general positions. Based on the results of the single crystal X-ray diffraction refinements combined with FTIR spectroscopy, three substitutions mechanisms are proposed: 2 Al3+ = Mg2+ + Si4+; ☐Mg2+ = 2H+; Si4+ = Al3+ + H+. Thus, in addition to the two general H positions, hydrogen is incorporated into the hydrous mineral via point defects. The elastic stiffness coefficients were measured for the Al-shy-B with 1.35 pfu Al by Brillouin scattering (BS). Al-bearing shy-B shows lower C11, higher C22 and similar C33 when compared to Al-free shy-B. The elastic anisotropy of Al-bearing shy-B is also higher than that of the Al-free composition. Such different elastic properties are due to the effect of lattice contraction as a whole and the specific chemical substitution mechanism that affect bonds strength. Al-bearing shy-B with lower velocity, higher anisotropy and wider thermodynamic stability can help to understand the low velocity zone and high anisotropy region in the subducted slab located in Tonga.

Description: Water reinjection into the formation is an indispensable operation in many energy engineering practices. This operation involves a complex hydromechanical (HM) coupling process and sometimes even causes unpredictable disasters, such as induced seismicity. It is acknowledged that the relative magnitude and direction of the principal stresses significantly influence the HM behaviors of rocks during injection. However, due to the limitations of current testing techniques, it is still difficult to comprehensively conduct laboratory injection tests under various stress conditions, such as in triaxial extension stress states. To this end, a numerical study of HM changes in rocks during injection under different stress states is conducted. In this model, the saturated rock is first loaded to the target stress state under drainage conditions, and then the stress state is maintained and water is injected from the top to simulate the formation injection operation. Particular attention is given to the difference in HM changes under triaxial compression and extension stresses. This includes the differences in the pore pressure propagation, mean effective stress, volumetric strain, and stress-induced permeability. The numerical results demonstrate that the differential stress will significantly affect the HM behaviors of rocks, but the degree of influence is different under the two triaxial stress states. The HM changes caused by the triaxial compression stress states are generally greater than those of extension, but the differences decrease with increasing differential stress, indicating that the increase in the differential stress will weaken the impact of the stress state on the HM response. In addition, the shear failure p otential of fracture planes with various inclination angles is analyzed and summarized under different stress states. It is recommended that engineers could design suitable injection schemes according to different tectonic stress fields versus fault occurrence to reduce the risk of injection-induced seismicity.