Similarity measurement on human mobility data with spatially weighted structural similarity index (SpSSIM)

Understanding diverse characteristics of human mobility provides profound knowledge of urban dynamics and complexity. Human movements are recorded in a variety of data sources and each describes unique mobility characteristics. Revealing similarity and difference in mobility data sources facilitates grasping comprehensive human mobility patterns. This study introduces a new method to measure similarities on two origin–destination (OD) matrices by spatially extending an image‐assessment tool, the structural similarity index (SSIM). The new measurement, spatially weighted SSIM (SpSSIM), utilizes weight matrices to overcome the SSIM sensitivity issue due to the ordering of OD pairs by explicitly defining spatial adjacency. To evaluate SpSSIM, we compared performances between SSIM and SpSSIM with resampling the orders of OD pairs and conducted bootstrapping to test the statistical significance of SpSSIM. As a case study, we compared OD matrices generated from three data sources in San Diego County, CA: U.S. Census‐based Longitudinal Employer–Household Dynamics Origin–Destination employment statistics, Twitter, and Instagram. The case study demonstrated that SpSSIM was able to capture similarities of mobility patterns between datasets that varied by distance. Some regions showed local dissimilarity while the global index indicated they were similar. The results enhance the understanding of complex mobility patterns from various datasets, including social media.     

Intra-oceanic island arc origin for Iratsu eclogites of the Sanbagawa belt,central Shikoku,southwest Japan

New geochemical and Sr–Nd isotopic data for the Iratsu eclogite and surrounding metamorphic rocks of the Sanbagawa belt, Japan, show that, while the protoliths of the metamorphic rocks formed in a variety of tectonic settings, the Iratsu body represents a deeply subducted and accreted island arc. The igneous protoliths of eclogites and garnet amphibolites were probably generated from a mantle source that had components of both a depleted mantle modified by slab-released fluid (as seen in a negative Nb anomaly) and an enriched mantle, similar to that of ocean island basalts (OIB). Fractional crystallization modeling indicates that the protoliths of some garnet clinopyroxenites from the Iratsu body are cumulates from a basaltic magma that crystallized under high O₂ and H₂O fugacities in the middle to lower crust. The source characteristics and crystallization conditions suggest that the protoliths of the Iratsu rocks formed in an oceanic island arc. Quartz eclogites from the marginal zone of the Iratsu body have geochemical signatures similar to turbidites from the Izu–Bonin island arc (as seen in a negative Nb anomaly and a concave REE pattern). The protoliths might be volcaniclastic turbidites that formed in a setting proximal to the oceanic island arc. Geochemical and isotopic signatures of the surrounding mafic schists are similar to normal (N-) and enriched (E-) mid-ocean-ridge basalt (MORB), and distinct from the rocks from the Iratsu body. The protoliths of the mafic schists likely formed in a plume-influenced mid-ocean ridge or back-arc basin. Pelitic schists from the surrounding rocks and pelitic gneisses from the marginal zone of the Iratsu body have evolved, continental geochemical signatures (as seen in a negative ε_Nd(t) value (~?5)), consistent with their origin as continent-derived trench-fill turbidites.     

Estimation of standing stock of chlorophylla and primary production from remote-sensed ocean color in the Oyashio region, the Western subarctic Pacific, during the spring bloom in 1997

In order to examine the applicability of remotely-sensed ocean color for the estimation of phytoplankton biomass and primary production in the Oyashio region, the western subarctic Pacific, vertical distributions of chlorophylla concentration and primary production were observed in April and May 1997. Spring bloom was observed in both April and May, and the surface concentration of chlorophylla exceeded 40 mg m⁻³. The relationship between the standing stocks of chlorophylla within the layer from the sea surface to one optical depth (0–1/k layer) and the surface chlorophylla concentration is expressed as a Michaelis-Menten equation. The mean ratio of the standing stock of chlorophylla in the euphotic layer to that in the 0–1/k layer was 4.41, this ratio did not significantly differ from 4.61 which was obtained at homogeneous distribution of chlorophylla within the euphotic layer. These facts suggest that the distribution of chlorophylla could be assumed to be homogeneous in the euphotic layer during the spring bloom. Results of primary production measurements by simulatedin situ method were compared with those by an algorithm with two variables; chlorphylla and non-spectral PAR. Daily primary production in the euphotic layer estimated by the algorithm varied in a range of 38–274% of that estimated by incubation, although the primary productions by the algorithm agreed with those by the incubation at a half of stations. Primary production within the euphotic layer calculated using simply the surface data was the same as that estimated using vertical distribution of chlorophylla. These results show that the primary production in the euphotic layer may be estimated from the remote sensed measurements during the spring bloom in the Oyashio region.     

Changes in the abundance and composition of picophytoplankton in relation to the occurrence of a Kyucho and a bottom intrusion in the Bungo Channel,Japan

The Bungo Channel in southwestern Japan receives both warm, called Kyucho, and cold deep-water intrusions (bottom intrusion) from the Pacific Ocean. Abundances of Prochlorococcus, Synechococcus, and eukaryotic picophytoplankton were monitored from 18 July to 17 August 2001 to clarify whether advected picophytoplankton from the Pacific Ocean can grow in the channel or not. Synechococcus cells were further discriminated into low- and high-PUB types according to their fluorescence property in flow cytometry. From 18 to 25 July, the water temperature decreased by 3 °C at a 5-m depth at all stations, indicating the occurrence of a bottom intrusion. From 25 July to 4 August, a Kyucho occurred and the water temperature rapidly increased. From 4 to 17 August, a bottom intrusion and a Kyucho both occurred twice, although the intensities were smaller than those occurring until 4 August. From 18 to 30 July, the abundance of both Prochlorococcus and a high-PUB type of Synechococcus drastically decreased because of a bottom intrusion; however, the abundances rapidly increased due to the advection by a Kyucho. These advected cells increased from 4 to 17 August in the channel and Kitanada Bay. Changes in the abundance of low-PUB type of Synechococcus and eukaryotic picophytoplankton were less noticeable than those in the abundance of Prochlorococcus and high-PUB type. The present study demonstrated that oceanic picophytoplankton advected by the Kyucho could grow in the channel. However, abundances of low-PUB type and eukaryotic picophytoplankton increased higher than those of Prochlorococcus and high-PUB type did. Thus, these oceanic phytoplankters will be excluded when Kyucho does not occur for a long time. The co-occurrence of various types of picophytoplankton found in the channel is probably achieved by both Kyucho event and their growth capability in the channel.     

Periodic Intrusion of Cold Water from the Pacific Ocean into the Bottom Layer of the Bungo Channel in Japan

A cold-water intrusion, called a “bottom intrusion”, occurs in the lower layer of the Bungo Channel in Japan. It is an intrusion from the shelf slope region of the Pacific Ocean margin in the south of the channel. In order to reveal the fundamental characteristics of the bottom intrusion, we conducted long-term observations of water temperature at the surface and bottom layers of the channel and 15-day current observations at the bottom of the shelf-break region. The long-term water temperature data indicated that the bottom intrusion occurs repeatedly between early summer and late autumn, and its reiteration between early and mid-summer causes a local minimum of water temperature in the lower layer in mid-summer. Moreover, the data revealed that most of the bottom intrusions occurred in neap tidal periods. The current meter recorded a bottom intrusion with a speed of approximately 15 cm⋅s⁻¹. The current meter also revealed that the intruded cold water slowly retreated back to the shelf slope region after the intrusion. This revised version was published online in July 2006 with corrections to the Cover Date.     

Genetic algorithm for fitting a mixed Bingham distribution to 3D orientations: a tool for the statistical and paleostress analyses of fracture orientations

The clustering of fracture orientations is important for tectonic studies and for geotechnical engineering. In this study, a real‐coded genetic algorithm was adopted to fitting a mixed Bingham distribution to orientation data by maximizing the log‐likelihood function of the distribution. The maximization is a difficult problem, because the function has multimodality and singularity. It was found that the algorithm was effective for this problem. Given the orientations of dilational fractures, the present method determines not only the stress axes and stress ratio of each of the fracture groups but also the maximum non‐dimensionalized fluid pressure at the time of their formation. In addition, the software calculates the 95 % error ellipses of the concentration axes. The present method found that the orientations of ore veins of the Akenobe Mine, SW Japan, should be partitioned into three clusters. It is shown that two of the groups had distinctive Zn and Sn contents, and that the ore fluids had overpressures only slightly greater than the minimum principal stress at the time of the deposition of Zn‐ and Sn‐rich veins.     

Solar Wind Speed and Expansion Rate of the Coronal Magnetic Field in Solar Maximum and Minimum Phases

Relationships between solar wind speed and expansion rate of the coronal magnetic field have been studied mainly by in-ecliptic observations of artificial satellites and some off-ecliptic data by Ulysses. In this paper, we use the solar wind speed estimated by interplanetary scintillation (IPS) observations in the whole heliosphere. Two synoptic maps of SWS estimated by IPS observations are constructed for two Carrington rotations CR 1830 and 1901; CR 1830 starting on the 11th of June, 1990 is in the maximum phase of solar activity cycle and CR 1901 starting on the 29th of September, 1995 is in the minimum phase. Each of the maps consist of 64800 (360×180) data points. Similar synoptic maps of expansion rate of the coronal magnetic field (RBR) calculated by the so-called potential model are also constructed under a radial field assumption for CR 1830 and CR1901. Highly significant correlation (r=–0.66) is found between the SWS and the RBR during CR1901 in the solar minimum phase; that is, high-speed winds emanate from photospheric areas corresponding to low expansion rate of the coronal magnetic field and low speed winds emanate from photospheric areas of high expansion rate. A similar result is found during CR 1830 in solar maximum phase, though the correlation is relatively low (r=–0.29). The correlation is improved when both the data during CR 1830 and CR 1901 are used together; the correlation coefficient becomes –0.67 in this case. These results suggest that the correlation analysis between the SWS and the RBR can be applied to estimate the solar wind speed from the expansion rate of the coronal magnetic field, though the correlation between them may depend on the solar activity cycle. We need further study of correlation analysis for the entire solar cycle to get an accurate empirical equation for the estimation of solar wind speed. If the solar wind speed is estimated successfully by an empirical equation, it can be used as an initial condition of a solar wind model for space weather forecasts.