地球电场（简称“地电场”主要由固体地球外部的电流体系在地球内部感应产生的分布于地表的快变化电场和地下介质中的物理、化学作用引起的慢变化电场两部分组成.地电场观测在资源勘探、地质调查等方面应用广泛，利用固定台站观测地电场的主要目的之一是进行地震等灾害事件的监测预报和地球电磁环境研究.20世纪初以来，国内外学者开展了固定台站的地电场观测和研究，对地电场时空变化的分布规律取得了重要认识^[1-8]，吉什和鲁尼在总结长期观测资料的基础上绘制了第一幅全球涡旋电流场分布图^[1] (1936年世界时18h).地电场观测在地震与火山等灾害事件监测预报实验中也得到广泛应用^[9-18].我国从“九五”后期开始了大空间范围、定点、连续、规范化的地电场观测^[19-20]，截止2011年我国已建成由100多个台站组成的地电场观测台网，其中部分台站持续观测时间已超过10年以上.国内进行的地电场观测是目前国内外地球电磁学观测的一大特色，目前关于地电场在地震监测预报实验中的研究工作较多，但关于地电场时空变化规律的研究甚至于国外都很薄弱.文献[21-23]应用我国部分地电场台站的观测数据研究了地电场日变化(以下简称“地电日变化”的时频分布特征，但数据应用不具中国大陆大区域的特点，研究成果有限、且认识有所不同，对于认识地电日变化的时空分布规律来讲是不够的.鉴于我国地电场台网在某些区域密度高、数据积累丰富，本文在文献[21-23]的研究基础上，应用中国大陆东部、西部、中低纬度、中纬度各两条经度链、纬度链的地电场台站观测数据，研究地电日变化的时空分布规律，期望能获得对中国大陆地电日变化的广域性、普遍性特点的进一步认识.

我国地电场台站观测装置是水平布设NS、EW两个正交观测方向再加一个斜交方向，长、短极距共6个测道^[20]，每测道采样率为1次/min，仪器分辨率为10 μV，24 h连续观测.为了研究地电场在中国大陆大空间区域的时、空变化规律，选取我国地电场台网中两条经度链、两条纬度链台站的观测数据.选取原则为：①台站分布范围尽量大，空间距离尽可能均匀；②尽可能选取经度链、纬度链上及近距离台站; ③台站观测数据质量合格.如图 1所示，沿118°E经线选14个台站组成经度链a，由南向北台站依次为:福建泉州、邵武，江苏海安、高邮，安徽嘉山，山东郯城、无棣，河北大柏舍、昌黎，北京通州，内蒙古宝昌，吉林榆树，黑龙江绥化、满洲里；沿99°E经线选12个台站组成经度链b，由南向北依次为：云南腾冲、洱源，四川泸沽湖、甘孜，青海大武、都兰，甘肃寺滩、古浪、山丹、高台、嘉峪关，宁夏石嘴山；沿40°N纬线选7个台站组成纬度链c，由西向东依次为：新疆乌什，甘肃瓜州、古浪、嘉峪关，山西大同，北京延庆、通州；沿34°N选12个台站组成链d，由西向东依次为:新疆和田，青海都兰、大武，甘肃陇南、静宁，陕西周至、乾陵，河南洛阳、周口，山东郯城，江苏高邮、海安.由于台站观测环境影响，某些近距离台站未人选.

图 1

Fig. 1

图 1 地电场观测台网中经纬链图示 Fig. 1 The longitudinal and latitudinal chains on geoelectric network

近年来，随着我国地电场台网的建设和台网规范化运行，执行了严格的观测质量监控程序，整体上数据质量有所提高，所以选取资料的时段主要为2010、2011年.个别台站受环境干扰，使用2008年以来其他年份的观测数据.

根据亥姆霍兹波茨坦中心（GFZ)公布的磁情指数，尽可能选取每月磁静日、且日期接近天数的观测数据.数据处理采用时序迭加方法，突出有规律性的日变化周期、抑制不规则噪声.时序迭加公式如下：

(1)

式(1)中X_n(t)表示所选时间范围内第n天的地电场数据，每天的数据长度为1440个(分钟值采样），k指第k天.

以经度链a上台站观测数据为主要研究对象. 图 2是链a中8个地电场台站每个台2010年2月3个磁静日（或一般磁静日）的观测数据时序迭加曲线.按纬度由低到高的顺序从上到下排列.由图明显可见：地电日变化主要表现为两峰一两谷变化波形，其中低纬度的泉州、海安、高邮台以及中低纬度的郯城、无棣和中高纬度的绥化台日变化波形表现为两峰一两谷波形；满洲里台表面上是中午前后1峰一1谷变化波形.

图 2

Fig. 2

图 2 2010年2月链a磁静日地电场数据时序迭加曲线 Fig. 2 The mean curves of geoelectric field on magnetically quiet day in Feb.2010

链a台站分布在25. 02°N-49. 57°N的纬度范围，由图 2似乎可见，总体上台站纬度高，日变化峰一谷幅度偏小(如绥化台纬度最高，北分量是0. 6 mV/km, 东分量0.4 mV/km); 纬度低幅度偏大(如泉州台纬度最低，北分量是20 mV/km, 东分量15 mV/km).且似乎午前峰一谷幅度略大于午后.

由图 2可见，链a各台站在午时（即12h左右）日变化出现波峰和波谷.部分台站之间存在少许相位差，是由于链a台站并不全在118°E经线上，位于此经线附近，经度偏东的台站相位早于经度偏西台站.例如，绥化台经度为126. 9°E、满洲里台117.43°E，绥化台波峰、波谷极值时间早于满洲里台，而满洲里、郯城、泉州台（自东向西顺序）经度接近，基本无相位差.

选取链d台站观测数据，并按劳埃德（Loyd)天文季节的D季节（11、12、1、2月）、季节（5、6、7、8月）和E季节（3、4、9、10月）进一步研究日变化随着地方时的变化（即“经度效应”）.在2011、2008、2010年中不同的天文季节中各选一个磁静日，比较链d台站在同一天的地电日变化波形（图 3).图中各台依从东到西的顺序上下排列，这些台站分布在79. 9°E-120.5°E之间，最大台间距超过3700 km.通过图中的辅助线可清晰地看到日变化相位有经度效应，随着经度偏西，日变化相位依次滞后.例如，2011年1月5日和2008年3月7日的日变化相位的经度效应最明显，2010年8月22日较明显，但不如前两个图（因数据质量影响）.从1月5日地电场日变化可明显看出：工苏海安、新疆和田台第二个波谷到达的时差约为2.7h，与两台站经度差（约40. 6°)对应的地方时差（约2. 7h)相符（见图中辅助线）；高邮一和田、郯城一和田波形相位约2. 6h，吻合台站当地的时差；其次，海安一周至、周至一都兰、都兰一和田台时差依次分别为0. 82、. 9、1h，日变化相位差依次分别为0. 82、. 08、1h.显然，地电场日变化两次起伏的相位存在明显的经度效应，相位差与台站之间当地的时差很吻合.

图 3

Fig. 3

图 3 链d台站地电场相位随经度的变化 Fig. 3 The phasic changes of geoelectrical field along with chain 'd'

同时还注意到，地电场快变化成分是相当同步的，无时差（见2011年1月5日图 0~240 min和1320 min前后、2010年8月22日图 0~120min).实际上，在数据处理中注意到此类快变化现象是普遍存在的，特别是在地电暴期间.说明地电日变化相位存在经度效应，而快变化同步，似乎反应了二者的“源”因素有区别.

分析了经度链b、纬度链c台站地电日变化波形、经度效应，其结果与经度链a、纬度链d的结果类似.另需说明，链a_d个别台站地电日变化波形、经度效应不很符合前面的结果，这种情况可能与台站观测环境、技术系统运行等因素有关; 个别台站日变化波形还表现为1峰一1谷的形态.其次，极个别台站观测数据不可用，文中未应用.

在地电场频率域变化的分析中，用快速傅里叶变换方法（FFT)处理了4条经、纬链上共14个台站2010-2011两年内116个月的地电场分钟值数据.由于日变化是频率极低的成分，为了客观获得日变化的周期成分，每个台站选取3~5个月的地电场北分量E_x和东分量E_y观测数据，且在选取数据时段的每个月各分量分钟值数据组成1个长度约为43200的时间序列，使时间序列长度内能足够包涵多个日变化周期成分.

图 4a是嘉山台2010年3月地电场E_x和E_y观测数据的FFT振幅谱.从图中可以清楚地看到, 优势周期的幅值从大到小排列依次为:12.4、8、24h、12、6、4.8h, 其中12.4h为最优势周期.在对其他台站数据的处理过程中, 较明显的优势周期中还出现了25.6、4h周期, 但其幅度明显小于前面几个周期成分的幅度.处理其他月份的观测数据, 其结果类似图 4a.图 4b为2010年全年嘉山台地电场E_x和E_y观测数据FFT振幅谱时频图.由图可见:① 6h、8、12.4、24h周期成分的谱幅度相对其他周期成分占优势, 其中12.4h周期成分的谱幅在全年显著, 24h周期成分在夏季比其他季节显著; ②夏季幅值大, 冬季幅值小, 特别是E_y分量更明显.

图 4

Fig. 4

图 4 嘉山台地电场日变化FFT振幅谱(上图:NS测道; 下图:EW测道) (a)2010年3月地电场FFT谱; (b)2010年地电场FFT谱时频图. Fig. 4 The amplitude spectrums of geoeletrical field on Jiashan station (the above: channel-NS, the below:channel-EW) (a) The results of the data in Mar.2010;(b) The time-frequency diagrams of the data in 2010.

表 1是2011年1、2、3、7、10共5个月, 链a的泉州、海安、高邮、嘉山4个台站两个测道的地电场数据FFT振幅谱统计的优势周期.除表 1中的4个台外, 还对链a的郯城、无棣、绥化、满洲里4个台观测数据做了谱分析(篇幅所限未列出).8个台每台每月各分量谱分析结果记为一组, 共计80组, 在振幅谱中可明显识别的优势周期共有333个.表 1中每组数据依谱幅从大到小排列.由表可见, 12.4h、12h以及24h周期为优势周期, 而12.4h周期为优势周期中出现频次最高的周期成分.其次8、6h等周期成分也多次出现, 还间或出现约25.6h的周期成分.同时注意到, 除泉州台1、2、3月外, 其余台站E_x、E_y两个分量优势周期相同或接近, 其次不同月份优势周期也相同.

表 1(Table 1)

表 1 链a四台站地电场数据振幅谱优势周期 Table 1 Preferred periods in amplitude spectrums of 4 geoeletrical field stations

表 1 链a四台站地电场数据振幅谱优势周期 Table 1 Preferred periods in amplitude spectrums of 4 geoeletrical field stations

表 2是对链a、链b台站地电场数据FFT振幅谱的统计(链b取腾冲、泸沽湖、古浪、石嘴山、山丹、高台6个台2011年2、8、10月的数据）.可以看出链a台站每月每台站出现8 h和12 h周期成分的频次最高，按频次统计组发生率分别为87. 5%和80%，6、24、12. 4h三种周期的组发生率相近，也都在60%以上，其次还有4.8、25. 6h和4h三种周期成分，但明显不及前5个周期成分出现的频次高.其他周期的组发生率的总和仅为13. 8%.因为地电场观测对环境干扰、电极稳定性、线路绝缘以及地下介质电性变化等因素敏感，易发生数据异常，少数不常见的周期成分可能是由于以上因素的影响.需指出，由于数据质量和谱估计误差，统计中提到的某一周期是指此周期±0.1 h范围(如12.4 h指12.4 h±0.1 h).与链a相比，链b台站每组FFT谱中可分辨的优势周期的个数平均为3. 11个，低于链a台站每组平均4. 16个，可能因为链b位于我国西部，台站环境艰苦、观测数据质量整体上不如链a台站.如云南腾冲、四川泸沽湖台存在观测场地环境干扰，可分辨的周期较少或无优势周期.链b台站主要周期按频次依次为12、8、24h，组发生率在60%以上，其中12h组发生率最高，为80. 6%;其余组发生率较高的周期成分为6、4 8h与链a不同的是，在链a台站出现12. 4h周期成分的频次高，组发生率高达62. 5%，而链b台站仅为2. 78%，25. 7h周期的频次也仅为2. 78%.

表 2(Table 2)

表 2 链a、链b、地电场优势周期统计表 Table 2 The statistical table of preferred periods of geoelectric field in chains‘a’ andchain ‘b’

表 2 链a、链b、地电场优势周期统计表 Table 2 The statistical table of preferred periods of geoelectric field in chains‘a’ andchain ‘b’

按最大FFT谱幅出现的频次统计（表 3)，在链a的80组数据（即80个月数据）中，12. 4h周期成分的频次最高，为42. 5%，其余依次为12、8、24h.说明在链a的8个台站中，地电日变化主要是由12. 4、12、8h和24h周期成分组成，其中12. 4h周期幅值最显著，虽然4. 8、25. 7、4h周期波也较常见，但非普遍存在、振幅不大；连b台站36组数据中最大谱幅的周期成分出现频次最高为12h，其次是8h和24h周期成分.

表 3(Table 3)

表 3 链a和链b地电场FFT最大振幅周期频次统计表 Table 3 The statistical table of the most preferred periods of geoelectric field in chains ‘a’ and chain ‘b’

表 3 链a和链b地电场FFT最大振幅周期频次统计表 Table 3 The statistical table of the most preferred periods of geoelectric field in chains ‘a’ and chain ‘b’

总之，在两个经度链上地电日变化可分辨的优势周期按周期成分的振幅大小主要有12. 4、12h、8、24、6、4. 8、4、25. 7h，但两条经度链上各自的优势周期成分又有区别.链a台站12. 4h周期最显著，链b台站12h周期成分最优显著.

对郯城2010年1月6-8日三个磁静日的地电场观测数据进行时序迭加处理，提取其中振幅最大的前5阶谐波进行反变换拟合.图 5中灰色曲线为原始曲线，黑色为拟合曲线，拟合曲线很好地近似了原始曲线.前5阶谐波振幅依次为1. 67、1. 01、0.88、0.33和0. 26 mV，其中12h和8h周期成分的振幅最大.

图 5

Fig. 5

图 5 郯城台地电场N测道2010年1月份3天静日平均后的谐波拟合 Fig. 5 The harmonic waves synthesis of 3-quiet-day average of Tancheng station channel Nin Jan.2010

文献[24]给出了用FFT处理的美国图森台1个月地磁场Sq日变化的主要周期，按幅值大小其主要周期成分依次为24 h、12 h和8 h周期等，说明地磁场Sq日变化主要周期成分为每日一次起伏.本文用每台、每测道1个月的数据组成长时间序列，其周期成分估算应准确，按幅值大小分别为12.4、12h、8、24 h，与文献[21-22]用最大熵方法估算的昌黎、宝坻、嘉峪关台地电日变化主要周期成分为12、24 h和8 h左右周期的结果吻合，说明地电日变化每日两次起伏为主，符合图 2-3中地电日变化波形.需说明：及个别台站日变化周期不是12. 4h/12h，例如满洲里台主要为8 h周期.

由于地电日变化主要为每日两次起伏，故对日变化幅度的分析针对12 h/12.4 h半日波.选取链a台站2010年每个天文季节中1个月磁静日3天数据统计日变化峰一谷幅度.选取2、8、10月代表D、J、E季节，并在图 6中分别用黑、红、蓝三种颜色分别表示（若某月出现数据异常，用相同天文季节内其他月份数据代替）.图 6a是链a的11个台E_x分量随纬度的变化，图 6b是E_x分量随纬度的变化，图 6c是各台站E_x、E_y分量分别平均的结果.由图可见，链a台站从约1~23 mV/km不等，且随纬度升高大致的逐渐减小，不过图 6a中E_x结果有些凌乱.根据电阻介质的欧姆定律J=σ · E(电导率σ=1/ρ)，电场强度与介质电阻率ρ有关.事实上，我国地电场观测的电极埋设在表层3 m左右深度，在欧姆定律中用地表视电阻率ρ_s代替ρ即可.文献[21]研究了“九五”几个地电场台站的与台址表层ρ_s的关系，得到：表层ρ_s大，则大.图 6中随纬度增高呈现略凌乱的现象是台址浅层电导率不同引起的.例如满洲里、泉州、宝昌台表层介质电测深视电阻率依次为150~200 Ωm，大于无棣、高邮、郯城、海安、通州、绥化、榆树、大柏舍等台站几Ωm~60 Ωm，宝昌、满洲里台较大、离散，无棣、高邮台小、离散.海安与郯城、通州与绥化、大柏舍与榆树台表层电阻率相差无几，则纬度高、日变幅度小.

图 6

Fig. 6

图 6 2010年经度链a地电场峰-谷幅度 Fig. 6 The changes of amplitude of geoelectrica lfield in chains "a" in 2010

同时，由图 6(a-b)得到，J季节的红色曲线高于E、D季节的绿色和黑色曲线，而D季节黑色曲线最低.分别统计图 6(a-b)展示的三个季节各台站的平均幅度得到：在J季节，E_x平均幅度为7. 65mV/km，E_y为6. 22 mV/km; 在E季节，E_x为6.28 mV/km，E_y为4. 92 mV/km; 在D季节，E_x为5.14 mV/km，E_y为4.15 mV/km..显然，在天文季节的夏季大，春秋季节其次，冬季最小.这个结果与文献[1]应用当时少数几个台站观测数据的研究结果一致.另外注意到：无论在哪个天文季节，都是地电场北向分量E_x日变幅度大于东向分量E_y的.

选取链b的12个台站观测数据，按链a方法统计的纬度效应得到，从约N35°开始随纬度增高、降低大致呈现逐步减小的现象(篇幅所限图未示）.据近年来的全国地电场台网观测资料质量评比成绩，链b台站观测质量整体上低于链a台站，且链b台站分布在南北地震带构造活动区，影响了对正常地电日变化的分析效果，与链a结果相比链b可应用数据的台站少，统计结果较为离散、凌乱，但总体来看链b台站从N35°开始随纬度增高、降低，而分别呈现逐步减小的现象.

应用纬度链d自西到东11个台站的地电场E_x、E_y，分量观测数据，按链a方法分析了各台站的，其中9个台站的集中在约1~11mV/km之间起伏变化，不呈现与经度有关的规律性变化.

吉什和鲁尼据全球地电场数据绘制了1936年世界时18时的全球涡旋地电流分布图的概貌^[1]，在北半球的中低纬度区、高纬度区各分布4个涡旋场，中低纬度涡旋中心在30°N附近，在高纬度区接近70°N附近，各涡旋分布相对于太阳的位置固定.其中，中低纬度涡旋中心大致在25°N、120°E和35°N、70°E的2个涡旋电流场覆盖或邻近中国大陆（在世界时18时），东部链a、西部链b两个子午链台站分布在25°N -49. 57°N的中低一中纬度范围，两个涡旋电流场覆盖链a、链b台站.由于各涡旋分布相对太阳的位置固定不变，在一天当中世界时的另外时刻，链a、链b台站还要经过另外两个中低纬度区的涡旋电流场•对于链a，泉州（25. 02°N、118. 51°E) -满洲里(49.57°N、117.43°E)共14个台站日变化幅度随纬度升高总体上逐步减小；对于链b，腾冲(25°N、98. 50°E)-嘉峪关（39. 80°N、98. 22°E)共12个台站从大致35°N开始随纬度增高和降低而逐步减小；且纬度链台站地电日变化的相位差相当对应当地的地方时.上述日变化幅度的表现似乎说明，位于中低纬度区的涡旋电流场分布与链a、链b台站的日变化幅度不无关系.普遍认为引起地电静日变化的电流系主要分布在中、低纬度区上空的电离层中，扰日变化与高纬的涡旋电流相对应^[7].在地电场变化的分析中，我们注意到地电日变化相位差明显对应地方时的时差，但快变化在纬度链上大空间范围（超过3700 km台站间距）同步变化，说明日变化和快变化似乎存在影响因素的差异，可能是在太阳风的作用下，电离层中的电流系引起了地面涡旋电流场电流强度的变化，造成了随纬度的变化.

地电场Sq和低磁情日变化主要为12. 4h (链a)、12h (链b)、8h和24h周期成分以及6、4.8h、4h等成分，其中12. 4h周期是主太阴半日分潮(M2波），12h周期是主太阳半日分潮（S2波），24h是太阳赤纬全日分潮（P1波），8、6、4. 8、4h也是潮汐调和分量.在这些周期成分中，24、12、8、6、4. 8 h成分与变化地磁场Sq谐波成分一致，文献[24]关于变化磁场中的上述谐波分量已做了详细研究，认为这些成分与中高层大气及电离层中的潮汐发电机理论息息相关.月球、太阳依次是地球潮汐现象的主要贡献者，链a台站地电日变化最主要是M2波周期，链b台站最主要是S2波周期，这些周期成分是主要的潮汐调和分量，因此认为地电日变化的主要周期与月日引潮力直接有关.半日潮每天两次高潮和两次低潮，全日潮每天一次循环，在引潮力作用下固体地球、大气层、电离层、海洋、湖泊等产生的潮汐现象引起了地面电流涡旋场强度发生周而复始的变化，引起了上述地电日变化的主要周期成分，即以12.4h、12、24h周期成分的幅度为最显著的日变化.

文献[23]研究了潮汐对岩石中裂隙走向的优势趋向作用和裂隙水渗流对地电场潮汐响应的影响.地壳介质中普遍存在裂隙诱发的各向异性^[25]，实际上该文讨论了介质中微裂隙走向优势趋向（或定向排列）构成的各向异性介质中传导电流引起的地电场变化的情况，导电机制主要是介质中的动电效应，导电溶液运移和微裂隙、孔隙是重要因素.地电场观测地表两个水平正交分量，本文应用上述微裂隙构成的方位各向异性介质进一步讨论月日引潮力对地电日变化的影响.在这种介质中欧姆定律是张量表示式，在微裂隙优势趋向方向及其垂直方向（即各向异性介质的两个水平电性主轴方向）的电场分量强度只与各自的电阻率ρ电流密度有关.地电场测量电极埋深在3 m左右的土壤介质中，欧姆定律中和.文献[26]给出了这种介质中视电阻率变化与饱水的固体介质骨架电阻率ρ₀、裂隙水电阻率ρ_f(一般情况下)和裂隙率ν及其变化Δν/ν的本构关系.据此本构关系，当无附加载荷时，介质无变形（即Δν/ν→0)，则介质视电阻率不变化（即Δρ_s/ν_s→0)，相应地介质中不变、亦不变；当月日正对地球某区域时引潮力最大，介质膨胀、微裂隙“张开”或产生新增微裂隙，即Δν/ν > 0, 介质中含电粒子、电荷等的水溶液等沿微裂隙优势取向方向运移，同时导电通道链通、Δρ_s/ν_s负向增大使ρ_s减小，介质中产生附加电流，此时电场强度将出现一次极值过程；在引潮力最小的区域，介质被压缩、微裂隙“闭合”，即Δν/ν < 0, 水溶液等不易运移，同时导电通道阻塞、Δρ_s/ρ_s正向增大使ρ_s增大，此时将出现一次相反的极值过程.于是，一天中顺潮、背朝的两次潮汐过程造成了涡旋电流强度周而复始的变化，造成了周期12.4、12h精确对应M2波、S2波的地电日变化半日波成分.

另一方面，空间电流系也是影响地电日变化的重要因素.地电日变化的午前午后表现类似于变化磁场Sq，据文献[24, 27]变化磁场Sq午前午后的不对称现象起因于低纬电流系和高纬电流系.其次，比较分析地电日变化波形与钻孔应变、静水位、重力等观测的半日潮波形可知，地电日变化与应变等半日波存在明显的相位差和视周期差异.这种差异表明电离层、大气层等的电动力学过程也是影响地电日变化的重要原因.因此，本文认为月日引潮力引起的地面涡旋电流场电流强度的变化和太阳风等引起的空间电磁活动共同产生了地电日变化的半日波成分.

链a台站一般位于我国东部海拔低于500 m的区域，且邻海/湖等水系居多，日变化半日波周期为12. 4h，而链b台站一般位于海拔约1500 m及以上的高原区，半日波周期为12h.其周期差异的原因可能是链a台站地下介质内部水网络发育，月球引潮力的作用更明显，对于链b台站空间电磁活动的作用更明显.

(1) 沿两个经度链、纬度链台站的地电日变化绝大多数表现为两峰一两谷的变化波形，波峰、波谷出现在当地时午前午后时段，一般午前日变化幅度略大午后.

(2) 地电日变化周期成分按最大幅度出现的频次排序依次为12. 4、12、8、24、6、4. 8、4、25. 7h周期等，这些周期是潮汐调和分量.其中在118°E链台站以主太阴半日波(M2波)12. 4 h为最主要周期，99°E链以主太阳半日分波（S2波）12h为最主要周期.

(3) 纬度链地电日变化半日波相位差与台站之间的地方时的时差一致，表现为经度偏东台站早于偏西台站的经度效应.

(4) 地电日变化半日波周期成分的幅度有纬度效应，总体上，118°E链从福建泉州一内蒙古满洲里台日变化幅度逐渐减小，99°E链从约35°N开始大致随纬度分别增高和降低而日变化幅度减小；Loyd天文季节的夏季日变化幅度最大，春秋季节其次，冬季最小；绝大多数台站地电场北分量艮日变化幅度大于东分量E_x；地电日变化幅度与台址浅层电阻率存在正比关系.

(5) 月日引潮力引起的地面涡旋电流场的电流强度变化和太阳风引起的大气层、电离层电磁活动共同产生了地电日变化的半日波周期成分.

十分感谢两位审稿专家和编辑部对本文付出的辛勤工作，审稿专家对文中逐句修改，并提出多条具体的修改意见和建议，体现了科学、严谨的学风.