Original Research ARTICLE
On Destruction of a Thermally Stable Layer by Compositional Convection in the Earth's Outer Core
- 1Research Institute for Mathematical Sciences, Kyoto University, Japan
- 2Department of Mathematics, Kyoto University, Japan
We discuss destruction of a thermally stable layer in the upper part of the Earth's outer core by compositional convection excited at the inner core boundary. We propose to use the radial distribution of power induced by thermal and compositional buoyancy (rate of kinetic energy production) as a measure of occurrence of thermal and compositional convection. The power consists of the terms proportional to convective entropy flux and convective compositional flux. In the region with positive power, convection is active because kinetic energy can be produced by buoyancy force, and a stably stratified layer could not be formed there. On the other hand, in the region with negative power, convection is suppressed and a stably stratified layer may be produced. Considering penetration effect of convection, we discuss possible maximum and minimum thicknesses of the stable layer based on the radial distribution of power and its radial integral, respectively. We construct a 1-dimensional thermal and compositional balance model of the Earth's core with a larger value of thermal conductivity recently suggested by high-pressure experiments and first principle calculations, and estimate radial distributions of power for various values of core mantle boundary (CMB) heat flux QCMB. When QCMB > QsCMB no thermally stable layer can exist, where QsCMB is the conductive heat flux along the adiabat at the CMB. On the other hand, when QCMB < QsCMB, formation of an upper thermally stable layer becomes possible, depending on the extent of penetration of compositional convection excited below. When QCMB is sufficiently lower than QsCMB, a thermally stable layer survives the maximum penetration of compositional convection. The results show that a thermally stable layer becomes effectively thinner when the effect of compositional convection is considered compared with the results of previous studies where the existence of a stable layer is evaluated based on the convective flux only.
It is widely thought that fluid motion driven in the Earth's outer core by cooling due to mantle convection generates and maintains the geomagnetic field through dynamo action. Although thermal convection would occur as a direct effect of cooling from the core mantle boundary (CMB) and as a result of the latent heat release at the inner core boundary (ICB) associated with growth of the inner core through solidification of Fe and Ni, the main energy source for geomagnetic field generation is considered to be supplied by compositional convection induced by the release of light elements at the ICB (e.g., Lister and Buffett, 1995). Estimation of the power available for magnetic energy as well as evolution of the structure of the Earth's core has been studied extensively using 1-dimensional thermal and compositional budget models (e.g., Loper, 1978; Gubbins et al., 1979; Stevenson, 1983).
Most of the models assumed vigorous convection mixes completely the whole fluid core, and entropy and composition are homogeneous. Alternative models considering formation of a stably stratified layer below the CMB were proposed (Labrosse et al., 1997; Lister and Buffett, 1998). They assumed that, when heat flow at the CMB QCMB determined by strength of mantle convection becomes smaller than upward conductive heat flow along the adiabatic temperature profile Qs, heat accumulates below the CMB, and a thermally stable layer starts to develop by conduction. Labrosse et al. (1997) assumed that erosion of the stable layer by mixing due to compositional convection generated from the deeper region does not occur, whereas Lister and Buffett (1998) considered accumulation of light elements due to inner core growth explicitly. They discussed the competition between light elements accumulating in the convecting region of the core and heat accumulating in the stratified layer due to a subadiabatic heat flux at the CMB.
Recently, it was suggested from high-pressure experiments and first principle calculations that the values of thermal conductivity under conditions of planetary cores are larger than those considered so far (Gomi et al., 2013; Pozzo et al., 2014). By using 1-dimensional thermal balance models with the updated values of thermal conductivity, generation and existence of a stably stratified layer in the Earth's outer core was discussed (Gomi et al., 2013; Labrosse, 2015). Their results showed that a stable layer with a thickness of O(1,000 km) could be produced when the heat flux across the CMB is small. They assumed that the region with negative convective heat flux is stably stratified. This assumption seems to be appropriate when convection is driven only by thermal effects, however, it is not correct when compositional convection occurs simultaneously. When compositional convection is sufficiently vigorous to overcome thermally stable stratification, it would mix up the stable layer and make it neutral.
Whether the stably stratified layer is formed and maintained below the CMB or it is destroyed by compositional convection from the deeper region is a significant issue for the evolution of the Earth's core structure. Formation of the stable layer affects the thermal history of the core, such as the age of the inner core and field intensity (e.g., Labrosse et al., 1997; Lister and Buffett, 1998). It is also related to geomagnetic secular variation observed at the surface of the Earth, since it could originate from the fluid motions in the stable layer if it exists (e.g., Braginsky, 1984, 1993; Buffett, 2014).
In the present study, we propose to use radial distribution of power induced by thermal and compositional buoyancy (rate of kinetic energy production) as a measure of occurrence of thermal and compositional convection. The power consists of the terms proportional to heat flux and compositional flux. In the region with positive power convection is active because kinetic energy can be produced by buoyancy force. On the other hand, in the region with negative power, convection is suppressed and the stably stratified layer may be produced.
In the following, formulation of a 1-dimensional thermal and compositional budget model is recalled, and a criterion for formation and destruction of the stable layer is proposed in section 2. Then we apply the model and the criterion to a simple core model, and illustrate formation of a stable layer and its thickness depending on the CMB heat flux in section 3. Summary and discussion are in the final section.
2. Model and Method
2.1. Estimation of Stratified Layer
The procedure for calculation of the power induced by thermal and compositional buoyancy (rate of kinetic energy production) in 1-dimensional thermal and compositional budget models was proposed by Lister and Buffett (1995). They did not apply it to the stable layer formation problem but estimated available energy for magnetic field generation using the integral value of the power in the whole domain. We propose to use the radial distribution of power as a criterion for the formation and destruction of the stable layer. The rationale for this criterion is as follows.
We assume that the whole outer core is well mixed by convection, so that concentrations of light elements and entropy are homogeneous in the radial direction. Temperature follows an adiabatic profile. Then the radial distribution of power by buoyancy force (kinetic energy production) wb(r) can be calculated from distributions of convective entropy flux and convective compositional flux, which are diagnostically obtained by thermal and compositional budgets as explained below.
The regions where kinetic energy production is positive (wb(r) > 0) can be mixed up by convection since it can emerge there, which is consistent with the assumption. Therefore, a stable layer is not formed there. On the other hand, in the regions where kinetic energy production is negative (wb(r) < 0) convection cannot be driven locally, and mixing may not occur there. However, the region would be mixed up even when wb(r) < 0 if convection driven in the lower part penetrates into there. Penetration could occur when total kinetic energy production in the convecting and penetrating regions is still positive. Then, possible penetration of convection can be evaluated by the radial integral of power , where c is the radius of ICB. We estimate the region with Wb(r)>0 which could be mixed up by penetrative convection (Figure 1).
Therefore, the maximum and minimum thicknesses of the stable layer are determined with wb(r) < 0 and Wb(r) < 0, respectively.
2.2. 1-dimensional Thermal and Compositional Balance Model
The equation for compositional budget is
where XC is averaged concentration increase of light elements, FC is convective compositional flux, and SC is light elements release at the ICB.
where r is radius, b and c are radii of the CMB and ICB, respectively, is the homogeneous concentration of light elements in the outer core, C′ is the concentration fluctuation from , ur the radial component of velocity, the convective compositional flux, ρa and ρc are the horizontally averaged density in the outer core and at the ICB, respectively. Here, the light elements are assumed to be perfectly incompatible in the inner core for simplicity. Total compositional budget is expressed by XC(b) = SC. Using (1) and the boundary condition FC(b) = 0, FC can be described alternatively as follows:
The equation for the thermal budget is
where Qconv, Qcooling, Ecomp, Qlatent and Qs are the convective heat flux, and its contributions by secular cooling, compositional energy, latent heat release, and thermal conduction, respectively, which are described as follows:
where Ta(r) is the adiabatic temperature profile in the outer core, Sa is entropy, which is assumed to be uniform, S′ is the entropy fluctuation from Sa, is the convective entropy flux, ΔS is the entropy jump associated with the crystallization of the inner core, TL is liquidus temperature, k(r) is thermal conductivity and μ′(r) is the chemical potential based on the value at the ICB. Following the discussion by Lister and Buffett (1995), we assume that Joule heating QJ, viscous dissipation Qvisc and the heat flow coming from the inner core QICB can be neglected in the thermal budget (7). From the total thermal budget in the whole core, we obtain
Previous studies showed that all the terms on the right hand side of (13) are proportional to dc/dt, and then, the radial distributions of all the terms are obtained for given CMB heat flow QCMB.
The convective entropy flux FS is expressed as
The power by buoyancy forces wb(r) is calculated with the convective fluxes as follows:
where and are thermal and compositional expansion coefficients, respectively.
2.3. Analytical Expressions of the Fluxes in a Simple Model
Here, we illustrate the result of a diagnosis on the formation of a stable layer in the present Earth's core with a simple setup similar to that used in Lister and Buffett (1995).
Density in the outer core is assumed to be constant. Gravity is given in the form of g(r) = go(r/b), where g0 is the gravitational acceleration at the CMB. Then adiabatic temperature Ta(r, t) and its radial gradient can be described as
where Ti(t) is temperature at the ICB and
Since for the values of the Earth's core, we approximate fa(r) as , and get
Accordingly, Qs(r) becomes
The present value of Ti is determined by setting the conductive heat flow along the adiabat at the CMB QsCMB = Qs(b).
On the other hand, the time derivative of Ti(t) should be evaluated to calculate Qcooling. From entropy variation with respect to temperature, pressure and composition variations
where Ti is temperature at the ICB, is heat of reaction, which is neglected for simplicity as well as contribution of pressure variation. Then, we have
Since the chemical potential is calculated as
its radial integral becomes
Then we obtain, from equation (10):
Here, we used the relation
derived by the total compositional budget XC(b) = SC.
The growth rate of the inner core dc/dt is determined by the global thermal budget Equation (13). By using the value of dc/dt, the convective compositional flux is calculated by (5), while the convective thermal flux is determined by
Finally, the convective entropy flux is estimated by (14), and the power by buoyancy forces is determined by (15).
The radial profile of thermal conductivity is approximated to a recent large estimated value (Gomi et al., 2013) with a quadratic function as , where kb, k0 are the values at the CMB and the center of the core, respectively. The adiabatic temperature at CMB is given so that Qs(b) = QsCMB = 9.3TW. The parameters used in the calculation are listed in Table 1.
Figure 2 shows the results in the case with QCMB = 12TW. The radial distributions of contributions to the thermal budget in the left panel of Figure 2 are qualitatively consistent with previous studies (e.g., Figure 9 in Gomi et al., 2013) although there are several differences. In this case, no stable layer is formed since Qconv is positive in the whole outer core. Kinetic energy production wb(r) shown in the center panel of Figure 2 is positive in the whole outer core. Note that both compositional and thermal contribution to wb(r) is positive in the whole outer core, indicating that thermal and compositional convection emerges there.
Figure 2. Result of 1-dimensional thermal and compositional budget model for QCMB = 12TW. (A) Radial distributions of heat flows. Black, blue, purple, red and green lines indicate Qconv, Qcooling, Ecomp, Qlatent, and Qs, respectively. (B) Radial distributions of work by buoyancy forces wb(r). Green and red lines indicate contributions from compositional and thermal buoyancies. Blue line shows total work. (C) Radially integrated work by buoyancy forces Wb(r).
Figure 3 shows the results in the case with QCMB = 8TW. In this case, a stable layer is possible since Qconv becomes negative. The region with wb(r) < 0 below the CMB appears in the center panel of Figure 3. The middle panel of Figure 3 shows the compositional contribution to wb(r) is positive in the whole outer core while the thermal contribution is negative at r≳2, 200 km, meaning that compositional convection is possible to emerge in the upper part of the outer core although thermal and compositional convection can be excited in the lower part. On the other hand, Wb(r) is positive in the whole outer core as shown in the right panel of Figure 3. This means that formation of a stable layer depends on the extent of penetration of compositional convection from below. A stable layer survives below the CMB when penetration is weak, whereas a stable layer is completely destroyed by compositional convection when its penetration is strong. Note that the stable layer is thinner than that estimated by the criterion with thermal convective heat flux only (Qconv < 0) used in the previous studies.
Figure 3. As in Figure 2 but for QCMB = 8TW.
Figure 4 shows the result in the case with QCMB = 4TW. In the right panel of Figure 4, there appears the region with Wb(r) < 0 below CMB, meaning that a stable layer is formed even when the extent of penetration of convection is maximum in this case. Note that the thermal contribution is negative in the almost whole outer core in this case.
Figure 4. As in Figure 2 but for QCMB = 4TW.
Figure 5 shows the summary of stable layer formation as a function of QCMB. The red and blue lines indicate the thickness of a stable layer in the case of minimum and maximum penetration, respectively. It is found that no stable layer is formed when QCMB>QsCMB = 9.3TW, whereas a stable layer is absolutely formed when QCMB < 4.4TW. We also show the bottom of a stable layer estimated by the criterion with thermal convective flux only by the black line for comparison. It is clear that destruction by compositional convection cannot be ignored in order to estimate the thickness of the stable layer.
Figure 5. Possible stable layer thickness as a function of CMB heat flow. The blue and red lines indicate the thickness of the minimum and maximum stable layers, respectively. The black line shows the thickness of a stable layer estimated by thermal convective flux only.
4. Conclusions and Discussion
We proposed to use radial distribution of power by thermal and compositional buoyancy (rate of kinetic energy production) as a new criterion for estimating formation and destruction of a stably stratified layer by thermal and compositional convection. The method for calculating power by buoyancy forces in a 1-dimensional thermal and compositional budget model was presented, and was applied to a simple 1-dimensional model for the Earth's core. The thickness of the stable layer formed below the CMB is estimated for various values of the CMB heat flow considering the effect of mixing by compositional convection in the deeper region. When QCMB>QsCMB(= 9.3TW), no stable layer is formed due to occurrence of compositional and thermal convection in the whole outer core. When QCMB<QsCMB(= 9.3TW), formation of an upper thermally stable layer becomes possible, depending on the extent of penetration of compositional convection excited at ICB. When QCMB < 4.4TW, a stable layer can be formed even when penetration of compositional convection is maximal. Note that this critical value of QCMB varies depending on the value of QsCMB and model configuration for Earth's core. For example, constant density distribution and degree 2 polynomial isentropic temperature used in this study would not be precise enough for estimation of stable layer formation, since some discrepancies between previous studies were found to originate from different levels of approximation in the isentropic temperature profile (e.g., Labrosse, 2015). It is necessary to use more precise Earth models with more realistic distributions of density, isentropic temperature and other physical properties for definitive conclusions about the Earth's core.
Since the present CMB heat flow is expected in the range 5–15 TW (e.g., Lay et al., 2008), our results also support possibility of a stable layer formation below the CMB as the previous studies. However our estimations show the thickness of the stable layer is significantly thinner than that estimated with the criterion only by thermal convective flux used in the previous studies due to mixing by compositional convection. For example, when QCMB = 8TW, the thickness of the stable layer can be 250 km at most, which is about 1,000 km thinner than the estimation by convective heat flux criterion (see Figure 5).
The present paper significantly contributes to thermal and compositional evolution problems in planetary cores by giving dynamically consistent estimation of stable layer formation. All the previous studies did not consider whether kinetic energy is adequately supplied to the convecting part or not, which should be checked by calculating kinetic energy production distribution.
We estimate emergence of a stable layer by assuming that the whole outer core is convecting. This procedure seems to be contradictory since existence of a stable layer breaks the assumption of the whole outer core convection. As a matter of fact, the solutions obtained by our method are first order approximations of quasi-steady states where concentrations of ight elements and temperature averaged in the whole outer core vary in time while their radial distributions are stationary. Spherically averaged conservation equation of light elements without assuming a homogeneous mixed state is described as follows:
where C(r, t) is concentration of light elements, FC(r, t) and FD(r, t) are convective and non-convective fluxes such as diffusion, respectively. Assuming quasi-steady state, , and integrating from r = c with the boundary condition FC(c)+FD(c) = SC, Equation (1) is modified as
This equation means that we can construct new solutions by replacing FC obtained with the assumption of whole outer core convection by FC+FD. For example, the flux distributions with FD = 0 (r<rt) and FC = 0 (r>rt) provide a quasi-steady solution for a given boundary r = rt. The upper part of C′(r) can be constructed from FD(r). This solution expresses the state where convection is active and mixes up the lower part r<rt, while convection is suppressed and diffusion is dominant in the upper part r>rt, forming stable stratification. Similar argument is possible for temperature by assuming T(r, t) = Ti(t)f(r), however, the thermal fluxes would be varied due to the difference of the temperature structure between the solutions with and without a stable layer. Then, in order to estimate thickness of a stable layer more precisely, after determining its bottom boundary with the criterion of kinetic energy production, we should give a certain temperature distribution in a stable layer and calculate distribution of kinetic energy production again. By repeating this procedure we can obtain more precise estimation of a thickness of a stable layer. For example, given a thermally diffusive temperature profile stable layer, re-calculations of wb lead to 0.07% and 1% differences of a stable layer thickness for QCMB = 8TW and 4TW (examples of Figures 3, 4), respectively.
Quasi-steady stable layer thickness obtained in this study does not consider time scale for its formation, which would be governed by thermal diffusion. In order to discuss stable layer formation in the Earth's core more precisely, we should compare the diffusion time of possible thickness of the stable layer with a possible age of the inner core. When the diffusion time is sufficiently short, estimation with our method would be effective. In contrast, when the diffusion time is longer, we should solve thermal and chemical evolution of the core as a time-dependent problem in order to estimate formation of the stable layer, for example, as performed by Lister and Buffett (1998). For example, when the thickness of the stable layer is 100km, and thermal diffusivity is 10−5 m2/s, the thermal diffusive time becomes 30Myr. This is presumably much less than the age of the inner core, since it was estimated as O(Gyr) by the previous studies without the effects of the stable layer and the stable layer would be so thin that it would not largely affect the age of the inner core. In contrast, when the thickness is 1000km, the thermal diffusive time becomes 3Gyr, which is probably larger than the age of the inner core if it is not affected by the presence of the stable layer. Note that in the case of destruction of a pre-existing stable layer due to thermal and compositional convection excited at the ICB, we should examine time scale for destruction (erosion rate), which would be determined not only by thermal diffusion but also by other physical properties, such as convective velocity, strength of stratification and rotation, and so on.
Our model assumes that thermal and/or compositional convection mixes entropy and composition uniformly because thermal and/or compositional Rayleigh numbers become so large that convective motion is turbulent due to the smallness of the molecular diffusivities of the outer core. However, there are possibilities that convection with coherent structures emerges. For example, large scale columnar convection affected by the Earth's rotation may coexist with small scale turbulent convective motion due to larger values of eddy diffusivities which may contribute to decrease the effective Rossby number for large scale fluid motions. The extent of penetration of columnar convection into the upper stable layer depends on the ratio between the Brunt-Väisälä frequency and the rotation frequency, and on the typical width scale of the convective columns in the linear theory (Takehiro and Lister, 2001). However, it is not known whether the penetration of columnar convection would be an effective mechanism to erode the stable layer. Furthermore, it has been proposed that the light elements released at the ICB may form chemical plumes and/or blobs rising through the outer core without mixing, which could create a chemically stratified layer at the top of the core (e.g., Braginsky, 1993; Shimizu and Loper, 1997). Recent numerical calculations by Manglik et al. (2010) show that chemical plumes can penetrate into the thermally stratified layer and keep rising to some extent. In order to clarify these possibilities, 3-dimensional numerical simulations are needed although it may be difficult to resolve fluid motions with wide spatial spectral ranges.
In the present paper, we illustrate formation of thermally stable layer below the CMB when the heat flow is drawn by mantle convection and its destruction by compositional convection excited at the ICB. The method presented here can also be applied to compositionally stable layer and its destruction by thermal convection. It has been proposed that a stable layer would be formed below the CMB by accumulation of light elements released at the ICB (e.g., Loper and Roberts, 1983) or by diffusion of light elements from the mantle (e.g., Buffett and Seagle, 2010, 2011). By introducing contribution terms of barodiffusion and injection of light elements from the CMB, we can model stable layer formation including the effects of mixing by convection in a similar way with radial distribution of power by buoyancy forces.
In the layer diagnosed as thermally stable stratification, radial temperature gradient contributes to stabilize but that of concentration to destabilize. This configuration permits emergence of finger-type double diffusive convection. However, since finger-type convection mainly transports light elements (Turner, 1967) and increase stability through weakening radial concentration gradient, the stable layer would not be destroyed by this effect. In the case of a compositional stable layer, diffusive-type convection may emerge which mainly transports heat (Turner, 1965). Therefore, the stable layer would be maintained even when double diffusive convection occurs, although double diffusive convection may break up into a series of alternating convecting layers in some specific conditions.
Note that our method can apply not only to planetary cores but also to other general situations where convective and stable layers coexist. In general, the stable layer thickness, which is expected to be between the maximum and minimum estimates evaluated by Wb(r) = 0 and wb(r) = 0, would be determined through complicated combined effects tied to the dynamics of chemical plumes and their possible interaction with the thermally stable layer, which are poorly understood. Nevertheless, in some extreme cases, the maximum or minimum estimates would give good stable layer thicknesses. For example, in planetary atmospheres where the effect of inertia is strong and advection of kinetic energy dominates local kinetic energy dissipation, the radially integrated kinetic energy production rate Wb(r) presumably gives a good estimation of stable layer thickness. In contrast, in the Earth's outer core, it is expected that the inertia terms in the equation of fluid motion can be neglected since the value of compositional diffusivity is significantly smaller than that of kinetic viscosity and the work done by convection is converted into magnetic energy. This means that advection of kinetic energy can be ignored and penetration of compositional convection scarcely occurs. Therefore, we can estimate thickness of a stable layer with local kinetic energy production rate wb(r) only. For example, the red line in Figure 5 shows that a stable layer with a thickness of 250 km is formed when QCMB = 8TW in the simple model presented in the previous section. We should verify our new criterion by performing numerical experiments of 2- or 3-dimensional hydrodynamic models, and check the expectation mentioned above. Effects of double diffusive convection possibly emerging in the stable layer should be examined as well. Moreover, by investigating parameter dependence of extent of penetration of convection, more precise estimation of stable layer formation will be possible.
ST did most of the theoretical formulation of this work. YS significantly contributed to 1D model development and figure generation.
This work was supported by JSPS (Japan Society for Promotion Science) KAKENHI Grant Numbers 16H01117 and 15H05834.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
We thank three reviewers and the editor for their valuable comments, which were very helpful for improving the manuscript. This work was supported by the Research Institute for Mathematical Sciences, a Joint Usage/Research Center located in Kyoto University.
Gomi, H., Ohta, K., Hirose, K., Labrosse, S., Caracas, R., Verstraete, M. J., et al. (2013). The high conductivity of iron and thermal evolution of the earth's core. Phys. Earth Planet. Inter. 224, 88–103. doi: 10.1016/j.pepi.2013.07.010
Loper, D. E., and Roberts, P. H. (1983). “Compositional convection and the gravitationally powered dynamo,” in Stellar and Planetary Magnetism, ed A. M. Soward (New York, NY: Gordon Bearch Science Publlication), 297–327.
Pozzo, M., Davies, C., Gubbins, D., and Alfe, D. (2014). Thermal and electrical conductivity of solid iron and iron-silicon mixtures at earth's core conditions. Earth Planet. Sci. Lett. 393, 159–164. doi: 10.1016/j.epsl.2014.02.047
Shimizu, H., and Loper, D. E. (1997). Time and length scales of buoyancy-driven flow structures in a rotating hydromagnetic fluid. Phys. Earth Planet. Inter. 104, 307–329. doi: 10.1016/S0031-9201(97)00035-6
Takehiro, S., and Lister, J. R. (2001). Penetration of columnar convection into an outer stably stratified layer in rapidly rotating spherical fluid shells. Earth Planet. Sci. Lett. 187, 357–366. doi: 10.1016/S0012-821X(01)00283-7
Keywords: thermal conductivity, heat flux, compositional flux, power by buoyancy forces, kinetic energy production, core dynamics, core stratification
Citation: Takehiro S and Sasaki Y (2018) On Destruction of a Thermally Stable Layer by Compositional Convection in the Earth's Outer Core. Front. Earth Sci. 6:192. doi: 10.3389/feart.2018.00192
Received: 28 May 2018; Accepted: 17 October 2018;
Published: 02 November 2018.
Edited by:Hagay Amit, University of Nantes, France
Reviewed by:Stéphane Labrosse, École Normale Supérieure de Lyon, France
Mathieu Bouffard, Max-Planck-Institut für Sonnensystemforschung, Germany
Copyright © 2018 Takehiro and Sasaki. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Shin-ichi Takehiro, email@example.com