Rypina Irina I.

No Thumbnail Available
Last Name
Rypina
First Name
Irina I.
ORCID
0000-0002-2115-4522

Search Results

Now showing 1 - 20 of 38
  • Article
    Investigating subsurface pathways of Fukushima cesium in the Northwest Pacific
    (American Geophysical Union, 2019-06-18) Cedarholm, Ella R. ; Rypina, Irina I. ; Macdonald, Alison M. ; Yoshida, Sachiko
    Advective pathways for Fukushima Daiichi Nuclear Power Plant (FDNPP)‐derived cesium observed in 2013 at 166°E south of the Kuroshio Extension (KE) at >500 m on the 26.5σθ isopycnal are investigated. Attention is paid to the KE's role in shaping these pathways. Using a high‐resolution model, particle trajectories were simulated backward and forward in time on 26.5σθ between the 2013 observations and the 2011 source. A large fraction of backtracked trajectories interacted with the mixed layer just offshore of the FDNPP. The likeliest pathway reaching the deepest 2013 observed cesium location runs along the KE out to ~165°E, where it turns sharply southward. Forward trajectory statistics suggest that for 26.5σθ waters originating north of the KE, this current acted as a permeable barrier west of 155–160°E. The deepest 2011 model mixed layers suggest that FDNPP‐derived radionuclides may have been present at 30°N in 2013 at greater depths and densities (700 m; 26.8σθ).
  • Article
    Chaotic advection in a steady, three-dimensional, Ekman-driven eddy
    (Cambridge University Press, 2013-12-05) Pratt, Lawrence J. ; Rypina, Irina I. ; Ozgokmen, Tamay M. ; Wang, P. ; Childs, H. ; Bebieva, Y.
    We investigate and quantify stirring due to chaotic advection within a steady, three-dimensional, Ekman-driven, rotating cylinder flow. The flow field has vertical overturning and horizontal swirling motion, and is an idealization of motion observed in some ocean eddies. The flow is characterized by strong background rotation, and we explore variations in Ekman and Rossby numbers, E and Ro, over ranges appropriate for the ocean mesoscale and submesoscale. A high-resolution spectral element model is used in conjunction with linear analytical theory, weakly nonlinear resonance analysis and a kinematic model in order to map out the barriers, manifolds, resonance layers and other objects that provide a template for chaotic stirring. As expected, chaos arises when a radially symmetric background state is perturbed by a symmetry-breaking disturbance. In the background state, each trajectory lives on a torus and some of the latter survive the perturbation and act as barriers to chaotic transport, a result consistent with an extension of the KAM theorem for three-dimensional, volume-preserving flow. For shallow eddies, where E is O(1), the flow is dominated by thin resonant layers sandwiched between KAM-type barriers, and the stirring rate is weak. On the other hand, eddies with moderately small E experience thicker resonant layers, wider-spread chaos and much more rapid stirring. This trend reverses for sufficiently small E, corresponding to deep eddies, where the vertical rigidity imposed by strong rotation limits the stirring. The bulk stirring rate, estimated from a passive tracer release, confirms the non-monotonic variation in stirring rate with E. This result is shown to be consistent with linear Ekman layer theory in conjunction with a resonant width calculation and the Taylor–Proudman theorem. The theory is able to roughly predict the value of E at which stirring is maximum. For large disturbances, the stirring rate becomes monotonic over the range of Ekman numbers explored. We also explore variation in the eddy aspect ratio.
  • Article
    The influence of an eddy in the success rates and distributions of passively advected or actively swimming biological organisms crossing the continental slope
    (American Meteorological Society, 2020-06-29) Rypina, Irina I. ; Pratt, Lawrence J. ; Entner, Samuel ; Anderson, Amanda ; Cherian, Deepak A.
    The Lagrangian characteristics of the surface flow field arising when an idealized, anticyclonic, mesoscale, isolated deep-ocean eddy collides with continental slope and shelf topography are explored. In addition to fluid parcel trajectories, we consider the trajectories of biological organisms that are able to navigate and swim, and for which shallow water is a destination. Of particular interest is the movement of organisms initially located in the offshore eddy, the manner in which the eddy influences the ability of the organisms to reach the shelf break, and the spatial and temporal distributions of organisms that do so. For nonswimmers or very slow swimmers, the organisms arrive at the shelf break in distinct pulses, with different pulses occurring at different locations along the shelf break. This phenomenon is closely related to the episodic formation of trailing vortices that are formed after the eddy collides with the continental slope, turns, and travels parallel to the coast. Analysis based on finite-time Lyapunov exponents reveals initial locations of all successful trajectories reaching the shoreline, and provides maps of the transport pathways showing that much of the cross-shelf-break transport occurs in the lee of the eddy as it moves parallel to the shore. The same analysis shows that the onshore transport is interrupted after a trailing vortex detaches. As the swimming speeds are increased, the organisms are influenced less by the eddy and tend to show up en mass and in a single pulse.
  • Article
    Search and rescue at sea aided by hidden flow structures
    (Nature Communications, 2020-05-26) Serra, Mattia ; Sathe, Pratik ; Rypina, Irina I. ; Kirincich, Anthony R. ; Ross, Shane D. ; Lermusiaux, Pierre F. J. ; Allen, Arthur ; Peacock, Thomas ; Haller, George
    Every year, hundreds of people die at sea because of vessel and airplane accidents. A key challenge in reducing the number of these fatalities is to make Search and Rescue (SAR) algorithms more efficient. Here, we address this challenge by uncovering hidden TRansient Attracting Profiles (TRAPs) in ocean-surface velocity data. Computable from a single velocity-field snapshot, TRAPs act as short-term attractors for all floating objects. In three different ocean field experiments, we show that TRAPs computed from measured as well as modeled velocities attract deployed drifters and manikins emulating people fallen in the water. TRAPs, which remain hidden to prior flow diagnostics, thus provide critical information for hazard responses, such as SAR and oil spill containment, and hence have the potential to save lives and limit environmental disasters.
  • Article
    Observing and quantifying ocean flow properties using drifters with drogues at different depths
    (American Meteorological Society, 2021-07-16) Rypina, Irina I. ; Getscher, Timothy ; Pratt, Lawrence J. ; Mourre, Baptiste
    This paper presents analyses of drifters with drogues at different depths—1, 10, 30, and 50 m—that were deployed in the Mediterranean Sea to investigate frontal subduction and upwelling. Drifter trajectories were used to estimate divergence, vorticity, vertical velocity, and finite-size Lyapunov exponents (FTLEs) and to investigate the balance of terms in the vorticity equation. The divergence and vorticity are O(f) and change sign along trajectories. Vertical velocity is O(1 mm s−1), increases with depth, indicates predominant upwelling with isolated downwelling events, and sometimes changes sign between 1 and 50 m. Vortex stretching is one of the significant terms, but not the only one, in the vorticity balance. Two-dimensional FTLEs are 2 × 10−5 s−1 after 1 day, 2 times as large as in a 400-m-resolution numerical model. Three-dimensional FTLEs are 50% larger than 2D FTLEs and are dominated by the vertical shear of horizontal velocity. Bootstrapping suggests uncertainty levels of ~10% of the time-mean absolute values for divergence and vorticity. Analysis of simulated drifters in a model suggests that drifter-based estimates of divergence and vorticity are close to the Eulerian model estimates, except when drifters get aligned into long filaments. Drifter-based vertical velocity is close to the Eulerian model estimates at 1 m but differs at deeper depths. The errors in the vertical velocity are largely due to the lateral separation between drifters at different depths and are partially due to only measuring at four depths. Overall, this paper demonstrates how drifters, heretofore restricted to 2D near-surface observations, can be used to learn about 3D flow properties throughout the upper layer of the water column.
  • Preprint
    Generalized Lagrangian coherent structures
    ( 2018-01-19) Balasuriya, Sanjeeva ; Ouellette, Nicholas T. ; Rypina, Irina I.
    The notion of a Lagrangian Coherent Structure (LCS) is by now well established as a way to capture transient coherent transport dynamics in unsteady and aperiodic fluid flows that are known over finite time. We show that the concept of an LCS can be generalized to capture coherence in other quantities of interest that are transported by, but not fully locked to, the fluid. Such quantities include those with dynamic, biological, chemical, or thermodynamic relevance, such as temperature, pollutant concentration, vorticity, kinetic energy, plankton density, and so on. We provide a conceptual framework for identifying the Generalized Lagrangian Coherent Structures (GLCSs) associated with such evolving quantities. We show how LCSs can be seen as a special case within this framework, and provide an overarching discussion of various methods for identifying LCSs. The utility of this more general viewpoint is highlighted through a variety of examples. We also show that although LCSs approximate GLCSs in certain limiting situations under restrictive assumptions on how the velocity field affects the additional quantities of interest, LCSs are not in general sufficient to describe their coherent transport.
  • Article
    Investigating the connection between complexity of isolated trajectories and Lagrangian coherent structures
    (Copernicus Publications on behalf of the European Geosciences Union and the American Geophysical Union, 2011-12-15) Rypina, Irina I. ; Scott, S. E. ; Pratt, Lawrence J. ; Brown, Michael G.
    It is argued that the complexity of fluid particle trajectories provides the basis for a new method, referred to as the Complexity Method (CM), for estimation of Lagrangian coherent structures in aperiodic flows that are measured over finite time intervals. The basic principles of the CM are explained and the CM is tested in a variety of examples, both idealized and realistic, and in different reference frames. Two measures of complexity are explored in detail: the correlation dimension of trajectory, and a new measure – the ergodicity defect. Both measures yield structures that strongly resemble Lagrangian coherent structures in all of the examples considered. Since the CM uses properties of individual trajectories, and not separation rates between closely spaced trajectories, it may have advantages for the analysis of ocean float and drifter data sets in which trajectories are typically widely and non-uniformly spaced.
  • Article
    Investigating the eddy diffusivity concept in the coastal ocean
    (American Meteorological Society, 2016-06-29) Rypina, Irina I. ; Kirincich, Anthony R. ; Lentz, Steven J. ; Sundermeyer, Miles A.
    This paper aims to test the validity, utility, and limitations of the lateral eddy diffusivity concept in a coastal environment through analyzing data from coupled drifter and dye releases within the footprint of a high-resolution (800 m) high-frequency radar south of Martha’s Vineyard, Massachusetts. Specifically, this study investigates how well a combination of radar-based velocities and drifter-derived diffusivities can reproduce observed dye spreading over an 8-h time interval. A drifter-based estimate of an anisotropic diffusivity tensor is used to parameterize small-scale motions that are unresolved and underresolved by the radar system. This leads to a significant improvement in the ability of the radar to reproduce the observed dye spreading.
  • Article
    Mode filters and energy conservation
    (Acoustical Society of America, 2010-04-07) Udovydchenkov, Ilya A. ; Rypina, Irina I. ; Brown, Michael G.
    The discrete form of the mode filtering problem is considered. The relevant equations constitute a linear inverse problem. Solutions to problems of this type are subject to a well-known trade-off between resolution and precision. But unlike the typical linear inverse problem, the correctly formulated mode filtering problem is subject to an energy conservation constraint. This letter focuses on the importance of satisfying, approximately at least, the energy conservation constraint when mode filtering is performed.
  • Article
    Near-surface transport pathways in the north Atlantic Ocean : looking for throughput from the subtropical to the subpolar gyre
    (American Meteorological Society, 2011-05) Rypina, Irina I. ; Pratt, Lawrence J. ; Lozier, M. Susan
    Motivated by discrepancies between Eulerian transport estimates and the behavior of Lagrangian surface drifters, near-surface transport pathways and processes in the North Atlantic are studied using a combination of data, altimetric surface heights, statistical analysis of trajectories, and dynamical systems techniques. Particular attention is paid to the issue of the subtropical-to-subpolar intergyre fluid exchange. The velocity field used in this study is composed of a steady drifter-derived background flow, upon which a time-dependent altimeter-based perturbation is superimposed. This analysis suggests that most of the fluid entering the subpolar gyre from the subtropical gyre within two years comes from a narrow region lying inshore of the Gulf Stream core, whereas fluid on the offshore side of the Gulf Stream is largely prevented from doing so by the Gulf Stream core, which acts as a strong transport barrier, in agreement with past studies. The transport barrier near the Gulf Stream core is robust and persistent from 1992 until 2008. The qualitative behavior is found to be largely independent of the Ekman drift.
  • Article
    Support for the Slope Sea as a major spawning ground for Atlantic bluefin tuna: evidence from larval abundance, growth rates, and particle-tracking simulations
    (Canadian Science Publishing, 2021-10-20) Hernández, Christina M. ; Richardson, David E. ; Rypina, Irina I. ; Chen, Ke ; Marancik, Katrin E. ; Shulzitsk, Kathryn ; Llopiz, Joel K.
    Atlantic bluefin tuna (Thunnus thynnus) are commercially and ecologically valuable, but management is complicated by their highly migratory lifestyle. Recent collections of bluefin tuna larvae in the Slope Sea off northeastern United States have opened questions about how this region contributes to population dynamics. We analyzed larvae collected in the Slope Sea and the Gulf of Mexico in 2016 to estimate larval abundance and growth rates and used a high-resolution regional ocean circulation model to estimate spawning locations and larval transport. We did not detect a regional difference in growth rates, but found that Slope Sea larvae were larger than Gulf of Mexico larvae prior to exogenous feeding. Slope Sea larvae generally backtracked to locations north of Cape Hatteras and would have been retained within the Slope Sea until the early juvenile stage. Overall, our results provide supporting evidence that the Slope Sea is a major spawning ground that is likely to be important for population dynamics. Further study of larvae and spawning adults in the region should be prioritized to support management decisions.
  • Article
    Radium-based estimates of cesium isotope transport and total direct ocean discharges from the Fukushima Nuclear Power Plant accident
    (Copernicus Publications on behalf of the European Geosciences Union, 2013-03-28) Charette, Matthew A. ; Breier, Crystaline F. ; Henderson, Paul B. ; Pike, Steven M. ; Rypina, Irina I. ; Jayne, Steven R. ; Buesseler, Ken O.
    Radium has four naturally occurring isotopes that have proven useful in constraining water mass source, age, and mixing rates in the coastal and open ocean. In this study, we used radium isotopes to determine the fate and flux of runoff-derived cesium from the Fukushima Dai-ichi Nuclear Power Plant (FNPP). During a June 2011 cruise, the highest cesium (Cs) concentrations were found along the eastern shelf of northern Japan, from Fukushima south, to the edge of the Kuroshio Current, and in an eddy ~ 130 km from the FNPP site. Locations with the highest cesium also had some of the highest radium activities, suggesting much of the direct ocean discharges of Cs remained in the coastal zone 2–3 months after the accident. We used a short-lived Ra isotope (223Ra, t1/2 = 11.4 d) to derive an average water mass age (Tr) in the coastal zone of 32 days. To ground-truth the Ra age model, we conducted a direct, station-by-station comparison of water mass ages with a numerical oceanographic model and found them to be in excellent agreement (model avg. Tr = 27 days). From these independent Tr values and the inventory of Cs within the water column at the time of our cruise, we were able to calculate an offshore 134Cs flux of 3.9–4.6 × 1013 Bq d−1. Radium-228 (t1/2 = 5.75 yr) was used to derive a vertical eddy diffusivity (Kz) of 0.7 m2 d−1 (0.1 cm2 s−1); from this Kz and 134Cs inventory, we estimated a 134Cs flux across the pycnocline of 1.8 × 104 Bq d−1 for the same time period. On average, our results show that horizontal mixing loss of Cs from the coastal zone was ~ 109 greater than vertical exchange below the surface mixed layer. Finally, a mixing/dilution model that utilized our Ra-based and oceanographic model water mass ages produced a direct ocean discharge of 134Cs from the FNPP of 11–16 PBq at the time of the peak release in early April 2011. Our results can be used to calculate discharge of other water-soluble radionuclides that were released to the ocean directly from the Fukushima NPP.
  • Article
    Eddy-induced particle dispersion in the near-surface North Atlantic
    (American Meteorological Society, 2012-12) Rypina, Irina I. ; Kamenkovich, Igor V. ; Berloff, Pavel S. ; Pratt, Lawrence J.
    This study investigates the anisotropic properties of the eddy-induced material transport in the near-surface North Atlantic from two independent datasets, one simulated from the sea surface height altimetry and one derived from real-ocean surface drifters, and systematically examines the interactions between the mean- and eddy-induced material transport in the region. The Lagrangian particle dispersion, which is widely used to characterize the eddy-induced tracer fluxes, is quantified by constructing the “spreading ellipses.” The analysis consistently demonstrates that this dispersion is spatially inhomogeneous and strongly anisotropic. The spreading is larger and more anisotropic in the subtropical than in the subpolar gyre, and the largest ellipses occur in the Gulf Stream vicinity. Even at times longer than half a year, the spreading exhibits significant nondiffusive behavior in some parts of the domain. The eddies in this study are defined as deviations from the long-term time-mean. The contributions from the climatological annual cycle, interannual, and subannual (shorter than one year) variability are investigated, and the latter is shown to have the strongest effect on the anisotropy of particle spreading. The influence of the mean advection on the eddy-induced particle spreading is investigated using the “eddy-following-full-trajectories” technique and is found to be significant. The role of the Ekman advection is, however, secondary. The pronounced anisotropy of particle dispersion is expected to have important implications for distributing oceanic tracers, and for parameterizing eddy-induced tracer transfer in non-eddy-resolving models.
  • Article
    Connection between encounter volume and diffusivity in geophysical flows
    (Copernicus Publications on behalf of the European Geosciences Union, 2018-04-04) Rypina, Irina I. ; Llewellyn Smith, Stefan ; Pratt, Lawrence J.
    Trajectory encounter volume – the volume of fluid that passes close to a reference fluid parcel over some time interval – has been recently introduced as a measure of mixing potential of a flow. Diffusivity is the most commonly used characteristic of turbulent diffusion. We derive the analytical relationship between the encounter volume and diffusivity under the assumption of an isotropic random walk, i.e., diffusive motion, in one and two dimensions. We apply the derived formulas to produce maps of encounter volume and the corresponding diffusivity in the Gulf Stream region of the North Atlantic based on satellite altimetry, and discuss the mixing properties of Gulf Stream rings. Advantages offered by the derived formula for estimating diffusivity from oceanographic data are discussed, as well as applications to other disciplines.
  • Article
    Competition between chaotic advection and diffusion: Stirring and mixing in a 3-D eddy model.
    (European Geosciences Union, 2019-04-05) Brett, Genevieve ; Pratt, Lawrence J. ; Rypina, Irina I. ; Wang, Peng
    The importance of chaotic advection relative to turbulent diffusion is investigated in an idealized model of a 3-D swirling and overturning ocean eddy. Various measures of stirring and mixing are examined in order to determine when and where chaotic advection is relevant. Turbulent diffusion is alternatively represented by (1) an explicit, observation-based, scale-dependent diffusivity, (2) stochastic noise, added to a deterministic velocity field, or (3) explicit and implicit diffusion in a spectral numerical model of the Navier–Stokes equations. Lagrangian chaos in our model occurs only within distinct regions of the eddy, including a large chaotic “sea” that fills much of the volume near the perimeter and central axis of the eddy and much smaller “resonant” bands. The size and distribution of these regions depend on factors such as the degree of axial asymmetry of the eddy and the Ekman number. The relative importance of chaotic advection and turbulent diffusion within the chaotic regions is quantified using three measures: the Lagrangian Batchelor scale, the rate of dispersal of closely spaced fluid parcels, and the Nakamura effective diffusivity. The role of chaotic advection in the stirring of a passive tracer is generally found to be most important within the larger chaotic seas, at intermediate times, with small diffusivities, and for eddies with strong asymmetry. In contrast, in thin chaotic regions, turbulent diffusion at oceanographically relevant rates is at least as important as chaotic advection. Future work should address anisotropic and spatially varying representations of turbulent diffusion for more realistic models.
  • Article
    Chaotic advection in an archipelago
    (American Meteorological Society, 2010-09) Rypina, Irina I. ; Pratt, Lawrence J. ; Pullen, Julie ; Levin, Julia C. ; Gordon, Arnold L.
    Techniques from dynamical systems theory have been applied to study horizontal stirring of fluid in the Philippine Archipelago. The authors’ analysis is based on velocity fields produced by two high-resolution (3 and 6 km) numerical models. Particular attention is paid to identifying robust surface flow patterns and associating them with dominant Lagrangian coherent structures (LCSs). A recurrent wind-driven dipole in the lee of the coastline is considered in detail. The associated LCSs form a template for stirring, exchange, and biological transport in and around the dipole. Chaotic advection is argued to provide a relevant framework for interpreting mesoscale horizontal stirring processes in an archipelago as a whole. Implications for the formation of filaments, the production of tracer variance, and the scale at which stirring leads to mixing are discussed in connection with an observed temperature record.
  • Article
    Trajectory encounter volume as a diagnostic of mixing potential in fluid flows
    (Copernicus Publications on behalf of the European Geosciences Union&the American Geophysical Union, 2017-05-03) Rypina, Irina I. ; Pratt, Lawrence J.
    Fluid parcels can exchange water properties when coming into contact with each other, leading to mixing. The trajectory encounter mass and a related simplified quantity, the encounter volume, are introduced as a measure of the mixing potential of a flow. The encounter volume quantifies the volume of fluid that passes close to a reference trajectory over a finite time interval. Regions characterized by a low encounter volume, such as the cores of coherent eddies, have a low mixing potential, whereas turbulent or chaotic regions characterized by a large encounter volume have a high mixing potential. The encounter volume diagnostic is used to characterize the mixing potential in three flows of increasing complexity: the Duffing oscillator, the Bickley jet and the altimetry-based velocity in the Gulf Stream extension region. An additional example is presented in which the encounter volume is combined with the u approach of Pratt et al. (2016) to characterize the mixing potential for a specific tracer distribution in the Bickley jet flow. Analytical relationships are derived that connect the encounter volume to the shear and strain rates for linear shear and linear strain flows, respectively. It is shown that in both flows the encounter volume is proportional to time.
  • Article
    Spreading pathways of Pilgrim Nuclear Power Station wastewater in and around Cape Cod Bay: Estimates from ocean drifter observations
    (Elsevier, 2022-10-19) Rypina, Irina I. ; Macdonald, Alison ; Yoshida, Sachiko ; Manning, James P. ; Gregory, Margaret ; Rozen, Nimrod ; Buesseler, Ken
    Near-surface drifter observations were used to study the spreading pathways in and around the Cape Cod Bay from a source region located just offshore of the Pilgrim Nuclear Power Station. The study was motivated by the recent closing of the power plant and a possible release of accumulated wastewater. The investigation applies several different techniques to the drifter data set to estimate and quantify various aspects of the circulation and spreading. Our goal was two-fold: first, to better understand and predict the fate of the Pilgrim wastewater should it be released; and second, to review, compare, and contrast several useful techniques that can be applied to drifter datasets in other parts of the global ocean. Our analysis suggests weaker spreading of the wastewater plume inside the Bay than outside, and sensitivity of the advection pathways to the location of the release. Statistical techniques predicted that part of the plume would likely be advected cyclonically around the inner coastline of the Bay towards the more quiescent eastern regions, while another part of the plume would likely pass close to the tip of Cape Cod and the beaches of the Outer Cape. For the soluble radionuclides, the levels observed in our statistical model offshore of Provincetown and Dennis/Brewster will be at least 100 times smaller than the initial concentrations.
  • Article
    Applying dynamical systems techniques to real ocean drifters
    (European Geosciences Union, 2022-10-07) Rypina, Irina I ; Getscher, Timothy ; Pratt, Lawrence J ; Ozgokmen, Tamay
    This paper presents the first comprehensive comparison of several different dynamical-systems-based measures of stirring and Lagrangian coherence, computed from real ocean drifters. Seven commonly used methods (finite-time Lyapunov exponent (FTLE), trajectory path length, trajectory correlation dimension, trajectory encounter volume, Lagrangian-averaged vorticity deviation, dilation, and spectral clustering) were applied to 144 surface drifters in the Gulf of Mexico in order to map out the dominant Lagrangian coherent structures. Among the detected structures were regions of hyperbolic nature resembling stable manifolds from classical examples, divergent and convergent zones, and groups of drifters that moved more coherently and stayed closer together than the rest of the drifters. Many methods highlighted the same structures, but there were differences too. Overall, five out of seven methods provided useful information about the geometry of transport within the domain spanned by the drifters, whereas the path length and correlation dimension methods were less useful than others.
  • Article
    Fukushima-derived radionuclides in the ocean and biota off Japan
    (National Academy of Sciences, 2012-04-02) Buesseler, Ken O. ; Jayne, Steven R. ; Fisher, Nicholas S. ; Rypina, Irina I. ; Baumann, Hannes ; Baumann, Zofia ; Breier, Crystaline F. ; Douglass, Elizabeth M. ; George, Jennifer ; Macdonald, Alison M. ; Miyamoto, Hiroomi ; Nishikawa, Jun ; Pike, Steven M. ; Yoshida, Sachiko