Revisiting ENSO impacts on the Indian Ocean SST based on a combined linear regression method

The El Ni?o-Southern Oscillation (ENSO) has great impacts on the Indian Ocean sea surface temperature (SST). In fact, two major modes of the Indian Ocean SST namely the Indian Ocean Basin (IOB) and the Indian Ocean Dipole (IOD) modes, exerting strong influences on the Indian Ocean rim countries, are both influenced by the ENSO. Based on a combined linear regression method, this study quantifies the ENSO impacts on the IOB and the IOD during ENSO concurrent, developing, and decaying stages. After removing the ENSO impacts, the spring peak of the IOB disappears along with significant decrease in number of events, while the number of events is only slightly reduced and the autumn peak remains for the IOD. By isolating the ENSO impacts during each stage, this study reveals that the leading impacts of ENSO contribute to the IOD development, while the delayed impacts facilitate the IOD phase switch and prompt the IOB development. Besides, the decadal variations of ENSO impacts are various during each stage and over different regions. These imply that merely removing the concurrent ENSO impacts would not be sufficient to investigate intrinsic climate variability of the Indian Ocean, and the present method may be useful to study climate variabilities independent of ENSO.     

Trace-element characteristics of east–west mantle geochemical hemispheres

Recent statistical analyses on the isotopic compositions of oceanic, arc, and continental basalts have revealed that the Earth's mantle is broadly divided into eastern and western hemispheres. The present study aimed to characterize the isotopically defined east–west geochemical hemispheres using trace-element concentrations. Basalt data with Rb, Sr, Nd, Sm, Pb, Th, and U in addition to the isotopic ratios ⁸⁷Sr/⁸⁶Sr, ¹⁴³Nd/¹⁴⁴Nd, ²⁰⁶Pb/²⁰⁴Pb, ²⁰⁷Pb/²⁰⁴Pb, and ²⁰⁸Pb/²⁰⁴Pb were selected mostly from the GEOROC and PetDB databases. A total of 4787 samples were used to investigate the global geochemical variations. The results show that the wide trace-element variations are broadly explained by the melting of melt-metasomatized and fluid-metasomatized mantle sources. The larger amount of the fluid component derived from subducted plates in the eastern hemisphere than that in the western hemisphere is inferred from the basalts. These characteristics support the hypothesis that focused subduction towards the supercontinent created the mantle geochemical hemispheres.     

Seismic Wavefields in the Deep Seafloor Area from a Submarine Landslide Source

We use the finite difference method to simulate seismic wavefields at broadband land and seafloor stations for a given terrestrial landslide source, where the seafloor stations are located at water depths of 1,900–4,300 m. Our simulation results for the landslide source explain observations well at the seafloor stations for a frequency range of 0.05–0.1 Hz. Assuming the epicenter to be located in the vicinity of a large submarine slump, we also model wavefields at the stations for a submarine landslide source. We detect propagation of the Airy phase with an apparent velocity of 0.7 km/s in association with the seawater layer and an accretionary prism for the vertical component of waveforms at the seafloor stations. This later phase is not detected when the structural model does not consider seawater. For the model incorporating the seawater, the amplitude of the vertical component at seafloor stations can be up to four times that for the model that excludes seawater; we attribute this to the effects of the seawater layer on the wavefields. We also find that the amplification of the waveform depends not only on the presence of the seawater layer but also on the thickness of the accretionary prism, indicating low amplitudes at the land stations and at seafloor stations located near the trough but high amplitudes at other stations, particularly those located above the thick prism off the trough. Ignoring these characteristic structures in the oceanic area and simply calculating the wavefields using the same structural model used for land areas would result in erroneous estimates of the size of the submarine landslide and the mechanisms underlying its generation. Our results highlight the importance of adopting a structural model that incorporates the 3D accretionary prism and seawater layer into the simulation in order to precisely evaluate seismic wavefields in seafloor areas.     

Physical properties of the shallow sediments in late Pleistocene formations,Ursa Basin,Gulf of Mexico,and their implications for generation and preservation of shallow overpressures

Understanding the evolution of abnormally high fluid pressures within sedimentary formations is critical for analysing hydrogeological processes and assessing drilling risks. We have constructed a two-dimensional basin model and have performed numerical simulations to increase the understanding of the history of fluid flow and shallow overpressures in the Pleistocene and Holocene formations in the Ursa basin, deepwater Gulf of Mexico. We measured physical properties of sediments, such as porosity and permeability, in the laboratory and estimated in situ pore pressures from preconsolidation pressures. We obtained porosity–effective stress relationships from measurements of bulk density, grain density and preconsolidation pressures in the laboratory. Porosity–effective stress relationships were also obtained from downhole density logs and measured pore pressures. The porosity–effective stress and porosity–permeability relationships obtained were applied in two-dimensional basin simulations. Results showed that high pore pressures developed shortly after sediment deposition. Peaks in pore pressure ratios were related to high sedimentation rates of mass transport deposits and the incision of the Ursa channel. Lateral flows from the area where the overburden is thick towards the area where it is thin have occurred at least since 30 ka. Present pore pressure and temperature distributions suggest that lateral flows play a role in re-distributing heat in the basin.     

Accurate ocean tide modeling in southeast Alaska and large tidal dissipation around Glacier Bay

An accurate prediction of ocean tides in southeast Alaska is developed using a regional, barotropic ocean model with a finite difference scheme. The model skill is verified by the observational tidal harmonics in southeast Alaska including Glacier Bay. The result is particularly improved in Glacier Bay compared to the previous model described by Foreman et al. (2000). The model bathymetry dominates the model skill. We re-estimate tidal energy dissipation in the Alaska Panhandle and suggest a value for tidal energy dissipation of 3.4 GW associated with the M₂ constituent which is 1.5 times the estimation of Foreman et al. (2000). A large portion of the M₂ energy budget entering through Chatham Strait is dissipated in the vicinity of Glacier Bay. Moreover, it is shown that the developed model has the potential to correct the ocean tide loading effect in geodetic data more efficiently than the model of Foreman et al. (2000), especially around Glacier Bay.     

东海风应力旋度驱动的21天周期的海底压力变化及其传播

在2015年6月至2017年6月期间,跨越庆良间水道的由2台加载压力传感器的倒置式回声仪（PIES）和5台加载压力传感器和海流计的倒置式回声仪（CPIES）组成的观测断面获得了近2年的海底压力时间序列。该时间序列中存在着显著的21天周期的振荡（P_bot21）,该信号在2016年7月至10月期间尤为强烈。P_bot21与东海陆架上的风应力旋度存在较显著的3天延迟相关,其相关系数达到0.65。本文采用正压海洋模式解释了这一信号的产生、传播以及耗散过程,模式结果显示东海陆架上的风应力旋度驱动产生P_bot21并向琉球岛链传播,而深海上的风应力旋度不能驱动产生这一信号。在陆架上,P_bot21伴随21天周期的风应力旋度由海岸向东南方向传播,但由于摩擦作用,信号在离开风场后几天内即耗散。断面能否观测到P_bot21与陆架上21天周期风应力旋度场的分布相关,长江口东南方向风应力旋度驱动的P_bot21能被观测到,而长江口东北方向产生的P_bot21不能被观测到。     

Predictability of the subtropical dipole modes in a coupled ocean–atmosphere model

Predictability of the subtropical dipole modes is assessed using the SINTEX-F coupled model. Despite the known difficulty in predicting subtropical climate due to large internal variability of the atmosphere and weak ocean–atmosphere coupling, it is shown for the first time that the coupled model can successfully predict the South Atlantic Subtropical Dipole (SASD) 1 season ahead, and the prediction skill is better than the persistence in all the 1–12 month lead hindcast experiments. There is a prediction barrier in austral winter due to the seasonal phase locking of the SASD to austral summer. The prediction skill is lower for the Indian Ocean Subtropical Dipole (IOSD) than for the SASD, and only slightly better than the persistence till 6-month lead because of the low predictability of the sea surface temperature anomaly in its southwestern pole. However, for some strong IOSD events in the last three decades, the model can predict them 1 season ahead. The co-occurrence of the negative SASD and IOSD in 1997/1998 austral summer can be predicted from July 1st of 1997. This is because the negative sea level pressure anomalies over the South Atlantic and the southern Indian Ocean in September–October (November–December) that trigger the occurrence of the negative SASD and IOSD are related to the well predicted tropical Indian Ocean Dipole (El Niño/Southern Oscillation). Owing to the overall good performances of the SINTEX-F model in predicting the SASD, some strong IOSD, and El Niño/Southern Oscillation, the prediction skill of the southern African summer precipitation is high in the SINTEX-F model.     

Locally and remotely forced atmospheric circulation anomalies of Ningaloo Niño/Niña

A recently identified climate mode called Ningaloo Niño (Niña) is associated with positive (negative) sea surface temperature (SST) anomalies off the west coast of Australia and negative (positive) sea level pressure (SLP) anomalies in the overlying atmosphere. By conducting a series of numerical experiments with an atmospheric general circulation model, generation mechanisms of the atmospheric circulation anomalies accompanied by Ningaloo Niño/Niña are examined. Even when SST is allowed to vary interannually only in the eastern South Indian Ocean, negative (positive) SLP anomalies are formed off the west coast of Australia in Ningaloo Niño (Niña) years, supporting the existence of local ocean–atmosphere interaction. When the model is forced by SST anomalies outside of the eastern South Indian Ocean, negative (positive) SLP anomalies are also generated in Ningaloo Niño (Niña) years owing to a Matsuno–Gill type response to atmospheric convection anomalies in the tropical Pacific. It is found that the latter impact is stronger in the current atmospheric general circulation model. Regarding climatic impacts, it is shown that Ningaloo Niño (Niña) induces wet (dry) anomalies over the northwestern part of Australia even when SST anomalies outside of the eastern South Indian Ocean are excluded from the SST forcing.     

Flux‐Free Fusion of Silicate Rock Preceding Acid Digestion for ICP‐MS Bulk Analysis

We present a new method for the decomposition of silicate rocks by flux‐free fusion in preparation for whole‐rock trace element determination (Sc, Rb, Sr, Y, Zr, Nb, Cs, Ba, rare earth elements and Hf) that is especially applicable to zircon‐bearing felsic rocks. The method was verified by analyses of RMs of mafic (JB‐1a, JB‐2, JGb‐1) and felsic rocks (JG‐3, JR‐3, JSd‐1, GSP‐2, G‐2). Pellets of powdered sample (up to 500 mg) without flux were weighed and placed in a clean platinum crucible. The samples were then fused in a Siliconit® tube furnace and quenched to room temperature. The optimum condition for the fusion of granitic rock was determined to be heating for 2–3 min at 1600 °C. The fused glass in the platinum crucible after heating was decomposed using HF and HClO₄ in a Teflon® beaker. Decomposed and diluted sample solutions were analysed using a quadrupole inductively coupled plasma‐mass spectrometer. Replicate analyses (n = 4 or 5) of the RMs revealed that analytical uncertainties were generally < 3% for all elements except Zr and Hf (～ 6%) in JG‐3. These higher uncertainties may be attributed to sample heterogeneity. Our analytical results for the RMs agreed well with recommended concentrations and recently published concentrations, indicating complete decomposition of our rock samples during fusion.     

Slab decarbonation and CO2 recycling in the Southwestern Colombian volcanic arc

The contribution of subducted carbonate sediments to the genesis of the Southwestern Colombian arc magmas was investigated using a comprehensive petrography and geochemical analysis, including determination of major and trace element contents and Sr, Nd, Hf and Pb isotope compositions. These data have been used to constrain the depth of decarbonation in the subducted slab, indicating that the decarbonation process continues into the sub-arc region, and ultimately becomes negligible in the rear arc. We propose on the basis of multi-isotope approach and mass balance calculations, that the most important mechanism to induce the slab decarbonation is the infiltration of chemically reactive aqueous fluids from the altered oceanic crust, which decreasingly metasomatize the mantle wedge, triggering the formation of isotopically different primary magmas from the volcanic front (VF) with relatively high ¹⁷⁶Hf/¹⁷⁷Hf, high ⁸⁷Sr/⁸⁶Sr, negative values of εNd and lower Pb isotopes compared to the rear arc (RA).The presence of more aqueous fluids at the volcanic front may increase the degree of decarbonation into carbonate-bearing lithologies. Moreover, with increasing pressure and temperature in the subduction system, the decrease in dehydration of the slab, leads to cessation of fluid-induced decarbonation reactions at the rear arc. This development allows the remaining carbonate materials to be recycled into the deep mantle.