Krishna Khanal
Current Organizations / Institutes
The organizations/institutes from the scholar's latest article
Research Subjects
Krishna Khanal has a total of 21 co-authors. Here are the top co-authors by number of shared publications:
- Gary P. Zank
University of Alabama in Huntsville
4 shared publications - Binod Adhikari
Department of Physics, St. Xavier’s College, Tribhuvan University, Kathmandu, Nepal
3 shared publications - Binod Bhattarai
10Department of Physics, University of California, Merced, 5200 North Lake Road, Merced, CA 95343, USA
3 shared publications Narayan P. Chapagain
3 shared publications- Laxman Adhikari
University of Alabama in Huntsville
2 shared publications Nirakar Sapkota
2 shared publications
Latest Publications
- Open Access
- 19 June 2026
- Published by American Astronomical Society in The Astrophysical Journal
- Vol. 1005 (1)
A Quasi-2D Turbulence Model of Coronal Heating Including Thermal Conduction and Radiative Cooling
One of the key mechanisms proposed for coronal heating is dissipation of quasi-two-dimensional (2D) turbulence generated above the photosphere by continual mixing and convective shuffling of open and closed magnetic field lines within the so-called magnetic carpet. This class of heating mechanism is modeled within the Nearly Incompressible—Magnetohydrodynamics (NI-MHD) framework and has previously been shown to supply sufficient energy to drive the solar wind to supersonic and super-Alfvénic speeds and heat the corona to temperatures exceeding one million kelvin. In addition to turbulent heating, the overall energy budget of the low corona and the chromospheric-coronal transition region is strongly influenced by loss mechanisms associated with thermal conduction and radiative cooling, which must therefore be included in realistic models. To that end, we present solutions to a coupled system of NI-MHD turbulence transport equations and solar wind evolution equations that self-consistently incorporates the effects of collisional thermal conduction and radiative cooling. Results show that an interplay between turbulence dissipation and heat-loss mechanisms can result in the expected solar-wind driving as well as coronal heating. The predicted density and thermal conduction flux at the inner boundary, as well as the overall wind speed, are broadly consistent with representative ranges of observed coronal and solar wind parameters. The model results are sensitive to the ratio of von Kármán–Taylor constant to correlation length, as well as to the assumed expansion factor of the associated flux tube. This sensitivity is explored through a parametric survey.
- 19 June 2026
- Published by American Astronomical Society in The Astrophysical Journal
- Vol. 1005 (1)
Open AccessA Quasi-2D Turbulence Model of Coronal Heating Including Thermal Conduction and Radiative Cooling
One of the key mechanisms proposed for coronal heating is dissipation of quasi-two-dimensional (2D) turbulence generated above the photosphere by continual mixing and convective shuffling of open and closed magnetic field lines within the so-called magnetic carpet. This class of heating mechanism is modeled within the Nearly Incompressible—Magnetohydrodynamics (NI-MHD) framework and has previously been shown to supply sufficient energy to drive the solar wind to supersonic and super-Alfvénic speeds and heat the corona to temperatures exceeding one million kelvin. In addition to turbulent heating, the overall energy budget of the low corona and the chromospheric-coronal transition region is strongly influenced by loss mechanisms associated with thermal conduction and radiative cooling, which must therefore be included in realistic models. To that end, we present solutions to a coupled system of NI-MHD turbulence transport equations and solar wind evolution equations that self-consistently incorporates the effects of collisional thermal conduction and radiative cooling. Results show that an interplay between turbulence dissipation and heat-loss mechanisms can result in the expected solar-wind driving as well as coronal heating. The predicted density and thermal conduction flux at the inner boundary, as well as the overall wind speed, are broadly consistent with representative ranges of observed coronal and solar wind parameters. The model results are sensitive to the ratio of von Kármán–Taylor constant to correlation length, as well as to the assumed expansion factor of the associated flux tube. This sensitivity is explored through a parametric survey.
- Open Access
- 10 June 2026
- Published by American Astronomical Society in The Astrophysical Journal
- Vol. 1004 (2) , 143
Turbulence and Its Potential Impact on Solar Chromospheric and Coronal Heating
Low-frequency turbulence in the solar chromosphere remains poorly understood. We address (1) the sources of low-frequency turbulence that potentially heat the chromosphere, and (2) how turbulence is transported and dissipated throughout the chromosphere and lower corona. We use particle-in-cell simulations to investigate mixed polarity magnetic fields corresponding to emergent magnetic carpet field in coronal holes or quiet Sun regions for strong (imbalanced) and weak (balanced) guide magnetic fields. The initial mixed polarity magnetic field transitions rapidly to a turbulent state dominated by advected small-scale nonlinear structures, with a minority slab turbulence population, and the emergent field is largely annihilated. Turbulence is anisotropic for an imbalanced magnetic field and more isotropic for balanced cases. We develop a transport model for turbulence advected and dissipated throughout the chromosphere by randomly distributed energy-containing scale dynamical flows described by log-normal statistics. We compute the expectations for the total energy per unit volume 〈 y 〉( h ) J m −3 , the Elsässer specific energy 〈 〈 Z ∞ 2 〉 〉 ( h ) m 2 s −2 , the heating rate 〈 H ̇ 〉 ( h ) J m −3 s −1 , and the correlation length 〈 λ 〉( h ) in kilometers as functions of height h above the photosphere. Turbulent energy is injected into the low corona by a random “patchwork” of sites across the transition region surface. The expected energy injection rates 〈 S ̇ 〉 J m −2 s −1 for the chromosphere and at the base of the corona exceed the estimated energy requirements needed to heat both the chromosphere and corona. Similarly, we show that spicules can be heated gradually with increasing height by entrained magnetic carpet and photospheric turbulence.
- 10 June 2026
- Published by American Astronomical Society in The Astrophysical Journal
- Vol. 1004 (2) , 143
Open AccessTurbulence and Its Potential Impact on Solar Chromospheric and Coronal Heating
Low-frequency turbulence in the solar chromosphere remains poorly understood. We address (1) the sources of low-frequency turbulence that potentially heat the chromosphere, and (2) how turbulence is transported and dissipated throughout the chromosphere and lower corona. We use particle-in-cell simulations to investigate mixed polarity magnetic fields corresponding to emergent magnetic carpet field in coronal holes or quiet Sun regions for strong (imbalanced) and weak (balanced) guide magnetic fields. The initial mixed polarity magnetic field transitions rapidly to a turbulent state dominated by advected small-scale nonlinear structures, with a minority slab turbulence population, and the emergent field is largely annihilated. Turbulence is anisotropic for an imbalanced magnetic field and more isotropic for balanced cases. We develop a transport model for turbulence advected and dissipated throughout the chromosphere by randomly distributed energy-containing scale dynamical flows described by log-normal statistics. We compute the expectations for the total energy per unit volume 〈 y 〉( h ) J m −3 , the Elsässer specific energy 〈 〈 Z ∞ 2 〉 〉 ( h ) m 2 s −2 , the heating rate 〈 H ̇ 〉 ( h ) J m −3 s −1 , and the correlation length 〈 λ 〉( h ) in kilometers as functions of height h above the photosphere. Turbulent energy is injected into the low corona by a random “patchwork” of sites across the transition region surface. The expected energy injection rates 〈 S ̇ 〉 J m −2 s −1 for the chromosphere and at the base of the corona exceed the estimated energy requirements needed to heat both the chromosphere and corona. Similarly, we show that spicules can be heated gradually with increasing height by entrained magnetic carpet and photospheric turbulence.
- Open Access
- 7 September 2025
- Published by American Geophysical Union (AGU) in Earth and Space Science
- Vol. 12 (9)
Suppression of Magnetopause Reconnection in the Presence of Cold Magnetospheric Plasma
The temporal variability of magnetopause reconnection is an important aspect of solar wind magnetosphere coupling. Even under stable solar wind driving, reconnection can be triggered, modulated, or suppressed because of magnetic field and plasma conditions near the magnetopause boundary. We analyze a unique event in which a THEMIS satellite crosses the subsolar magnetopause three times within a 5 min interval in the presence of a cold‐ion population on the magnetospheric side of the boundary. During the first crossing, the satellite detects reconnection outflow and a D‐ shaped ion velocity distribution earthward from the boundary, indicating an active reconnection. The signatures disappear during the second crossing when the magnetospheric cold‐ion density increases significantly and reappear during the third crossing when the magnetospheric density drops to a level comparable to that of the first crossing. The solar wind and magnetosheath conditions do not change much during the interval. The magnetospheric population is evidently associated with a plasmaspheric plume with considerable variation in density. According to the theory of mass loading, the presence of such a plume population results in the local Alfvén speed at the second crossing being 40% smaller compared to the first and third crossings. However, the theory itself does not suggest suppression. We discuss possible suppression mechanisms considering the additional effects of the prevailing solar wind and local magnetopause conditions.
- 7 September 2025
- Published by American Geophysical Union (AGU) in Earth and Space Science
- Vol. 12 (9)
Open AccessSuppression of Magnetopause Reconnection in the Presence of Cold Magnetospheric Plasma
The temporal variability of magnetopause reconnection is an important aspect of solar wind magnetosphere coupling. Even under stable solar wind driving, reconnection can be triggered, modulated, or suppressed because of magnetic field and plasma conditions near the magnetopause boundary. We analyze a unique event in which a THEMIS satellite crosses the subsolar magnetopause three times within a 5 min interval in the presence of a cold‐ion population on the magnetospheric side of the boundary. During the first crossing, the satellite detects reconnection outflow and a D‐ shaped ion velocity distribution earthward from the boundary, indicating an active reconnection. The signatures disappear during the second crossing when the magnetospheric cold‐ion density increases significantly and reappear during the third crossing when the magnetospheric density drops to a level comparable to that of the first crossing. The solar wind and magnetosheath conditions do not change much during the interval. The magnetospheric population is evidently associated with a plasmaspheric plume with considerable variation in density. According to the theory of mass loading, the presence of such a plume population results in the local Alfvén speed at the second crossing being 40% smaller compared to the first and third crossings. However, the theory itself does not suggest suppression. We discuss possible suppression mechanisms considering the additional effects of the prevailing solar wind and local magnetopause conditions.
- Open Access
- 7 March 2025
- Published by American Geophysical Union (AGU) in Journal of Geophysical Research: Space Physics
- Vol. 130 (3)
Controlling Factors of the Local Time Extent of Magnetopause Reconnection: A Statistical Study
Magnetopause reconnection is the dominant mechanism for transporting solar wind energy and momentum into the magnetosphere‐ionosphere system. Magnetopause reconnection can occur along X‐lines of variable extent in the direction perpendicular to the reconnection plane. Identifying the spatial extent of X‐lines using satellite observations has critical limitations. However, we can infer the azimuthal extent of the X‐lines by probing the ionospheric signature of reconnection, the antisunward flow channels across the ionospheric Open‐Closed Field Line Boundary (OCB). We study 39 dayside magnetopause reconnection events using conjugate in situ and ionospheric observations to investigate the variability and controlling factors of the spatial extent of reconnection. We use spacecraft data from Time History of Events and Macroscale Interactions during Substorms (THEMIS) to identify in situ reconnection events. The width of the antisunward flow channels across the OCB is measured using the concurrent measurements from Super Dual Auroral Radar Network (SuperDARN). Also, the X‐line lengths are estimated by tracing the magnetic field lines from the ionospheric flow boundaries to the magnetopause. The solar wind driving conditions upstream of the bow shock are studied using solar wind monitors located at the L1 point. Results show that the magnetopause reconnection X‐lines can extend from a few Earth Radii (RE) to at least 22 RE in the GSM‐Y direction. Furthermore, the magnetopause reconnection tends to be spatially limited during high solar wind speed conditions.
- 7 March 2025
- Published by American Geophysical Union (AGU) in Journal of Geophysical Research: Space Physics
- Vol. 130 (3)
Open AccessControlling Factors of the Local Time Extent of Magnetopause Reconnection: A Statistical Study
Magnetopause reconnection is the dominant mechanism for transporting solar wind energy and momentum into the magnetosphere‐ionosphere system. Magnetopause reconnection can occur along X‐lines of variable extent in the direction perpendicular to the reconnection plane. Identifying the spatial extent of X‐lines using satellite observations has critical limitations. However, we can infer the azimuthal extent of the X‐lines by probing the ionospheric signature of reconnection, the antisunward flow channels across the ionospheric Open‐Closed Field Line Boundary (OCB). We study 39 dayside magnetopause reconnection events using conjugate in situ and ionospheric observations to investigate the variability and controlling factors of the spatial extent of reconnection. We use spacecraft data from Time History of Events and Macroscale Interactions during Substorms (THEMIS) to identify in situ reconnection events. The width of the antisunward flow channels across the OCB is measured using the concurrent measurements from Super Dual Auroral Radar Network (SuperDARN). Also, the X‐line lengths are estimated by tracing the magnetic field lines from the ionospheric flow boundaries to the magnetopause. The solar wind driving conditions upstream of the bow shock are studied using solar wind monitors located at the L1 point. Results show that the magnetopause reconnection X‐lines can extend from a few Earth Radii (RE) to at least 22 RE in the GSM‐Y direction. Furthermore, the magnetopause reconnection tends to be spatially limited during high solar wind speed conditions.
- 1 October 2019
- Published by Elsevier in Journal of Atmospheric and Solar-Terrestrial Physics
- Vol. 192, 104777
Spectral characteristic of geomagnetically induced current during geomagnetic storms by wavelet techniques
- 1 October 2019
- Published by Elsevier in Journal of Atmospheric and Solar-Terrestrial Physics
- Vol. 192, 104777
Spectral characteristic of geomagnetically induced current during geomagnetic storms by wavelet techniques
- Open Access
- 7 February 2019
- Published by American Geophysical Union (AGU) in Space Weather
- Vol. 17 (2) , 238-251
HILDCAA‐Related GIC and Possible Corrosion Hazard in Underground Pipelines: A Comparison Based on Wavelet Transform
High‐intensity long‐duration continuous auroral electrojet activities (HILDCAAs) are recognized as the continuous (duration >2 days) auroral electrojet activities (AE > 1,000 nT) caused by interplanetary Alfvén waves driving magnetic reconnection at the Earth's magnetopause. This paper focuses on the study of geoeffectiveness of HILDCAA on the basis of associated geomagnetically induced currents (GICs) and pipe‐to‐soil voltage (PSV) observed in an underground pipeline at Finnish Natural Gas Station, Mäntsälä. In this study, 113 HILDCAAs with different interplanetary sources (72 corotating interaction region‐storm preceded, 29 interplanetary coronal mass ejection‐storm preceded, and 12 nonstorm) have been selected and their possible contributions in underground pipeline corrosion have been quantitatively compared by assessing GIC/PSV profiles for each of these events. In addition, the spectral characteristics of GIC during these events have been studied using continuous wavelet transforms. The Morlet wavelet has been applied to 10‐s modeled GIC data corresponding to each event to explore the main band structures and the periodicities found in GIC spectrum. GIC/PSV modeling is based on plane wave method and distributed source transmission line analogy between pipelines and electric circuits. HILDCAAs are found to drive a small‐amplitude, but continuous, fluctuation in GIC throughout the event duration of several days. This makes the cumulative effect of HILDCAAs in pipeline corrosion noteworthy. The spectral analysis shows that GIC possesses both short‐term, as well as continuous, power distribution with different periodicities. CIR‐preceded HILDCAAs are found to be more geoeffective while taking associated GIC and PSV into account. Possible physical explanations supporting the results have been presented.
- 7 February 2019
- Published by American Geophysical Union (AGU) in Space Weather
- Vol. 17 (2) , 238-251
Open AccessHILDCAA‐Related GIC and Possible Corrosion Hazard in Underground Pipelines: A Comparison Based on Wavelet Transform
High‐intensity long‐duration continuous auroral electrojet activities (HILDCAAs) are recognized as the continuous (duration >2 days) auroral electrojet activities (AE > 1,000 nT) caused by interplanetary Alfvén waves driving magnetic reconnection at the Earth's magnetopause. This paper focuses on the study of geoeffectiveness of HILDCAA on the basis of associated geomagnetically induced currents (GICs) and pipe‐to‐soil voltage (PSV) observed in an underground pipeline at Finnish Natural Gas Station, Mäntsälä. In this study, 113 HILDCAAs with different interplanetary sources (72 corotating interaction region‐storm preceded, 29 interplanetary coronal mass ejection‐storm preceded, and 12 nonstorm) have been selected and their possible contributions in underground pipeline corrosion have been quantitatively compared by assessing GIC/PSV profiles for each of these events. In addition, the spectral characteristics of GIC during these events have been studied using continuous wavelet transforms. The Morlet wavelet has been applied to 10‐s modeled GIC data corresponding to each event to explore the main band structures and the periodicities found in GIC spectrum. GIC/PSV modeling is based on plane wave method and distributed source transmission line analogy between pipelines and electric circuits. HILDCAAs are found to drive a small‐amplitude, but continuous, fluctuation in GIC throughout the event duration of several days. This makes the cumulative effect of HILDCAAs in pipeline corrosion noteworthy. The spectral analysis shows that GIC possesses both short‐term, as well as continuous, power distribution with different periodicities. CIR‐preceded HILDCAAs are found to be more geoeffective while taking associated GIC and PSV into account. Possible physical explanations supporting the results have been presented.
- Open Access
- 7 September 2018
- Published by American Geophysical Union (AGU) in Earth and Space Science
- Vol. 5 (9) , 440-455
Field‐Aligned Current and Polar Cap Potential and Geomagnetic Disturbances: A Review of Cross‐Correlation Analysis
Polar cap potential (PCV) is an important parameter used for determining what kind of interaction takes place between solar wind and magnetosphere. Highly energetic particles from Sun driven by solar wind constantly bombard with Earth's magnetosphere–ionosphere system that results into a phenomenon like auroras, and major geomagnetic disturbances. Solar wind electron deposition determines the magnitude of field‐aligned current (FAC) and ultimately leads to PCV variation. Several studies found that increase in magnitude of IMF‐Bz causes an electric field of cross magnetosphere to increase, and it leads to increase in magnitude of ionospheric cross‐polar cap potential (PCV). Moreover, PCV was found to be a linear function of Vsw. In this research, we aim to study how field‐aligned current (FAC), for example, region 1 current and PCV, is related during different forms of geomagnetic disturbances. In all events, FAC and PCV are found to have corresponding fluctuations—especially at times of significant variation of IMF‐Bz (negative Bz interval) following the linearity of equation suggested by Moon in Moon (2012, https://doi.org/10.5140/JASS.2012.29.3.259). We found one‐to‐one correspondence between FAC and PCV. We did CWT analysis and found that FAC and PCV have more or less same spectral behaviors for each event considered. The cross‐correlation analysis shows a high and positive correlation between FAC and PCV at 0‐min time lag for all geomagnetic activity. The CWT analysis clearly supports the result of cross correlation between FAC and PCV. We found that FAC and Vsw, FAC‐B, and FAC and AE are also positively correlated with high‐correlation coefficient at lag 0 min for all geomagnetic storm. However, FAC‐Bz, FAC‐By, and FAC‐SYM (H) have varying correlation in different events. For a particular storm and substorm, the parameters Bz and By may not necessarily be varied with FAC in regular sequence but IMF (B) always show positive correlation with FAC for all geomagnetic activity. This paper presents a clear relation between FAC and PCV. This result will help to identify some of the outstanding issues in determining the causal mechanism of PCV variation, a crucial thing to understanding the coupling between the solar wind and M‐I system.
- 7 September 2018
- Published by American Geophysical Union (AGU) in Earth and Space Science
- Vol. 5 (9) , 440-455
Open AccessField‐Aligned Current and Polar Cap Potential and Geomagnetic Disturbances: A Review of Cross‐Correlation Analysis
Polar cap potential (PCV) is an important parameter used for determining what kind of interaction takes place between solar wind and magnetosphere. Highly energetic particles from Sun driven by solar wind constantly bombard with Earth's magnetosphere–ionosphere system that results into a phenomenon like auroras, and major geomagnetic disturbances. Solar wind electron deposition determines the magnitude of field‐aligned current (FAC) and ultimately leads to PCV variation. Several studies found that increase in magnitude of IMF‐Bz causes an electric field of cross magnetosphere to increase, and it leads to increase in magnitude of ionospheric cross‐polar cap potential (PCV). Moreover, PCV was found to be a linear function of Vsw. In this research, we aim to study how field‐aligned current (FAC), for example, region 1 current and PCV, is related during different forms of geomagnetic disturbances. In all events, FAC and PCV are found to have corresponding fluctuations—especially at times of significant variation of IMF‐Bz (negative Bz interval) following the linearity of equation suggested by Moon in Moon (2012, https://doi.org/10.5140/JASS.2012.29.3.259). We found one‐to‐one correspondence between FAC and PCV. We did CWT analysis and found that FAC and PCV have more or less same spectral behaviors for each event considered. The cross‐correlation analysis shows a high and positive correlation between FAC and PCV at 0‐min time lag for all geomagnetic activity. The CWT analysis clearly supports the result of cross correlation between FAC and PCV. We found that FAC and Vsw, FAC‐B, and FAC and AE are also positively correlated with high‐correlation coefficient at lag 0 min for all geomagnetic storm. However, FAC‐Bz, FAC‐By, and FAC‐SYM (H) have varying correlation in different events. For a particular storm and substorm, the parameters Bz and By may not necessarily be varied with FAC in regular sequence but IMF (B) always show positive correlation with FAC for all geomagnetic activity. This paper presents a clear relation between FAC and PCV. This result will help to identify some of the outstanding issues in determining the causal mechanism of PCV variation, a crucial thing to understanding the coupling between the solar wind and M‐I system.