Swath Mapping Research Articles

Kaiwi shoreline basalts fed by the west rift zone of East Molokaʻi

Bathymetry and acoustic imagery swath mapping, along with observations and samples from four manned submersible and four ROV dives, confirm that a seafloor slope break on the northern approaches to Kaiwi Channel, between the islands of Oʻahu and Molokaʻi, Hawaiʻi is a former shoreline, now submerged ∼800 m below present sea level. Subaerially emplaced, low-relief basaltic lavas above the slope break transition to submarine morphologies below. The entire region has been tilted about 1° to the SSE (150°), and is cut by an 8–15 m-high, north-facing scarp, 100–400 m south of the slope break. The distribution of platy, table-top, and rarer mounded branching corals indicates the former presence of fringing reefs around low-relief paleo-islands. We infer that the regional tilt resulted from loading by younger Hawaiian volcanoes, compounded by flexural uplift and back tilting away from the unloaded footwall of a flank landslide to the north.Basalt samples collected from both above and below the slope break have petrography, chemical composition, and age (1.64–1.80 Ma) indicating correlation with the (late-shield) Lower Member of the East Molokaʻi Volcanics, rather than with the more proximal volcano of West Molokaʻi. The most likely source of the Kaiwi basalts is a submarine ridge (rift zone) that extends northwest away from ʻĪlio Point on West Molokaʻi. Although the submarine ridge was previously assumed to be an extension of West Molokaʻi's northwest rift, we conclude that regional bathymetry and gravity are consistent with this feature being an extension of the west rift of East Molokaʻi. A corallary of this interpretation is that the shoreline slope break (SSB 7 of Taylor, 2019) in this area is distinct from and younger than the southern SSB 7 formed on West Molokaʻi volcano (∼1.65 Ma vs. ∼1.8 Ma).

Marine geological investigation of Edward VIII Gulf, Kemp Coast, East Antarctica

AbstractA physical oceanographic, geophysical and marine geological survey of Edward VIII Gulf, Kemp Coast, collected data from conductivity–temperature–depth casts, multi-beam bathymetric swath mapping and 3.5 kHz sub-bottom surveying. Modified circumpolar deep water (mCDW) is observed in Edward VIII Gulf, as well as notable bathymetric features including mega-scale glacial lineations and a 1750 m-deep trough. Sedimentological, geochemical, rock-magnetic and micropalaeontological analysis of two kasten cores document regional palaeoclimate and palaeo-oceanographic conditions over the past 8000 years, with a warm period occurring fromc.8 to 4 ka and a shift to cooler conditions beginning atc.4 ka and persisting until at least 0.9 ka. Sediment packages > 40 m thick within deep troughs in Edward VIII Gulf present potential targets for higher-resolution Holocene and deglacial climate studies. Despite the presence of mCDW on the shelf, inland bed topography consisting of highland terrain suggests the likelihood of relative stability of this sector of the East Antarctic Ice Sheet.

Bare Earth DEM Generation for Large Floodplains Using Image Classification in High-Resolution Single-Pass InSAR

The increasing frequency and severity of flood events demand improved accuracy of hydrodynamic models to better mitigate the societal and economic consequences of these disasters. Ideally, these models derive elevation data from high-accuracy LiDAR, which often captures the subtle variations in elevation and microtopography—elements that are critical to high-resolution hydrodynamic modeling. Due largely to the cost of acquisitions, these data, however, are not widely available at high spatial resolutions for large-scale areas, i.e. floodplains and coastal regions, especially outside of the U.S. and Western Europe, where flood-prone communities already lack the infrastructure and resources to manage these disasters. High-resolution interferometric synthetic aperture radar (InSAR) systems may offer a viable and cost-effective alternative to LiDAR, with comparable spatial resolutions and vertical accuracies. Like any swath mapping sensor, systematic errors in calibration knowledge increase as the lever-arm increases when viewing further off vertical. InSAR-derived digital elevation models (DEMs) are often corrected using LiDAR or ground control points, although this limits the application of InSAR to those areas where these data exist. In this paper, we present a novel approach for using near-range InSAR data to correct inherent systematic bias that propagates across adjacent, overlapping swaths for the same airborne acquisition. This eliminates the need for LiDAR in generating InSAR DEMs, while maintaining a vertical error of approximately 0.26 m relative to LiDAR, at a spatial resolution of 30 m. Data was acquired using the NASA/JPL airborne, single-pass, Ka-band GLISTIN-A (Glacier and Ice Surface Topography Interferometer) InSAR over the Red River Basin (RRB), ND. The accuracy of the final DEM is evaluated using National Geodetic Survey (NGS) GPS and a high-resolution LiDAR DEM.

Submarine platform development by erosion of a Surtseyan cone at Capelinhos, Faial Island, Azores

AbstractErosion of volcanic islands ultimately creates shallow banks and guyots, but the ways in which erosion proceeds to create them over time and how the coastline retreat rate relates to wave conditions, rock mass strength and other factors are unclear. The Capelinhos volcano was formed in 1957/58 during a Surtseyan and partly effusive eruption that added an ~2.5 km2 tephra and lava promontory to the western end of Faial Island (Azores, central North Atlantic). Subsequent coastal and submarine erosion has reduced the subaerial area of the promontory and created a submarine platform. This study uses historical information, photos and marine geophysical data collected around the promontory to characterize how the submarine platform developed following the eruption. Historical coastline positions are supplemented with coastlines interpreted from 2004 and 2014 Google Earth images in order to work out the progression of coastline retreat rate and retreat distance for lava‐ and tephra‐dominated cliffs. Data from swath mapping sonars are used to characterize the submarine geometry of the resulting platform (position of the platform edge, gradient and morphology of the platform surface). Photographs collected during SCUBA and ROV dives on the submarine platform reveal a rugged surface now covered with boulders. The results show that coastal retreat rates decreased rapidly with time after the eruption and approximately follow an inverse power‐law relationship with coastal retreat distance. We develop a finite‐difference model for wave attenuation over dipping surfaces to predict how increasing wave attenuation contributed to this trend. The model is verified by reproducing the wave height variation over dipping rock platforms in the UK (platform gradient 1.2° to 1.8°) and Ireland (1.8°). Applying the model to the dipping platform around Capelinhos, using a diversity of cliff resistance predicted from known lithologies, we are able to predict erosion rate trends for some sectors of the edifice. We also explore wider implications of these results, such as how erosion creates shallow banks and guyots in reef‐less mid‐oceanic archipelagos like the Azores. © 2019 John Wiley & Sons, Ltd. © 2019 John Wiley & Sons, Ltd.

Validation of Glacier Topographic Acquisitions from an Airborne Single-Pass Interferometer.

The airborne glacier and ice surface topography interferometer (GLISTIN-A) is a single-pass radar interferometer developed for accurate high-resolution swath mapping of dynamic ice surfaces. We present the first validation results of the operational sensor, collected in 2013 over glaciers in Alaska and followed by more exhaustive collections from Greenland in 2016 and 2017. In Alaska, overlapping flight-tracks were mosaicked to mitigate potential residual trends across-track and the resultant maps are validated with lidar. Furthermore, repeat acquisitions of Columbia Glacier collected with a three day separation indicate excellent stability and repeatability. Commencing 2016, GLISTIN-A has circumnavigated Greenland for 4 consecutive years. Due to flight hour limitations, overlapping swaths were not flown. In 2016, comparison with airborne lidar data finds that residual systematic errors exhibit evenly distributed small slopes (all less than 10 millidegrees) and nadir biases were typically less than 1 m. Similarly 2017 data exhibited up to meter-scale nadir biases and evenly distributed residual slopes with a standard deviation of ~10 millidegrees). All satisfied the science accuracy requirements of the Greenland campaigns (3 m accuracy across an 8 km swath).

Satellite gravity – enhancements from new satellites and new altimeter technology

ABSTRACTThis paper reviews the impacts of new satellite altimeter data sets and new technology on the production of satellite gravity. It considers the contribution of the increased data volume, the application of new altimeter acquisition technology and the potential for future developments. Satellite altimeter derived gravity has provided gravity maps of the world's seas since the 1980s, but, from 1995 to 2010, virtually all improvements were in the processing as there were no new satellite data with closely spaced tracks. In recent years, new data from CryoSat‐2 (launched in 2010) and the geodetic mission of Jason‐1 (2012–2013) have provided a wealth of additional coverage and new technology allows further improvements. The synthetic aperture radar mode of CryoSat‐2 uses a scanning approach to limit the size of the altimeter sea surface footprint in the along‐track direction. Tests indicate that this allows reliable data to be acquired closer to coastlines. The synthetic aperture radar interferometric mode of CryoSat‐2 uses two altimeters to locate sea‐surface reflection points laterally away from the satellite track. In a study to generate gravity for freshwater lakes, this mode is found to be valuable in extending the available satellite coverage. The AltiKa altimeter uses higher frequency radar to provide less noisy sea‐surface signals and its new orbit mode gives potential for further improvements in satellite gravity. Future developments include the potential for swath mapping to provide further gravity improvements.

Introduction to the Special Issue “Seafloor backscatter data from swath mapping echosounders: from technological development to novel applications”

Introduction to the Special Issue “Seafloor backscatter data from swath mapping echosounders: from technological development to novel applications”

A framework to quantify uncertainties of seafloor backscatter from swath mapping echosounders

Multibeam echosounders (MBES) have become a widely used acoustic remote sensing tool to map and study the seafloor, providing co-located bathymetry and seafloor backscatter. Although the uncertainty associated with MBES-derived bathymetric data has been studied extensively, the question of backscatter uncertainty has been addressed only minimally and hinders the quantitative use of MBES seafloor backscatter. This paper explores approaches to identifying uncertainty sources associated with MBES-derived backscatter measurements. The major sources of uncertainty are catalogued and the magnitudes of their relative contributions to the backscatter uncertainty budget are evaluated. These major uncertainty sources include seafloor insonified area (1–3 dB), absorption coefficient (up to > 6 dB), random fluctuations in echo level (5.5 dB for a Rayleigh distribution), and sonar calibration (device dependent). The magnitudes of these uncertainty sources vary based on how these effects are compensated for during data acquisition and processing. Various cases (no compensation, partial compensation and full compensation) for seafloor insonified area, transmission losses and random fluctuations were modeled to estimate their uncertainties in different scenarios. Uncertainty related to the seafloor insonified area can be reduced significantly by accounting for seafloor slope during backscatter processing while transmission losses can be constrained by collecting full water column absorption coefficient profiles (temperature and salinity profiles). To reduce random fluctuations to below 1 dB, at least 20 samples are recommended to be used while computing mean values. The estimation of uncertainty in backscatter measurements is constrained by the fact that not all instrumental components are characterized and documented sufficiently for commercially available MBES. Further involvement from manufacturers in providing this essential information is critically required.

Marsili and Cefalù basins: The evolution of a rift system in the southern Tyrrhenian Sea (Central Mediterranean)

Abstract The Marsili Basin (in the southern Tyrrhenian Sea), whose mode of extension is still a controversial issue, is the youngest bathyal basin of the Central Mediterranean. A thin sedimentary cover in the basin permits to image basement fabric and structure by swath mapping and seismic reflection data. We investigate the crustal structure of the southern Tyrrhenian Sea, extending from the bathyal Marsili Basin to the adjacent Sicily continental margin. Interpretation of seismic reflection profiles (calibrated by well and dredge data) and crustal cross-sections were used to identify the stratigraphic infill, structural pattern and large-scale crustal features of the region. We recognized three basins in the southern Tyrrhenian Sea: (1) the Termini basin that is an overfilled sedimentary basin on the continental shelf area; (2) the Cefalu basin, located on the continental slope area and filled by thick deep water turbidite deposits; (3) the distal Marsili basin, filled by hemipelagic and thin distal turbidite deposits. The sequence stratigraphy interpretation permitted us to recognize fourth-order depositional sequences and the stratigraphic signature of the rift stages. An important increase in the sedimentary supply from the continental shelf to the bathyal basin occurred approximately over the last 0.5 Ma and is related to the uplift of the coastal area. The stratigraphic constrains indicate a Lower Pleistocene age for the opening of the southern Tyrrhenian Sea basin. The structural map reveals a complex fault pattern, in which coeval normal faults disclose a triangular basin from the Marsili bathyal basin to the Sicily continental margin, associated to an Euler pole in the northern Sicily. Tacking in account the faults pattern that developed in the whole Southern Tyrrhenian Sea during the Lower Pleistocene, we reconstructed two opposite triangular basins separated by a perpendicular rift. For analogy with the contiguous Vavilov basin, we propose that the extension in the Marsili basin reached the mantle exhumation stage.

Remotely sensed habitat variables are poor surrogates for functional traits of rocky reef fish assemblages

SUMMARYConservation planning often relies on measures such as species richness and abundance to prioritize areas for protection. Nonetheless, alternative metrics such as functional traits have recently been shown to be useful complementary measures for detecting biological change. Timely conservation planning often precludes the collection of such detailed biological data relying instead on remotely-sensed habitat mapping as a surrogate for diversity. While there is evidence that habitat maps may predict taxonomic species richness and diversity in some coastal ecosystems, it is unknown whether similar strong relationships exist for functional traits and functional multimetrics. We compared the performance of physical habitat structural complexity obtained from high definition swath mapping in explaining variation in traditional taxonomic metrics as well as functional traits (e.g., maximum length, trophic level, gregariousness) and functional multimetrics (e.g., functional richness, dispersion) of fish assemblages. Reef complexity measures were good surrogates for fish species richness and abundance but not for functional traits or multimetrics, except functional richness at the scale of 1 m. Remotely sensed habitat maps may not be a good surrogate for predicting functional traits and multimetrics of fish assemblages, and must be used with caution when maximizing such aspects of assemblages is a priority for conservation planning.

The present-day number of tectonic plates

The number of tectonic plates on Earth described in the literature has expanded greatly since the start of the plate tectonic era, when only about a dozen plates were considered in global models of present-day plate motions. With new techniques of more accurate earthquake epicenter locations, modern ways of measuring ocean bathymetry using swath mapping, and the use of space based geodetic techniques, there has been a huge growth in the number of plates thought to exist. The study by Bird (2003) proposed 52 plates, many of which were delineated on the basis of earthquake locations. Because of the pattern of areas of these plates, he suggested that there should be more small plates than he could identify. In this paper, I gather together publications that have proposed a total of 107 new plates, giving 159 plates in all. The largest plate (Pacific) is about 20 % of the Earth’s area or 104 Mm2, and the smallest of which (Plate number 5 from Hammond et al. 2011) is only 273 km2 in area. Sorting the plates by size allows us to investigate how size varies as a function of order. There are several changes of slope in the plots of plate number organized by size against plate size order which are discussed. The sizes of the largest seven plates is constrained by the area of the Earth. A middle set of 73 plates down to an area of 97,563 km2 (the Danakil plate at number 80, is the plate of median size) follows a fairly regular pattern of plate size as a function of plate number. For smaller plates, there is a break in the slope of the plate size/plate number plot and the next 32 plates follow a pattern of plate size proposed by the models of Koehn et al. (2008) down to an area of 11,638 km2 (West Mojave plate # 112). Smaller plates do not follow any regular pattern of area as a function of plate number, probably because we have not sampled enough of these very small plates to reveal any clear pattern.

Origin of washboard moraines of the Des Moines Lobe: Spatial analyses with LiDAR data

The Des Moines Lobe (DML)—the largest lobe along the southern margin of the Laurentide ice sheet—deposited low-relief (1–5m) transverse moraine ridges, herein called washboard moraines, that are most prominent in central Iowa. Their origin as either subglacial or ice-marginal features bears on the dynamics of the lobe, geomorphic reconstructions of it, and interpretations of similar ridges elsewhere. Data from airborne laser swath mapping provide the first digital elevation models with sufficient spatial resolution to study in detail the geometric attributes of these topographically subtle features. Spectral analysis of profiles perpendicular to ridge crests indicates that most of them are spaced with statistically significant periodicities, with dominant wavelengths of 70–150m. Normalizing and stacking these profiles indicate that, on average, they display no systematic asymmetry. Locally, washboard moraines are intersected by discontinuous longitudinal ridge segments. Trends of the crests of moraines are scalloped, with cusps that point upglacier and are coincident with outwash trains, which were overridden by the lobe. Our interpretation of these data is that the moraines formed subglacially as crevasse-squeeze ridges. Overridden outwash trains likely supported low basal water pressures and resulted in anomalously slow basal slip, causing upglacier deflection of crevasses that extended to the bed. This crevasse-squeeze interpretation reinforces evidence that the DML surged to its maximum position. Bering Glacier, as a temperate, surge-type glacier that deposited crevasse-squeeze ridges similar to the washboard moraines of the DML, may be a good analog for the lobe.

Quantifying the effects of vegetation structure on snow accumulation and ablation in mixed‐conifer forests

AbstractSnowmelt from forested, mountainous environments in the western United States is a critical regional water resource for streamflow and ecological productivity. These landscapes are undergoing rapid changes from the combined effects of forest fires, insect infestation and climate change. Numerous observational studies demonstrate that trees control snowpack accumulation and ablation over scales of tens of metres. Representing forest heterogeneity in models is important for understanding how changes in climate and vegetation cover affect the snowpack; yet, many snow models simplify a forest into two categories: canopy covered and non‐canopy covered. We combine existing parameterizations of mass and energy fluxes within a new three‐dimensional framework informed by Airborne Laser Swath Mapping (ALSM)‐derived canopy maps and evaluated with ALSM‐derived snow depth maps to explicitly simulate snow cover in relation to heterogeneous canopy. Model results capture much of the observed snow variability depicted in the 1‐m ALSM‐derived snow depth maps. Observations and modelled results identify open areas <15 m from tree canopies as having more snow and more snow variability than areas >15 m from tree canopies, and modelled results predict that open areas <15 m from tree canopies have 30–40% more net snow water input than areas that are underneath tree canopies and 10–25% more net snow water input than areas that are >15 m from tree canopies. Furthermore, 1‐m simulations give higher estimates for net snow water input than coarser resolution simulations, mainly in areas with fewer trees. These results suggest the importance of explicitly representing canopy edges in snow models. Copyright © 2014 John Wiley & Sons, Ltd.

Deviations from self-similarity in barchan form and flux: The case of the Salton Sea dunes, California

[1] Barchans are the type of aeolian dune associated with a relatively uniform wind direction, incomplete sand coverage of the substrate, and low vegetation cover. Here I present an analysis of the morphology and migration rates of 40 dunes in the Salton Sea dune field using historical aerial orthophotographs, airborne laser swath mapping, terrestrial laser scanning, and measurements of the aerodynamic roughness length derived from wind velocity profiles. The data demonstrate that the Salton Sea dunes deviate from self-similarity such that smaller dunes have a lower ratio of slip face height to crest height and a lower slope, on average, compared with larger dunes and that smaller dunes migrate more slowly than would be predicted based on an inverse relationship between migration rate and dune height. The lack of self-similarity in barchans has been attributed to the dependence of speed-up ratios on dune size and the presence of a finite saturation length in the physics of aeolian transport. Here I argue that deviations from self-similarity at this study site are more likely due to the systematic decrease in aerodynamic roughness length with increasing elevation on stoss slopes. The data set I developed should prove useful to the aeolian geomorphic community for the further testing of models for barchan evolution.

Impact of flow routing on catchment area calculations, slope estimates, and numerical simulations of landscape development

[1] Flow routing across real or modeled topography determines the modeled discharge and wetness index and thus plays a central role in predicting surface lowering rate, runoff generation, likelihood of slope failure, and transition from hillslope to channel forming processes. In this contribution, we compare commonly used flow-routing rules as well as a new routing rule, to commonly used benchmarks. We also compare results for different routing rules using Airborne Laser Swath Mapping (ALSM) topography to explore the impact of different flow-routing schemes on inferring the generation of saturation overland flow and the transition between hillslope to channel forming processes, as well as on location of saturation overland flow. Finally, we examined the impact of flow-routing and slope-calculation rules on modeled topography produced by Geomorphic Transport Law (GTL)-based simulations. We found that different rules produce substantive differences in the structure of the modeled topography and flow patterns over ALSM data. Our results highlight the impact of flow-routing and slope-calculation rules on modeled topography, as well as on calculated geomorphic metrics across real landscapes. As such, studies that use a variety of routing rules to analyze and simulate topography are necessary to determine those aspects that most strongly depend on a chosen routing rule.

In situ calibration of seafloor acoustic backscatter data from swath mapping sonars

Geometric and radiometric corrections are necessary to convert raw seafloor acoustic backscatter data to imagery in which the intensity variations are due to the geoacoustical properties of the bottom and the angular dependence of the backscattering process. Radiometric corrections include removal of the combined effects of the sonar's transmit and receive beam patterns, and removal of the area ensonified. To this end, a method is presented to estimate the transmit/receive beam pattern of a multibeam swath mapping sonar using seafloor acoustic backscatter data collected over an entire survey area. This involves estimating the scattering area at each sounding based on local seafloor slopes. For a given sonar installation, results show that such beam pattern estimates remain stable to within ± 0.5 dB over different survey areas.

Active tectonics and the LiDAR revolution

A revolution in remote sensing, light detection and ranging (LiDAR) laser altimetry swath mapping, reveals the details of topographic features at such high resolution that they have transformed our understanding of tectonic forcing of the shape of the Earth’s surface. Meter-scale DEMs (digital elevation models) capture fault offsets, fault zone structure, off-fault deformation, and landscape properties at microgeomorphic scale, highlighting that the surface faithfully records the complexity and sensitivity of deformation in detail.

Recognition of surface flow processes influenced by roads and trails in mountain areas using high-resolution topography

Road networks in mountainous forest landscapes have the potential to increase the susceptibility to erosion and shallow landsliding. The same issue is observed also for minor trail networks, with evidences of surface erosion due to surface flow redistribution. This could be a problem in regions such as the Italian Alps where forestry and tourist activities are a relevant part of the local economy. This is just one among the several effects of modern anthropogenic forcing: it is now well accepted by the scientific community that we are living in a new era where human activities may leave a significant signature on the Earth, by altering its morphology, and significantly affecting the related surface processes. In this work, we proposed a methodology for the automatic recognition of roads and trails induced flow direction changes. The algorithm is based on the calculation of the drainage area variation in the presence, or in the absence of anthropic features such as roads and trails on hillslopes. To simulate the absence of alteration, the surface was smoothed considering moving windows of varying size. In the analysis, we used a 1 and 0.5 m Airborne Laser Swath Mapping technology (ALSM), using LiDAR (Light Detection And Ranging), and 0.2 m Terrestrial Laser Scanner (TLS) derived Digital Terrain Models (DTMs). The aim of the work is to underline the effectiveness of the proposed method based on high resolution topography in the detailed recognition of surface flow direction alteration due to roads, but also trail networks. We propose an automatic method to map at a large scale such alterations, also in areas where it is difficult to recognize them without a trail network surveyed in the field. This methodology could be considered as a support for modeling (i.e., terrain stability and erosion models), and it can be used to interactively assist the design of new infrastructure to reduce their effects on surface instabilities. The reported methodology could also have a role in risk management and environmental planning for mountain areas where tourism and the related economic activities are critical, and where also trails deserve attention due to induced slope instabilities.

Technical Note: Detection of gas bubble leakage via correlation of water column multibeam images

Abstract. Hydroacoustic detection of natural gas release from the seafloor has been conducted in the past by using singlebeam echosounders. In contrast, modern multibeam swath mapping systems allow much wider coverage, higher resolution, and offer 3-D spatial correlation. Up to the present, the extremely high data rate hampers water column backscatter investigations and more sophisticated visualization and processing techniques are needed. Here, we present water column backscatter data acquired with a 50 kHz prototype multibeam system over a period of 75 seconds. Display types are of swath-images as well as of a "re-sorted" singlebeam presentation. Thus, individual and/or groups of gas bubbles rising from the 24 m deep seafloor clearly emerge in the acoustic images, making it possible to estimate rise velocities. A sophisticated processing scheme is introduced to identify those rising gas bubbles in the hydroacoustic data. We apply a cross-correlation technique adapted from particle imaging velocimetry (PIV) to the acoustic backscatter images. Temporal and spatial drift patterns of the bubbles are assessed and are shown to match very well to measured and theoretical rise patterns. The application of this processing to our field data gives clear results with respect to unambiguous bubble detection and remote bubble rise velocimetry. The method can identify and exclude the main source of misinterpretations, i.e. fish-mediated echoes. Although image-based cross-correlation techniques are well known in the field of fluid mechanics for high resolution and non-inversive current flow field analysis, we present the first application of this technique as an acoustic bubble detector.

Geomorphological variations at hydrothermal sites in the southern Mariana Trough: Relationship between hydrothermal activity and topographic characteristics

This study presents the first detailed geomorphological characterization of field-scale geological features associated with hydrothermal systems in the southern Mariana Trough, using near-bottom swath mapping data collected by the autonomous underwater vehicle (AUV) Urashima during cruise YK09-08 and dive observation data acquired by the submersible Shinkai6500 during cruise YK10-11. The motivation of this study is to examine the relationship between geomorphological characteristics and hydrothermal activity, and to examine the nature of tectonic and volcanic controls on the hydrothermal system in this area. Two of the hydrothermal sites in the study area (near 12°57′N, 143°37′E) are located on the active backarc spreading axis (the Snail and Yamanaka sites), one is located at the eastern foot of the axial high (the Archean site), and two are located on an off-axis knoll about 5km from the spreading axis (the Pika and Urashima sites). The on-axis area is divided into tectonically dominant and volcanically dominant zones; volcanically dominant zones are characterized by mounds (height, 5–30m; diameter, 250–320m) cut by fissures. The Snail and Yamanaka sites are located adjacent to these fissures, and are possibly represented local activity associated with a 4th order segment-scale diking event (on the basis of comparisons with previously studied cases on the East Pacific Rise with similar on-axis geological characteristics). In contrast to the on-axis sites, the off-axis sites show no evidence of faulting. The Archean site at the foot of the axial high is characterized by a single mound (height, 50–100; diameter, 250–300m), pronounced off-axis lava flows, and the presence of high-amplitude rugged seafloor features; the site is located at the top of the mound. Numerous ridge lines (height, mainly 2–6m) extend radially from the top of the mound, and several chimney-like structures (up to approximately 6m high) occur on the top and slopes of the mound. The Pika site is located on the western peak of an off-axis knoll, and the newly discovered Urashima site is located at the northern foot of the western peak of the same knoll. The western peak is characterized by bumpy seabed textures formed by numerous smaller-scale mounds and ridge lines; however, the eastern peak has a very smooth top and slope, and shows no signs of hydrothermal activity. Numerous mounds (heights, 5–75m; diameters, 50–350m) are developed on the comparatively gentle slope of the knoll, in contrast to the numerous ridge lines (height, mainly 1–6m) developed on the relatively steep slopes of the knoll. On the basis of the associated geomorphological features, the three off-axis sites (Archean, Pika, and Urashima) were identified as localities created by relatively long-term large-scale hydrothermal activity, as compared with sites in the on-axis area. The sustained activity at off-axis sites appears closely related to an off-axis upwelling magma system. The three off-axis hydrothermal sites are composed mainly of breccia assemblages that probably originated from hydrothermal activity with black smoker venting. These areas are characterized by numerous ridge lines, conical mounds, and bumpy seabed texture, whereas the on-axis sites are characterized by the absence of ridge lines, and the presence of white smoker and shimmering observed on dome-shaped pillow mounds. Hence, the distribution of ridge lines, mound morphology, and bumpy seabed texture is likely to correlate with hydrothermal activity.