ABSTRACT
We investigate the rising flux tube and the formation of sunspots in an unprecedentedly deep computational domain that covers the whole convection zone with a radiative magnetohydrodynamics ...simulation. Previous calculations had shallow computational boxes (<30 Mm) and convection zones at a depth of 200 Mm. By using our new numerical code Radition and RSST for Deep Dynamics(r2d2), we succeed in covering the whole convection zone and reproduce the formation of the sunspot from a simple horizontal flux tube because of the turbulent thermal convection. The main findings are as follows. (1) The rising speed of the flux tube is larger than the upward convection velocity because of the low density caused by the magnetic pressure and the suppression of the mixing. (2) The rising speed of the flux tube exceeds 250 m s−1 at a depth of 18 Mm, while we do not see any clear evidence of the divergent flow 3 h before the emergence at the solar surface. (3) Initially, the root of the flux tube is filled with the downflows, and then the upflow fills the centre of the flux tube during the formation of the sunspot. (4) The essential mechanisms for the formation of the sunspot are the coherent inflow and the turbulent transport. (5) The low-temperature region is extended to a depth of at least 40 Mm in the matured sunspot, with the high-temperature region in the centre of the flux tube. Some of the findings indicate the importance of the deep computational domain for the flux emergence simulations.
This paper presents a three-dimensional simulation of chromospheric jets with twisted magnetic field lines. Detailed treatments of the photospheric radiative transfer and the equations of state allow ...us to model realistic thermal convection near the solar surface, which excites various MHD waves and produces chromospheric jets in the simulation. A tall chromospheric jet with a maximum height of 10-11 Mm and lifetime of 8-10 minutes is formed above a strong magnetic field concentration. The magnetic field lines are strongly entangled in the chromosphere, which helps the chromospheric jet to be driven by the Lorentz force. The jet exhibits oscillatory motion as a natural consequence of its generation mechanism. We also find that the produced chromospheric jet forms a cluster with a diameter of several Mm with finer strands. These results imply a close relationship between the simulated jet and solar spicules.
ABSTRACT We investigate the effect of coronal temperature on the formation process of solar chromospheric jets using two-dimensional magnetohydrodynamic simulations of the region from the upper ...convection zone to the lower corona. We develop a new radiative magnetohydrodynamic code for the dynamic modeling of the solar atmosphere, employing an LTE equation of state, optically thick radiative loss in the photosphere, optically thin radiative loss in the chromosphere and the corona, and thermal conduction along the magnetic field lines. Many chromospheric jets are produced in the simulations by shock waves passing through the transition region. We find that these jets are projected farther outward when the coronal temperature is lower (similar to that in coronal holes) and shorter when the coronal temperature is higher (similar to that in active regions). When the coronal temperature is high, the deceleration of the chromospheric jets is consistent with the model in which deceleration is determined by the periodic chromospheric shock waves. However, when the coronal temperature is low, the gravitational deceleration becomes more important and the chromospheric jets approach ballistic motion.
Wildlife living around urbanized areas is often a cause of crucial issues such as zoonosis and wildlife-vehicle collisions. Despite this, residents hold positive views on the presence of urban ...wildlife primarily due to aesthetic reasons. This accepting attitude towards our coexistence with urban wildlife has made it difficult for wildlife managers to come to a consensus concerning the importance of human-urban wildlife conflicts. Although countermeasures such as lethal force and/or fencing are commonly used to control human-wildlife conflicts, these approaches are rarely applied in the case of urban wildlife. It is essential to recognize the gap between the current state of urban wildlife management and advanced scientific knowledge of urban wildlife behavior in order to mitigate urban deer conflicts. Fortunately, behavioral ecologists have been attempting to apply the perspective of individual differences, such as animal personality, to wildlife management. Studies have shown how the personalities of wildlife contribute to their adaptation to urban habitats. In order to prevent human-urban wildlife conflicts, recognizing the personalities of wildlife and selective culling of bold individuals should be conducted for deliberate selection for shyness when developing wildlife management plans. Making wildlife shy away from humans is essential to urban wildlife management. The aim of this study is to review observed measures against human-urban wildlife conflicts in Japan and to propose a new direction for innovative and effective approaches that takes animals personality into account to mitigate urban-wildlife conflicts. For this review we will target deer as a model species because deer are among the most serious of problem-causing urban wildlife.
Display omitted
•Urban deer increase the risk of zoonosis and wildlife–vehicle collisions.•No effective measures against the urban deer problem have been implemented to date.•Bolder deer prefer forest edge and tend to be urban deer.•Culling in forest edge area reduces boldness of the population.•Managing the boldness of deer populations could effectively resolve the problem.
Aims. We report the small temporal variation of the axial dipole moment near the solar minimum and its application to the solar-cycle prediction by the surface flux transport (SFT) model. Methods. We ...measure the axial dipole moment using the photospheric synoptic magnetogram observed by the Wilcox Solar Observatory (WSO), the ESA/NASA Solar and Heliospheric Observatory Michelson Doppler Imager (MDI), and the NASA Solar Dynamics Observatory Helioseismic and Magnetic Imager (HMI). We also use the SFT model for the interpretation and prediction of the observed axial dipole moment. Results. We find that the observed axial dipole moment becomes approximately constant during the period of several years before each cycle minimum, which we call the axial dipole moment plateau. The cross-equatorial magnetic flux transport is found to be small during the period, although a significant number of sunspots are still emerging. The results indicate that the newly emerged magnetic flux does not contribute to the build up of the axial dipole moment near the end of each cycle. This is confirmed by showing that the time variation of the observed axial dipole moment agrees well with that predicted by the SFT model without introducing new emergence of magnetic flux. These results allow us to predict the axial dipole moment at the Cycle 24/25 minimum using the SFT model without introducing new flux emergence. The predicted axial dipole moment at the Cycle 24/25 minimum is 60–80 percent of Cycle 23/24 minimum, which suggests the amplitude of Cycle 25 is even weaker than the current Cycle 24. Conclusions. The plateau of the solar axial dipole moment is an important feature for the longer-term prediction of the solar cycle based on the SFT model.
The morphological asymmetry of leading and following sunspots is a well-known characteristic of the solar surface. In the context of the large-scale evolution of the surface magnetic field, the ...asymmetry has been assumed to have only a negligible effect. Using the surface flux transport (SFT) model, we show that the morphological asymmetry of leading and following sunspots has a significant impact on the evolution of the large-scale magnetic field on the solar surface. By evaluating the effect of the morphological asymmetry of each bipolar magnetic region (BMR), we observe that the introduction of asymmetry to the BMR model significantly reduces the contribution to the polar magnetic field, especially for large and high-latitude BMRs. Strongly asymmetric BMRs can even reverse regular polar field formation. The SFT simulations based on the observed sunspot record show that the introduction of morphological asymmetry reduces the root-mean-square difference from the observed axial dipole strength by 30%-40%. These results indicate that the morphological asymmetry of leading and following sunspots has a significant effect on the solar cycle prediction.
The solar convection zone is filled with turbulent convection in highly stratified plasma. Several theoretical and observational studies suggest that the numerical calculations overestimate the ...convection velocity. Since all deep convection zone calculations exclude the solar surface due to substantial temporal and spatial scale separations, the solar surface, which drives the thermal convection with efficient radiative cooling, has been thought to be the key to solve this discrepancy. Thanks to the recent development in massive supercomputers, we are successful in performing the comprehensive calculation covering the whole solar convection zone. We compare the results with and without the solar surface in the local domain and without the surface in the full sphere. The calculations do not include the rotation and the magnetic field. The surface region has an unexpectedly weak influence on the deep convection zone. We find that just including the solar surface cannot solve the problem.
This study examined 1) the clinical relevance of trunk movement asymmetry, which was evaluated using a trunk-mounted inertial measurement unit (IMU), and 2) the relationship between trunk movement ...asymmetry and lower limb muscle strength asymmetry in individuals with knee osteoarthritis (OA).
One-hundred-thirty-one participants (mean age, 74.2 years; 71.8% female; Kellgren and Lawrence K&L grade ≥1) underwent gait analysis at their preferred pace for IMU-based measurement of trunk movement asymmetry (harmonic ratio HR and improved HR). The isometric strength of quadriceps and hip abductors was evaluated using a hand-held dynamometer. Pain and disability level were evaluated using a validated self-reported questionnaire. Multiple regression analyses with covariate adjustment were performed to examine the relationship between trunk movement asymmetry (independent variable) and pain, disability level, or muscle strength asymmetry (dependent variables).
Individuals with severe knee OA (K&L grade ≥3) had increased trunk movement asymmetry in the medio-lateral axis compared to those with a K&L grade of 1. Increased trunk movement asymmetry was associated with a greater knee pain and disability. The increased trunk movement asymmetry was significantly associated with an increase in the asymmetry of quadriceps strength, but not with asymmetry in the strength of hip abductor.
Our findings indicate that increased medio-lateral trunk movement asymmetry may be an indicator of impairment, rather than adaptation, in individuals with knee OA. This preliminary finding warrants validation by future study. Paying close attention to medio-lateral trunk movement asymmetry may be key to our understanding of OA-related pain and disability.
To test the hypothesis that the interaction between low back pain (LBP) and knee pain intensity contributes to the disability level of individuals with knee osteoarthritis (OA).
Community-dwelling ...participants with knee OA (Kellgren/Lawrence K/L grade ≥1) were enrolled. LBP and its severity were identified using questionnaires. Knee pain severity and disability level were evaluated using the Japanese Knee Osteoarthritis Measure (JKOM) subscale. Multiple linear regression analyses were performed to examine the effect of the LBP–knee pain interaction, an independent variable, on disability, a dependent variable.
A total of 260 participants (age, 48–88 years; 77.7% women) were included. Of them, 151 (58.1%) had LBP. The LBP–knee pain interaction was significantly associated with disability after the adjustment for covariates. A post-hoc subgroup analysis revealed that the relationship between knee pain intensity and disability level was higher in individuals with LBP (beta: 0.621 points; 95% confidence interval CI: 0.511 to 0.731 points) than in those without LBP (beta: 0.402 points; 95% CI: 0.316 to 0.487 points).
LBP interacts with knee pain intensity and contributes to disability level in individuals with knee OA. Coexisting LBP and knee pain had a stronger impact on disability level than LBP or knee pain alone. These findings highlight the potential deteriorative effects of the LBP–knee interaction on disability. Maximal treatment effects for disability might be achieved when LBP and knee pain are targeted simultaneously, rather than separately.
Context. The reduced speed of sound technique (RSST) has been used for efficient simulation of low Mach number flows in solar and stellar convection zones. The basic RSST equations are hyperbolic and ...are suitable for parallel computation by domain decomposition. The application of RSST is limited to cases in which density perturbations are much smaller than the background density. In addition, nonconservative variables are required to be evolved using this method, which is not suitable in cases where discontinuities such as shock waves coexist in a single numerical domain. Aims. In this study, we suggest a new semiconservative formulation of the RSST that can be applied to low Mach number flows with large density variations. Methods. We derive the wave speed of the original and newly suggested methods to clarify that these methods can reduce the speed of sound without affecting the entropy wave. The equations are implemented using the finite volume method. Several numerical tests are carried out to verify the suggested methods. Results. The analysis and numerical results show that the original RSST is not applicable when mass density variations are large. In contrast, the newly suggested methods are found to be efficient in such cases. We also suggest variants of the RSST that conserve momentum in the machine precision. The newly suggested variants are formulated as semiconservative equations, which reduce to the conservative form of the Euler equations when the speed of sound is not reduced. This property is advantageous when both high and low Mach number regions are included in the numerical domain. Conclusions. The newly suggested forms of RSST can be applied to a wider range of low Mach number flows.