پیاده‌سازی سه مجموعه روابط پارامتری انتقال بار الکتریکی در مدل یک‌بُعدی ابر

نویسندگان

1 استادیار، پژوهشکده علوم جوّی، پژوهشگاه ملی اقیانوس‌شناسی و علوم جوّی، تهران، ایران

2 استادیار، گروه فیزیک فضا، موسسه ژئوفیزیک دانشگاه تهران، ایران

چکیده

از آنجا که نتایج شبیه‌سازی میدان الکتریکی درون ابر به مقدار زیادی به روابط مورد استفاده در پارامترسازی انتقال بار الکتریکی حساس است، در این پژوهش سه مجموعه روابط پارامتری پیشنهاد شده تاکاهاشی در 1978 و 1984 (TAK)، جایاراتنه و همکاران در 1983، گاردینر و همکاران در 1985 و زیگلر و همکاران در 1991 (JGZ) و سندرز و همکاران در 1991 (SAN) کُدنویسی و در مدل یک‌بُعدی ابر قائم پیاده‌سازی شد.
نتایج حاصل از بررسی متوسط بار الکتریکی شبیه‌سازی شده با کاربست هر سه مجموعه روابط پارامتری نشان داد که الگوی دوقطبی بار شبیه‌سازی شده برای TAK و  JGZبه‌صورت دوقطبی مثبت و برای SAN به‌صورت دوقطبی منفی است. همچنین میدان الکتریکی حاصل از اجرای طرحواره TAK بین دقیقه‌‌‌‌های 25 تا 56، برای طرحواره JGZ بین دقیقه‌‌‌‌های 19 تا 34 و برای طرحواره SAN بین دقیقه‌‌‌‌های 27 تا 57 شکل گرفته است. مقایسه زمان و ارتفاع‌‌‌‌های وقوع بیشینه میدان الکتریکی مثبت و منفی شبیه‌سازی شده روشن ساخت که برای مجموعه روابط پارامتری TAK و SAN این مقادیر نزدیک به هم است، درحالی‌که برای JGZ سه دقیقه اختلاف زمانی و 5/12 کیلومتر اختلاف ارتفاع مشاهده شد. در ادامه با مقایسه بین شدت میدان الکتریکی شبیه‌سازی شده و میدان الکتریکی آستانه، مشخص شد که از بین سه مجموعه روابط پارامتری، TAK بیشترین و JGZ کمترین رویداد آذرخش را شبیه‌سازی کرده است.

کلیدواژه‌ها


عنوان مقاله [English]

Implementation of three sets of electric charge transfer parameterization in a one-dimensional cloud model

نویسندگان [English]

  • Pegahfar N 1
  • Gharaylou M 2
چکیده [English]

 Results of numerical simulation of intra-cloud electrification depend on a mechanism that determines the sign and magnitude of charge transferred to hydrometeors (including graupels and ice crystals), through their collision (Mansell et al., 2005). Some of the microphysical processes play an important role in the above mechanism. In order to estimate the amount and sign of the transmitted charge per collision for numerical purpose, the results of laboratory researches are commonly used. Two kinds of inductive and non-inductive (used in the current paper) mechanisms could be applied. The research studies conducted for the non-inductive one can be determined based on the liquid water content (LWC), temperature (T), ice accretion rate and the particle size spectrum (Takahashi, 1978; Jayaratne et al., 1983; Gardiner et al., 1985; Saunders et al., 1991; Ziegler et al., 1991; Saunders and Peck, 1998; Pereyra and Avila, 2002).
According to the importance of this issue, in this research three sets of relations resulted from laboratory studies have been used. These sets were proposed by Takahashi (TAK, 1978 and 1984), Jayaratne / Gardiner / Ziegler (JGZ, Jayaratne and colleagues 1983; Gardiner et al. 1985; Ziegler et al. 1991) and Sanders et al. (SAN, Sanders et al. 1991). These parameterizations relate the mean charge transferred per collision to liquid water content and temperature. Following these studies, the prepared schemes of three sets have been implemented in an explicit time-dependent one-dimensional cloud model (ETM), based on Chen and Sun (2002). In the 1-D cloud model entrainment-detrainment and eddy diffusion processes have been considered. Also microphysical processes have been parameterized using Lin et al. (1983) and Rutledge and Hobbs (1984) schemes. The convection is initiated using potential temperature perturbation, defined by Chen and Sun (2004). The input data for the simulation of vertical cloud is from an idealized sounding including pressure, temperature and water vapor mixing ratio. This cloud model simulates vertical velocity (w), equivalent ice potential temperature (θei), water vapor mixing ratio (qv), cloud water mixing ratio (qc), ice mixing ratio (qi), rain water mixing ratio (qr), snow mixing ratio (qs) and graupel mixing ratio (qg). The cloud model was set up with 1 second time step, 70-minutes simulation duration and 250 m spatial resolution in the vertical direction up to a height of 15 km.  The initial radius for the cloud column was considered as 3000 m.
The results of simulated mean charge transferred per collision using three sets of parameterizations (TAK, JGZ and SAN) show that their dipole pattern outputs are not the same. Simulations based on TAK and JGZ relations produced positive dipole (positive charge distribution lies over negative charge one), while simulation using SAN parameterization produced negative dipole pattern (the negative charge distribution on the top of the positive charge one). The simulation results show that the electric field was produced between 25-56, 19-34 and 27-57 minutes using TAK, JGZ and SAN parametric relations, respectively. It is noteworthy that the maximum values of positive and negative intra-cloud electric fields were obtained when applying TAK relations in the charge-transfer simulations. While, simulations using JGZ and SAN parameterizations led to the minimum values for positive and negative intra-cloud electric fields.      
The time and height of positive and negative electric field occurrences based on three sets of applied parameterizations were also compared. The results of comparisons demonstrated that the values acquired from TAK and SAN parameterizations were close. However, the values for simulation using JGZ indicated three minutes time discrepancy and 12.5 km height difference.
Finally, the simulated intra-cloud electric fields using three TAK, JGZ and SAN parameterization sets were compared with the threshold electric field, defined by Marshall et al. (1995), to extract the number of lightning occurrences. Our findings show that the maximum and minimum values of lightning events were seen in simulations using TAK and JGZ parametric relations, respectively. The number of lightning occurrences was 40, 12 and 30 for simulations using TAK, JGZ and SAN parameterizations respectively.    

کلیدواژه‌ها [English]

  • One dimensional cloud model
  • Intra-cloud charge transfer schemes
  • Lightning
  • The intra-cloud electric field
  • TAK
  • JGZ
  • SAN