Characteristic features in Venus’ nightside cloud-top temperature obtained by Akatsuki/LIR

Mid-infrared images of almost the entire Venus nightside hemisphere obtained by the Longwave Infrared Camera (LIR) onboard Akatsuki on December 9 and 10, 2010 reveal that the brightness temperature of the cloud-top ranges from 237 K in the cold polar collars to 243 K in the equatorial region, significantly higher than the values obtained by Venera 15. Other characteristic features of the temperature distributions observed are zonal belt structures seen in the middle and low latitudes and patchy temperature structures or quasi-periodic streaks extending in a north–south direction in the northern middle latitudes and southern low latitudes.     

Re-examination of the Source Mechanism of the 1998 Papua New Guinea Earthquake and Tsunami

— Simulation of tsunami propagation and runup of the 1998 Papua New Guinea (PNG) earthquake tsunami using the detailed bathymetry measured by JAMSTEC and adding bathymetric data at depths less than 60 m is carried out, reproducing the tsunami energy focus into Warapu and Arop along the Sissano Lagoon. However, the computed runup heights in the lagoon are still lower than those measured. Even if the error in estimating the fault parameters is taken into consideration, computational results are similar. Analysis by the wave ray method using several scenarios of the source size of the tsunami and location by the wave ray method suggests that a source characterized by small size in water 1,000-m deep approximately 25 km offshore the lagoon, best fits the arrival determined from the interviews with eyewitnesses. A two-layer numerical model simulating the interaction of the tsunami with a landslide is employed to study the behavior of a landslide-generated tsunami with different size sand depths of the initial slide just outside the lagoon. A landslide model with a volume of 4–8 × 10⁹ m³ is selected as the best in order to reproduce the distribution of the measured tsunami runup in the lagoon. The simulation of a tsunami generated in two stages, fault and landslide, could show good agreement with the runup heights and distribution of the arrival time, but a time gap of around 10 minutes remains, suggesting that a tsunami generated by the mainshock at 6:49 PM local time is too small for people to notice, and the following tsunami triggered by landslide or mass movement near the lagoon about ten minutes after the mainshock attacked the coast and caused the huge damage.     

Damage to coastal villages due to the 1992 Flores Island earthquake tsunami

A field survey of the 1992 Flores Island earthquake tsunami was conducted during December 29, 1992 to January 5, 1993 along the north coast of the eastern part of Flores Island. We visited over 40 villages, measured tsunami heights, and interviewed the inhabitants. It was clarified that the first wave attacked the coast within five minutes at most of the surveyed villages. The crust was uplifted west of the Cape of Batumanuk, and subsided east of it. In the residential area of Wuring, which is located on a sand spit with ground height of 2 meters, most wooden houses built on stilts collapsed and 87 people were killed even though the tsunami height reached only 3.2 meters. In the two villages on Babi Island, the tsunami swept away all wooden houses and killed 263 of 1,093 inhabitants. Tsunami height at Riang-Kroko village on the northeastern end of Flores Island reached 26.2 meters and 137 of the 406 inhabitants were killed by the tsumani. Evidence of landslides was detected at a few points on the coast of Hading Bay, and the huge tsunami was probably formed by earthquake-induced landslides. The relationship between tsunami height and mortality was checked for seven villages. The efficiencies of trees arranged in front of coastal villages, and coral reefs in dissipating the tsunami energy are discussed.     

Numerical identification of tsunami boulders and estimation of local tsunami size at Ibaruma reef of Ishigaki Island,Japan

Tsunami boulders deposited along the coast constitute important geological evidence for paleotsunami activity. However, boulders can also be deposited by large storm waves. Although several sedimentological and theoretical methods have been proposed to differentiate tsunami and storm wave affected boulders, no appropriate numerical method exists for their differentiation. Therefore, we developed a new numerical scheme to differentiate tsunami and storm wave boulders for coastal boulders on Ishigaki Island, Japan. In this area, tsunami and storm waves have emplaced numerous boulders on the reef and the coast. By conducting numerical calculations of storm waves in this region, we estimated the size of a storm wave that can explain the maximum clast size distribution of boulders on the reef. Consequently, we showed that a wave with a combination of 8 m in initial wave height and 10 s period can satisfy the above conditions when we assume mean sea level. In contrast to the boulders on the reef, all boulders deposited along the shore are heavier than the calculated possible maximum clast size distribution by the storm wave. Therefore, we confirmed these boulders as being of tsunami origin. Results of previous studies showed that they were most likely deposited or reworked by the 1771 Meiwa tsunami. Then, using the tsunami boulders, we numerically estimated the wave period and amplitude of the 1771 Meiwa tsunami, which should have had a 4–5 min period and 5.6–5.9, 6.3–7.0 m amplitude, respectively. Using the proposed scheme, it is possible to differentiate tsunami and storm wave boulders and estimate the size of past storm waves and tsunami waves, although it is noteworthy that there are exceptions for which the scheme cannot be applied.     

Modeling boulder transport by coastal waves on cliff topography: Case study at Hachijo Island,Japan

Clifftop coastal boulders transported by storm waves or tsunamis have been reported around the world. Although numerical calculation of boulder transport is a strong tool for the identification of tsunami or storm boulders, and for estimation of the wave size emplacing boulders, models which can reasonably solve boulder transport from below a cliff or from a cliff-edge onto a cliff-top do not yet exist. In this study, we developed a new numerical formulation for cliff-top deposition of boulders from the cliff edge or below the cliff, with validation from laboratory tests. We then applied the model using storm and tsunami wave forcing to simulate the observed boulder deposits at the northwest coast of Hachijo Island, Japan. Using the model, the actual distribution of boulders was explained well using a reasonable storm wave height without assumption of anomalously high-water level by storm surge. Results show that boulder transport from the cliff edge or under the cliff onto the cliff-top was possible from a tsunami with periods of 5~10 min or storm waves with no storm surge. However, the actual distribution of boulders on the cliff was explained only by storm waves, but not by tsunami. Therefore, the boulders distributed at this site are likely of storm wave origin. Our developed model for the boulder transport calculation can be useful for identifying a boulder's origin and can reasonably calculate cliff-top deposition of boulders by tsunami and storm waves. © 2019 John Wiley & Sons, Ltd. © 2019 John Wiley & Sons, Ltd.     

The reduction effects of mangrove forest on a tsunami based on field surveys at Pakarang Cape,Thailand and numerical analysis

Using an integrated approach including satellite imagery analysis, field measurements, and numerical modeling, we investigated the damage to mangroves caused by the 2004 Indian Ocean tsunami at Pakarang Cape in Pang Nga Province, Thailand. Comparing pre- and post-tsunami satellite imagery of the study area, we found that approximately 70% of the mangrove forest was destroyed by the tsunami. Based on field observations, we found that the survival rate of mangroves increased with increasing stem diameter. Specifically, we found that 72% of Rhizophora trees with a 25–30 cm stem diameter survived the tsunami impact, whereas only 19% with a 15–20 cm stem diameter survived. We simulated the 2004 Indian Ocean tsunami using the nonlinear shallow-water wave theory to reproduce the tsunami inundation flow and investigated the bending moment acting on the mangrove trees. Results of the numerical model showed that the tsunami inundated areas along the mangrove creeks, and its current velocity reached 5.0 m s⁻¹. Based on the field measurements and numerical results, we proposed a fragility function for mangroves, which is the relationship between the probability of damage and the bending stress caused by the maximum bending moment. We refined the numerical model to include the damage probability of mangrove forests using the obtained fragility function to investigate the tsunami reduction effect of mangrove forest. Under simple numerical conditions related to the mangrove forest, ground level, and incident wave, the model showed that a mangrove forest of Rhizophora sp. with a density of 0.2 trees m⁻² and a stem diameter of 15 cm in a 400 m wide area can reduce the tsunami inundation depth by 30% when the incident wave is assumed to have a 3.0 m inundation depth and a wave period of 30 min at the shoreline. However, 50% of the mangrove forest is destroyed by a 4.5 m tsunami inundation depth, and most of the mangrove forest is destroyed by a tsunami inundation depth greater than 6 m. The reduction effect of tsunami inundation depth decreased when the tsunami inundation depth exceeded 3 m, and was mostly lost when the tsunami inundation depth exceeded 6 m.     

Throughfall partitioning by trees

Although we know that rainfall interception (the rain caught, stored, and evaporated from aboveground vegetative surfaces and ground litter) is affected by rain and throughfall drop size, what was unknown until now is the relative proportion of each throughfall type (free throughfall, splash throughfall, canopy drip) beneath coniferous and broadleaved trees. Based on a multinational data set of >120 million throughfall drops, we found that the type, number, and volume of throughfall drops are different between coniferous and broadleaved tree species, leaf states, and timing within rain events. Compared with leafed broadleaved trees, conifers had a lower percentage of canopy drip (51% vs. 69% with respect to total throughfall volume) and slightly smaller diameter splash throughfall and canopy drip. Canopy drip from leafless broadleaved trees consisted of fewer and smaller diameter drops (D_{50_DR}, 50th cumulative drop volume percentile for canopy drip, of 2.24 mm) than leafed broadleaved trees (D_{50_DR} of 4.32 mm). Canopy drip was much larger in diameter under woody drip points (D_{50_DR} of 5.92 mm) than leafed broadleaved trees. Based on throughfall volume, the percentage of canopy drip was significantly different between conifers, leafed broadleaved trees, leafless broadleaved trees, and woody surface drip points (p ranged from <0.001 to 0.005). These findings are partly attributable to differences in canopy structure and plant surface characteristics between plant functional types and canopy state (leaf, leafless), among other factors. Hence, our results demonstrating the importance of drop‐size‐dependent partitioning between coniferous and broadleaved tree species could be useful to those requiring more detailed information on throughfall fluxes to the forest floor.     

Multi-decadal changes in the relationships between rainfall characteristics and debris-flow occurrences in response to gully evolution after the 1990–1995 Mount Unzen eruptions

Explosive volcanic eruptions can cause long-term landscape change, leading to increased sediment discharge that continues after the cessation of the eruptions. During the period 1990–1995, eruptions of Mount Unzen, Japan, generated large amounts of pyroclastic material, resulting in 57 debris-flow events during 1991–2018. To investigate changes in the relationships between rainfall characteristics and debris-flow occurrence, we conducted the following: geometric analysis of two gullies (i.e., debris-flow initiation zones) using LiDAR (light detection and ranging)-generated 1 m DEMs (digital elevation models); rainfall analysis, based on the relationship between rainfall duration and mean intensity (i.e., considering the intensity–duration, or ID, threshold); and debris-flow monitoring during 2016–2018. Since 1991, rainfall runoff has caused erosion of the supplied pyroclastic material, generating a channel network consisting of incised gullies. With sufficient rainfall, debris flows formed, accompanied by further gully erosion; this resulted in both vertical and lateral adjustments of the cross-sectional geometry. In the two decades since the eruptions ceased, readily mobilized pyroclastic material has become scarce as the gullies have adjusted to local hydrographic conditions. At the same time, the infiltration capacity of the volcanic flank has increased, reducing the capacity for overland flow. As a result, since 2000, rainfall events with intensities above the ID threshold have occurred; however, the lack of sediment supplied by the gullies appears to have hindered the occurrence and development of debris flows. This suggests that debris flows in volcanically perturbed landscapes may occur at lower rainfall thresholds as long as the corresponding upland channels are evolving as a result of intense overland flow. However, as such channels evolve towards equilibrium geometries, the frequency of debris flows decreases in response to the reduction in sediment availability.