Diverse tsunamigenesis triggered by the Hunga Tonga-Hunga Ha'apai eruption.

On the evening of 15 January 2022, the Hunga Tonga-Hunga Ha'apai volcano1 unleashed a violent underwater eruption, blanketing the surrounding land masses in ash and debris2,3. The eruption generated tsunamis observed around the world. An event of this type last occurred in 1883 during the eruption of Krakatau4, and thus we have the first observations of a tsunami from a large emergent volcanic eruption captured with modern instrumentation. Here we show that the explosive eruption generated waves through multiple mechanisms, including: (1) air-sea coupling with the initial and powerful shock wave radiating out from the explosion in the immediate vicinity of the eruption; (2) collapse of the water cavity created by the underwater explosion; and (3) air-sea coupling with the air-pressure pulse that circled the Earth several times, leading to a global tsunami. In the near field, tsunami impacts are strongly controlled by the water-cavity source whereas the far-field tsunami, which was unusually persistent, can be largely described by the air-pressure pulse mechanism. Catastrophic damage in some harbours in the far field was averted by just tens of centimetres, implying that a modest sea level rise combined with a future, similar event would lead to a step-function increase in impacts on infrastructure. Piecing together the complexity of this event has broad implications for coastal hazards in similar geophysical settings, suggesting a currently neglected source of global tsunamis.

Tsunami Runup and Inundation in Tonga from the January 2022 Eruption of Hunga Volcano.

On January 15th, 2022, at approximately 4:47 pm local time (0347 UTC), several weeks of heightened activity at the Hunga volcano 65 km northwest of Tongatapu, culminated in an 11-h long violent eruption which generated a significant near-field tsunami. Although the Kingdom of Tonga lies astride a large and tsunamigenic subduction zone, it has relatively few records of significant tsunami. Assessment activities took place both remotely and locally. Between March and June 2022, a field team quantified tsunami runup and inundation on the main populated islands Tongatapu and Eua, along with several smaller islands to the north, including the Ha'apai Group. Peak tsunami heights were ~ 19 m in western Tongatapu, ~ 20 m on south-eastern Nomuka Iki island and ~ 20 m on southern Tofua, located ~ 65 km S and E and 90 km N from Hunga volcano, respectively. In western Tongatapu, the largest tsunami surge overtopped a 13-15 m-high ridge along the narrow Hihifo peninsula in several locations. Analysis of tide gauge records from Nukualofa (which lag western Tongatapu arrivals by ~ 18-20 min), suggest that initial tsunami surges were generated prior to the largest volcanic explosions at ~ 0415 UTC. Further waves were generated by ~ 0426 UTC explosions that were accompanied by air-pressure waves. Efforts to model this event are unable to reproduce the timing of the large tsunami wave that toppled a weather station and communication tower on a 13 m-high ridge on western Tongatapu after 0500 UTC. Smaller tsunami waves continued until ~ 0900, coincident with a second energetic phase of eruption, and noted by eyewitnesses on Tungua and Mango Islands. Despite an extreme level of destruction caused by this tsunami, the death toll was extraordinarily low (4 victims). Interviews with witnesses and analysis of videos posted on social media suggest that this can be attributed to the arrival of smaller 'pre tsunami' waves that prompted evacuations, heightened tsunami awareness due to tsunami activity and advisories on the day before, the absence of tourists and ongoing tsunami education efforts since the 2009 Niuatoputapu, Tonga tsunami. This event highlights an unexpectedly great hazard from volcanic tsunami worldwide, which in Tonga's case overprints an already extreme level of tectonic tsunami hazard. Education and outreach efforts should continue to emphasize the 'natural warning signs' of strong ground shaking and unusual wave and current action, and the importance of self-evacuation from coastal areas of low-lying islands. The stories of survival from this event can be used as global best practice for personal survival strategies from future tsunami. Supplementary Information: The online version contains supplementary material available at 10.1007/s00024-022-03215-5.

Swept away: ocean currents and seascape features influence genetic structure across the 18,000 Km Indo-Pacific distribution of a marine invertebrate, the black-lip pearl oyster Pinctada margaritifera.

BACKGROUND: Genetic structure in many widely-distributed broadcast spawning marine invertebrates remains poorly understood, posing substantial challenges for their fishery management, conservation and aquaculture. Under the Core-Periphery Hypothesis (CPH), genetic diversity is expected to be highest at the centre of a species' distribution, progressively decreasing with increased differentiation towards outer range limits, as populations become increasingly isolated, fragmented and locally adapted. The unique life history characteristics of many marine invertebrates such as high dispersal rates, stochastic survival and variable recruitment are also likely to influence how populations are organised. To examine the microevolutionary forces influencing population structure, connectivity and adaptive variation in a highly-dispersive bivalve, populations of the black-lip pearl oyster Pinctada margaritifera were examined across its ~18,000 km Indo-Pacific distribution. RESULTS: Analyses utilising 9,624 genome-wide SNPs and 580 oysters, discovered differing patterns of significant and substantial broad-scale genetic structure between the Indian and Pacific Ocean basins. Indian Ocean populations were markedly divergent (F st = 0.2534-0.4177, p < 0.001), compared to Pacific Ocean oysters, where basin-wide gene flow was much higher (F st = 0.0007-0.1090, p < 0.001). Partitioning of genetic diversity (hierarchical AMOVA) attributed 18.1% of variance between ocean basins, whereas greater proportions were resolved within samples and populations (45.8% and 35.7% respectively). Visualisation of population structure at selectively neutral loci resolved three and five discrete genetic clusters for the Indian and Pacific Oceans respectively. Evaluation of genetic structure at adaptive loci for Pacific populations (89 SNPs under directional selection; F st = 0.1012-0.4371, FDR = 0.05), revealed five clusters identical to those detected at neutral SNPs, suggesting environmental heterogeneity within the Pacific. Patterns of structure and connectivity were supported by Mantel tests of isolation by distance (IBD) and independent hydrodynamic particle dispersal simulations. CONCLUSIONS: It is evident that genetic structure and connectivity across the natural range of P. margaritifera is highly complex, and produced by the interaction of ocean currents, IBD and seascape features at a broad scale, together with habitat geomorphology and local adaptation at regional levels. Overall population organisation is far more elaborate than generalised CPH predictions, however valuable insights for regional fishery management, and a greater understanding of range-wide genetic structure in a highly-dispersive marine invertebrate have been gained.

Understanding marine larval dispersal in a broadcast-spawning invertebrate: A dispersal modelling approach for optimising spat collection of the Fijian black-lip pearl oyster Pinctada margaritifera.

Fisheries and aquaculture industries worldwide remain reliant on seed supply from wild populations, with their success and sustainability dependent on consistent larval recruitment. Larval dispersal and recruitment in the marine environment are complex processes, influenced by a multitude of physical and biological factors. Biophysical modelling has increasingly been used to investigate dispersal and recruitment dynamics, for optimising management of fisheries and aquaculture resources. In the Fiji Islands, culture of the black-lip pearl oyster (Pinctada margaritifera) is almost exclusively reliant on wild-caught juvenile oysters (spat), through a national spat collection programme. This study used a simple Lagrangian particle dispersal model to investigate current-driven larval dispersal patterns, identify potential larval settlement areas and compare simulated with physical spat-fall, to inform targeted spat collection efforts. Simulations successfully identified country-wide patterns of potential larval dispersal and settlement from 2012-2015, with east-west variations between bi-annual spawning peaks and circulation associated with El Niño Southern Oscillation. Localised regions of larval aggregation were also identified and compared to physical spat-fall recorded at 28 spat collector deployment locations. Significant and positive correlations at these sites across three separate spawning seasons (r(26) = 0.435; r(26) = 0.438; r(26) = 0.428 respectively, p = 0.02), suggest high utility of the model despite its simplicity, for informing future spat collector deployment. Simulation results will further optimise black-lip pearl oyster spat collection activity in Fiji by informing targeted collector deployments, while the model provides a versatile and highly informative toolset for the fishery management and aquaculture of other marine taxa with similar life histories.

Keeping It Local: Dispersal Limitations of Coral Larvae to the High Latitude Coral Reefs of the Houtman Abrolhos Islands.

In 2011 the first recorded bleaching event for the high latitude Houtman Abrolhos Islands (HAI) coral communities was documented. This bleaching event highlighted the question of whether a supply of 'heat tolerant' coral recruits from the tropical north would be sufficient to provide a level of resistance for these reefs to future warming events. Using Lagrangian modelling we showed that due to its regional isolation, large-scale larval input from potential tropical northern source populations to the HAI is unlikely, despite the southward flowing Leeuwin current. Successful recruitment to artificial substrates was recorded following the bleaching event. However, this was negligible (0.4 ± 0.1 recruits per tile) compared to 2013 post impact recruitment (128.8 ± 15.8 recruits per tile). Our data therefore provides preliminary evidence suggesting that the connectivity of the HAI with coral communities in the north is limited, and population maintenance and recovery is likely driven primarily by self-recruitment. Given the low thermal tolerance of the HAI coral communities, the dominance of Acropora, and the apparent reliance on self-recruitment, an increased frequency of thermally anomalous conditions at the HAI (such as experienced in 2011) has the potential to reduce the long-term stability of the HAI coral populations and species that depend upon them.

A Parallel Population Genomic and Hydrodynamic Approach to Fishery Management of Highly-Dispersive Marine Invertebrates: The Case of the Fijian Black-Lip Pearl Oyster Pinctada margaritifera.

Fishery management and conservation of marine species increasingly relies on genetic data to delineate biologically relevant stock boundaries. Unfortunately for high gene flow species which may display low, but statistically significant population structure, there is no clear consensus on the level of differentiation required to resolve distinct stocks. The use of fine-scale neutral and adaptive variation, considered together with environmental data can offer additional insights to this problem. Genome-wide genetic data (4,123 SNPs), together with an independent hydrodynamic particle dispersal model were used to inform farm and fishery management in the Fijian black-lip pearl oyster Pinctada margaritifera, where comprehensive fishery management is lacking, and the sustainability of exploitation uncertain. Weak fine-scale patterns of population structure were detected, indicative of broad-scale panmixia among wild oysters, while a hatchery-sourced farmed population exhibited a higher degree of genetic divergence (Fst = 0.0850-0.102). This hatchery-produced population had also experienced a bottleneck (NeLD = 5.1; 95% C.I. = [5.1-5.3]); compared to infinite NeLD estimates for all wild oysters. Simulation of larval transport pathways confirmed the existence of broad-scale mixture by surface ocean currents, correlating well with fine-scale patterns of population structuring. Fst outlier tests failed to detect large numbers of loci supportive of selection, with 2-5 directional outlier SNPs identified (average Fst = 0.116). The lack of biologically significant population genetic structure, absence of evidence for local adaptation and larval dispersal simulation, all indicate the existence of a single genetic stock of P. margaritifera in the Fiji Islands. This approach using independent genomic and oceanographic tools has allowed fundamental insights into stock structure in this species, with transferability to other highly-dispersive marine taxa for their conservation and management.