Prediction of earthquake energy input from smoothed fourier amplitude spectrum

This paper attempts to show analytically that the energy-input spectra of damped SDOF systems and undamped MDOF systems excited by an earthquake motion can be predicted by smoothing the Fourier amplitude spectrum of the base acceleration. The spectral window for smoothing in the frequency domain for a damped SDOF system is identical with the probability density function of the time-variant or instantaneous vibration frequency resulting from non-linear hysteresis. The spectral window for an undamped MDOF system is identical with the set of squared participation factors associated with vibration modes. It was found that the increase in damping factor and the increase in participation of higher modes provide wider spectral windows, resulting in more flattened or unaltered energy-input spectra due to enhanced smoothing effects.     

Long-term high-frequency measurements of dibromomethane in the atmosphere at algae-rich and algae-poor coastal sites

Dibromomethane (CH₂Br₂), a natural stratospheric ozone depleting substance, is mostly emitted from the ocean, but the relative importance of coastal (or macroalgae) and open ocean emissions is unknown. We made long-term high-frequency measurements of CH₂Br₂ concentrations at two remote coastal sites in Japan, on the subtropical Hateruma Island (poor in macroalgae) and at Cape Ochiishi (rich in macroalgae). CH₂Br₂ concentrations at Hateruma showed prominent seasonal variation, being lower in summer (around 0.94 ppt) than in winter (around 1.23 ppt). In contrast, CH₂Br₂ concentrations at Ochiishi were highly variable, often exceeding 2 ppt in the summer but with minimum baseline concentrations close to those from Hateruma; in the winter the concentrations were almost constant at about 1.3 ppt. Analysis of the data suggested that (1) emissions from macroalgae were not likely to extend offshore, but instead were localized near the shore, (2) strong macroalgal emissions of CH₂Br₂ were almost limited to the summer, but it was not reflected in the seasonality of the baseline concentrations of CH₂Br₂ in the atmosphere, and therefore (3) macroalgal or coastal emissions of CH₂Br₂ in the temperate zone might have a rather limited contribution to the global CH₂Br₂ sources. These findings are especially important for the understanding of the tropospheric and stratospheric bromine budget.     

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.     

Potential of submarine-cave sediments and oxygen isotope composition of cavernicolous micro-bivalve as a late Holocene paleoenvironmental record

A sediment layer (43 cm thick) and surface sediments (5 cm thick) in a submarine limestone cave (31 m water depth) on the fore-reef slope of Ie Island, off Okinawa mainland, Japan, were examined by visual, mineralogical and geochemical means. Oxygen isotope analysis was performed on the cavernicolous micro-bivalve Carditella iejimensis from both cored sediments and surface sediments, and the water temperature within the cave was recorded for nearly one year. These data show that: (1) water temperature within the cave is equal to that at 30 m deep in the open sea; (2) the biotic and non-biotic environments within the cave have persisted for the past 2000 years; (3) mud-size carbonate detritus is a major constituent of the submarine-cave deposit, and may have come mainly from the suspended carbonate mud produced on the emergent Holocene reef flat over the past two millennia; (4) the δ¹⁸O-derived temperature (T_δ18_O) of C. iejimensis suggests that the species grows between April and July; (5) the T_δ18_O of C. iejimensis from cored sediments implies that there were two warmer intervals, at AD 340 ± 40 and AD 1000 ± 40, which correspond to the Roman Warm Period and Medieval Warm Period, respectively. These suggest that submarine-cave sediments provide unique information for Holocene reef development. In addition, oxygen isotope records of cavernicolous C. iejimensis are a useful tool to reconstruct century-scale climatic variability for the Okinawa Islands during the Holocene.     

Shock-wave heating model for chondrule formation: Hydrodynamic simulation of molten droplets exposed to gas flows

Millimeter-sized, spherical silicate grains abundant in chondritic meteorites, which are called as chondrules, are considered to be a strong evidence of the melting event of the dust particles in the protoplanetary disk. One of the most plausible scenarios is that the chondrule precursor dust particles are heated and melt in the high-velocity gas flow (shock-wave heating model). We developed the non-linear, time-dependent, and three-dimensional hydrodynamic simulation code for analyzing the dynamics of molten droplets exposed to the gas flow. We confirmed that our simulation results showed a good agreement in a linear regime with the linear solution analytically derived by Sekyia et al. [Sekyia, M., Uesugi, M., Nakamoto, T., 2003. Prog. Theor. Phys. 109, 717-728]. We found that the non-linear terms in the hydrodynamical equations neglected by Sekiya et al. [Sekiya, M., Uesugi, M., Nakamoto, T., 2003. Prog. Theor. Phys. 109, 717-728] can cause the cavitation by producing negative pressure in the droplets. We discussed that the fragmentation through the cavitation is a new mechanism to determine the upper limit of chondrule sizes. We also succeeded to reproduce the fragmentation of droplets when the gas ram pressure is stronger than the effect of the surface tension. Finally, we compared the deformation of droplets in the shock-wave heating with the measured data of chondrules and suggested the importance of other effects to deform droplets, for example, the rotation of droplets. We believe that our new code is a very powerful tool to investigate the hydrodynamics of molten droplets in the framework of the shock-wave heating model and has many potentials to be applied to various problems.     

Geochemical relation between gas hydrates and water venting in the seismic blanking zone on the northern Cascadia continental margin offshore Vancouver Island, Canada

The isotopic composition (δD and δ¹⁸O) and chloride concentration (Cl⁻) of pore waters from the northern Cascadia continental margin offshore Vancouver Island were measured to characterize the relations between the water flow regime and the distribution, formation and dissociation of gas hydrates. The δD values of pore waters in gas hydrate-bearing sediments are slightly higher ( 1‰) than those of seawater as the result of gas hydrate dissociation during core recovery and handling. Within the seismic blanking zone, the δD values were slightly lower (− 1‰) than values measured from sites outside the blanking area (0‰). We attribute these differences to 1) distillation of D-rich water during hydrate formation in the center of the blanking zone and 2) limited migration of pore water between inside and outside of the blanking zone due to different fluid fluxes. In contrast, the δ¹⁸O values and Cl⁻ concentrations do not show significant spatial variation due to decreased isotopic fractionation of oxygen and small fraction of chloride relative to hydrogen isotope during gas hydrate formation. The δD value of pore water, therefore, appears to be a sensitive indicator of gas hydrate occurrence. We estimate that gas hydrate occupied at least 2.0 to 6.3% of sediment pore space using δD distribution in this area.     

Effect of bed roughness on turbulent boundary layer and net sediment transport under asymmetric waves

This paper presents an investigation of the roughness effects in the turbulent boundary layer for asymmetric waves by using the baseline (BSL) k–ω model. This model is validated by a set of the experimental data with different wave non-linearity index, N_i (namely, N_i = 0.67, N_i = 0.60 and N_i = 0.58). It is further used to simulate asymmetric wave velocity flows with several values of the roughness parameter (a_m/k_s) which increase gradually, namely from a_m/k_s = 35 to a_m/k_s = 963. The effect of the roughness tends to increase the turbulent kinetic energy and to decrease the mean velocity distribution in the inner boundary layer, whereas in the outer boundary layer, the roughness alters the turbulent kinetic energy and the mean velocity distribution is relatively unaffected. A new simple calculation method of bottom shear stress based on incorporating velocity and acceleration terms is proposed and applied into the calculation of the rate of bed-load transport induced by asymmetric waves. And further, the effect of bed roughness on the bottom shear stress and bed-load sediment transport under asymmetric waves is examined with the turbulent model, the newly proposed method, and the existing calculation method. It is found that the higher roughness elements increase the magnitude of bottom shear stress along a wave cycle and consequently, the potential net sediment transport rate. Moreover, the wave non-linearity also shows a big impact on the bottom shear stress and the net sediment transport.     

Triplications for the converted wave in transversely isotropic media with a tilted symmetry axis

In seismic data processing, serious problems could be caused by the existence of triplication and need to be treated properly for tomography and other inversion methods. The triplication in transversely isotropic medium with a vertical symmetry axis has been well studied and concluded that the triplicated traveltime only occurs for S wave and there is no triplication for P and converted PS waves since the P wave convexity slowness always compensates the S wave slowness concavity. Compared with the vertical symmetry axis model, the research of the triplication in transversely isotropic medium with a tilted symmetry axis is still keeping blank. In order to analyse the triplication for the converted wave in the tilted symmetry axis model, we examine the traveltime of the triplication from the curvature of averaged P and S wave slowness. Three models are defined and tested in the numerical examples to illustrate the behaviour of the tilted symmetry axis model for the triplicated traveltime with the change of the rotation angle. Since the orientation of an interface is related to the orientation of the symmetry axis, the triplicated traveltime is encountered for the converted wave in the tilted symmetry axis model assuming interfaces to be planar and horizontal. The triplicated region is influenced by the place and level of the concave curvature of the P and S wave slowness.     

Similarity solution for the interaction of the stellar wind with surrounding interstellar medium

Similarity solution for the interaction of stellar wind with surrounding interstellar medium is obtained on the base of momentum conservation. The conversion efficiency of kinetic energy of stellar wind into the kinetic energy of expanding shell is derived. The results are compared with the observations on the ring nebulae associated with the Wolf-Rayet star and it is shown that the observed values of energy conversion are in good agreement with the momentum conserving model.     

Empirical formula to relate the auroral electrojet intensity with interplanetary parameters

An empirical formula has been constructed using the results of correlative analyses to determine in what form the AL index, as a measure of the intensity of the westward auroral electrojet, depends on interplanetary parameters. The formula thus obtained shows that AL is mainly determined by B_sV² where B_s is the southward component of the IMF and V is the solar wind velocity, and is modulated in a characteristic way by the combined effect of the east-west component of the IMF and the tilt angle of the Earth's dipole axis toward the Sun-Earth line. In contrast, effects of the solar wind density and the IMF variability were found to be insignificant.Implications of the empirical formula are discussed mainly in relation to the problem of the location in the dayside magnetosphere of the region where the reconnection process to initiate the substorm takes place.