青藏高原那曲地区雷电特征初步分析

通过对2002年夏季青藏高原那曲地区雷暴过程及闪电观测资料的初步分析,发现该地区雷暴电荷结构具有多样性和复杂性,地闪明显偏少. 对高原地闪的一些基本特征参量的统计分析表明,无论正地闪还是负地闪梯级先导前都具有持续时间较长的云内放电过程,地闪以单次回击为主. 与中低纬度地区相比,高原地闪中正地闪比例明显要高,为33髎;负地闪为67髎;正、负地闪回击后常常伴随短时间的连续电流.     

Spatiotemporal characteristics of positive cloud-to-ground lightning discharges and bidirectional leader of the lightning

The negative CG lightning discharges neutralizing negative charges in cloud usually dominate for most of thunderstorms. However, a lot of positive CG light-ning discharges often occur in the disappearing stage of thunderstorms, in the stratiform region of mesoscale convective systems and some supercells producing hail and tornado. Because the positive CG lightning discharges produce larger current of the return stroke and neutralize more charges due to the continuing currents with longer las…     

Ten-year study of cloud-to-ground lightning activity in the Iberian Peninsula

The characteristics (annual and diurnal cycle, polarity, multiplicity and first stroke peak current) of ∼4.3×10⁶ cloud-to-ground flashes recorded in the Iberian Peninsula during the first decade of measurements of the lightning detection network installed in Spain are analyzed. The mean monthly variation shows maximum lightning activity between May and September, while minimum values are observed in January and February. The mean diurnal cycle shows maximum values at 1700 LT and minimum values at around 1000 LT. The average maximum flash density (not corrected for detection efficiency) is 2.1 flashes km⁻² year⁻¹. Maximum lightning activity is associated with mountainous areas. The effect of the Mediterranean Sea is also seen. The percentage of positive flashes is 9%, although this changes over the year from 6.5% in June to 22.6% in January. The average multiplicity is found to be 2.0 for the negative flashes and 1.1 for the positive flashes, and the percentages of single-stroke flashes are 53.6% and 89%, respectively. The monthly distribution of multiplicity for negative flashes peaks in the summer and minimum values are found in the winter. The multiplicity of the positive flashes does not seem to be function of the month. The median (mean) first stroke peak currents are found to be 23.5 kA (27.3 kA) for the negative flashes and 35.3 kA (47.1 kA) for the positive flashes. For both polarities, the peak current is higher in the summer than in the winter. The percentage of positive flashes and the mean peak currents for negative flashes are higher over the sea areas than over land.     

Some observations of bipolar flashes during summer thunderstorms near Warsaw

Lightning discharges monitored by the SAFIR network system in Poland have been additionally identified over the 100×100 km area near Warsaw by single-point independent recordings of electric field and Maxwell current rapid changes. The data collected in summer thunderstorm days of 2002 showed some untypical properties of the lightning discharges which are rarely observed. Especially remarkable was a number of ground multi-stroke flashes with the return strokes (RS) which transported to the earth charges of opposite signs. Bipolar flashes (BF) of this kind were mostly involved in the events in which the nearby intracloud (ic) and cloud-to-ground (c-g) discharges were very closely associated in time. Events of such a close collocation of two different types of lightning discharges, previously called the complex lightning discharge events (CLDE), were quite often observed during summer thunderstorms in Poland. The events of this kind, i.e. 8 flashes, identified by the SAFIR detection system as BF’s present the multiple stroke flashes of the mean horizontal separation distance between striking points of particular RS equal to (2.8 ± 2.1) km and of the mean time interval between strokes of (46.8 ± 74.4) ms. The time separation between the observed BF and the adjacent ic flashes was from 0.1 to 335 ms, and horizontal separation distance between them ranged from 1.8 to 14.5 km. The multiplicity of the recorded BF’s ranged from 2 to 4 strokes. Four of these BF’s followed the ic discharge, but the other three preceded the ic and one was alone with no close ic.     

不同电结构雷暴中的地闪特征

对观测资料的分析表明,甘肃地区和海南岛的地闪特征具有很大差异。海南岛负地闪（P型）首次回击前的电场变化波形与国外的报道类似,而甘肃地区仅有30％左右的负地闪属P型,约70％的负地闪首次回击前的电场变化波形具有明显的云闪特征（C型负地闪）。海南岛没有发现正地闪,甘肃地区有正地闪。该两地区存在两类电结构不同的雷暴,两类雷暴中地闪放电过程及特征的很大差异,说明我国的防雷规范及措施很可能需要考虑雷暴的地区差异性。     

Low-frequency E-field Detection Array(LFEDA)—Construction and preliminary results

In recent years, locating total lightning at the VLF/LF band has become one of the most important directions in lightning detection. The Low-frequency E-field Detection Array(LFEDA) consisting of nine fast antennas was developed by the Chinese Academy of Meteorological Sciences in Guangzhou between 2014 and 2015. This paper documents the composition of the LFEDA and a lightning-locating algorithm that applies to the low-frequency electric field radiated by lightning pulse discharge events(LPDEs). Theoretical simulation and objective assessment of the accuracy and detection efficiency of LFEDA have been done using Monte Carlo simulation and artificial triggered lightning experiment, respectively. The former results show that having a station in the network with a comparatively long baseline improves both the horizontal location accuracy in the direction perpendicular to the baseline and the vertical location accuracy along the baseline. The latter results show that detection efficiencies for triggered lightning flashes and return strokes are 100% and 95%, respectively. The average planar location error for return strokes of triggered lightning flashes is 102 m. By locating LPDEs in thunderstorms, we find that LPDEs are consistent with convective regions as indicated by strong reflectivity columns, and present a reasonable distribution in the vertical direction.In addition, the LFEDA can reveal an image of lightning development through mapping the channels of lightning. Based on three-dimensional locations, the vertical propagation speed of the preliminary breakdown and the changing trend of the leader's speed in an intra-cloud and a cloud-to-ground flash are investigated. The research results show that the LFEDA has the capability for three-dimensional location of lightning, which provides a new technique for researching lightning development characteristics and thunderstorm electricity.     

On the NOx production by laboratory electrical discharges and lightning

Different approaches are used in estimating the global production of NO_x by lightning flashes, including field measurements carried out during thunderstorm conditions, theoretical studies combining the physics and chemistry of the electrical discharges, and measurements of NO_x yield in laboratory sparks with subsequent extrapolation to lightning. In the latter procedure, laboratory data are extrapolated to lightning using the energy as the scaling quantity. Further, in these studies only the return strokes are considered assuming that contributions from other processes such as leaders, continuing currents, M components, and K processes are negligible. In this paper, we argue that the use of energy as the scaling quantity and omission of all lightning processes other than return strokes are not justified. First, a theory which can be used to evaluate the NO_x production by electrical discharges, if the current flowing in the discharge is known, is presented. The results obtained from theory are compared with the available experimental data and a reasonable agreement is found. Numerical experiments suggest that the NO_x production efficiency of electrical discharges depends not only on the energy dissipated in the discharge, but also on the shape of current waveform. Thus, the current signature, can influence extrapolation of laboratory data to lightning flashes. Second, an estimation of the NO_x yield per lightning flash is made by treating the lightning flash as a composite event consisting of several discharge processes. We show that the NO_x production takes place mainly in slow discharge processes such as leaders, M components, and continuing currents, with return strokes contributing only a small fraction of the total NO_x. The results also show that cloud flashes are as efficient as ground flashes in NO_x generation. In estimating the global NO_x production by lightning flashes the most influencing parameter is the length of the lightning discharge channel inside the cloud. For the total length of channels inside the cloud of a typical ground flash of about 45 km, we estimate that the global annual production of NO_x is about 4 Tg(N).     

Observation and study on the whole process of cloud-to-ground lightning using narrowband radio interferometer

A narrowband radio interferometer has been developed and used to locate the entire sources of VHF radiations from a negative cloud-to-ground (CG) lightning discharge which contains 19 strokes. This system uses five antennas to form an array consisting of short- and long-baselines along two or- thogonal directions. The system error which comes from frequency conversion is reduced by phase detection through direct high frequency amplifying. An interactive graphic analysis procedure is used to remove the fringe ambiguities which exist inherently in interferometry and to determine the direction of lightning radiation sources in two dimensions (azimuth and elevation) as a function of time at a time resolution of microsecond orders. With the developed system, the whole progression process in time and space of a lightning flash can be reconstructed. In this paper, combining the synchronous data of electric filed change and VHF radiation, the whole processes of an example negative CG flash have been studied in detail. It is found that the preliminary breakdown event of the CG flash started from negative charge region and exhibited firstly a downward pregression and then an upward propagation. There were very intense and continuous radiations during stepped leaders which became much stronger when the first return stroke began. In contrast, there were less and only discrete radiations during dart leaders. Stepped leader and dart leader may transform to each other depending on the state of the ionization of the path. The progression speed of initial stepped leaders was about 105 ms?1, while that was about 4.1×106 and 6.0×106 ms?1 for dart leaders and dart-stepped leaders, respectively. M events produced hook-shaped field changes accompanied by active burst of radiations at their begin- nings. Followed these active radiation processes, M events appeared to contact finally into conducting main discharge channels. The mean progression speed of M events was about 7×107 ms?1, greater than that of the dart leaders and dart-step leaders. K events and attempted leaders were essentially the same as dart leaders except that they could not reach the ground and initiate return strokes.     

Development characteristics of cloud-to-ground lightning with multiple grounding points

Using the optical images of a cloud-to-ground lightning flash with multiple grounding points obtained by a highspeed video system in the Qinghai Province of China along with synchronous radiated electric field information, the propagation characteristic and the electric field change features of the leaders and the grounding behavior of discharge channels are analyzed.In addition, the two-dimensional velocity of the leader was estimated and its correlation with the time interval of the corresponding subsequent return stroke, and that with the peak current of return stroke are investigated. The results show that the average distance between the three obvious grounded points of the first return stroke channel is about 512.7 m, and the average time interval between the pulses of the corresponding electric field fast changes is 3.8 μs. Further, the average time interval between electric field pulses from the stepped leader is smaller than that of normal single grounding lightning. The observed lightning in our study has two main channels, namely the left and right channels. Based on our observations, it is clear that the dart leader comes close to the ground in case of the left channel after the first return stroke, but it fails to form a return stroke.However, the right channel exhibits a relatively rare phenomenon in that the subsequent return stroke R2 occurred about 2.1 ms after the dart leader arrived at the ground, which was unusually long; this phenomenon might be attributed to the strong discharge of the first return stroke and insufficient charge accumulation near the grounded point in a timely manner. The two-dimensional velocities for the stepped leader of the two main channels are about 1.23×105 and 1.16×105 m s-1, respectively. A sub-branch of stepped leader for the left channel fails to reach the ground and develops into an attempt leader eventually; this might be attributed to the fact that the main branch connects considerably many sub-branches, which leads to the instantaneous decline of the potential difference between the sub-branch and ground. Furthermore, it might also be because the propagation direction of this sub-branch is almost perpendicular to the atmospheric electric field direction, which is not conducive to charge transfer. The two-dimensional velocities for the dart leaders of five subsequent return strokes are all in the normal range, and they positively correlate with the peak current of the subsequent return stroke.     

Field-aligned currents in the winter and summer hemispheres caused by IMF <Emphasis Type="Italic">By</Emphasis>

Structures controlled by the IMF By sign and season of the year have been detected based on the decomposition of field-aligned current maps constructed using magnetic field measurements on polar low-orbiting satellites. It has been indicated that field-aligned currents have identical structures, composed of the main polar circular current and the return current at the polar cap dayside boundary, at any By sign in the summer hemisphere. Two different types of structures are implemented under winter conditions depending on the By sign. For the northern winter, it is the polar circular current and the return current at the polar cap nightside boundary at By < 0; current sheets are strongly stretched along latitudes below 80° MLat, and only small part of the current is in the noon sector of the polar cap. For the summer winter, the corresponding structures are implemented at opposite By signs. The intensities of the field-aligned currents, originating as a result of the interhemispheric asymmetry and flowing along closed geomagnetic field lines near the polar cap boundary, have been estimated. The maximum of the interhemispheric current density is 0.25 μA m⁻² in the summer and 0.1 μA m⁻² in the winter; the total current is 5 × 10⁵ and 5 × 10⁴ A, respectively.     

Charge analysis on lightning discharges to the ground in Chinese inland plateau (close to Tibet)

Since the summer of 1996, scientists from China and Japan have conducted a joint observation of natural cloud-to-ground lightning discharges in the Zhongchuan area that is located close to Qinghai-Xizang (Tibet) Plateau, China. It has been found that the long-duration of intracloud discharge processes, just before the first return stroke, lasted more than 120 ms for 85% of cloud-to-ground flashes in this area, with a mean duration of 189.7 ms and a maximum of 300 ms. We present the results of charge sources neutralized by four ground flashes and two intracloud discharge processes, just before the first return stroke, by using the data from a 5-site slow antenna network synchronized by GPS with 1 s time resolution. The result shows that the altitudes of the neutralized negative charge for three negative ground flashes were between 2.7 to 5.4 km above the ground, while that of neutralized positive charges for one positive ground flash and one continuing current process were at about 2.0 km above the ground. The comparison with radar echo showed that the negative discharges initiated in the region greater than 20 dBZ or near the edge of the region with intense echoes greater than 40 dBZ, while positive discharge initiated in the weak echo region.     

Belowground carbon balance and carbon accumulation rate in the successional series of monsoon evergreen broad-leaved forest

The belowground part of terrestrial ecosystem is a huge carbon pool. It is believed that of the total 2500Gt carbon stored in global terrestrial ecosystem, soil carbon storage within the 1 m surface layer ac- counts for 2000Gt, which is 4-fold of vegetation car- bon storage[1,2]. Compared with the carbon in the vegetation, carbon in the deep soil layers is much more stable, and it will stay in soil profile permanentlyunless geological vicissitude occurs. Essentially, forest restoration is the…     

SO2 Emission from Active Volcanoes Measured Simultaneously by COSPEC and mini-DOAS

We measured SO₂ emission rate from six volcanoes in Latin America (Santa Ana, El Salvador; San Cristóbal and Masaya, Nicaragua; Arenal and Poás, Costa Rica; Tungurahua and Sierra Negra, Ecuador) and from Mt. Etna, Italy, using two different remote sensing techniques: COSPEC (COrrelation SPECtrometer) and miniDOAS (miniaturized Differential Optical Absorption Spectroscopy). One of the goals of this study was to evaluate the differences in SO₂ emission rates obtained by these two methods. The observed average SO₂ emission rates measured during this study were 2688 t·d⁻¹ from Tungurahua in July 2006, 2375 t·d⁻¹ in September 2005 and 480 t·d⁻¹ in February 2006 from Santa Ana, 1200 t·d⁻¹ in May 2005 from Etna, 955 t·d⁻¹ in March 2006 and 1165 t·d⁻¹ in December 2006 from Masaya, 5400 t·d⁻¹ of March 7, 2006 and 265 t·d⁻¹ in March 2006 from San Cristobal, 113 t·d⁻¹ in April 2006 from Arenal, 104 t·d⁻¹ in April 2006 from Poás and 11 t·d⁻¹ in July 2006 from Sierra Negra volcano. Most of the observed relative differences of SO₂ emission measurements from COSPEC and miniDOAS were lower than 10%.     

Uncertainty analysis of CO<Subscript>2</Subscript> flux components in subtropical evergreen coniferous plantation

We present an uncertainty analysis of ecological process parameters and CO₂ flux components (R _eco, NEE and gross ecosystem exchange (GEE)) derived from 3 years’ continuous eddy covariance measurements of CO₂ fluxes at subtropical evergreen coniferous plantation, Qianyanzhou of ChinaFlux. Daily-differencing approach was used to analyze the random error of CO₂ fluxes measurements and bootstrapping method was used to quantify the uncertainties of three CO₂ flux components. In addition, we evaluated different models and optimization methods in influencing estimation of key parameters and CO₂ flux components. The results show that: (1) Random flux error more closely follows a double-exponential (Laplace), rather than a normal (Gaussian) distribution. (2) Different optimization methods result in different estimates of model parameters. Uncertainties of parameters estimated by the maximum likelihood estimation (MLE) are lower than those derived from ordinary least square method (OLS). (3) The differences between simulated R_eco, NEE and GEE derived from MLE and those derived from OLS are 12.18% (176 g C·m⁻²·a⁻¹), 34.33% (79 g C·m⁻²·a⁻¹) and 5.4% (92 g C·m⁻²·a⁻¹). However, for a given parameter optimization method, a temperature-dependent model (T_model) and the models derived from a temperature and water-dependent model (TW_model) are 1.31% (17.8 g C·m⁻²·a⁻¹), 2.1% (5.7 g C·m⁻²·a⁻¹), and 0.26% (4.3 g C·m⁻²·a⁻¹), respectively, which suggested that the optimization methods are more important than the ecological models in influencing uncertainty in estimated carbon fluxes. (4) The relative uncertainty of CO₂ flux derived from OLS is higher than that from MLE, and the uncertainty is related to timescale, that is, the larger the timescale, the smaller the uncertainty. The relative uncertainties of R_eco, NEE and GEE are 4%−8%, 7%−22% and 2%−4% respectively at annual timescale. Supported by the National Natural Science Foundation of China (Grant No. 30570347), Innovative Research International Partnership Project of the Chinese Academy of Sciences (Grant No. CXTD-Z2005-1) and National Basic Research Program of China (Grant No. 2002CB412502)     

Study of charge structure and radiation characteristic of intracloud discharge in thunderstorms of Qinghai-Tibet Plateau

The comprehensive observations on lightning discharges were conducted in Naqu area of Qinghai-Tibet Plateau in summer of 2002. The electric structures of thunderstorms and the characteristics of lightning discharges at initial stage were analyzed by using the observation data. The results show that most of intracloud (IC) lightning flashes were polarities inverted in thunderstorms with tripole electric charge structure and occurred between negative charge region located in the middle of the thunderstorm and positive charge region located at the bottom of the thunderstorm. The radiation characteristics of discharge processes in cloud with longer lasting time involved in Cloud-to-Ground (CG) lightning flashes were similar to that of IC discharges.A lot of radiation pulses were produced in these discharge processes. Because the IC discharges took place at the bottom of thundercloud and were near the ground, they may produce more serious damage to equipment on the ground therefore should not be neglected in lightning protection.     

Observation and modeling of NPP for Pinus elliottii plantation in subtropical China

Based on the stem analysis of 59 individuals of Pinus elliottii in combination with tree biomass models, we calculated annual biomass increment of forest plots at Qianyanzhou Ecological Station, Chinese Academy of Sciences in subtropical China. In addition, canopy layer and community NPP were calculated based on 12 years’ litter fall data. NPP of the 21-year-old forest was estimated by using the BIOME BGC model; and both measured NPP and estimated NPP were compared with flux data. Community biomass was 10574 g · m⁻²; its distribution patterns in tree layer, shrub layer, herbaceous layer, tree root, herbaceous and shrub roots and fine roots were 7542, 480, 239, 1810, 230, 274 and 239 g · m⁻², respectively. From 1999 to 2004, the average annual growth rate and litter fall were 741 g · m⁻² · a⁻¹ (381.31 gC · m⁻² · a⁻¹) and 849 g · m⁻² · a⁻¹ (463 gC · m⁻² · a⁻¹), respectively. There was a significant correlation between annual litter fall and annual biomass increment; and the litter fall was 1.19 times the biomass increment of living trees. From 1985 to 2005, average NPP and GPP values based on BGC modeling were 630.88 (343.31–906.42 gC · m⁻² · a⁻¹) and 1 800 gC · m⁻² · a⁻¹ (1351.62–2318.26 gC · m⁻² · a⁻¹). Regression analysis showed a linear relationship (R ²=0.48) between the measured and simulated tree layer NPP values. NPP accounted for 30.2% (25.6%–32.9%) of GPP, while NEP accounted for 57.5% (48.1%–66.5%) of tree-layer NPP and 41.74% (37%–52%) of stand NPP. Soil respiration accounted for 77.0% of measured tree NPP and 55.9% of the measured stand NPP. NEE based on eddy covariance method was 12.97% higher than the observed NEP. Supported by the National Key Basic Research Special Foundation of China (Grant No. 2002CB4125), International Joint Research Project under Ministry of Science and Technology of China (Grant No. 2006DFB91920)     

Study of charge structure and radiation characteristic of intracloud discharge in thunderstorms of Qinghai-Tibet Plateau

The comprehensive observations on lightning discharges were conducted in Naqu area of Qinghai-Tibet Plateau in summer of 2002. The electric structures of thunderstorms and the characteristics of lightning discharges at initial stage were analyzed by using the observation data. The results show that most of intracloud (IC) lightning flashes were polarities inverted in thunderstorms with tripole electric charge structure and occurred between negative charge region located in the middle of the thunderstorm and positive charge region located at the bottom of the thunderstorm. The radiation characteristics of discharge processes in cloud with longer lasting time involved in Cloud-to-Ground (CG) lightning flashes were similar to that of IC discharges. A lot of radiation pulses were produced in these discharge processes. Because the IC discharges took place at the bottom of thundercloud and were near the ground, they may produce more serious damage to equipment on the ground therefore should not be neglected in lightning protection.     

He,Ne and Ar isotopic composition of Fe-Mn crusts from the western and central Pacific Ocean and implications for their genesis

The noble gas nuclide abundances and isotopic ratios of the upmost layer of Fe-Mn crusts from the western and central Pacific Ocean have been determined. The results indicate that the He and Ar nu- clide abundances and isotopic ratios can be classified into two types: low 3He/4He type and high 3He/4He type. The low 3He/4He type is characterized by high 4He abundances of 191×10-9 cm3·STP·g-1 on average, with variable 4He, 20Ne and 40Ar abundances in the range (42.8―421)×10-9 cm3·STP·g-1, (5.40―141)×10-9 cm3·STP·g-1, and (773―10976)×10-9 cm3·STP·g-1, respectively. The high 3He/4He samples are characterized by low 4He abundances of 11.7×10-9 cm3·STP·g-1 on average, with 4He, 20Ne and 40Ar abundances in the range of (7.57―17.4)×10-9 cm3·STP·g-1, (10.4―25.5)×10-9 cm3·STP·g-1 and (5354―9050)×10-9 cm3·STP·g-1, respectively. The low 3He/4He samples have 3He/4He ratios (with R/RA ratios of 2.04―2.92) which are lower than those of MORB (R/RA=8±1) and 40Ar/36Ar ratios (447―543) which are higher than those of air (295.5). The high 3He/4He samples have 3He/4He ratios (with R/RA ratios of 10.4―12.0) slightly higher than those of MORB (R/RA=8±1) and 40Ar/36Ar ratios (293―299) very similar to those of air (295.5). The Ne isotopic ratios (20Ne/22Ne and 21Ne/22Ne ratios of 10.3―10.9 and 0.02774―0.03039, respectively) and the 38Ar/36Ar ratios (0.1886―0.1963) have narrow ranges which are very similar to those of air (the 20Ne/22Ne, 21Ne/22Ne, 38Ar/36Ar ratios of 9.80, 0.029 and 0.187, respectively), and cannot be differentiated into different groups. The noble gas nuclide abundances and isotopic ratios, together with their regional variability, suggest that the noble gases in the Fe-Mn crusts originate primarily from the lower mantle. The low 3He/4He type and high 3He/4He type samples have noble gas characteristics similar to those of HIMU (High U/Pb Mantle)- and EM (Enriched Mantle)-type mantle material, respectively. The low 3He/4He type samples with HIMU-type noble gas isotopic ratios occur in the Magellan Seamounts, Marcus-Wake Seamounts, Marshall Island Chain and the Mid-Pacific Sea- mounts whereas the high 3He/4He type samples with EM-type noble gas isotopic ratios occur in the Line Island Chain. This difference in noble gas characteristics of these crust types implies that the MagellanSeamounts, Marcus-Wake Seamounts, Marshall Is- land Chain, and the Mid-Pacific Seamounts originated from HIMU-type lower mantle material whereas the Line Island Chain originated from EM-type lower mantle material. This finding is consistent with varia- tions in the Pb-isotope and trace element signatures in the seamount lavas. Differences in the mantlesource may therefore be responsible for variations in the noble gas abundances and isotopic ratios in the Fe-Mn crusts. Mantle degassing appears to be the principal factor controlling noble gas isotopic abundances in Fe-Mn crusts. Decay of radioactive isotopes has a negligible influence on the nuclide abundances and isotopic ratios of noble gases in these crusts on the timescale of their formation.     

Disturbance of the ionospheric <Emphasis Type="Italic">D</Emphasis> region by the electric current of the atmospheric-ionospheric electric circuit

The mechanism by which electron and ion densities change in the ionospheric D region due to the electric current flowing in the atmospheric-ionospheric electric circuit is studied. The current disturbance in this circuit exists over the regions of increased seismic, meteorological, and thunderstorm activity. In the framework of the model, the influence of the electron and ion transportation under the action of the electric field on the formation of a disturbance in the D region and heating of the plasma electron component by the field are considered. The calculation results show that the densities of electrons and ions can change by an order of magnitude at an increase in the current density up to ∼(10⁻⁹–10⁻⁸) A m⁻², the sign of the disturbance depending on the current direction.     

Numerical simulation of the effect of lower positive charge region in thunderstorms on different types of lightning

Combined with the existing stochastic lightning parameterization scheme, a classic tripole charge structure in thunderstorms is assumed in the paper, and then 2-dimensional fine-resolution lighting discharge simulations are performed to quantitatively investigate the effect of lower positive charge (LPC) on different types of lightning. The results show: (1) The LPC plays a key role in generating negative cloud-to-ground (CG) flashes and inverted intra-cloud (IC) lightning, and with the increase of charge density or distribution range of LPC region, lightning type changes from positive polarity IC lightning to negative CG flashes and then to inverted IC lightning. (2) Relative to distribution range of charge regions, the magnitude of charge density of the LPC region plays a dominant role in lightning type. Only when the maximal charge density value of LPC region is within a certain range, can negative CG flashes occur, and the occurrence probability is relatively fixed. (3) In this range, the charge density and distribution range of LPC region jointly determine the occurrence of negative CG flashes, which has a linear boundary with the trigger condition of IC lightning. (4) The common effect of charge density and distribution range of the LPC region is to change the distribution of positive potential well of bottom part of thunderstorms, and inverted IC lightning occurs when the initial reference potential is close to 0 MV, and negative CG flashes occur when the initial reference potential is far less than 0 MV.