Singh 8.0_REVFC_Layout 6 ANNALS OF GEOPHYSICS, 57, 6, 2014, A0657; doi:10.4401/ag-6596 A0657 Shift of effective lightning areas during pre to post period of solar cycle minimum of 2008-2009 as determined from Schumann resonance studies at Agra, India Birbal Singh*, Devbrat Pundhir RBS Engineering Technical Campus, Department of Electronics and Communication Engineering, Bichpuri, Agra, India ABSTRACT Employing a set of 3-component search coil magnetometer, Schumann resonance studies have been in progress at Agra (Geograph. lat. 27.2°N, long. 78°E), India since 01 April, 2007. We have analysed the data for two periods; first from 01 April, 2007 to 31 March, 2008 (period-I), and then from 01 March, 2011 to 29 February, 2012 (period-II) which cor- respond to pre and post periods of solar cycle minimum of 2008-2009. From the diurnal variation of first mode intensity and frequency, we study the seasonal variations of global thunderstorm activity, effective source distance and level of lightning during both the periods. We show that world thunderstorm activity shifts to summer in the northern hemisphere as the effective source distance approaches close to the ob- server, and the level of intense lightning shifts from the month of July, 2007 in period-I to August, 2011 in period-II. This is supported by Lightning Imaging Sensor (LIS) satellite data also. A possible expla- nation in terms of increasing solar activity is suggested. 1. Introduction The phenomenon of Schumann resonance (SR) oc- curring in the earth-ionosphere cavity as a result of reso- nance between direct and round the globe propagation of extremely low frequency (ELF) waves radiated from lightning discharges was first predicted by Schumann [1952] and experimentally verified by Balser and Wagner [1960]. The SR appears in the form of standing waves at the frequencies of 8, 14, 20, 26... Hz which are known as SR modes. Owing to its widespread applications in the studies of global thunderstorm activities, ground surface temperature, lower region of the ionosphere, and fore- cast of monsoon etc. extensive research work has been done in this field and a detailed description of the early work is presented in an excellent monopgraph by Nick- olaenko and Hayakawa [2002]. Some very interesting re- cent results which have come out of the extensive mor- phological and varying geophysical condition studies have been presented by many workers. For example, Price and Melnikov [2004] have studied the diurnal, sea- sonal, and inter-annual variations in the SR phenomenon and have found three dominant maxima in the diurnal cycle related to lightning activity in south-east Asia, Africa, and South America. Further, the largest global lightning activity occurs during the northern hemisphere summer ( JJA), with the southern hemisphere summer (DJF) having the least lightning around the globe. Green- berg and Price [2007] have compared the ELF data on the ground with the Optical Transient Detector (OTD) data on satellite and found discrepancies in lightning activities between the Asian and South American centers. To un- derstand the relative strength of the different source re- gions they calculated the attenuation rate for the electromagnetic waves in SR band in both the North- South and East-West direction reaching their station. They found that the simulated magnetic field relative am- plitude based on OTD observations showed a peak at 2000 UT corresponding to American source 40% larger than African and Asian sources, whereas in the East-West direction the African peak at 1400 UT was 80% larger than the American peak and Asian peak was not visible. Hayakawa et al. [2005] and Nickolaenko et al. [2006] have reported anomalies in SR bands as they found increase in the intensity of fourth harmonic due to earthquake which were interpreted in terms of superposition of di- rect ELF signals from a distant thunderstorm source in America with those scattered from a conducting distur- bance in atmosphere over the epicenter of earthquake. Williams and Satori [2007] have studied the solar radia- Article history Received May 30, 2014; accepted December 26, 2014. Subject classification: Schumann resonance, Thunderstorm activity, Effective lightning areas. tion-induced changes in ionospheric height and SR wave- guide on different time scales. They find that, in general, order of magnitude changes in radiation are needed to cause relative changes in ionospheric height as large as 10% as in the case on both the diurnal and 11-year time scales. On day-night asymmetry of the earth-ionosphere cavity and solar terminator effects Nickolaenko and Hayakawa [2002] concluded theoretically that the ampli- tude variations to be expected at the sunrise/sunset ter- mination line are practically undetectable. Figure 7 of Nickolaenko et al [2006] demonstrates that the lens-like “terminator” pattern is pertinent to the global thunder- storm activity itself. This is why a uniform earth-ionos- phere cavity provides the “terminator effect” similar to observations (in the UT) at an orbitary observatory. The confusion with the ionosphere local height occurred due to closeness of the start-stop times of the planetary light- ning activity with the sunrise and sunset moments at the middle East and central Europe. Similar inference was drawn also by Penchony et al. [2007] who observed that the SR field variations are governed primarily by the vari- ations in source intensity and source-receiver geometry and that the effect of the day-night asymmetry in the ion- osphere is secondary. In contrast to the above, Melnikov et al. [2004] presented observational evidence for the sun- rise/sunset terminator effect based on multistation SR observations with hourly time resolution and found steep increase/decrease of SR amplitude during terminator times. A similar result was found also by Ondraskova et al. [2007] at Modra observatory, Slovakia for the first four modes of electric field component. The later results were supported also by the observations of Satori et al. [2007a] who found jump-like increase of SR amplitude between the local ionosphere and surface sunrise times and sharp decreases between the local surface and ionospheric sun- set times. The SR phenomenon is due to contribution from both the global and local thunderstorm activities. In order to separate the local contribution from the global one Sentman and Fraser [1991] have suggested a tech- nique which has been applied by many workers [Pechony and Price 2006, Nickolaenko and Hayakawa 2008, Nick- olaenko et al 2011]. It is shown that local modulation function depends primarily on source-receiver distance geometry and it can be separated with the technique sug- gested by Sentman and Fraser [1991]. Ondraskova et al. [2009, 2011] have studied the decrease of SR frequencies and changes in effective lightning areas toward the solar cycle minimum of 2008-2009 using the data from Modra observatory Slovakia. The extraordinary fall in the first mode frequency was interpreted in the light of the sug- gestion made by Satori et al [2005] that the reduction in X-ray radiation causing decrease in the ionospheric con- ductivity may be attributed responsible for the observed effect. They further calculated the effective thunderstorm areas from the monthly mean diurnal frequency range (DFR) of the electric field component. They found that the difference in the northern and southern hemisphere summer areas not only declined with the decreases in solar activity but almost vanished during the deep solar minimum of 2008-2009. Semi-annual variation in the areas dominates in the years of the deep solar minimum. In our recent paper [Tyagi et al. 2013, hereafter paper I] we have studied the characteristics of global thunder- storm activities extracted from SR data for the period of 01 April, 2007-31 March, 2008. In the present paper we extend some of the studies made earlier in paper I espe- cially those related to variation of global thunderstorm activities and compare the results. Then, we consider the first mode frequency and intensity of the magnetic field data to study the effective source distance and level of lightning for both the periods 01 April, 2007-31 March, 2008 (period-I) and 01 March, 2011-29 February, 2012 (pe- riod-II) which are pre and post periods of solar cycle min- imum of 2008-2009 and compare the results. 2. Experimental set up and method of data processing The details of experimental set up and method of data processing are not available in paper I. However, these details are available in our another paper [Singh et al. 2014]. Here, we mention the details briefly as fol- lows;We are using a 3-component search coil magne- tometer (LEMI-30) which has been procured from Lviv center of Institute of Space Research, Ukraine. The search coils are buried under-ground in a relatively noise free area in Agriculture fields of the Bichpuri campus of R.B.S. college, Agra which is located 12 Km west of Agra city in rural area. The three sensors are oriented in geographical North-South (X-component), East-West (Y-component), and vertical (Z-component) directions and each of them has following specifications; The LEMI-30 system has sampling rate of 256 Hz and sends these samples to dedicated PC. The LEMI-30i soft- SINGH ET AL. 2 Dimension : 870 mm (length), 85 mm (diameter) Frequency range : 0.001-30 Hz Measuring range : ± 200 nT Transformation factor 0.001-1Hz : 20×f mV/nT 1-30 Hz : 20 mV/nT Auxiliary output gain : 20 dB Magnetic noise level at 0.01-10 Hz : ≤ 20pT/Hz1/2-≤0.04 pT/Hz1/2 Mains interference rejection : > 60 dB 3 ware in the PC takes average of each four samples sim- ply by summing and dividing so that a binary file con- taining data at the sampling rate of 64 Hz is stored in the PC. The recorded data on PC in amplitude-time may be seen in frequency-time (dynamic spectra) by performing spectral analysis using FFT available in MATLAB with 1024 words of data length (temporal resolution=16 sec, frequency resolution=0.06 Hz) at a time. The power spectral density (PSD) of the input signal is evaluated using Welch spectral technique [Welch 1967] which uses averaged modified periodograms. The PSD are prepared for each one hour data (230400 data points) and a Ham- ming window of 1024 data points with sliding of half the window is used to compute the modified peri- odogram of each segment. In this method the spectrum obtained per hour is the average of 450 spectra with 1024 data points. The Welch method is closely related to the method of complex demodulation described by Bing- ham et al [1967] and hence a separate method to extract SR modal frequencies in the analysis of data is not ap- plied. An example of dynamic spectrum and correspon- ding PSD obtained from the data recorded by the X-component of the sensors on 18 February, 2012 is pre- sented in Figure 1a and 1b respectively. The PSD shows amplitude in dB which is converted in intensity (in nT2/Hz) and the results corresponding to Figure 1b is shown in Figure 1c. In Figure 1a the time is shown in UT which is related with local time LT=UT+5.5hrs. Here, it may be noted that Figure 1a resembles in many respect with Figure 1a of paper I especially the distribution of spikes in the coloured figure. However, these two figures are different because they were observed on different days, and the differences are indicated by PSDs also. Since the time of measurement is same during which a passenger train passes every day, similar spikes appear in both the figures. The detailed examination of the qual- ity of the SR data shows that due to local noises includ- ing the movement of trains on a railway track (about 300 m away in East-West direction) the X-component is af- fected least, the Y-component moderately, but the Z- component severely. Hence, we use X-component data mostly in our analysis. However, we use Y-component data also wherever necessary by applying careful selec- tion of the data. In order to further process the data we write a program in MATLAB to extract hourly informa- tion about the peak frequencies and intensities (deter- mined by the program itself ) of three modes and arrange them in separate columns in EXCEL. In this way, we prepare a time series for the frequencies and intensi- ties of the SR data for 24 hours each day and then for 12 months from 01 March, 2011 to 29 February, 2012. This time series is used to deduce various characteristics of SR phenomenon observed at Agra station. 3. Lightning Imaging Sensor (LIS) data While electromagnetic radiations emitted by light- ning form an important tool for the study of global lightning, another very useful advanced space based tool is the Optical Transient Detector (OTD) which al- lows us to understand the spatial and temporal distri- bution of lightning around the globe, Christian et al. [2003] have given a detailed description of instruments and their characteristics and measurements used to construct lightning climatology maps demonstrating the geographical and seasonal distribution of lightning activity for the globe. However, since OTD satellite stopped working prior to the beginning of measure- ments at Agra, we have used data for this paper from Lightning Image Sensor (LIS) available on the website: http://thunder.nsstc.nasa.gov/data/query/distribu- tion/html collected on Tropical Rainfall Measurement Mission satellite (TRMM) with the following coverage; The data have been used for both the periods of 01 April, 2007 to March 31, 2008 and 01 March, 2011 to 29 February, 2012. It may be noted here that both the elec- tromagnetic (SR) and OTD/LIS measurements suffer SHIFT OF EFFECTIVE LIGHTNING AREAS Figure 1. (a) Frequency-time spectrogram of the ULF/ELF data recorded on 18 February, 2012 at Agra showing four Schumann res- onance lines (b) Power spectrum density (PSD) of the data in Fig- ure 1a and (c) corresponding intensity-frequency variation. Location : Sub-Tropical Latitude range : ± 40° Longitude range : ± 180° Spatial resolution : 3-6 Km Temporal resolution : ORBIT due to some drawbacks. The problem with SR meas- urements lie in a too high global rate of strokes of about 100 events per second. With the basic 8 Hz frequency of SR we have about 12 independent pulses arriving at an observer during the single period of this frequency. It is obvious that signals will overlap severely making impos- sible their separate processing. The main problem in OTD observations is that an observer does not know what is going on outside the frame of optical device at any particular time. Furthermore, there are ‘invisible’ strokes, i.e. the discharges whose optical emissions do not reach the top of the cloud [Nickolaenko and Hayakawa 2002]. In contrast to the above, a space sen- sor capable of mapping both inter cloud and cloud-to- cloud lightning discharges from geostationary orbit during day and night with a spatial resolution of 10 Km and a detection efficiency of 90% has been developed. In addition, this sensor, which is called the Lightning Map- ping Sensor (LMS) can monitor storms on a continual basis. The combination of modern solid state mosaic focal planes with extensive on-board signal processing in the LMS provides a powerful technique for the detection of weak background-contaminated signals and permits the detection of lightning during the day [Christian et al. 1989]. Such data could be very helpful in defining the source diameter of prominent lightning regions. How- ever, the capabilities of lightning detection from geosta- tionary orbit are limited by several factors which include visible area limited to a nominal field of view, resolution depends on geographical position and drops with dis- tance from nadir, earth observation is interrupted and detection efficiency reduces if a spinning satellite is used, light scattering in clouds limits the location to 8-10 Km etc. [Finke and Hauf 2002]. Hence, we have employed the LIS data only in this paper. As far as SR observations are concerned we have mentioned in section 2 how movements of trains on nearby tracks produce spikes in our data whose intensities are not as significant as to ap- pear in the PSDs of the SR data. Especially, the use of North-South (X-component) component of the data per- mits us to process the SR characteristics without much problem. In general, the X-components contain more power than Y-components due to the fact that energies generated by lightning from nearby thunderstorm cen- ter (Asian source) and those from American source trans- mitted along East-West direction are higher than those transmitted in the North-South direction. The raw tem- poral data show sporadic spikes in amplitude due to pas- sage of trains but no saturation effect is observed. The noises arising from other sources (electrical and atmos- pheric noises) are sometimes serious and do affect the quality of data. However, they are occasional and due precaution is taken while processing the data. 4. Observational results 4.1 Monthly distribution of thunderstorm activities In paper I, we have made a correlative study be- tween the first mode SR intensity observed at Agra and Lightning Imaging Sensor (LIS) data for the period 01 April, 2007-31 March, 2008 (period I). Here, we make a similar study for the period 01 March, 2011-29 February, 2012 (period II). Both the results are shown in Figure2a,b. The dark histograms show the LIS data whereas the light dark histograms show the SR intensity. The intensity variations shown in both parts of the figure are inte- grated intensity of the first three modes because a modal SR intensity depends not only on the activities of the source but also on source-observer geometries. From a general scanning of the data we find that the intensity of the second and third modes constitutes approximately 40 percent of the total intensity of the three modes. The consideration of integrated intensity may effectively neu- tralize the geometrical factors. It may be noted from Fig- ure 2 that peak thunderstorm activity occurred in the month of July during period I whereas the same oc- curred in the month of august during period II.The cor- relation coefficients between the SR intensity and LIS data for the two periods are 0.81 and 0.70 respectively. In section 3, we have indicated the short falls in the meas- urements of both the SR and OTD data which may be responsible for the low correlation coefficient. The cor- relation coefficients are low but not poor possibly be- cause both the SR and LIS data are influenced more by south-Asian thunderstorm center (which is just 3000 Km from the observing stations) in spite of the fact that both the measurements respond to global phenomena. Fur- ther, looking at the two panels of Figure 2 we find that the LIS data in the upper panel corresponding to period I are 3-5 times larger than those in the lower panel cor- responding to period II, whereas the magnitude of SR intensities during the two periods are opposite to each other i.e. the intensities are 2-3 times larger in period II than these in period I. These anomalies in SR intensity and LIS data will be discussed later. 4.2 Effective source distance and level of lightening The dependence of frequency range, i.e. the max- imum value minus its minimum value (dF1=f1 max– f1 min) has been used for establishing the effective size of zone occupied by the world wide thunderstorm activ- ity from the SR records [Nickolaenko and Rabinovicz 1995, Nickolaenko et al. 1998]. The diurnal variation in dF1 is related to movement of the lightening region along the equator from one continent to another around the globe. The physical background is based on the dis- continuity of f1 when the point source crosses the nodal SINGH ET AL. 4 5 line. When the source becomes wider the discontinuity turns into the abrupt transition from small to great fre- quency values. The variation in df1 vanish when the lightening sources uniformly cover the globe. A de- tailed description of these results are given in section 4.1.5 and 4.2.2 of Nickolaenko and Hayakawa [2002]. Here, it may be emphasized that dF1 is connetcted with the source size only when measured in the Ez field com- ponent. It reflects the daily alternations of the source distance when measured in the horizontal magnetic field. However, the above results are based on meas- urement of vertical electric field component and one has to measure the peak frequency in the <|E(f )|2> spectra to evaluate the source width [Nickolaenko and Hayakawa 2002, Satori et al. 2009, Ondraskova et al. 2011]. The reason for why electric field measurement is more important than magnetic field measurement for the above purpose may be understood from Figure 3 in which power spectra for both E and H fields in the vicinity of the first mode over the distance-frequency plane is presented (the point source is considered). This figure has been received with the courtesy of Prof. A.P. Nickolaenko. Although the figure is in black and white the results are clearly evident. The left frame shows the power spectra of the vertical electrical field Ez, and the right frame shows the dynamic spectrum of the hori- zontal magnetic field Hφ. The obvious distinction of these spectra is the node around the 10 Mm distance in the Ez-field. The horizontal magnetic field component has the “tilted” maximum here. The tilt is conditioned by the influence of higher resonance modes, and the same influence makes the nodal line asymmetric. As a result, the f1 frequency in the Ez power spectrum has the discontinuity at d=10 Mm. Since we do not have electric filed measurement setup at Agra station we do not use dF1 to evaluate source width during SHIFT OF EFFECTIVE LIGHTNING AREAS Figure 2. (a) Correlation between the global distribution of thunderstorm activities deduced from SR data at Agra for the period of 01 April, 2007 - 31 March, 2008 with Lightning Image Sensor (LIS) data on the satellite (b) the same as above but for the period of 01 March, 2011 - 29 February, 2012. the two periods under consideration using magnetic field data. However, some of the parameters of SR magnetic field data such as f1 and intensity I1 can be used to derive effective source distance (from f1) and, after words, the effective level of lightning activity (from I1, and effective source distance D). For this purpose we follow the equation [Nickolaenko and Hayakawa 2002], f1(d)=ao+a1d+a2d 2+a3d 3+……………………...........(1) Here d is the source-observer distance in Mm, ao= 6.71336, a1=0.181732, a2=-0.008612582, a3=-0.0001477123. This equation can be used to evaluate peak frequency from the known source distance. If we rotate the coor- dinate axes of the graph obtained from above and plot the same data in the form d(f1) the polynomial fit pro- vides [see Nickolaenko 2014] d(f1)=bo+b1f1+b2f1 2+b3f1 3+b4f1 4 ……………………(2) where bo=207049, b1=-110796, b2=22229.7, b3=-1982.03, b4=66.2696. We can use this equation to evaluate the source distance from the experimental peak frequency. Since the real source width is non-zero, this distance should be regarded as effective one. The results of these calculations are presented in graphical forms for the two periods in Figure. 4a to 4c. In the two panels of Fig- ure 4a are shown the seasonal variation of first mode SR frequency of the North-South component (X- com- ponent) of the horizontal magnetic field and corre- sponding source-observer distance derived from equation (2) for period-I (upper panel) and period-II (lower panel) respectively. The range of frequencies shown on the ordinate of the panel of period-I is higher than that on the ordinate of panel-II. An explanations for the decrease in the first mode frequency in period- II and increasing trend with the months will be given in discussion section. Here, we see that in both the pe- riods the summer thunderstorms shift towards the northern hemisphere (NH) summer (May-August) as the minimum source-observer distances reach around 12.4 Mm and 9 Mm respectively. Annual drifts of these sources in the two periods are approximately 4 Mm and SINGH ET AL. 6 Figure 4(a). Seasonal variation of first mode SR frequency of the North-South component (X- component) of the horizontal magnetic field for period-I (upper panel) and period-II (lower panel). The ordi- nate on the right shows the corresponding source-observer distance. Figure 4(b). Seasonal variation of source-observer distance derived from observed first mode frequency for the two periods. Figure 3. Power spectra of the electric filed E (left) and magnetic field H (right) in the vicinity of the first mode SR over the distance-fre- quency plane (received with the courtesy of Prof. A.P. Nickolaenko).. 7 1.3 Mm corresponding to 36° and 12° respectively. The seasonal variations of the effective source-observer dis- tance derived from the first mode SR frequency for the two periods are separately presented in Figure 4b. It should be noted that the obtained estimates of the dis- tances to the source in the two periods are more than 10 Mm (for period-I) and less than 10 Mm (for period- II) respectively. Clearly visible are the seasonal drift of world thunderstorms to the NH summer. Satori [2003] has carried out long-term monitoring of the Ez component of the SR parameters (peak amplitude and peak frequencies) and provides very useful informa- tion on seasonal motion of the global lightning activ- ity [see also Hayakawa et al 2010]. The variations of average intensity of the first mode SR corresponding to North-South components for the two periods are shown in the upper and lower panels of Figure 4c. The intensity data are the same as shown in Figure 2. How- ever, they are shown here to be examined in the light of the variation of source-observer distances pre- sented in Figure 4b. Here, we see that the intensities are higher in the summer months (May-August) when the source-observer distances are lowest indicating that nearer the sources larger are the intensities. An interesting point to be noted here is that the level of maximum intensity corresponding to intense light- ning exists in the months of July in period-I whereas the same exists in the months of August in period-II. These months of largest SR intensities are supported by LIS data also (see Figure 2). Hence, we can say that there is a shift of maximum lightening areas from July in period-I to August in period-II. There is an increase in the intensities in the month of October in period-I also but it is not supported by LIS data which are less than those in the month of July. The intensities fall in the southern hemisphere (SH) summer (November- January) because the thunderstorms move away from the observer as shown by the increasing trend in source-observer distance. 5. Discussion As per NASA report the solar cycle minimum of 2008-2009 was extraordinarily deep and the sun was quietest ever seen in almost a century. The X-ray radia- tion, which ionizes the earth’s atmosphere and is re- sponsible for the conductivity profile of the ionosphere, was by an order of magnitude lower than the solar cycle minimum studied previously [Onraskova et al. 2011]. In Figure 5 we show the variation of total num- ber of M and X classes of flares occurred between the years 2006 and 2012. These flares are shown by hatched and open portions of histograms respectively. As it may be seen from the figure, the number of the two classes of flares decreases towards the solar cycle minimum such that there is no flare of either class in 2009. Then the number increases as the solar cycle increases. The solar flare data are obtained from NOAA National Geophysical Data Center on the website: www.ngdc.noaa.gov. It is well known that according to classical Schumann formula, the eigen-frequencies are given by the expression fn=7.5 √n (n+1) which are 10.6, 18.4, 26… Hz for a spherical Earth-ionosphere cavity. However, the real (experi- mental) values are 8, 14, 20… Hz. The deviations are caused by the influence of the finite conductivity of the ionosphere [Nickolaenko and Hayakawa 2002]. From Figure 1 it may be seen that peak frequencies of the dif- ferent modes do not coincide with their corresponding customary frequencies but are lower. For example, the fundamental mode looks to be near 7.5 Hz (instead of closer to 8 Hz), the third at 19 Hz rather than 20 Hz, and the fourth at 25 Hz (instead of 26 Hz). A possible reason for the low frequencies seems to be the fact that since the period of observation corresponds to the in- creasing phase of solar cycle minimum of 2008-2009, the effect of deep minimum activity still continues SHIFT OF EFFECTIVE LIGHTNING AREAS Figure 4(c). Seasonal variation of first mode intensity for the two periods. Figure 5. Variation of M and X classes of solar flares between the years 2006 and 2012. under which the modal frequencies are influenced. It may be mentioned here that Ondraskova et al [2011] have reported similar decrease in the frequencies of all the four modes from their observations at Modra ob- servatory from solar cycle maximum of 2001-2002 to solar cycle minimum of 2008-2009. Satori et al [2005] have also reported a similar decrease in the frequencies. From Figure 4a we see that the peak frequencies of the first SR mode are much lower in period-II than those in period-I. The peak frequencies are lower than that of customary frequency of 8 Hz in period-I also. Reduc- tions are also seen in peak frequencies of SR data ob- tained at Shillong station in India for the months of July and August 2007 and March and May, 2008 (Data ob- tained with the courtesy of Dr. B.M. Pathan, Indian In- stitute of Geomagnetism, Mumbai, May, 2010). Recently, Ondraskova et al. [2009, 2011] have reported a significant decrease of about 0.30 Hz in the funda- mental SR frequency from observations at Modra ob- servatory in Slovak Republic during the latest solar cycle having minimum around 2008-2009. This ex- traordinary fall of the fundamental mode frequency is attributed to the unprecedented drop in the ionizing ra- diation in X-ray frequency band, although the patterns of the daily and seasonal variation remain unchanged. Similar decreases in the first mode SR frequency have also been reported by Satori et al. [2005] for the previ- ous 11 year solar cycle. Further decreases in frequency in period-II but having an increasing trend with the month may be interpreted in terms of the minimum solar cycle effect and increasing period of current solar cycle (2009-2020). According to Satori et al. [2005], changes of the X-ray radiation dominate the variations in the conductivity profile within the upper character- istic layer (90-100 Km portion of E-region). The de- crease of this conductivity by up to one order of magnitude over the solar cycle is responsible for the ob- served SR frequency decrease by several tenths of Hz. Although, there is a controversy that this is the inter- mediate region where electron conductivity is substi- tuted by the ionic one. Therefore knee models of the conductivity profiles are used while modelling SR [e.g. Mustak and Williams 2002]. In general two possibilities have been suggested for the reduction in frequencies of the first mode; (i) variation in median distance between the source and the observer, and (ii) ionospheirc modi- fication. Peak frequency changes to the extent of 10-20 percent have been observed which are attributed to the uncertainties arising from spatial distribution of light- ning sources exciting the SR modes [De 2007]. De et al [2009] have also observed large diurnal variation in the peak frequency of the first mode and found similar vari- ation in the data obtained at Moshiri in Japan. Model calculations have shown that a uniform decrease in the reference ionospheirc height reduces the resonant fre- quencies and Q factors (i.e. it simultaneously increases the wave attenuation). Specifically, it is found that a shift from 8 to 7.6 Hz in the peak frequency of the first mode is caused by a global reduction of ionospheirc height by approximately 5 Km. However, such reduction will be observed in all the modes, and the effect will be ob- served globally. Since measurements differ from place to place it is concluded that the deviations observed experimentally are conditioned by the displacement of the global thunderstorms rather than by modifica- tion of the global ionospheirc profiles [Nickolaenko and Hayakawa 2002]. The long-term observations of SR carried out by Ondraskova et al [2007] have also confirmed that variation in peak frequency of the lower SR modes can be attributed mainly to the source-observer distance effect. The intensity varia- tion during the two periods (see Figure 2) is influenced by solar activities. Since the period I corresponds to close to solar cycle minimum of 2008-2009, the inten- sities are lower. In contrast, the period II corresponds to increasing solar cycle, the intensities are higher. These results are very well supported by the results of earlier workers. For example, Kulak et al. [2003] have analyzed six years of SR data at a mid-latitude station which included the periods of minimum to maximum (1995-2001) of solar cycle 23. The results of the analy- sis and modelling have shown that first SR frequency increases from 7.75 Hz at minimum to about 7.95 Hz at maximum while the global mean attenuation rate at 8 Hz varies from 0.31 dB/Mm at minimum to about 0.26 dB/Mm at maximum. Recently, Koloskov et al. [2013] have reported the results of long period (2002- 2012) horizontal magnetic field SR observations con- ducted at Ukrainian Antarctic station “Academician Vernadsky” and shown that both the SR first mode frequency and horizontal component’s intensity fol- low the trend in variation of solar activity. The LIS data relate to cloud formation and lightning flashes near the earth’s atmosphere and are also influenced by solar activities. A detail discussion on this topic will follow later. A possible explanation for shifting of the peak thunderstorm activity from the month of July, 2007 in period-I to the month of August, 2011 in pe- riod-II may be given in terms of galactic cosmic ray (GCR) variation with solar activity. Many workers [Tinsley and Yu 2004, Siingh and Singh 2010, Siingh et al. 2011, Williams and Mareev 2014] have reviewed the cosmic ray- mediated cloud microphysics on the 11 year solar cycle. It has been mentioned that in all spe- cific mechanisms, increased GCRs during solar mini- mum are linked with increased cloud condensation SINGH ET AL. 8 9 nuclei (CCN) and ultimately with increased cloudi- ness, and (in the atmospheric electrical context) with increased cloud water in the mixed phase region of moist convection [Williams et al. 2002]. This has been supported by other workers also [Pierce and Adams 2009, Erlykin and Wolfandale 2011].It is now a well- established fact that GCR flux is reduced as solar ac- tivity increases. This is because the increased solar magnetic flux associated with enhanced solar wind during high solar activity periods deflects the GCRs away from the heliosphere and hence from entering the earth’s atmosphere. The reduced levels of GCRs modulated by increasing solar wind can significantly decrease cosmogenic ionization resulting in decreased conductivity and thus increase electric potential. Any significant change in electric field in the region of cloud formation can affect cloud formation, cloud movement etc. through a number of processes not yet understood [Lakshmi 2008]. The possible long term correlation between GCRs and low cloud cover (LCC) could be explained considering the influence of solar activity [Erlykin et al. 2009]. The enhanced solar irra- diance causes a rise of the mean surface temperature and results in enhanced vertical convection. From below 3 Km warm air rises to greater heights and causes the LCC to decrease and median cloud cover (MCC) to increase setting a change in vertical gradient of the electric field. Thus, during high solar activity periods the decrease in GCRs causes a decrease in LCC, although based on this argument no casual con- nection between GCRs and LCC could be established. Keeping the above arguments in view, it is possible that the intense level of lightning in the month of July in period-I is influenced by the increasing solar activ- ity in period-II and shifted to the month of August in period-II. Since the increased solar activity during the period-II is responsible for decreased cloud formation, this would result in reduction in lightning activities and hence a reduction in LIS data as found in Figure 2 during period-II. However, this explanation is tenta- tive and a detailed explanation for the fall of LIS data during period-II is under study. In all the results pre- sented in Figure 4, we have adopted a single source model due to Nickolaenko and Hayakawa [2002] as- suming that the global lightening can be approxi- mated by a single source migrated in longitude around the earth. This is not a good model for the real global lightening source, which exists in three well rec- ognized continental zones, none of which ever really goes to sleep over the course of 24 UT hours. Hence, a model incorporating the three sources is better and this has been recognized by Nickolaenko et al. [1998]. However, it is also far from the ideal because the source parameters depend (although slightly) on the particular field site. Applications of OTD/LIS data also does not provide one-to-one correspondence to observations. And also we do not know spatial distri- butions of strokes [A.P. Nickolaenko, personal com- munication, September 2014]. 6. Conclusion Employing a set of 3-component search coil mag- netometer and analysing the X-component data, we study the variation of first mode SR intensity and fre- quency during two periods of 01 April, 2007 to 31 March, 2008 (Period I) and 01 March, 2011 to 29 Febru- ary, 2012 (Period II) which correspond to pre and post periods of solar cycle minimum of 2008-2009. Using the intensity data, we first study the monthly distribution of global thunderstorm activity and then using both the frequency and intensity data we study the effective source distance and level of lightning during both the periods. We show that the world thunderstorm activity shifts to summer month in the northern hemisphere summer as the effective source distance approaches close to the observer, by 12 and 9 Mm respectively in both the periods. The level of lightening in different months of the two periods are supported by the LIS data with correlation coefficients of 0.81 and 0.70 re- spectively. The level of intense lightning lies in the month of July, 2007 in period-I which shifts to August, 2011 in period-II. This shifting is attributed to increasing solar activity in period-II due to which cloud formation and movement are reduced and hence the effective lightning also is delayed. Acknowledgements. The authors are highly grateful to the Department of Science and Technology, Government of India, New Delhi for the financial support in the form of major research project. The authors also thank Mr. Ajay Singh for typing the manuscript and help with the figures. References Balser, M. and C. A. Wagner (1960). Observation of Earth-ionosphere cavity resonances, Nature, 188, 638-641. Bingham, C., M.D. Godfrey and J.W. Turkey (1967). Modern techniques of power spectrum estimation, IEEE trans. Audio and Electroacoust, AU-15, 56-66. Christian, H.J., R.J. Blakeslee and S.J. Goodman (1989). The detection of lightening from geostationary orbit, J. Geophys. Res: Atmospheres (1984-2012), 94, 13329-13337. Christian, H.J., R.J. Blakeslee, D.J. Boccippio, W.L. Boeck, D.E. Buechler, K.T. Driscoll, S.J Goodman, J.M. Hall, W.J., Koshak, D.M. Mach and M.F. Stewart (2003). Global frequency and distribution of light- SHIFT OF EFFECTIVE LIGHTNING AREAS ening as observed from space by the Optical Tran- sient Detector, J. Geophys. Res., 108, doi: 10.1029/ 2002 JD002347. De, S. S. (2007), Guest Editorial, Schumann resonance. Its different aspects and latest wonders, Indian J. Radio & Space Phys., 36, 359-361. De, S.S., B. K. De, B. K. Sarkar, B. Bandyopadhyay, D. K. Halder, P. Suman and S. Barui (2009), Analysis of Schumann resonance spectra from Kolkata and their possible inter pretations, Indian J. Radio & Space Phys., 38, 208-214. Erlykin, A.D., T. Sloan and A.W. Wolfendale (2009). The search of cosmic ray effects in the clouds, J. Atmos. Solar-Terr. Phys., 71, 955-958. Erlykin, A.D. and A.W. Wolfendale (2011). Cosmic rays effects on the cloud cover and their relevance to cloud change, J. Atmos. Solar-Terr. Phys., 73 (13), 1681-1686. Finke, U. and T. Hauf (2002). Detect and locate lighten- ing events from geostationary satellite observations, report part II: Feasibility of lightening location from a geostationary orbit, No. EUM/CO/02/1016/SAT, Institute for Meteorologie und Krimatologie, Uni- versity Hannover. Greenberg, E. and C. Price (2007). Diurnal variation of ELF transients and background noise in the Schu- mann resonance band, Radio Sci., 42, RS2S08, doi: 10.1029/2006RS003477. Hayakawa, M., K. Ohta, A. P. Nickolaenko and Y. Ando (2005). Anomalous effect in Schumann resonance phenomena observed in Japan, possibly associated with the Chi-Chi earthquake in Taiwan, Ann. Geo- phys, 23, 1335-1346. Hayakawa, M., A.V. Shvets, Y. Hobara and A.P. Nicko- laenko (2010). Recent studies of Schumann reso- nance and ELF transients, In lightening: properties, Formation and Types, chapter 3, Ed. (Mathew D. Wood), Nova Science Pub., London, pp: 1-33. Koloskov, A.V., A.V. Spikelets, N.A. Baru, O. Budanov, V.G. Rutlers, B.J. Gavriliouk, A.V. Paznuhov and Y.M. Yampolsky (2013). Diagnosis of global light- ening activity on the basis of long-term monitoring signals Schumann resonance on the UAS, Ukrainian Antarctic Journal, N.12, 170-176. Kulak, A., J. Kubisz, A. Michalec, S. Zieba and Z. Nieckarcz (2003). Solar variation in extremely low fre- quency propagation parameters: 2. Observation of Schumann resonances and compulation of the ELF attenuation parameters, J. Geophys. Res., 108, 1271, doi 10.1029/2002 J/A009305. Lakshmi, D. R. (2008). A brief review of solar activity modulations of global electric circuit, in the pro- ceedings of DST sponsored workshop on Electro- dynamical coupling of atmospheric regions, organ- ized by Indian Institute of Geomagnetism, Navi Mumbai, 25-26, November, pp. 6-8. Melnikov, A., C. Price and G. Satori (2004). Influence of solar terminator passages on Schumann reso- nance parameters, J. Atmos Solar-Terr. Phys., 66, 1187-1194. Mustak, V.C. and E.R. Williams (2002). ELF propaga- tion parameters for uniform models of the earth- ionosphere wave guides, J. Atmos. Solar-Terr. Phys., 64, 1989-2001. Nickolaenko, A. P. and L. M. Rabinowicz (1995). Study of annual changes of global lightning distribution and frequency variations of the first Schumann resonance mode, J. Atmos. Solar-Terr. Phys., 57, 1345-1348. Nickolaenko, A.P., G. Satori, V. Ziegler, L.M. Rabinow- icz and L.G. Kudintseva (1998). Parameters of global thunderstorm activity deduced from the long-term Schumann resonance records, J. Atmos. Solar-Terr. Phys., 60, 387-399. Nickolaenko, A. P. and M. Hayakawa (2002). Reso- nances in the earth-ionosphere cavity, Kluwer Aca- demic Publishers, Dordrecht, the Netherlands. Nickolaenko, A. P., M. Hayakawa, M. Sekiguchi, Y. Ando and K. Ohta (2006). Model modification in Schumann resonance intensity caused by a localized ionosphere disturbance over the earthquake epi- center, Ann. Geophys., 24, 567-575. Nickolaenko A. P., O. Penchony and C. Price (2006), Model variations of Schumann resonance based on Optical Transient Detector maps of global lightning activity, J. Geophys. Res. 111, D23102, doi: 10.1029/ 2005JD006844. Nickolaenko, A. P. and M. Hayakawa (2008). Universal and local time components in Schumann resonance intensity, Ann. Geophys., 26, 813-822. Nickolaenko, A. P., E. I. Yatsevich, A. V. Shvets, M. Hayakawa and Y. Hobara, (2011). Universal and local time variations deduced from simultaneous Schumann resonance records at three widely sepa- rated observatories, Radio Sci., 46, 5, October, doi: 10.1029/2011RS004663. Nickolaenko, A. P. (2014), Deducing the world thun- derstorm activity from the Schumann resonance record (once again about point source model in Schumann studies), Radio Phys. and Electronics, 5(10), 33-41 (in Russian). Ondraskova, A., P. Kostecky, S. Sevcik and L. Rosenberg (2007). Long term observation of Schumann reso- nance at Modra observatory, Radio Sc., 42, RS2509, doi: 10.1029/2006RS003478. Ondraskova, A., S. Sevcik and P. Kostecky (2009). A sig- nificant decrease of the fundamental Schumann res- SINGH ET AL. 10 11 onance frequency during the solar cycle minimum of 2008-9 as observed at Modra Observatory, Con- tributions to Geophysics and Geodesy, 39/4, 345-354. Ondraskova, A., S. Sevcik and P. Kostecky (2011). De- crease of Schumann resonance frequencies and changes in the effective lightening areas towards the solar cycle minimum of 2008-2009, J. Atmos. Solar - Terr. Phys., 73, 534-543. Pechony, O. and C. Price (2006). Schumann resonances: Interpretation of local diurnal intensity modula- tions, Radio Sci., 41, 2, doi: 10.1029/2006RS003455. Pechony, O., Price C. and A. P. Nickolaenko (2007). Rel- ative importance of the day-night asymmetry in Schumann resonance amplitude records, Radio Sci- ence, doi: 1029/2006RS003456. Pierce, J.R. and P.S. Adams (2009). Can cosmic rays affect cloud condensation nuclei by altering new particle formation rates? Geophys. Res. Lett., 36 (LO 9820). Price, C. and A. Melnikov (2004). Diurnal, seasonal and inter-annual variations of the Schumann resonance parameters. J. Atmos. Solar-Terr. Phys., 66, 1179-1185. Satori, G. (2003). On the dynamics of the north-south seasonal migration of global lightning, Pro. of 12th ICAE, Versailles, France, global lightening and cli- mate, pp. 1-4. Satori, G., E. Williams and V. Mushtak (2005). Response of the Earth-ionosphere cavity resonator to the 11- year solar cycle in X-radiation, J. Atmos. Solar-Terr. Phys., 67, 553-562. Satori, G., M. Neska, E. Williams and J. Szendroi (2007a). Signatures of the day-night asymmetry of the Earth- ionosphere cavity in high time resolution Schumann resonance records, Radio Sci., 42, RS2510, doi: 10.1029/2006RS103483. Satori, G., I. Lemperger and J. Bor (2007b). Modulation of annual and semiannual areal variations of global lightning on the 11-year solar cycle, In: Second In- ternational Symposium on lightening physics and effects, Vienna, April 19 and 20. Satori, G., V. Moshtak and E.R. William (2009). Schu- mann resonance signatures of global lightning ac- tivity, In Lightning: Principles, instruments and applications [eds] Betz, H.D., Schumann, U., Laroche, P., Springer, pp. 641. Schumann, W. O. (1952). On the free oscillations of a conducting sphere which is surrounded by an air layer and an ionosphere shell, Z. Naturforschaftung, 7a, 149-154. Sentman, D. D. and B. J. Fraser (1991). Simultaneous observations of Schumann resonances in California and Australia: evidence for intensity modulation by the local height of the D region. J. Geophys. Res., 96 (A9), 15973-15984. Siingh, D. and R.P. Singh (2010). The role of cosmic rays in the earth’s atmospheric processes, Pramana, 74, 153-168. Siingh, D., R.P. Singh, A.K. Singh, M.K. Kulkarni, A.S. Gautam and A.K. Singh (2011). Solar Activity, Light- ening and climate, Article No. 912, Springer, doi: 10.1007/s 10712-011-9127-1. Singh, B., R. Tyagi, Y. Hobara and M. Hayakawa (2014). X-rays and solar proton event induced changes in the first mode Schumann resonance frequency observed at a low latitude station Agra, India, J. Atmos. Solar- Terr. Phys., 113C, 1-9, doi: 10.1016/j.jstp.2014.02010. Tinsley, B.A. and F. Yu (2004). Atmospheric ionization and clouds as link between solar activity and cli- mate, In Solar variability and its effects on climate, Geophysical monograph No. 141, A.G.U., doi: 10.1029/141GM22. Tyagi, R., B. Singh, O.P. Singh and R.L. Dua (2013). Characteristics of global thunderstorm extracted from Schumann resonance at Agra, Int. J. Adv. Res. in Electrical Electronics and Instrumentations Engg., 2(3), 1129-1137. Welch, P. D. (1967). The use of Fast Faurier transform for the estimation of power spectra: A method passed on time averaging over short modified peri- odgrams, IEEE Trans. Audio and Electroacoust., AU-16, 70-73. Williams, E.R. and Coauthors (2002). Contrasting con- nective regims over the Amazon: Implication to cloud electrification, J. Geophys. Res., 107, http://dx. doi.org/10.1029/2001JD000380 (LBA special issue, D20, 8081). Williams, E.R. and G. Satori (2007). Solar radiation- in- duced changes in ionospheric height and the Schu- mann resonance waveguide on different timescales, Radio Sci., 42, RS2S11, doi: 1029/2006RS003494. Williams, E. and E. Mareev (2014). Recent progress on the global electric circuit, Atmos. Res. (135-136), 208-227, doi.org/10.1016/j.armores. 2013.05.015. *Corresponding author: Birbal Singh, RBS Engineering Technical Campus, Department of Electronics and Communication Engineering, Bichpuri, Agra, India; email: bbsagra@gmail.com. © 2014 by the Istituto Nazionale di Geofisica e Vulcanologia. All rights reserved. SHIFT OF EFFECTIVE LIGHTNING AREAS << /ASCII85EncodePages false /AllowTransparency false /AutoPositionEPSFiles false /AutoRotatePages /None /Binding /Left /CalGrayProfile (Dot Gain 20%) /CalRGBProfile (sRGB IEC61966-2.1) /CalCMYKProfile (U.S. Web Coated \050SWOP\051 v2) /sRGBProfile (sRGB IEC61966-2.1) /CannotEmbedFontPolicy /Warning /CompatibilityLevel 1.3 /CompressObjects /Tags /CompressPages true /ConvertImagesToIndexed true /PassThroughJPEGImages true /CreateJobTicket false /DefaultRenderingIntent /Default /DetectBlends true /DetectCurves 0.1000 /ColorConversionStrategy /LeaveColorUnchanged /DoThumbnails false /EmbedAllFonts true /EmbedOpenType false /ParseICCProfilesInComments true /EmbedJobOptions true /DSCReportingLevel 0 /EmitDSCWarnings false /EndPage -1 /ImageMemory 1048576 /LockDistillerParams true /MaxSubsetPct 100 /Optimize false /OPM 1 /ParseDSCComments true /ParseDSCCommentsForDocInfo true /PreserveCopyPage true /PreserveDICMYKValues true /PreserveEPSInfo true /PreserveFlatness true /PreserveHalftoneInfo false /PreserveOPIComments false /PreserveOverprintSettings true /StartPage 1 /SubsetFonts true /TransferFunctionInfo /Apply /UCRandBGInfo /Preserve /UsePrologue false /ColorSettingsFile (None) /AlwaysEmbed [ true /AndaleMono /Apple-Chancery /Arial-Black /Arial-BoldItalicMT /Arial-BoldMT /Arial-ItalicMT /ArialMT /CapitalsRegular /Charcoal /Chicago /ComicSansMS /ComicSansMS-Bold /Courier /Courier-Bold /CourierNewPS-BoldItalicMT /CourierNewPS-BoldMT /CourierNewPS-ItalicMT /CourierNewPSMT /GadgetRegular /Geneva /Georgia /Georgia-Bold /Georgia-BoldItalic /Georgia-Italic /Helvetica /Helvetica-Bold /HelveticaInserat-Roman /HoeflerText-Black /HoeflerText-BlackItalic /HoeflerText-Italic /HoeflerText-Ornaments /HoeflerText-Regular /Impact /Monaco /NewYork /Palatino-Bold /Palatino-BoldItalic /Palatino-Italic /Palatino-Roman /SandRegular /Skia-Regular /Symbol /TechnoRegular /TextileRegular /Times-Bold /Times-BoldItalic /Times-Italic /Times-Roman /TimesNewRomanPS-BoldItalicMT /TimesNewRomanPS-BoldMT /TimesNewRomanPS-ItalicMT /TimesNewRomanPSMT /Trebuchet-BoldItalic /TrebuchetMS /TrebuchetMS-Bold /TrebuchetMS-Italic /Verdana /Verdana-Bold /Verdana-BoldItalic /Verdana-Italic /Webdings ] /NeverEmbed [ true ] /AntiAliasColorImages false /CropColorImages true /ColorImageMinResolution 150 /ColorImageMinResolutionPolicy /OK /DownsampleColorImages true /ColorImageDownsampleType /Bicubic /ColorImageResolution 300 /ColorImageDepth -1 /ColorImageMinDownsampleDepth 1 /ColorImageDownsampleThreshold 1.10000 /EncodeColorImages true /ColorImageFilter /DCTEncode /AutoFilterColorImages true /ColorImageAutoFilterStrategy /JPEG /ColorACSImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /ColorImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /JPEG2000ColorACSImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /JPEG2000ColorImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /AntiAliasGrayImages false /CropGrayImages true /GrayImageMinResolution 150 /GrayImageMinResolutionPolicy /OK /DownsampleGrayImages true /GrayImageDownsampleType /Bicubic /GrayImageResolution 300 /GrayImageDepth -1 /GrayImageMinDownsampleDepth 2 /GrayImageDownsampleThreshold 1.10000 /EncodeGrayImages true /GrayImageFilter /DCTEncode /AutoFilterGrayImages true /GrayImageAutoFilterStrategy /JPEG /GrayACSImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /GrayImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /JPEG2000GrayACSImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /JPEG2000GrayImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /AntiAliasMonoImages false /CropMonoImages true /MonoImageMinResolution 1200 /MonoImageMinResolutionPolicy /OK /DownsampleMonoImages true /MonoImageDownsampleType /Bicubic /MonoImageResolution 1200 /MonoImageDepth -1 /MonoImageDownsampleThreshold 1.08250 /EncodeMonoImages true /MonoImageFilter /CCITTFaxEncode /MonoImageDict << /K -1 >> /AllowPSXObjects false /CheckCompliance [ /None ] /PDFX1aCheck false /PDFX3Check false /PDFXCompliantPDFOnly false /PDFXNoTrimBoxError true /PDFXTrimBoxToMediaBoxOffset [ 0.00000 0.00000 0.00000 0.00000 ] /PDFXSetBleedBoxToMediaBox true /PDFXBleedBoxToTrimBoxOffset [ 0.00000 0.00000 0.00000 0.00000 ] /PDFXOutputIntentProfile (None) /PDFXOutputConditionIdentifier () /PDFXOutputCondition () /PDFXRegistryName (http://www.color.org) /PDFXTrapped /Unknown /CreateJDFFile false /SyntheticBoldness 1.000000 /Description << /ENU (Use these settings to create PDF documents with higher image resolution for high quality pre-press printing. The PDF documents can be opened with Acrobat and Reader 5.0 and later. These settings require font embedding.) /JPN /FRA /DEU /PTB /DAN /NLD /ESP /SUO /NOR /SVE /KOR /CHS /CHT /ITA >> >> setdistillerparams << /HWResolution [2400 2400] /PageSize [595.000 842.000] >> setpagedevice