Adjoint waveform tomography, which is an emerging seismic imaging method for the crust- and global-scale problems, has gained popularity in the past and present decade. This study, for first time, applies adjoint waveform tomography to the large volume of seismic data recorded by the densely spaced, permanent monitoring network that covers the entirety of Japan. We develop a heterogeneous shear-wave velocity model of central Japan that agrees with the geology and lithology. The results reduce the time-frequency phase misfit by 16.4% in the 0.02–0.05 Hz frequency band and 6.7% in the 0.033–0.1 Hz band, respectively. We infer that some velocity anomalies resolved in this work would reflect the subsurface structures such as the volcanic fluids, dehydration of the subducted crust, and sedimentary basin. In addition, dense distributions of deep earthquakes are visible beneath the high-velocity blocks estimated in this study. The results of this study suggest the possibility of imaging large scale heterogeneous subsurface structures using waveform tomography with a densely distributed network of permanent seismometers.

The relative roles of parameters governing relative permeability, a crucial property for two-phase fluid flows, are incompletely known. To characterize the influence of viscosity ratio (M) and capillary number (Ca), we calculated relative permeabilities of nonwetting fluids (knw) and wetting fluids (kw) in a 3D model of Berea sandstone under steady-state condition using the lattice Boltzmann method. We show that knw increases and kw decreases as M increases due to the lubricating effect, locally occurred pore-filling behavior, and instability at fluid interfaces. We also show that knw decreases markedly at low Ca (log Ca < −1.25), whereas kw undergoes negligible change with changing Ca. An M–Ca–knw correlation diagram, displaying the simultaneous effects of M and Ca, shows that they cause knw to vary by an order of magnitude. The color map produced is useful to provide accurate estimates of knw in reservoir-scale simulations and to help identify the optimum properties of the immiscible fluids to be used in a geologic reservoir.

Monitoring the hydraulic properties within subsurface fractures is vitally important in the contexts of geoengineering developments and earthquakes. Geophysical observations are promising tools for remote determination of subsurface hydraulic properties; however, quantitative interpretations are hampered by the paucity of relevant geophysical data for fractured rock masses. This study explored simultaneous changes in hydraulic and geophysical properties of natural rock fractures with increasing normal stress and correlated these property changes through coupling experiments and digital fracture simulations. We show that electrical resistivity is linked with permeability and flow area regardless of fracture roughness, whereas elastic wave velocity is roughness dependent. We also are able to categorize fracture flow patterns as aperture-dependent, aperture-independent, or disconnected flows, with transitions at specific stress levels. Elastic wave velocity offers potential for detecting the transition between aperture-dependent flow and aperture-independent flow, and resistivity is sensitive to detect the connection/disconnection of the fracture flow.

We applied a polarization analysis of InSight seismic data to estimate the temporal variation and frequency dependence of the Martian ambient noise field. Low-frequency (<1 Hz) P-waves show a diurnal variation in their dominant back-azimuths that are apparently related to wind and the direction of sunlight in a distant area. Low-frequency Rayleigh waves (0.25–1 Hz) show diurnal variations and a dominant back-azimuth related to the wind direction in a nearby area. Low-frequency signals that are derived mainly from wind may be sensitive to subsurface structure deeper than the lithological boundary derived from an autocorrelation analysis. On the other hand, dominant back-azimuths of high-frequency (>1 Hz) waves point toward the InSight lander, especially in daytime, indicating that wind-induced lander noise is dominant at high frequencies. These results point to the presence of several ambient noise sources as well as geologic structure at the landing site.

Research interest in the Kinki region, southwestern Japan, has been aroused by the frequent occurrence of microearthquake activity that do not always coincide with documented active fault locations. Previous studies in the Kinki region focused mainly on deep, large-scale structures and could not efficiently resolve fine-scale (~10 km) shallow crustal structures. Hence, characterization of the upper crustal structure of this region at an improved spatial resolution is required. From the cross-correlation of the vertical components of the ambient seismic noise data recorded by a densely-distributed seismic array, we estimated Rayleigh wave phase velocities using a frequency domain method. Then, we applied a direct surface wave tomographic method for the measured phase velocity dispersion data to obtain the 3D S-wave velocity model of the Kinki region. The estimated velocity model reveals a NE-SW trending low-velocity structure coinciding with the Niigata-Kobe Tectonic Zone (NKTZ) and the active Biwako-seigan Fault Zone (BSFZ). Also, we identified fine-scale low-velocity structures coinciding with known active faults on the eastern side of the NKTZ, as well as sets of low-velocity structures across the Tanba region, that may be attributable to the weathering effects or activity of unidentified concealed fault zones. Furthermore, sedimentary basins manifest as low-velocity zones extending to depths ranging from ~1.5 to 2 km, correlating with those reported in previous studies. Our results therefore contribute towards fundamental understanding of earthquake faulting as well as tectonic boundary and will be useful for hazard assessment and disaster mitigation.

Geothermal systems consisting of fractures in impermeable rocks are difficult to characterize by in situ methods. In an effort to link characteristics of small-scale and large-scale fractures, this study investigated possible relations between their geophysical parameters. We upscaled the relationship between fracture permeability and formation factor in a laboratory specimen to larger fracture dimensions. Microscopic flow characteristics indicate that this relationship is related to the tortuosity of flow paths. We derived an empirical formula that directly predicts changes in fracture permeability from changes in formation factor. This relation may make it possible to monitor subsurface hydraulic activities through resistivity observations.

Active-source seismic surveys have resolved the fine-scale P-wave velocity (Vp) of the subsurface structure in subduction forearcs. In contrast, the S-wave velocity (Vs) structure is poorly resolved despite its usefulness in understanding rock properties. This study estimates Vp and Vs structures of the Nankai Trough forearc, by applying transdimensional inversion to high-frequency receiver function waveforms. As a result, a thin (~1 km) low-velocity zone (LVZ) is evident at ~6 km depth beneath the sea level, which locates ~3 km seaward from the outer ridge. Based on its high Vp/Vs ratio (~2.5) and comparison to an existing seismic reflection profile, we conclude that this LVZ reflects a high pore pressure zone at the upper portion of the underthrust sediment. We infer that this overpressured underthrust sediment hosts slow earthquake activities, and that accompanied strain release helps impede coseismic rupture propagation further updip.

In an effort to reveal the subsurface hydraulic changes in fractures by seismic monitoring, aperture-related velocity changes need to be investigated. We developed a numerical approach for calculating changes in elastic wave velocity with fracture aperture opening by determining the internal energy of a digitized fracture model based on natural rough surfaces. The simulated local elastic energy revealed that the interaction energy converged within 1.5 mm of the mean fracture position, and was insignificant unless the fractures intersected. This energetic approach clarified the aperture–velocity relationship and reproduced the experimental results. Further calculations using digital fractures with various sizes and density demonstrated that the velocity can be accounted for by the superposition of a linear function of fracture density and quadratic function of aperture, and is insensitive to the fracture size. Although the relationship between fracture permeability and elastic wave velocity (i.e., the k-V relationship) depends on the fracture density, the offset-normalized k-V relationship shows clear linearity with the fracture density. The proposed k-V relationship as a function of the aperture and fracture density indicates that laboratory-scale fracture properties of a single fracture can be applied to multiple fractures on a larger scale. Our findings can be used to interpret temporal changes in seismic observations and to monitor fluid flow in fractures.