Hyporheic exchange flows induced by constructed riffles and steps in lowland streams in southern Ontario,Canada

Stream–subsurface water interaction induced by natural riffles and constructed riffles/steps was examined in lowland streams in southern Ontario, Canada. The penetration of stream water into the subsurface was analysed using hydrometric data, and the zone of > 10% stream water was calculated from a chemical mixing equation using tracer injection of bromide and background chloride concentrations. The constructed riffles studied induced more extensive hyporheic exchange than the natural riffles because of their steeper longitudinal hydraulic head gradients and coarser streambed sediments. The depth of > 10% stream water zone in a small and a large constructed riffle extended to > 0·2 m and > 1·4 m depths respectively. Flux and residence time distribution of hyporheic exchange were simulated in constructed riffles using MODFLOW, a finite‐difference groundwater flow model. Hyporheic flux and residence time distribution varied along the riffles, and the exchange occurring upstream from the riffle crest was small in flux and had a long residence time. In contrast, hyporheic exchange occurring downstream from the riffle crest had a relatively short residence time and accounted for 83% and 70% of total hyporheic exchange flow in a small and large riffle respectively. Although stream restoration projects have not considered the hyporheic zone, our data indicate that constructed riffles and steps can promote vertical hydrologic exchange and increase the groundwater–surface water linkage in degraded lowland streams. Copyright © 2006 John Wiley & Sons, Ltd.     

Seismological monitoring on the 2003 Tokachi-oki earthquake, derived from off Kushiro permanent cabled OBSs and land-based observations

Japan Marine Science and Technology Center installed a cabled geophysical observatory system off Kushiro, Hokkaido Island in July 1999. This observatory system comprises three ocean bottom seismographs (OBSs), two tsunami gauges, and a geophysical/geochemical monitoring system. 4 years and 2 months after the installation, a megathrust earthquake (the 2003 Tokachi-Oki earthquake, 26th September in Japan Standard Time (JST), M_JMA 8.0) occurred along a plate boundary underneath a forearc basin where the system is located. The system recorded clear unsaturated seismograms just at 28.6 km from the epicenter. This paper demonstrates advantages brought by the cabled observatory to record the megathrust earthquake showing how earthquake detectability is improved dramatically combining permanent OBS and land-based observations around the region, and importance of the in situ monitoring on the seismogenic zone. In the present study, processing OBSs and land-based network together, and comparing magnitudes of common observed earthquakes with national authorized network, event detection level improved down to M 1.5, which is much lower than the previously designed as down to M  2. Comparing detection level before and after installing OBSs, we found dramatic improvement of the earthquake detection level in the interesting region. Real-time continuous observations of microearthquakes since 1999 have brought us tremendous findings. First, a seismic quiescence started about 10 days before the 2003 Tokachi-Oki earthquake. Second, aftershock distribution is not uniform over the focal area and can be divided into several sub-regions, which might indicate an existence of several asperities. We think that the geophysical observations helped to understand the initiation process of the rupture of the 2003 Tokachi-Oki earthquake and that observations including seismological, geodynamic, hydrogeological, and the other multidisciplinary observations would provide a clue to future understanding of seismogenic processes at subduction zones.     

Parameters of split transverse waves from local earthquakes along eastern Hokkaido recorded before the Tokachi-Oki September 26, 2003, earthquake (M = 8.0)

The parameters of split S waves from local weak earthquakes along eastern Hokkaido Island are studied over the period of 2003, including the strong Tokachi-oki September 26, 2003 earthquake (M = 8.0). Earthquake records of five stations belonging to the ISV seismological network were used. The studies of the split S wave parameters showed that they vary in space and time along Hokkaido Island. The zones of the Hidaka Mountains (ERM, MYR), Tokachi Plain (IWN, URH), and Kushiro Plain (AKK) are distinguished along Hokkaido. The anisotropy coefficients beneath the ERM, MYR, IWN, URH, and AKK stations attain 10.5, 10, 5, 3.5, and 6.5%, respectively. Beneath ERM, azimuths of the fast S wave (?) are predominantly in the N-S direction until July and in the E-W direction from July (parallel and normal to the Japan trench strike). By the time of the Tokachi-oki earthquake, the ? directions were oriented SE in agreement with the direction of the Pacific plate motion. The ? directions on the northern side of the Hidaka Range (MYR) are predominantly orthogonal to those beneath ERM, which can be evidence for differences in the direction of deformations on opposite sides of the range. Higher seismicity, the variation of S wave parameters, and a high anisotropy of the medium point to an intense development of deformation (dilatancy) processes in the area of the Hidaka Mountains. The fast wave azimuths beneath AKK are predominantly 50°–70°, and this orientation is consistent with the direction of migration of the Kurile arc front along the trench. Beneath IWN, the azimuths ? are oriented along the N-NE directions, and beneath URH, along the direction of the Pacific plate motion (100°–150°). Strengthening of mechanical properties of the medium and development and accumulation of shear deformations in a subhorizontal plane are supposed to take place in the Tokachi Plain area.     

Robust and exploratory analysis of active mesoscale tectonic zones in Japan utilizing the nationwide GPS array

A monitoring GPS array recently developed in Japan can yield nationwide maps of active inland tectonic zones (ATZs) on a mesoscale, approximately 70 to several hundred kilometers in lateral extent. But it has been difficult to characterize ATZs in Japan, as they are in fact operational on multiple scales and our efforts are often hindered by various irregularities in the data. The key to overcoming these problems would be to gain an insight into the available data before any precise kinematic modeling is performed with indefinite assumptions. In this study, horizontal velocity fields, deduced from the nationwide GPS array, were treated with a set of techniques in robust smoothing and exploratory data analysis that brought out exceptionally powerful mesoscale ATZs, and made them easier to characterize. The resolved ATZs were then retrospectively monitored to study their regional and temporal variations, using a set of approx. 840 observation stations, about 30 km apart, for a 4-year series of fixed observation time-intervals, 810 days each. The smoothing operation involved three steps: (1) imputation of the velocity fields for the purpose of anti-aliasing, (2) robust smoothing of the velocity fields with the median operative, and (3) visualization of deformation-rate distributions in several coordinate independent parameters, and post-filtering. The geometrical resolvability of mesoscale ATZs was confirmed by calibrating the smoothing scheme against synthetic tectonic boundary models before it was applied to the case study in Japan. ATZs in Japan, which are essentially visible as systematic deviations in the velocity fields on the International Terrestrial Reference Frame (ITRF) and as strain rate anomalies, were highlighted sharply along some known tectonic zones, chains of active volcanoes, and areas above low seismic velocity anomalies in the crust and upper mantle, all of which generally paralleled the offshore trench axes. The geometrical agreements among the mapped ATZs and the physical anomalies in the crust are presumably due to their common structural weakness on the mesoscale. In the four main islands of Japan, all but 30–40% of the strain rate anomalies persisted during the entire 6 years of the case study period, while the rest sporadically appeared or disappeared in a period from several months to a few years. The transient shifts in the deformation rates were remarkably synchronous with some nearby major tectonic episodes: large earthquakes and slow events. Differential plate coupling strengths along the subduction zones can also be inferred from the persistent pattern of rotational strain rate anomalies forming clockwise and counterclockwise pairs along the Pacific. Our empirical observations suggest that the first-order features of interseismic crustal deformations in Japan can be characterized as collateral processes behaving in response to fluctuations of the tectonic stresses on multiple scales, likely influenced by changes of plate coupling strengths on the contiguous subduction faults.     

An impulsive solar X-ray burst in the energy range 90–400 keV

A hard solar X-ray burst in the energy range 90–400 keV was observed at about 0528:30 UT on 13 April 1974. This event is peculiar in the sense that the X-ray flux is considerably high whereas the corresponding radio and optical flares are of moderate sizes, that the duration of the X-ray burst is very short, and that the X-ray peak is delayed by about 2.5 min relative to the first radio peak.     

Normal faulting of the Daiichi-Kashima Seamount in the Japan Trench revealed by the Kaiko I cruise, Leg 3

A detailed topographic and geophysical survey of the Daiichi-Kashima Seamount area in the southern Japan Trench, northwestern Pacific margin, clearly defines a high-angle normal fault which splits the seamount into two halves. A fan-shaped zone was investigated along 2–4 km spaced, 100 km long subparallel tracks using narrow multi-beam (Seabeam) echo-sounder with simultaneous measurements of gravity, magnetic total field and single-channel seismic reflection records. Vertical displacement of the inboard half was clearly mapped and its normal fault origin was supported. The northern and southern extensions of the normal fault beyond the flank of the seamount were delineated. Materials on the landward trench slope are displaced upward and to sideways away from the colliding seamount. Canyons observed in the upper landward slope terminate at the mid-slope terrace which has been uplifted since start of subduction of the seamount. Most of the landward slope except for the landward walls aside the seamount comprises only a landslide topography in a manner similar to the northern Japan Trench wall. This survey was conducted on R/V “Jean Charcot” as a part of the Kaiko I cruise, Leg 3, in July–August 1984 under the auspices of the French-Japanese scientific cooperative program.     

Paleotectonic structures and their influence on recent seismo-tectonics in the south Kuril subduction zone

Abstract Bathymetric data from south of Hokkaido obtained during a cruise of R/V Hakuho-Maru are summarized, and their correlation with earthquake occurrence is discussed. There are structural lineations on the seaward slope of the Kuril Trench, oblique to the Kuril Trench axis and parallel to the magnetic lineations in the Pacific plate. The structural lineations comprise horst-grabens generated by normal faulting. This suggests that Cretaceous tectonic structures originating at the spreading centre affect present seismotectonics around the trench axis. The structural-magnetic relation is compared to the case of the Japan Trench. North-east of the surveyed area, there are two major fracture zones (Nosappu Fracture Zone and Iturup Fracture Zone) that divide the oceanic plate into three segments. If the fracture zones (FZ) and the zone of paleo-mechanical weakness, represented by magnetic lineations, can control the direction of normal faults at a trench, the extent of the resulting topographic roughness on the seaward slope of the trench would be different across an FZ because of the differences in ages. By studying recent large earthquakes occurring in the south Kuril region, it is shown that several main-aftershock distributions for large earthquakes in this region are bounded by the Nosappu FZ and the Iturup FZ. Two models (Barrier model and Rebound model) are presented to interpret earthquake occurrence near the south Kuril Islands. The Barrier model explains seismic boundaries seen in several examples for earthquake occurrence in the south Kuril regions. The fracture zone forming the boundary of two segments with different magnetic lineations is also the boundary of two different normal fault systems on their ocean bottom, and the difference in sea-bottom roughness between two normal fault systems should affect the seismic coupling at a plate interface. Due to the difference of seismic coupling, earthquake occurrence is controlled by an FZ and then the FZ acts as a seismic boundary (Barrier model). Existing normal faults created by plate bending of subducting oceanic plate should rebound after its subduction (Rebound model). This rebound of normal faults may cause intraplate earthquakes with a high-angle reverse-fault mechanism such as the 1994 Shikotan Earthquake. The energy released by an intraplate earthquake generated by normal-fault rebounding is not directly related to that of interplate earthquakes such as low-angle thrust earthquakes. It is a reason why large earthquakes occurred in the same region during a relatively short period.     

Three-dimensional P- and S-wave velocity structures beneath the Japan Islands obtained by high-density seismic stations by seismic tomography

We construct fine-scale 3D P- and S-wave velocity structures of the crust and upper mantle beneath the whole Japan Islands with a unified resolution, where the Pacific (PAC) and Philippine Sea (PHS) plates subduct beneath the Eurasian (EUR) plate. We can detect the low-velocity (low-V) oceanic crust of the PAC and PHS plates at their uppermost part beneath almost all the Japan Islands. The depth limit of the imaged oceanic crust varies with the regions. High-V_P/V_S zones are widely distributed in the lower crust especially beneath the volcanic front, and the high strain rate zones are located at the edge of the extremely high-V_P/V_S zone; however, V_P/V_S at the top of the mantle wedge is not so high. Beneath northern Japan, we can image the high-V subducting PAC plate using the tomographic method without any assumption of velocity discontinuities. We also imaged the heterogeneous structure in the PAC plate, such as the low-V zone considered as the old seamount or the highly seismic zone within the double seismic zone where the seismic fault ruptured by the earthquake connects the upper and lower layer of the double seismic zone. Beneath central Japan, thrust-type small repeating earthquakes occur at the boundary between the EUR and PHS plates and are located at the upper part of the low-V layer that is considered to be the oceanic crust of the PHS plate. In addition to the low-V oceanic crust, the subducting high-V PAC plate is clearly imaged to depths of approximately 250 km and the subducting high-V PHS zone to depths of approximately 180 km is considered to be the PHS plate. Beneath southwestern Japan, the iso-depth lines of the Moho discontinuity in the PHS plate derived by the receiver function method divide the upper low-V layer and lower high-V layer of our model at depths of 30–50 km. Beneath Kyushu, the steeply subducting PHS plate is clearly imaged to depths of approximately 250 km with high velocities. The high-V_P/V_S zone is considered as the lower crust of the EUR plate or the oceanic crust of the PHS plate at depths of 25–35 km and the partially serpentinized mantle wedge of the EUR plate at depths of 30–45 km beneath southwestern Japan. The deep low-frequency nonvolcanic tremors occur at all parts of the high-V_P/V_S zone—within the zone, the seaward side, and the landward side where the PHS plate encounters the mantle wedge of the EUR plate. We prove that we can objectively obtain the fine-scale 3D structure with simple constraints such as only 1D initial velocity model with no velocity discontinuity.     

On the correlation between ionospheric perturbations as detected by subionospheric VLF/LF signals and earthquakes as characterized by seismic intensity

The long-term data during seven years from January 2001 to December 2007, as observed by the Pacific VLF/LF network consisting of several Japanese stations and one station in Kamchatka, are extensively utilized to perform a statistical correlation study between the lower ionospheric perturbations as detected by subionospheric propagation and earthquakes (EQs). In this paper, we adopt a physical parameter, the maximum seismicity intensity observed (I) to define the strength of an EQ unlike the previously and conventionally used EQ magnitude and depth, which is a combined effect of EQ magnitude and depth, together with the Earth’s surface information and geological condition around the EQ epicenter. After considering EQs only take place on the land because of the use of seismicity and by using the superimposed epoch analysis, it is found for the larger EQs with I from 5 to 7 (we feel serious trembling and we expect serious damage) that the most important VLF/LF parameter, trend (nighttime average amplitude), shows a definite decrease about 10 days before the EQ by exceeding 2σ (σ: standard deviation) criterion; the dispersion shows a maximum about 10 days before the EQ but not exceeding 2σ line and finally the nighttime fluctuation shows an enhancement about 10 days before the EQ (with exceeding the 2σ level). A definite statistical correlation is confirmed between the ionospheric perturbations and I when I is strong enough in a range from 5 to 7. Whereas, there is no significant correlation between the two when I is in a range from 3 to 4. Finally, together with the corresponding results for EQs in the sea, but close to the land, these results are discussed in the light of lithosphere–ionosphere coupling mechanism.     

Spatio-temporal streamflow generation in a small,steep headwater catchment in western Japan

Headwaters contribute a substantial part of the flow in river networks. However, spatial variations of streamflow generation processes in steep headwaters have not been well studied. In this study, we examined the spatio-temporal variation of streamflow generation processes in a steep 2.98-ha headwater catchment. The time when baseflow of the upstream section exceeded that downstream was coincident with the time when the riparian groundwater switched from downwelling to upwelling. This suggests that upwelling of the riparian groundwater increased considerably in the upstream section during the wet period, producing a shift in the relative size of baseflow between the upstream and downstream sections. The timing of fluctuations among hillslope soil moisture, hillslope groundwater and streamflow reveals that the hillslope contributed to storm flow, but this contribution was limited to the wet period. Overall, these results suggest that streamflow generation has strong spatial variations, even in small, steep headwater catchments.

EDITOR A. Castellarin ASSOCIATE EDITOR X. Chen