本篇論文模擬電離層E層的電子濃度與因風切效應所引發的散塊E層發生機制。為了得知E層的電子濃度，我們利用E層中存在的五個離子(NO+、O2+、N2+、O+、Fe+)濃度加總來推算電子濃度。一開始單純只考慮經由光化學作用(產生作用與消失作用)產生的背景電子濃度分佈，則模擬的方法為：把隨高度、時間變化的離子產生率與離子消失率代入離子的連續方程式，以Runge-Kutta Method來解此連續方程式(一階微分方程式)，如此可得到各個離子濃度隨著高度與時間變化的分佈。
接著我們模擬潮汐風場造成散塊E層的發生現象。在E層，因離子碰撞頻率大於離子旋繞頻率，緯向潮汐風可藉由U×B drift導致離子在高度上有上下移動，當上層緯向風為西風、下層緯向風為東風，則會使離子在垂直漂移速度的輻合區域內累積，形成散塊E層。依上述想法，潮汐風造成的離子垂直漂移速度亦即為連續方程式中離子的傳輸項。接著把傳輸項代回連續方程式，以Crank-Nicholson Method來解此連續方程式(二維一階微分方程式)，如此可見因風切效應所造成的散塊E層發生現象。

In this study, we not only simulate height-time variations of background electron density of E region ionosphere but also study the formation of Sporadic E layers that is highly associated with the neutral wind shear. In order to model the electron density of E region ,we calculate the densities of NO+、O2+、N2+、O+、Fe+ that are assumed to consist in the ions in E region. The initial condition of the simulation is the photochemical equilibrium in E region, we calculate theoretical production rate and loss rates of the ions and substitude them into continuity equation to solve the first-order partial differential equation in accordance with Runge-Kutta method.
We then model the tidal wind that is responsible for the formation of Sporadic E layers. In the lower E region, the zonal neutral wind is the primary factor causing vertical ion drifts, which is the physical process of U×B dynamo in E region. Therefore, it requires that the zonal neutral wind with weatward in upper and eastward in the lower E region will result in ion convergence in northern hemisphere. The Sporadic E layers will occur at the region where the ion converges. With this configuration, the vertical ion drift velocity induced by tidal wind consititudes the transport term in continuity equation. Once the transport term is considered in the continuity equation, we use the Crank-Nicholson method to solve the resulting partial differential equation. The result indicates that a sporadic E layers indeed occurs at the regions where the ion are converged in the node of the neutral wind shear.

﹝1﹞ A. J. Carrasco, I.S. Batista, and M. A. Abdu, “Simulation of the sporadic E layer response to prereversal associated evening vertical electric field enhancement near dip equator”, Journal Of Geophysical Research, Vol. 112, 2007
﹝2﹞ K. Davies, Ionospheric Radio, 1989
﹝3﹞ Kelly, M. C., The Earth’s Ionosphere: Plasma Physics and Electrodynamics, 6pp., Academic, San Diego, 1989
﹝4﹞ R. W. Schunk, A. E. Nagy, Ionospheres, 356pp., Cambridge, 2000
﹝5﹞ Chapman, S., “The absorption and dissociative or ionizing effects of monochromatic radiation in an atmosphere on a rotating earth”, 2, Grazing incidence, Proc, Phys. Soc. London, 43,483, 1931.
﹝6﹞ F. L. Smith, Cody Smith, “Numerical evaluation of chapman’s grazing incidence integral ch(ηp, χp) ”,Journal of geophysical research,vol.77, 1973
﹝7﹞ Swider, W., Jr., The determination of the optical depth at large solar zenith distances, Planet, Space Sci., 12, 761, 1964
﹝8﹞ Rishbeth, H., and O. K. Gariott, Introduction to Ionospheric Physics, p.96, Academic, New Your, 1969
﹝9﹞ Banks P. M., and G. Kockarts, Aeronomy, Academic, San Diego, Calif., 1973
﹝10﹞ J. E. Titheridge, Model results for the ionospheric E region: solar and seasonal changes, Ann. Geophysicae 15, 63-78, 1996
﹝11﹞ Fesen, C. G., D. W. Rusch, and J.-C. Gerard, The latitude gradient of the NO peak density, J. Geophys. Res., 95, 19052-19059,1990.
﹝12﹞ Magdalena Helmer and John M.C. Plane, “A model of meteoric iron in the upper atmosphere. ”, Journal Of Geophysical Research, Vol.103,1998
﹝13﹞ A. D. Danilov and U. A. Kalgin, Seasonal and latitudinal variations of eddy diffusion coefficient in the mesosphere and lower thermosphere.,Journal of Atmospheric and Terrestrial Physics, Vol. 54 1992
﹝14﹞ Theodore L. Brown, The Chemistry of Metallic Elements in the Ionosphere and Mesosphere,1973
﹝15﹞ L. N. Carter , J. M. Forbes, Global transport and localized layering of metallic ions in the upper atmosphere. Ann. Geophysicae, 1999
﹝16﹞ D. M. Hunten, R. P. Turco, O. B. Toon, Smoke and Dust Particles of Meteoric Origin In the Mesosphere and Stratosphere. Atmosphere Science, 1980
﹝17﹞ T. J. Kane, C. S. Gardner, Structure and Seasonal Variability of the Nighttime Mesospheric Fe Layer at Midlatitudes,JGR,1993
﹝18﹞ Chen, W. M., and R. D. Harris, An ionospheric E-region nighttime model, J. Atmos. Terr. Phys., 1971
﹝19﹞ M. K. Rees, Physics and chemistry of the upper atmosphere, Cambridge, 1989
﹝20﹞ E. Turunen, H. Matveinen, J. Tolvanen and H. Ranta, D-region Ion Chemistry Model,
﹝21﹞ D. F. Strobel and T. R. Young, The Nighttime Ionosphere: E Region and Lower F Region
﹝22﹞ Abdu, M. A., J. Macdougall, I. S. Batista, J. H. A. Sobral, and P. T. Jayachandran, Equatorial evening prereversal electric field enhancement and sporadic E layer disruption: A Manifestation of E and F region coupling, J. Geophys. Res., 2003
﹝23﹞ Morton, Y. T., J. D. Mathews, and Q. Zhou, Further evidence for a 6-h tide above Arecibo, J. Atmos. Terr. Phys.,1993
﹝24﹞ Steven C. Chapra 著, 應用數值方法, 王晉中譯, 美商麥格羅-希爾國際股份有限中司, 2005