DOI QR코드

DOI QR Code

A New Optimized Localized Technique of CG Return Stroke Lightning Channel in Forest

  • Kabir, Homayun (Faculty of Engineering, Department of Electrical Engineering, University of Malaya) ;
  • Kanesan, Jeevan (Faculty of Engineering, Department of Electrical Engineering, University of Malaya) ;
  • Reza, Ahmed Wasif (Faculty of Engineering, Department of Electrical Engineering, University of Malaya) ;
  • Ramiah, Harikrishnan (Faculty of Engineering, Department of Electrical Engineering, University of Malaya) ;
  • Dimyati, Kaharudin (Department of Electrical and Electronics Engineering, Faculty of Engineering, National Defence University of Malaysia)
  • 투고 : 2015.01.02
  • 심사 : 2015.06.19
  • 발행 : 2015.11.01

초록

Localization of lightning strike point (LSP) in the forest is modeled to mitigate the forest fire damage. Though forest fire ignited by lightning rarely happens, its damage on the forest is grievousness. Therefore, predicting accurate location of LSP becomes crucial in order to control the forest fire. In this paper, we proposed a new hybrid localization algorithm by combining the received signal strength (RSS) and the received signal strength ratio (RSSR) to improve the accuracy by mitigating the environmental effect of lightning strike location in the forest. The proposed hybrid algorithm employs antenna theory (AT) model of cloud-to-ground (CG) return stroke lightning channel to forecast the location of the lightning strike. The obtained results show that the proposed hybrid algorithm achieves better location accuracy compared to the existing RSS method for predicting the lightning strike location considering additive white Gaussian noise (AWGN) environment.

키워드

1. Introduction

Forest is a prime facet of bionomical research. Forest fire extenuation has become a vital issue at the present time to minimize the mass ecological catastrophe. Thousands of square kilometers are cremated by forest fires with an immense amount of damage. More than 100 thousand forest fires occurred in the last decade in the whole world. The main cause of forest fire includes lightning, volcanic eruption, sparks from rock falls, and human activities. Though only one third of forest fires are enkindled by lightning, more than 85% of the total burned area is caused by lightning in Canada [1]. In Austria, 15% of forest fire is due to lightning and this figure rises to 40% during summer season [2].

Cloud-to-ground (CG) lightning is mainly responsible for forest fire. In nature, 90% of CG lightning is negative and the rest is positive. For the negative return stroke, peak current range, rate of change of the current range, and stroke duration are (4.6-80) kA, (5.5-32) kA/µs, and (30-200) µs, respectively. On the contrary, for the positive return stroke, peak current range, rate of change of the current range, and stroke duration are (4.6-250) kA, (5.5-32) kA/µs, and (25-2000) µs, respectively [3]. High temperature, which is generated by current and sufficient duration time of return stroke, is necessary to ignite woody fuels. Hence, the positive lightning is the prime candidate rather than negative lightning for setting off the forest fire [3]. To control lightning ignited forest fire, detection of igniting strike point in the forest (mainly positive lightning) is therefore essential [4]. Several techniques are studied to locate the lightning strike point (LSP). Most of the techniques are adopting time of arrival (TOA), time difference of arrival (TDOA), received signal strength (RSS) of receivers, and magnetic direction finder (MDF) [5-10].

Improving the location accuracy is the main challenge in localizing the lightning system to mitigate the forest fire effect. Cost and complexity of the algorithm are also vital issues for the localization. In TOA and TDOA techniques, time synchronization between global positioning systems (GPS) is required [11]. Nevertheless, it is complicated and costly [12]. In contrast, the RSS based localization is more straightforward and cost effective. In this study, the direct localization method is simulated first using the received electrical signal strength in a 2D environment. This position estimation is used as the boundary limits of genetic algorithm (GA) search space. The ratio of the received electrical signal of different couples of sensors is embedded in the fitness function of the GA. In this paper, the composite model of RSS and the received signal strength ratio (RSSR) is combined with the GA to estimate better localization accuracy of the LSP that will be compared with the RSS method.

The whole paper has been arranged in the following manner. Section 2 explains the propagation of lightning signal. In Section 3, we propose the hybrid localization algorithm of RSS and RSSR. In Section 4, we present the results analysis and performance evaluation. Finally, Section 5 provides the concluding remarks.

 

2. Lightning Signal Propagation

The CG lightning is a natural electrical spark, which is almost two or three kilometers in length. In this study, the CG return stroke lightning channel radiation propagation model is investigated for better representation of real life CG lightning strike. From the previous study, Bruce-Golde (BG) (equal current in whole return stoke channel), traveling current source (TCS), transmission line (TL), modified transmission line with linear current decay (MTLL), and modified transmission line with exponential current decay (MTLE) are developed to model the lightning return-stroke channel [13]. Those engineering models are grouped into two categories: the transmission line type models and traveling-current type models. The first category contains the TL, MTLL and MTLE models that are viewed as incorporating a current source at the channel as injecting a specified current wave which propagates upward into the channel without any distortion. In the traveling current type models, to which belongs the TCS and the BG models, the return current is viewed as being generated at the upward moving return stoke front and then propagating downward.

In the present decade, researchers have started to develop a new model based on antenna theory (AT). Researches proved that AT is one of the best models of electromagnetic radiation of CG return stroke lightning channel [14], [15]. The diameter of lightning channel is considered around 0.01 to 0.18 m [16]. Hence, lightning flash is a finite length of current carrying thin conductor, which radiates electromagnetic wave [17], [18].

The electromagnetic waves that are emitted from CG return stroke in forest propagate throughout the environment, are received by the sensors, as shown in Fig. 1.

Fig. 1.Lightning strike in forest and measurement system of electromagnetic signal

Generally, the lightning current is higher in the channel closer to the ground. Therefore, the maximum field radiates from the lower section of lightning channel [18]. Hence, the effective part of the lightning channel can be treated as a dipole antenna [19]. From the channel, three types of electromagnetic signal, namely, radioactive, inductive, and electrostatic are radiated. In the far field, the receiver only picks up the vertical polarized electrical signal [20]. The vertical polarized electric signal is given in (1) [19].

where,

and

Here, E0 is the peak electric signal strength of lightning channel, μ0 is the permittivity of free space, ω is the angular velocity, Ei is the RSS at receiver, (x, y) is the unknown coordinate of the electromagnetic source of LSP, ( xi, yi ) is the known coordinate of the receivers, Zi is the height of the receivers from the ground, z1 is the height of the lower part of the source (LSP) from the ground, z2 is the height of the upper part of the source (LSP) from the ground, ri is the horizontal distance between the source and the receivers, i is the index of the receivers 1, 2, 3, 4, …, N, Ia is the pick channel based current, I(z1,t) is the standard channel based current of lightning return stroke, I0 is the amplitude of the channel based current, τ1 is the front time constant, τ2 is the decay time constant and n is the exponent having the value of 2 to 10.

Most of the theoretical studies are dealing with the electromagnetic field radiation from lightning strikes to the forest, assuming the ground to be a perfectly-conducting (P-C) plane. However, there is experimental evidence that the finite ground conductivity affects the radiated electromagnetic fields [21]. More precisely, the propagation of lightning electromagnetic fields along the earth surface results in the attenuation of the radiated components of the field spectrum associated with the initial rise to the peak of the electric. The equation (1a) is obtained after considering the ground conductivity of σ.

where is the attenuation function for the case of homogeneous ground [22].

Due to high amplitude of lightning radiated signals, the noise is considered based on additive white Gaussian noise (AWGN) only [23]. Thus, the measurement field strength can be expressed as follows:

where Ei = [E1, E2,..., EN] is the RSS of the receivers without environmental effect and n(An,σ) is AWGN [24]. Here, σ is the variance and An is the mean amplitude of the AWGN noise.

 

3. Proposed Hybrid Localization Algorithm

In this section, we discuss the proposed hybrid localization algorithm of RSSR and RSS, which is based on AT of CG return strokes lightning channel, followed by GA.

The sensor and LSP geometry of RSS method is shown in Fig. 2, where minimum three sensors are needed to estimate the position of LSP. Three circles are created by the three sensors, which are intersected in one point and that point is the estimated position of LSP. In the RSSR method, minimum three sensors are also required. The received signal ratio of the two sensors is formed as a hyperbola, which has two focus points. Two hyperbolas are intersected in one point, which is the estimated position of LSP, as shown in Fig. 3.

Fig. 2.Localization geometry of RSS method

Fig. 3.Localization geometry of RSSR method

The equation of i-th sensor is obtained from (2) as follows:

where rei = ri + distance error for AWGN .

In RSSR method, one reference sensor is needed, which is denoted by S1 . The received signal strength ratio of i-th sensor with respect to reference sensor is,

To estimate the optimum horizontal distances re1, re2, re3,…,rei between LSP and different sensors from (4), we use here GA that is a search heuristic, which mimics the process of natural evolution. This heuristic is routinely used to generate useful solutions to optimization and search problems. GA belongs to the larger class of evolutionary algorithms, which generate solutions to optimization problems using techniques inspired by natural evolution, such as inheritance, mutation, selection and crossover [25]. Hence, GA is one of the best methods to do the global optimization [25]. Thereby, we have applied the GA to find out the best approximation of auxiliary vector. A simple GA procedure has been presented in Fig. 4.

Fig. 4.Block diagram of simple GA

3.1 Initial population

A GA is a probabilistic system that has established the standards of genetics. In this system, the first style is to depict the length of the genetic string named Chromosome. These Chromosomes having a quality of the destination role are termed as Chromosome’s fitness. After selecting the length of Chromosome, the initial samples of the Chromosome relies on the arbitrary collection of a number of chromosomes.

In this suggested issue, to determine the values of the horizontal distance, GA is applied, which permitted maximizing the value of (4). GA is designed as follows: the construction of the Chromosome is 30 bits separated into ith genes G1, G2,…,Gi presenting the horizontal distance of source and sensors. Those chromosomes are obtained by rearranging (3) as follows,

Squaring both left and right sides of (5) yields,

As , we can derive (6) as follows:

where Ai = Zi - z1, Fi = Zi - z2, , .

From (7), four roots can be derived. For the boundary condition of GA, it is very important to choose a correct root from those four roots. In the following way, the correct root can be found. If only one root is positive, then it can be used for the boundary condition of the GA. If all roots are positive, it selects one of which gives the almost theoretical value of Ei(r) in (1). If all are negative or imaginary, then the absolute

As an example, two binary strings of the Chromosomes are B1 and B2, where

Based on the fitness of individual population, a probabilistic choice is made at this stage. The program is stopped when the fitness value of any individual population reaches the expected value. Otherwise, it goes through the next step called crossover after selecting the best fit population.

3.2 Crossover

The genetic operator is utilized after the determination operation for generating new results depending on the current results in the sample. Two sorts of essential specialists can be found, one is a crossover and another is mutation. The input of the crossover operation is a couple of Chromosome called Guardians and the output gives a reproduced Chromosome known as Children. Sets associated with individuals were picked up arbitrarily from the mating pool according to the probability of crossover. Provided that the parallel letter set was utilized within GA then every spot was given transformation probability from 0 to 1. The variation in the quality is generally low.

The crossover is done in two steps for selected two Chromosomes: (1) the samples are represented by the binary strings and (2) according to the probability distribution, two Chromosomes do crossover in which binary strings can be switched. Therefore, single point crossover and two point crossover are used for these Chromosomes. By using uniform probability distribution, the crossover point is chosen. As an example, if 5 is the crossover site, the aforementioned B1 and B2 is changed after first step crossover and A1 and A2 will be

The sample uses the uniform crossover in the first step. In the step, binary strings are created and the common terms are searched in both parents, to find out the common mark in the common terms. In the string template, the binary bit will be set to 1 when there is matching term in both parents, otherwise, it will be zero. Again from the population, two numbers are selected randomly and exchanged with each other. As an example, if M1 and M2 are the sets of random number as

And, these M1 and M2 are made crossover with B1 and B2, the two parents Chromosomes are changed after the first step crossover. The new children will be like as

Second step crossover is done to increase the diversity and the new Chromosomes becomes

Within GA to reach the very least or perhaps the greatest position in search engine optimization, progress as well as procedure of selection is usually duplicated. When the fitness of the best Chromosome does not change significantly, the GA operation is completed.

3.3 Mutation

Changing a small number of bases randomly with small probability is known as mutation. It is performed for maintaining the diversity in the population. To avoid the parameter convergence, this operation is essential. The proper convergence of the GA is dependent upon its configuration, which is a very important issue. This convergence is made by defining the control parameters, like the population size, the Chromosome size, mutation rate and crossover point. Unfortunately, it is a very difficult task to determine these parameters and in most of the cases, they are defined empirically. From the population, each selected Chromosome was mutated separately, like crossover. The operator exchanges with a randomly selected term in the corresponding complimentary subset of the string. In this case, a simple point mutation is used. As an example, after mutation T1 will be as follows:

Parameters for the genetic algorithm configuration are:

● Size of the population: 700

● Rate of mutation: 0.35%

● Stopping criteria: 200

This whole process runs until reaching the best optimum value of the fitness function. The best population (best horizontal distance from the LSP to receivers) is obtained with respect to the best fitness value of the fitness function.

The obtained distance data vector from GA is,

And, the true distance vector is,

Using the maximum likelihood (ML) estimation, (8) is obtained where θ = [ x y ]T.

where L is the likelihood function and PDF is the probability density function. After some steps, (9) is obtained from (8),

From (9), we obtain the estimated position of the electromagnetic source. Let the actual coordinate of electromagnetic source be ( Xa,Ya ) and the number of independent trials in the simulation is 2500, then we obtain the root mean square (RMS) error of lightning location accuracy (LLA) by using (10).

 

4. Results and Discussion

In this section, we represent the results analysis and performance comparison between the proposed method and the RSS method based on AT model of CG return stroke lightning channel. In this simulation, we assume the receiver with high pass filter for picking up 1 MHz center frequency of the electric field of lightning channel [3], 125 km baseline of the receivers and the average height of the forest is from 10 m to 13 m [26]. We have simulated the proposed model and the existing RSS model of lightning at various peak current (2-20) kA of lightning channel under AWGN. For simulation, MATLAB is used in core i7 processor with 8 GB RAM. The simulation procedure is given below:

1) First step:

1.1 Select the initial value of horizontal distances between the LSP and receivers by using (7).

2) Second step:

2.1 Initial values are used as the boundary condition of GA.

2.2 Find out the best fit horizontal distances between the LSP and receivers from GA by using (4) as the fitness function of GA.

3) Third step:

3.1 Apply the ML into the best fit horizontal distances, as shown in (9).

3.2 Find out the estimated position of lightning strike.

4) Fourth step:

4.1 Error calculation by following (10).

The obtained results are illustrated in Fig. 5 after completing all the steps of the simulation.

Fig. 5.RMS positioning error (m) of LSP with respect to SNR (dB) at peak current values of lightning channel 2 kA, 10 kA, and 20 kA, respectively

In Fig. 5, the simulated results are shown at peak current 2 kA, 10 kA, and 20 kA, respectively, and at a frequency of 1 MHz of lightning channel. For the first scenario, the obtained error of LLA are 2.765 km, 2.605 km, and 2.489 km for the RSS method and 1.305 km, 1.256 km, and 1.185 km for the proposed method at peak current of lightning channel 2 kA, 10 kA, and 20 kA, respectively when the SNR is 10 dB. Next, the error of LLA becomes 1.120 km, 1.072 km, and 0.970 km for the RSS and 0.407 km, 0.3627 km, and 0.330 km for the proposed method when the SNR increases to 30 dB. Finally, the error of LLA are found to be 0.124 km, 0.108 km, and 0.097 km for the RSS and 0.0448 km, 0.0372 km, and 0.0331 km for the proposed method when the SNR is 50 dB at peak current 2 kA, 10 kA, and 20 kA, respectively.

The RMS position error comparison of LSP estimation for the proposed RSSR method and existing RSS method at 30 dB noise is represented in Table 1, where the numbers of sensors of network are varied. At 4 sensors, the error of LLA are observed 630.89 m, 557.45 m, and 523.51 m for the RSS method and 245.12 m, 232.05 m, and 201.47 m for the proposed method at peak current of lightning channel 2 kA, 10 kA, and 20 kA, respectively when the SNR is 30 dB. In addition, the error of LLA is reduced with the increment of network sensors due to the inconsistence reduction of non-linear equation.

Table 1.RMS positioning error (m) of LSP with respect to different number of sensors at 30 dB noise

The RMS comparison of LSP error (m) for the proposed method and RSS at noise 30 dB is represented in Table 2, where LSP is considered at the top of the forest and ground. Table 2 shows that, the RMS positioning error (m) of the LSP at ground is significantly greater than LSP at the top of the forest for the same pick channel based current. It is happening because the electrical signal strength of the LSP at ground is lower than the LSP at the tree.

Table 2.RMS positioning error (m) of LSP with respect to lightning strike at tree and ground for 30 dB noise

The RMS comparison of LSP error (m) for the proposed method and RSS at noise 30 dB for different ground conductivity is represented in Table 3. The RMS error increases with decreasing ground conductivity. It is happening because of decreasing amplitude of the field and increasing its rise time with increasing propagation distance and/or decreasing ground conductivity. The equation (1a) is utilized for the replacement of (1) to estimate the ground conductivity effects for locating the LSP.

Table 3.RMS positioning error (m) of LSP with respect to different ground conductivity for 30 dB noise

The recorded simulation time of the proposed method and RSS at noise 40 dB and 20 dB is represented in Table 4. The processing time of the proposed method is slightly greater than the RSS method due to little bit complex equation of the proposed method than that of the RSS. In addition, the simulation time is increased with the increment of noise because of high inconsistence data for high noise.

Table 4.Recorded simulation time of the proposed method and RSS at 40 dB and 20 dB noise

Finally, the proposed RSSR method to locate the position of the lightning is better than the existing RSS method, which is clearly justified in Fig. 5, Table 1 and Table 2 at low to high noise level and different sensors geometry. In practical field, the received field strength is increased with the increase of the peak current of CG lightning return stroke. For this reason, the noise effect of higher peak current signal is less than the low peak current signal. Therefore, the error of localization is less when the peak current is high. When the received signal is taken during a short period for analysis, the peak amplitude of the received signal must be greater than the peak value of theoretical received signal. After taking the ratio of practical received signal between two sensors, the noise effect must be mitigated, which is confirmed by our obtained results.

 

5. Conclusion

A hybrid model based on AT model of CG return stroke lightning channel in the central analyzer has been proposed in this paper. In this study, RSSR and RSS of receiver network are combinedly used for mitigating the noise effect in LLA. The proposed method is viable and achieves superior performance compared to the existing RSS technique. The proposed model can be easily implemented by changing the central analyzer program of the RSS based localization algorithm of CG return stroke lightning channel in the forest or other applications.

참고문헌

  1. Y. Wang, K. R. Anderson, “An evaluation of spatial and temporal patterns of lightning-and human-caused forest fires in Alberta, Canada, 1980-2007,” International Journal of Wildland Fire, vol. 19, no. 8, pp. 1059-1072, 2011. https://doi.org/10.1071/WF09085
  2. M. M. Muller, H. Vacik, G. Diendorfer, A. AepaciI, F. Herbert, H. Gossow, “Analysis of lightning-induced forest fires in Austria,” Theoretical and Applied Climatology, pp. 1-11, 2012.
  3. D. Fuquay, D. M. Fuquay, A. R. Taylor, R. G. Hawe, C. W. Scemid, “Lightning discharges that caused forest fires,” Journal of Geophysical Research, vol. 77, no. 12, pp. 2156-2158, 1972. https://doi.org/10.1029/JC077i012p02156
  4. M. A. Uman, The Art and Science of Lightning Protection, 1st ed., Cambridge University Press, New York, USA, 2010.
  5. C. Migli, D. Tao, D. Yaping, “Performance of TOA/DF Lightning Location Network in China-Site errors and detection efficiency,” in 7th Asia-Pacific International Conference on Lightning, 2011, pp. 48-54.
  6. K. L. Cummins, E. P. Krider, M. D. Malone, “The US National Lightning Detection Network and applications of cloud-to-ground lightning data by electric power utilities,” IEEE Transactions on Electromagnetic Compatibility, vol. 40, no. 4, pp. 465-480, 1998. https://doi.org/10.1109/15.736207
  7. H. Zhixiang, W. Yinping, Z. Wenguag, Z. Hongping, “Accuracy Analysis of the TDOA Method in a Lightning Location System,” International Conference on Management and Service Science, 2009, pp. 1-4.
  8. S. Xiubin, X. Zhang, X. Mou, Y. Li, “Five-channel VHF lightning detection method for cloud flashes,” in 7th Asia-Pacific International Conference on Lightning, 2011, pp. 878-881.
  9. H. Kwon, D. J. Pack, “Cooperative target localization by multiple unmanned aircraft systems using sensor fusion quality,” Optimization Letters, vol. 6, no. 8, pp. 1707-1717, 2012. https://doi.org/10.1007/s11590-011-0360-9
  10. P. Moravek, D. Komsny, M. Simek, D. Girbau, A. Lazabro, “Energy analysis of received signal strength localization in wireless sensor networks,” Radio engineering, vol. 20, no. 4, pp. 937-945, 2011.
  11. K. L. Cummins, M. J. Murphy, “An Overview of Lightning Locating Systems: History, Techniques, and Data Uses, With an In-Depth Look at the U.S. NLDN,” IEEE Transactions on Electromagnetic Compatibility, vol. 51, no. 3, pp. 499-518, 2009. https://doi.org/10.1109/TEMC.2009.2023450
  12. N. Bulusu, J. Heidemann, D. Estrin, “GPS-less low-cost outdoor localization for very small devices,” IEEE Personal Communications, vol. 7, no. 5, pp. 28-34, 2000.
  13. D. Pavaello, "Electromagnetic radiation from lightning return strokes to tall structures," Ph.D Thesis, Swiss Federal Institute of Technology, 2007.
  14. Y. Baba, V. A. Rakov, “On the interpretation of ground reflections observed in small-scale experiments simulating lightning strikes to towers,” IEEE Transactions on Electromagnetic Compatibility, vol. 47, no. 3, pp. 533-542, 2005. https://doi.org/10.1109/TEMC.2005.852220
  15. S. Bonyadi-Ram, R. Moini, S. H. H. Sadeghi, V. R. Rakov, “On representation of lightning return stroke as a lossy monopole antenna with inductive loading,” IEEE Transactions on Electromagnetic Compatibility, vol. 50, no. 1, pp. 118-127, 2008. https://doi.org/10.1109/TEMC.2007.913221
  16. R. Kumar, S. Ranade, B. Gowda, “An effective localization algorithm based on received signal strength,” In 2010 IEEE Aerospace Conference, 2010, pp. 1-8.
  17. A. S. Shoory, R. Moini, S. H. H. Sadeghi, V. R. Rakov, “Analysis of lightning-radiated electromagnetic fields in the vicinity of lossy ground,” IEEE Transactions on Electromagnetic Compatibility, vol. 47, no. 1, pp. 131-145, 2005. https://doi.org/10.1109/TEMC.2004.842104
  18. S. H. S. Moosavi, R. Moini, S. H. H. Sadeghi, “Representation of a lightning return-stroke channel as a nonlinearly loaded thin-wire antenna,” IEEE Transactions on Electromagnetic Compatibility, vol. 51, no. 3, pp. 488-498, 2009. https://doi.org/10.1109/TEMC.2009.2025655
  19. X. Wang, P. Hoole, E. Gunawan, “An electromagnetic time delay method for determining the positions and velocities of mobile stations in a GSM network,” Progress in Electromagnetics Research, vol. 23, pp.165-186, 1999. https://doi.org/10.2528/PIER98102603
  20. P. R. P. Hoole, S. R. H. Hoole, “Simulation of lightning attachment to open ground, tall towers and aircraft,” IEEE Transactions on Power Delivery, vol. 8, no. 2, pp. 732-740, 1993. https://doi.org/10.1109/61.216882
  21. V. Cooray, The Lightning Flash, 2nd Edition, IET, London, UK, 2014.
  22. J. R. Wait, “Transient fields of a vertical dipole over homogeneous curved ground,” Can. J. Phys., vol. 36, pp. 9-17, 1956.
  23. T. G. Wood, "Geo-location of individual lightning discharges using impulsive VLF electromagnetic wave forms," Ph.D Thesis, Stanford University, 2005.
  24. S. J. Orfanidis, Optimum Signal Processing: An introduction, New York: Macmillan, London: Collier Macmillan, 1988.
  25. P. Bajpai, M. Kumar, “Genetic Algorithm - An Approach to Solve Global Optimization Problems,” Indian Journal of Computer Science and Engineering, vol. 1, no. 3, pp. 199-206, 2010.
  26. S. C. Popescu, R. H. Wynne, “Seeing the trees in the forest: using lidar and multispectral data fusion with local filtering and variable window size for estimating tree height,” Photogrammetric Engineering and Remote Sensing, vol. 70, no. 5, pp. 589-604, 2004. https://doi.org/10.14358/PERS.70.5.589