Next Article in Journal
Prevalence and Risk Factors Associated with Musculoskeletal Discomfort in Spay and Neuter Veterinarians
Next Article in Special Issue
Unusual Childhood Waking as a Possible Precursor of the 1995 Kobe Earthquake
Previous Article in Journal
Modelling Niche Differentiation of Co-Existing, Elusive and Morphologically Similar Species: A Case Study of Four Macaque Species in Nakai-Nam Theun National Protected Area, Laos
Previous Article in Special Issue
Possible Electromagnetic Effects on Abnormal Animal Behavior Before an Earthquake
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Early Results of Three-Year Monitoring of Red Wood Ants’ Behavioral Changes and Their Possible Correlation with Earthquake Events

1
Department of Geology, Faculty of Biology, University Duisburg-Essen, 45141 Essen, Universitätsstr. 5, Germany
2
Am Plexer 7, 50374 Erftstadt, Germany
3
Image Analysis Group, Faculty of Electrical Engineering and Information Technology, TU Dortmund, Otto-Hahn-Straße 4, 44227 Dortmund, Germany
*
Author to whom correspondence should be addressed.
Animals 2013, 3(1), 63-84; https://doi.org/10.3390/ani3010063
Submission received: 14 December 2012 / Revised: 25 January 2013 / Accepted: 28 January 2013 / Published: 4 February 2013
(This article belongs to the Special Issue Biological Anomalies Prior to Earthquakes)

Abstract

:

Simple Summary

For three years (2009–2012), two red wood ant mounds (Formica rufa-group), located at the seismically active Neuwied Basin (Eifel, Germany), have been monitored 24/7 by high-resolution cameras. Early results show that ants have a well-identifiable standard daily routine. Correlation with local seismic events suggests changes in the ants’ behavior hours before the earthquake: the nocturnal rest phase and daily activity are suppressed, and standard daily routine does not resume until the next day. At present, an automated image evaluation routine is being applied to the video streams. Based on this automated approach, a statistical analysis of the ant behavior will be carried out.

Abstract

Short-term earthquake predictions with an advance warning of several hours or days are currently not possible due to both incomplete understanding of the complex tectonic processes and inadequate observations. Abnormal animal behaviors before earthquakes have been reported previously, but create problems in monitoring and reliability. The situation is different with red wood ants (RWA; Formica rufa-group (Hymenoptera: Formicidae)). They have stationary mounds on tectonically active, gas-bearing fault systems. These faults may be potential earthquake areas. For three years (2009–2012), two red wood ant mounds (Formica rufa-group), located at the seismically active Neuwied Basin (Eifel, Germany), have been monitored 24/7 by high-resolution cameras with both a color and an infrared sensor. Early results show that ants have a well-identifiable standard daily routine. Correlation with local seismic events suggests changes in the ants’ behavior hours before the earthquake: the nocturnal rest phase and daily activity are suppressed, and standard daily routine does not resume until the next day. At present, an automated image evaluation routine is being applied to the more than 45,000 hours of video streams. Based on this automated approach, a statistical analysis of the ants’ behavior will be carried out. In addition, other parameters (climate, geotectonic and biological), which may influence behavior, will be included in the analysis.

1. Introduction

One of the major challenges of seismology is to predict earthquakes. However, despite applying modern space and/or ground based technologies, a reliable method for short-term earthquake forecasting with reliable warnings several hours or days in advance is currently not possible and remains limited to only a few minutes before the quake happens. During the past decade, earthquake precursors, like ionospheric perturbations, electromagnetic phenomena, ground heating and composition changes of soil-efflux (radon, CO2, etc.), have been scientifically investigated. However, they are far from being successfully applicable, because these precursors are highly regionally variable depending not only on the large-scale tectonic main stress regime, but also on the small-scale tectonic regime [1,2,3].
Abnormal animal behaviors have been reported previously, but create problems in monitoring and reliability [4,5,6,7]. However, long-term in situ data collection, statistically significant validation of these precursors in focus regions over a long period of time and plausible scenarios explaining the evolution of such behaviors are still missing.
The situation is different with red wood ants (RWA; Formica rufa-group (Hymenoptera: Formicidae). A particular advantage of RWA is their high sensitivity to environmental changes. Forced by some million years of evolutionary selection, they may have developed anticipatory mechanisms [8,9]. Besides an extremely strong temperature sensitivity (0.25 K) [10], they have chemo-receptors for the detection of CO2-concentrations [11,12] and an electromagnetic field sensitivity [13,14,15]. Recent research in the West Eifel (West Germany), combined with geostatistical analyses, have demonstrated the correlation of soil gas anomalies, especially helium anomalies, as outstanding fault zone tracers and spatial distribution of red wood ant (RWA) mounds along tectonically active, gas-permeable faults (cf. Figure 1) [16,17].
Figure 1. GeoBioScience: relation of red wood ant (RWA) mound sites and gas-permeable tectonic fault systems. Along the open systems, geogenic gas concentrations exceed atmospheric standards (He > 5,220 ppb) and background values (CO2 > 2 Vol%; Rn up to 110 Bq/L). Linear alignment of RWA mounds occurs if the gas-permeable fault is running open for a very long distance, e.g., 3 km at the Oberehe site (West Eifel, cf. Figure 2(a)). Here, mound numbers of 120 mounds in a line can be mapped. At intersections of two or more tectonic faults, clusters of RWA mounds occur. Solitary mounds can be found mostly at unusual locations with spotty degassing [16].
Figure 1. GeoBioScience: relation of red wood ant (RWA) mound sites and gas-permeable tectonic fault systems. Along the open systems, geogenic gas concentrations exceed atmospheric standards (He > 5,220 ppb) and background values (CO2 > 2 Vol%; Rn up to 110 Bq/L). Linear alignment of RWA mounds occurs if the gas-permeable fault is running open for a very long distance, e.g., 3 km at the Oberehe site (West Eifel, cf. Figure 2(a)). Here, mound numbers of 120 mounds in a line can be mapped. At intersections of two or more tectonic faults, clusters of RWA mounds occur. Solitary mounds can be found mostly at unusual locations with spotty degassing [16].
Animals 03 00063 g001
At the reference location Oberehe (cf. Figure 2(a)), a site with small soil/sediment covers, a direct degassing of geogenic gases into the atmosphere occurs. The linear alignment of the RWA mounds directly corresponds to that of the helium anomalies. The Hough transform (a well-established algorithm for identifying linear structures in sets of points [18]) applied on the spatial distribution of RWA mounds (cf. Figure 2(b)) in the West Eifel (1,140 km²) and the directions of fault zones published in [19,20] (cf. Figure 2(c)) shows that the centers of the modes, i.e., the local maxima in the histograms, denote the preferential alignment directions, while their widths indicate the corresponding directional variations. In the RWA histogram, the local maxima are separated by vertical black lines. For each mode, the average value is indicated by a solid vertical red line and the ±1 standard deviation interval surrounded by two dashed vertical red lines. More specifically, the azimuthal directions in which the RWA mound positions are aligned correspond to those of the main fault directions. The comparison of both the spatial RWA mound distribution with the significantly increased gas anomalies confirms that the RWA colony is situated in a tectonically sheared region. Although the reference site is surrounded by various forest stands, no further RWA mounds were found in those forest stands during the area-wide mapping. The statistical analyses show that the spatial distribution of RWA mounds maps tectonically active, gas-permeable strike-slip faults (cf. Figure 2) and that RWA mounds can be used as biological indicators of tectonically active, gas-permeable strike-slip faults [16,17,21]. This is especially useful when information about the active tectonic regime is incomplete or the resolution by technical means is insufficient [16,17]. These faults may be potential earthquake areas and are simultaneously information channels deeply reaching into the crust. Recent research shows that, among others, electromagnetic field changes emitted from the pre-focal area several days before the upcoming earthquake are discussed as one trigger mechanism for bioanomalies prior to large earthquakes [22,23,24,25]. Also mid-term to short-lived thermal infrared anomalies from the rock surface in the 8–12 μm region along thrust or strike-slip faults prior to major earthquakes are reported [26,27,28]. Consequently, an integrated approach in understanding geophysical phenomena requires consideration also of the interaction between the abiotic and biotic environment.
Figure 2. Relationship between helium anomalies (He > 5,220 ppb), fault zones and RWA mounds for the Oberehe site in the West Eifel (approx. 40 km distance from Koblenz, East Eifel); (a) showing helium anomalies (orange dots), RWA mounds (red dots), mineral springs (blue dots) and volcanoes (black dots); (b) showing alignment of the RWA mounds (degrees from N) and (c) fault zone directions [19,20](degrees from N). The comparison of both the spatial RWA mound distribution with the significantly increased gas anomalies confirms that the RWA colony is situated in a tectonically sheared region. Although the reference site is surrounded by various forest stands, no further RWA mounds were found in those forest stands during the area-wide mapping. These results suggest that RWA mounds can be used as biological indicators of active, gas-permeable faults [16,17].
Figure 2. Relationship between helium anomalies (He > 5,220 ppb), fault zones and RWA mounds for the Oberehe site in the West Eifel (approx. 40 km distance from Koblenz, East Eifel); (a) showing helium anomalies (orange dots), RWA mounds (red dots), mineral springs (blue dots) and volcanoes (black dots); (b) showing alignment of the RWA mounds (degrees from N) and (c) fault zone directions [19,20](degrees from N). The comparison of both the spatial RWA mound distribution with the significantly increased gas anomalies confirms that the RWA colony is situated in a tectonically sheared region. Although the reference site is surrounded by various forest stands, no further RWA mounds were found in those forest stands during the area-wide mapping. These results suggest that RWA mounds can be used as biological indicators of active, gas-permeable faults [16,17].
Animals 03 00063 g002
For a detailed and statistically well-based study, it is necessary to investigate whether changes in the ants’ standard daily routine can be correlated with local seismic events. The seismically active study area at the Neuwied Basin (Eifel, Germany) is located in an area with a complex tectonic history. Commencing during the Neogene and persisting during the Quaternary, the uplift of the Rhenohercynian Zone affects the dynamics of the study area. The currently NW-SE striking main stress direction opens pathways for geogenic gases and potential magmas [29,30,31]. Concurrently, a conjugate dextral wrench fault system (striking in 100°–110° WNW-ESE representing the extensional regime of the recent compressional stress field and a conjugated shear system trending in 140°–150° NNW-SSE) in combination with a clockwise block rotation and relatively small kinematic WSW-movements exists [29,32,33,34]. At the same time, Variscan faults as part of a conjugated shear system are reactivated [35,36]. This situation has led to a complex tectonic break clod. Due to the ongoing extension of the Cenozoic rift system in the western part of the Eurasian plate, geotectonic processes still influence the active regional tectonics and lead to a large number of low-magnitude seismic events [37,38]. In the study area, an average of about 100 earthquakes per year with magnitudes up to M 3.9 occur located on different tectonic fault regimes (strike-slip faults and/or normal or thrust faults) [29,39,40].
For three years, from 2009 to 2012, we have monitored two RWA mounds (Formica rufa-group; F. pratensis and F. polyctena, 24/7 by high-resolution cameras. Both mounds are located at two different fault regimes (on a WNW-ESE and on a NNW-SSE trending fault, respectively) approximately 22 km apart. The location of the two cameras is shown in Figure 3. We refer to the cameras as AntCams. During the observation period, ten magnitude > 2.0 earthquakes in a 40 km radius occurred (cf. Table 1).
Table 1. Magnitude > 2.0 earthquakes events (from 2009 to 2012) in a 40 km radius around the location of both cameras (AntCam 1 and 2) within the tectonically active Neuwied basin [39,40].
Table 1. Magnitude > 2.0 earthquakes events (from 2009 to 2012) in a 40 km radius around the location of both cameras (AntCam 1 and 2) within the tectonically active Neuwied basin [39,40].
EpicenterDateTime (UTC)Focal Depth (km)Magnitude (M)Longitude (°)Latitude (°)Linear distance (km) of AntCam 1 and 2 to epicenter
AntCam 1AntCam 2
Bad Ems09/10/200912:36:21.705 *3.27.7250.404022
Kruft04/11/201011:16:16.5010 *2.97.3450.401410
Ochtendung08/05/201109:29:37.0510 *2.37.3850.36158
Plaidt02/21/201202:34:57.62112.57.3850.38176
Kobern03/12/201219:00:42.66122.77.4350.33221
Ochtendung03/29/201201:23:41.0610 *2.27.4150.34202
Plaidt05/22/201200:08:47.9810 *2.37.3750.37176
Lonnig11/15/201216:40:24.6210 *2.47.4150.32214
Lonnig11/22/201215:43:15.9010 *2.87.4150.32214
Winningen12/16/201205:54:35.406.92.27.5350.30308
* Focal depth manually determined by [40].
Figure 3. Location of both cameras (AntCam 1 and 2) within the tectonically active Neuwied basin. Shown are the East Eifel volcanic field (black dots), earthquake events (green dots) 2000–2012 [39], mineral springs and mofettes (blue dots) and a schematic depiction of the tectonically active fault systems (colored lines) [41]. Although within this region the RWA mounds have not been completely mapped so far, it is apparent that RWA mounds indicate tectonically active faults.
Figure 3. Location of both cameras (AntCam 1 and 2) within the tectonically active Neuwied basin. Shown are the East Eifel volcanic field (black dots), earthquake events (green dots) 2000–2012 [39], mineral springs and mofettes (blue dots) and a schematic depiction of the tectonically active fault systems (colored lines) [41]. Although within this region the RWA mounds have not been completely mapped so far, it is apparent that RWA mounds indicate tectonically active faults.
Animals 03 00063 g003

2. Methods

2.1. Technical Equipment of the Monitoring Stations (AntCams)

Both AntCams are high-resolution cameras (Mobotix MX-M12D-Sec-DNight-D135N135; 1,280 × 960 pixels) equipped with both a color (day) and an infrared sensor (night). The video streams of the ants’ activities were recorded with 12 frames per second and annotated with a time stamp (UTC). The network compatible AntCams were connected to a NAS (network attached storage; 1TB) via a POE-supply (power over Ethernet). Because there is a vast amount of data (in summer up to 50 GB per day), the video streams were stored temporarily on the NAS. Every 10 to 14 days, the NAS is exchanged. Since September 2009 (AntCam 1) and January 2010 (AntCam 2), more than 45,000 hours of video streams have been recorded and stored in a video archive.

2.2. Manual Analyses of Video Streams

In the first project phase, the video streams were inspected by three independent operators who were not informed about the occurrence of earthquakes in the regarded period of time. To classify the ant swarm behavior on top of the mounds by a category C(t), eight activity categories were fixed: activity category 0 = no ants on top of the mound, activity category 6 = crowded top of the mound and activity category 7 = sunbathing (cf. Figure 4). The sunbathing in early spring, which takes place only for a few hours at daytime, was classified as a special activity category. The ants are very active on top of the mound within this short period, lasting approximately two weeks starting in March. Within this period of time, a large amount of the entire population is resting on top of the mound, collecting heat with their bodies and transferring the heat into the mound.
Figure 4. AntCam 1: Definition of activity indices with the help of difference images extracted from the master video streams to classify the ant swarm behavior on top of the mounds. Examples: Image 1a (near-infrared image acquired at night-time) shows activity index 0, meaning that there are no moving ants on top of the mound. Image 1b is the corresponding difference image between two subsequent video frames acquired at a time interval of 10 s. Because there are no moving ants on top of the mound, the difference image is black. Image 1a (color image acquired at daytime) shows activity index 1, meaning that there are only a few moving ants on top of the mound. Image 2b is the corresponding difference image. Moving ants appear as small bright dots. It is clearly apparent from the difference images that the number of RWAs on top of the mound is increasing with increasing activity index.
Figure 4. AntCam 1: Definition of activity indices with the help of difference images extracted from the master video streams to classify the ant swarm behavior on top of the mounds. Examples: Image 1a (near-infrared image acquired at night-time) shows activity index 0, meaning that there are no moving ants on top of the mound. Image 1b is the corresponding difference image between two subsequent video frames acquired at a time interval of 10 s. Because there are no moving ants on top of the mound, the difference image is black. Image 1a (color image acquired at daytime) shows activity index 1, meaning that there are only a few moving ants on top of the mound. Image 2b is the corresponding difference image. Moving ants appear as small bright dots. It is clearly apparent from the difference images that the number of RWAs on top of the mound is increasing with increasing activity index.
Animals 03 00063 g004
Figure 5. Results of manually inspected video streams (mean values generated from the three single curves) for the time segment September 08–30, 2009 (AntCam 1), in combination with local earthquakes (magnitude and depth of epicenters, cf. also Table 2), indicates changes in the ants’ behavior hours before the earthquake took place: e.g., the nocturnal rest phase and daily activity were suppressed September 10, 2009 and standard daily routine did not resume until the next day (cf. Figure 5(a) and Figure 8(b)). The spectrum of the tidal potential for the mentioned time segment is given in Figure 5(b)) and the weather data in Figure 5(c–f).
Figure 5. Results of manually inspected video streams (mean values generated from the three single curves) for the time segment September 08–30, 2009 (AntCam 1), in combination with local earthquakes (magnitude and depth of epicenters, cf. also Table 2), indicates changes in the ants’ behavior hours before the earthquake took place: e.g., the nocturnal rest phase and daily activity were suppressed September 10, 2009 and standard daily routine did not resume until the next day (cf. Figure 5(a) and Figure 8(b)). The spectrum of the tidal potential for the mentioned time segment is given in Figure 5(b)) and the weather data in Figure 5(c–f).
Animals 03 00063 g005
Table 2. Earthquake events in September 2009 in a 40 km radius of the location of AntCam 1 within the tectonically active Neuwied basin [39].
Table 2. Earthquake events in September 2009 in a 40 km radius of the location of AntCam 1 within the tectonically active Neuwied basin [39].
EpicenterDateTime (UTC)Focal Depth (km)Magnitude (M)Longitude (°)Latitude (°)Linear distance (km) of AntCam 1 to epicenter
Bad Ems09/10/200912:36:21.705.0 *3.27.7250.4040
Polch09/12/200904:26:57.83.20.97.3350.3117
Nickenich09/12/200922:57:19.22.80.47.3250.4213
Ochtendung09/15/200923:06:23.03.00.67.4150.3620
Plaidt09/29/200904:38:26.42.10.57.3850.3817
* Focal depth manually determined by [39].
Every 10 minutes of each day (24 hours), a one minute sequence was evaluated at 8‐fold speed by all three operators. At this frame rate, the motion behavior of the ants is well discernible for the human eye. They scored 144 discrete numerical values as activity indices per day according to the annotated UTC time stamp. To ensure the same basis of RWA activity for the eight activity categories, the three operators were provided with master video streams showing clear examples of the eight activity categories. The results were processed in Excel, and daily activity curves (rounded mean values generated from the three single curves) connecting the discrete numerical values (activity indices) were generated. For example, if no ants were present on top of the mound, the 3 observers scored 0, showing a green curve touching the x-axis (cf. Figure 5(a) and Figure 8(a)). If the rounded mean value was, e.g., 4, the curve raised to the left-handed activity index 4 on the y-axis (cf. Figure 5(a) and Figure 8(a)). An example of the manual analysis is presented in Figure 5(a)). The Ants’ activity shows an abnormal behavior, specifically before an earthquake with a magnitude of 3.2 (Bad Ems earthquake in 40 km distance from AntCam 1 [40]). Although this magnitude can be addressed as low compared to large earthquake events, e.g., in Italy, the tectonic regime and the nature of a fault has to be taken into account. This earthquake may have triggered the tectonically active and gas-permeable “Laacher See Strike-slip Fault”, which can be tracked from the AntCam 1 position in the West, Wehrer Kessel, Laacher See, Plaidt to Bad Ems in the East [41]. In this case, bioanomalies prior to M > 2 earthquakes should be monitored in contrast to normal or reverse faults that do not show any gas-permeabilities, due to its nature. The results will be discussed in detail in Section 3.
Additionally, all abiotic (duration of sunshine, rain, snow, wind, aerating of the mound etc.) and biotic parameters, as the appearance and duration of predators (woodpecker, magpie, blackbird), and/or disturbing factors (boars, badger, “anting” of birds etc.) on top of the mound and various activities of ants (repeated nuptial flights per month, transport of larval and pupa, mound construction activities etc.) that might have an effect on the ants’ daily routine were read on a per-second basis and correlated with the analyzed daily activity to document the possible influence of these factors on the ants’ activity (cf. Figure 5(a), e.g., nuptial flights). At intervals, nest air measurements (CO2, helium, radon, H2S and CH4) were performed.
Earth tide data provided by the University of Bern [42] were evaluated (cf. Figure 5(b)). The earth’s response to the gravitational influence of both the Sun and the Moon results in two slight lunar and two solar tidal bulges (earth tides). The two bulges occur at the surface part of the earth that approximately faces the Moon and at the opposite surface part while the Earth rotates around its axis. Earth tides can be correlated with volcanic eruptions [43]. Whether earth tides also trigger earthquakes is still under discussion [44]. Another question is if earth tides might also trigger the ants’ behavior with regard to their daily routine or nuptial flights (cf. Figure 12(a)).
Additionally, diurnal, monthly and mean values of a number of meteorological elements (minimum, maximum and mean temperature [°C], sunshine duration [hours], wind force [Bft] and velocity [m/s], amount of precipitation [mm] and relative humidity [%]) of the meteorological station Nuerburg-Barweiler provided by the German Meteorological Service (DWD) were considered [45].
Meanwhile, more than 15,000 of hours of video data of both AntCams were analyzed manually.

2.3. Automated Image Analysis

The key question was whether the ants’ behavioral changes and their correlation with earthquake events are statistically significant and if detection by an automated system is possible. Based on the manually gained results, an automated image evaluation routine was developed with the software Matlab and applied to the video streams to gain an objective interpretation of the ant behavior. In our image analysis system, an activity index is computed automatically in a first step based on a difference image technique. Based on a regression approach, this activity index is mapped to the range of the manually estimated activity indices to allow direct comparisons.
The basis of the routine is a comparison of the ants’ activity on top of the mound every 10 seconds, i.e., a video image is acquired every 10 seconds. In the style of the concept proposed in [46] (cf. [47] for further discussion of their method), for each pixel Iuv(t) with coordinates (u,v) of the video image acquired at time t, the system determines the pixel-wise absolute difference Duv(t) with respect to the previous image Iuv(t ‒ ∆t) acquired at the previous time step according to:
Duv(t) = |Iuv(t) ‒ Iuv(t‒∆t)|
In the absolute difference image Duv(t), moving ants appear as bright spots (cf. Figure 6). The well-known disadvantage of the concept of subtracting subsequent images of a sequence, namely that objects tend to be only partially visible in the difference image [46,47], is not relevant in our scenario, due to the small size of the ants in the image of at most only a few pixels (to overcome this problem for large objects, a statistical method to generate a specific “background image” is proposed in [47]). To reduce negative influences caused by, e.g., moving blades of grass, a mask was used covering only the visible top of the mound. To compensate slight movements of the camera, e.g., due to wind, before the determination of the absolute difference image, an image registration of the current image on the previous image is performed based on normalized cross-correlation of well-distinguishable image points [48].
Figure 6. Near-infrared video image of top of the RWA mound monitored by AntCam 1 with active near-infrared illumination (a). Grass containing chlorophyll appears bright, as it strongly reflects near-infrared light. The ants are visible as very small dark spots. Change detection between two video frames acquired at nighttime at a time interval of 10 s (b). Moving ants are apparent as bright spots.
Figure 6. Near-infrared video image of top of the RWA mound monitored by AntCam 1 with active near-infrared illumination (a). Grass containing chlorophyll appears bright, as it strongly reflects near-infrared light. The ants are visible as very small dark spots. Change detection between two video frames acquired at nighttime at a time interval of 10 s (b). Moving ants are apparent as bright spots.
Animals 03 00063 g006
In our system, the pixel intensities of Duv(t), which exceed three times the pixel noise σ, are summed over the image region, M, corresponding to the top of the mound, where the ant movements are visible:
Animals 03 00063 i001
The sum value A(t) is regarded as a measure for the ants’ average activity, but it cannot be compared directly to the category C(t) defined in Section 2.2 based on manual inspection of the video sequences. To derive an activity index from the image analysis system, which is directly comparable to these categories, it has to be taken into account that images are acquired at daytime with a CMOS RGB color camera, while for images acquired at night time, a CMOS near-infrared camera sensor with active illumination is used, leading to a systematically different appearance of the images acquired at daytime and night time, respectively. Hence, the activities A(t) are subdivided into daytime values, Ad(t), and night time values, An(t). Over time intervals T of about 2–3 weeks duration, two sets of values:
Animals 03 00063 i002
are defined. To the daytime set, Sd, a first-order polynomial regression function, pd(x), and to the nighttime set, Sn, a second-order polynomial regression function, pn(x), is fitted in such a way that the root-mean-squared difference between C(t) and pd(Ad(t)) and pn(An(t)), respectively, is minimized. Hence, the function values, pd(Ad(t)) and pn(An(t)), of the polynomials correspond to automatically determined activity indices, which can be directly compared to the manually inferred categories, C(t). To suppress the influence of image noise, a temporal median filter with a window length of 10 minutes is applied to the activity indices, pd(Ad(t)) and pn(An(t)).
The polynomials and the manual analysis are in good agreement (Figure 7(a–c)) and show the same abnormal behavior as the manual analysis. The result will be discussed in detail in Section 3.
Figure 7 indicates changes in the ants’ behavior hours before the earthquake took place: the nocturnal rest phase and daily activity were suppressed on September 10, 2009, and standard daily routine did not resume until the next day (Figure 5(a–c) and Figure 8(b)). A comparison of the earth tides for the specific time interval as provided by the University of Bern (Figure 7(d–f)) and the earthquake events (Figure 7(a–c)) for the given time period suggest no possible correlation. For the given time period, ventilation of the mound (Figure 7(a–c)) takes place only at positive vertical displacements (Figure 7(d–f)). Weather conditions in September 2009 (Figure 7(g–j)) do not have a great influence on the RWA behavior. In the night from 19 to 20 September, a mismatch between the manual and the automatic analyses is apparent (Figure 7(b)). A possible explanation is the change in brightness of the material on top the mound at night. It is necessary to further analyze this issue in long-term studies.
Figure 7. Detailed presentation of the activity index (blue: automatically determined; green: manually determined (ac)) for the time interval September 08–30, 2009 (AntCam1) in correlation with local earthquakes (magnitude and depth of epicenters), earth tide displacement (df), and meteorological data (gj).
Figure 7. Detailed presentation of the activity index (blue: automatically determined; green: manually determined (ac)) for the time interval September 08–30, 2009 (AntCam1) in correlation with local earthquakes (magnitude and depth of epicenters), earth tide displacement (df), and meteorological data (gj).
Animals 03 00063 g007

3. Results and Discussion

Early results of the manually analyzed long-term observations show that ants have a well-identifiable standard daily routine that can be best characterized by an M-shaped curve: activity rises from dawn until approximately 11 A.M. Activity peaks can be observed at midday and in the late afternoon (cf. Figure 8(a)). During the evening, the activity is reduced again. This approach is contradictory to the observation documented so far [49]. It shows that it is possible to obtain an objective and testable null hypothesis and detailed data by a long-term, 24/7 in situ monitoring of the ants’ behavior using suitable technical equipment and experimental setup. A previous study [49] merely considers a serendipitous event. Additionally, no further information is provided on whether the acquired very short (30 s) single video streams have been objectively analyzed by an automated routine.
Figure 8. RWA have a well-identifiable standard daily routine that can be best characterized by an M-shaped curve (a). The correlation between the video observations and the times of occurrence of local earthquakes indicates changes in the ants’ behavior hours before the earthquake took place (b); the nocturnal rest phase and daily activity were suppressed and standard daily routine did not resume until the next day.
Figure 8. RWA have a well-identifiable standard daily routine that can be best characterized by an M-shaped curve (a). The correlation between the video observations and the times of occurrence of local earthquakes indicates changes in the ants’ behavior hours before the earthquake took place (b); the nocturnal rest phase and daily activity were suppressed and standard daily routine did not resume until the next day.
Animals 03 00063 g008
The video observations and the times of occurrence of local earthquakes indicate changes in the ants’ behavior hours before the earthquake took place: the nocturnal rest phase and daily activity were suppressed, and standard daily routine did not resume until the next day (cf. Figure 8(b)). The results show that it is essential to monitor the ants’ behavior before the earthquake strikes a region. The eight minutes observation period in the run-up of the Landers M 7.4 earthquake in the Mojave Desert is too short and cannot produce statistically significant results, as it cannot be determined for such a short period of time whether the observed ants’ behavior is normal or abnormal. To define the “normal” state, it would have been necessary in that scenario to observe the ants’ behavior over longer periods of time during earthquake-free intervals in order to obtain a statistically relevant reference for observations during earthquakes. Furthermore, our observations suggest that a documentation of the aftershock activity alone is insufficient and that the conclusion that ants do not show any reaction on earthquakes, as it was made in the literature so far, is not justified based on the limited amount of data examined in previous studies [49]. Another prerequisite is to have detailed information on the active tectonic regime in the corresponding region and an understanding of geotectonic processes for interpretation of geo-biological processes, i.e., RWA as bioindicators for the identification of active tectonic systems. Observations some 100 km away from the epicenter should be interpreted very carefully, since tectonic fault types and systems can change strongly within this distance. Additionally, it has not been examined so far, e.g., by geochemical analyses, if Messor pergandei has a close relation to tectonic fault systems.
One of the major challenges of seismology is to determine whether fault zones generate any indication of impending earthquakes. The evolutionary mechanism of exaptation, i.e., a vibration triggered early warning response or animals’ perception of these stimuli to P-waves, is reported in [5]. Precursors, as electromagnetic phenomena or the change of the gas content and concentration, might be sensed by the ants due to their electromagnetic sensitivity and ability to detect CO2 [16,17,18]. The objective of this methodological approach was to examine if ants show reactions prior to earthquakes. The specific mechanisms for how and why they show such reactions need to be researched in in situ experiments and in close cooperation with biologists.
In the night from September 19 to 20, a mismatch between the manual and the automatic analyses is apparent (Figure 7(b)). Analyses of the corresponding video material showed that the algorithm works correctly and clearly detects and counts the ants’ movements. A possible explanation is the change in brightness of the material on top of the mound at night that does not play a role for the manual analyses, but for the automated ones. It is necessary to further analyze this issue in long-term studies.
Although we have monitored two different species (F. pratensis and F. polyctena), their daily activity patterns do not show systematic differences at days without any earthquake event and can be considered the same. According to the correlation coefficients (R between 0.82 and 0.97), the statistical correlation can be described as strong. This highly similar behavior of the two different species, as presented in Figure 9, is contradictory to previous studies and the common tenet [50,51].
External influence factors, like predators (woodpeckers, blackbirds, etc.) or occasionally “anting” of birds to banish parasites with formic acid, which only last for some seconds to some minutes at maximum, create only a slight increase of the peak activity, but do not have any influence on the ants’ daily routine. In practice, such influences can be avoided by fencing the mound.
Figure 9. Comparison of the ants’ activity of F. pratensis (AntCam 1, blue curve) and F. polyctena (AntCam 2, red curve) at given days in April 2010 and their strong statistical correlation (R between 0.82 and 0.97).
Figure 9. Comparison of the ants’ activity of F. pratensis (AntCam 1, blue curve) and F. polyctena (AntCam 2, red curve) at given days in April 2010 and their strong statistical correlation (R between 0.82 and 0.97).
Animals 03 00063 g009
Figure 10. On April 20, 2010, a badger appeared on top of the mound (21:50:56 to 21:51:30 (UTC)). As can be seen, there is only a small peak activity caused by the badger, because it hits the main exit at this time interval. After this short event, the ants’ activity is calmed down to the normal activity level.
Figure 10. On April 20, 2010, a badger appeared on top of the mound (21:50:56 to 21:51:30 (UTC)). As can be seen, there is only a small peak activity caused by the badger, because it hits the main exit at this time interval. After this short event, the ants’ activity is calmed down to the normal activity level.
Animals 03 00063 g010
The external influence of larger animals, like badgers or boars, does not show raised activity peaks, unless a part of the mound is affected by digging (badger) or destruction (woodpecker, black bird etc.). Then, the ants’ activity rises to protect the mound and the population. This can be documented by a peak activity (e.g., reconstruction of the mound), which amounts to only a small fraction of the daily routine and is negligible, for the daily routines have to be considered across the whole day. Nevertheless, this extraordinary situation is documented (cf. Figure 10).
During the winter season, the ants’ activity is—as it is well known—very small to non-existent. However, we have observed ants to be visible on top of the mound when an earthquake with a magnitude of about 2 occurs, despite snow cover and freezing temperatures and, moreover, with delayed activity compared to the spring-autumn seasons. The question is whether the ants’ activity will rise to a higher level for earthquakes with higher magnitudes.
The automatically analyzed activity index (“purple curve” in Figure 11) shows that our image analysis system works independently of given manual results. Furthermore, a comparison of the manually inspected video observations (green curve) and the “raw” activity indices, Ad(t) and An(t) (prior to the mapping to the manually inferred categories C(t) based on the polynomial regression functions) (blue curve), shows that both activity curves show the same trends in the daily ants’ routine.
Figure 11. Automatically analyzed raw activity index A(t) (before polynomial regression) in comparison with the manually analyzed data. A qualitatively similar behavior is clearly apparent.
Figure 11. Automatically analyzed raw activity index A(t) (before polynomial regression) in comparison with the manually analyzed data. A qualitatively similar behavior is clearly apparent.
Animals 03 00063 g011
According to our first analyses, clearly anomalous ants behavior may occur for the local earthquake with a magnitude > 2 (Figure 8(b)). At lower magnitudes (cf. Table 3), an anomalous behavior may possibly be present as well (Figure 12). The question is: how RWA would react in regions with higher magnitude (M ≥ 4) earthquakes and if there is a habituation effect with regard to lower magnitudes. Possibly, the ants show reactions only for higher magnitudes. A comparison of the earth tides and the earthquake events suggest a possible correlation for the given time interval. An accumulation of earthquakes (June 09 (day 161); June 12–14 (days 163–165) and June 26–29 (days 177–180)) occurs only when vertical displacements are positive and at their maximum. Also, nuptial flights occur only when vertical displacements are positive at their maximum in the given time interval (June 22–30 (days 174–181)). It is necessary to further analyze these relationships in long-term studies.
Possibly, the ants’ behavior may also depend on the specific fault type. Strike-slip faults create deep-reaching tectonic fractures that might be gas-permeable and show mineralizations. Earthquakes occurring on these faults would trigger gas anomalies prior to the event that can propagate over a certain distance from the future quake epicenter, whereas along normal or reverse faults, due to their physical characteristics, such anomalies would be smaller.
Figure 12. The automatically determined activity index (blue) for the time interval June 08–30, 2010 (AntCam 1) combined with the manually analyzed activities (green) in correlation with local earthquakes with a magnitude of <2. At lower magnitudes, an anomalous behavior may possibly be present as well. A comparison of the earth tides and the earthquake events suggest a possible correlation for the given time interval. An accumulation of earthquakes (June 09 (day 161); June 12–14 (days 163–165) and June 26–29 (days 177–180)) occurs only when the absolute values of the vertical displacements are at their maximum. A similar, but less pronounced, relation can be observed for nuptial flights in the given time interval (June 22–30, days 174–181). It is necessary to further analyze these relationships in long-term studies.
Figure 12. The automatically determined activity index (blue) for the time interval June 08–30, 2010 (AntCam 1) combined with the manually analyzed activities (green) in correlation with local earthquakes with a magnitude of <2. At lower magnitudes, an anomalous behavior may possibly be present as well. A comparison of the earth tides and the earthquake events suggest a possible correlation for the given time interval. An accumulation of earthquakes (June 09 (day 161); June 12–14 (days 163–165) and June 26–29 (days 177–180)) occurs only when the absolute values of the vertical displacements are at their maximum. A similar, but less pronounced, relation can be observed for nuptial flights in the given time interval (June 22–30, days 174–181). It is necessary to further analyze these relationships in long-term studies.
Animals 03 00063 g012
Table 3. Earthquake events in June 2010 in a 40 km radius around the location of AntCam1 within the tectonically active Neuwied basin [39].
Table 3. Earthquake events in June 2010 in a 40 km radius around the location of AntCam1 within the tectonically active Neuwied basin [39].
EpicenterDateTime (UTC)Focal Depth (km)Magnitude (M)Longitude (°)Latitude (°)Linear distance (km) of AntCam 1 to epicenter
Polch06/09/201003:31:25.03.40.17.3350.3018
Koblenz06/13/201014:46:34.75.0 *0.87.5650.3430
Nickenich06/13/201020:40:15.215.50.77.3250.3918
Plaidt06/14/201012:32:30.96.40.97.3450.4014
Plaidt06/14/201020:41:47.45.30.87.3550.3915
Kruft06/26/201005:23:34.729.80.97.3650.3915
S. Nickenich06/26/201005:24:18.615.20.67.3250.3913
Winningen06/26/201008:15:56.32.10.67.5150.3228
N. Ochtendung06/27/201014:27:27.12.10.67.3950.3718
E. Kruft06/27/201015:04:37.72.10.67.3550.3815
Boppard06/29/201005:59:34.02.10.27.6250.2339
* Focal depth could not be determined exactly by [39].

4. Conclusions

RWA bioanomalies (suppression of the nocturnal rest phase and daily activity, continuation of the standard daily routine not before the next day) have been recorded and analyzed for several earthquakes with magnitudes of up to about 3. The automated image analyses routine has been shown to provide a valuable tool to objectively identify and classify the ants’ activity on top of the mounds and an examination of its correlation with earthquakes. The results show that the automated system works and can be utilized in practice.
According to our first analyses, clearly anomalous ant behavior occurs for the local earthquake on September 10, 2009, with a magnitude of 3.2. At lower magnitudes, an anomalous behavior may possibly be present as well. The specific mechanisms for how and why they show such reactions need to be researched in in situ experiments and in close cooperation with biologists. An open question is whether RWA would react in the same way on earthquakes with higher magnitudes (M > 4), which have not occurred during our period of observation. This will be one topic of future research. In the next step, the analyses routine will be integrated into the AntCams to transmit the activity indices online and in real-time. After a test phase of about half of a year, we plan to install our monitoring system in tectonically more active regions.
Although the investigation and results presented here are promising, they are only a first step towards a completely new research complex. Long-term studies have to show whether confounding factors and climatic influences can clearly be distinguished. However, our early results suggest that it makes sense to consolidate and extend the research to determine a pattern for exceptional activity situations. These studies are not provided for earthquake prediction, but are an important step towards the understanding of geobiological processes.

Acknowledgments

We would like to thank Elisabeth and Rainer Degen (Weibern) and Bernhard Mauel, First District Councillor, Landkreis Mayen-Koblenz, as well as Hans-Toni Dickers (Treasurer Goloring e.V.) for their hospitality and for supporting our research project throughout the years by granting permission for installing our AntCams on their properties.

Conflict of Interest

The authors declare no conflict of interest.

References and Notes

  1. Fletcher, L.E.; Christensen, I.A.; Liberda, J.J.; Rojas, J.I.; Borrero del Pino, C. Methods of Multivariable Earthquake Precursor Analysis and a Proposed Prototype Earthquake Early Warning System. In Proceeding of the International Space University Summer Session Program, Beijing, China, 25 June–24 August 2007; pp. 1–11.
  2. van Dam, J.W.; Horton, W.; Tsintsadze, N.L.; Kaladze, T.D.; Garner, T.W.; Tsamalashvili, L.V. Some physical mechanisms of precursors to earthquakes. J. Plasma Fusion Res. 2009, 8, 199–202. [Google Scholar]
  3. Yasuokaa, Y.; Igarashib, G.; Ishikawac, T.; Tokonamic, S.; Shinogia, M. Evidence of precursor phenomena in the Kobe earthquake obtained from atmospheric radon concentration. Appl. Geochem. 2006, 21, 1064–1072. [Google Scholar] [CrossRef]
  4. Schaal, R.B. An evaluation of the animal-behavior theory for earthquake prediction. Calif. Geol. 1988, 41, 41–45. [Google Scholar]
  5. Kirschvink, J.L. Earthquake prediction by animals: Evolution and sensory perception. Bull. Seism. Soc. Am. 2000, 90, 312–323. [Google Scholar] [CrossRef]
  6. Ikeya, M.; Yamanaka, C.; Mattsuda, T.; Sasaoka, H.; Ochiai, H.; Huang, Q.H.; Ohtani, N.; Komuranani, T.; Ohta, M.; Ohno, Y.; Nakagawa, T. Electromagnetic pulses generated by compression of granitic rocks and animal behaviour. Episodes 2000, 23, 262–265. [Google Scholar]
  7. Bhargava, N.; Katiyar, V.K.; Sharma, M.L.; Pradhan, P. Earthquake prediction through animal behavior: A review. Indian J. Biomech. 2009, 7-8, 159–165. [Google Scholar]
  8. Schreiber, U.; Brennholt, N.; Simon, J. Gas permeable deep reaching fracture zones encourage site selection of ants. Ecol. Ind. 2009, 9, 508–517. [Google Scholar] [CrossRef]
  9. Berberich, G.; Schreiber, U. First Results of 1 Year Monitoring of Red Wood Ant Behaviour as Short-Term (>1 h) Indicators for Earthquake Prediction. In Proceedings of the EGU General Assembly, Vienna, Austria, 4–8 April 2011.
  10. Erpenbeck, A.; Kirchner, W. Zur Kälteresistenz der Kleinen Roten Waldameise Formica polyctena Foerst. (Hymenoptera, Formicidae). Z Angew. Entomol. 1983, 96, 271–281. [Google Scholar] [CrossRef]
  11. Hetz, S.K.; Bradley, T.J. Insects breathe discontinuously to avoid oxygen toxicity. Nature 2005, 433, 516–519. [Google Scholar]
  12. Hölldobler, B.; Wilson, E.O. Der Superorganismus—Der Erfolg von Ameisen, Bienen, Wespen und Termiten; Springer: Berlin, Germany, 2010. [Google Scholar]
  13. Camlitepe, Y.; Stradling, D.J. Wood ants orient to magnetic fields. Proc. Roy. Soc. Lond. 1995, 261, 37–41. [Google Scholar] [CrossRef]
  14. Camlitepe, Y.; Aksoy, V.; Uren, N.; Yilmaz, A.; Becenen, I. An experimental analysis on the magnetic field sensitivity of the black-meadow ant formica pratensis retzius (Hymenoptera: Formicidae). Acta Biol. Hungarica 2005, 56, 215–224. [Google Scholar] [CrossRef]
  15. de Oliveira, J.F.; Wajnberg, E.; de Souza Esquivel, D.M.; Weinkauf, S.; Winklhofer, M.; Hanzlik, M. Ant antennae: Are they sites for magnetoreception? J. Roy. Soc. Interface 2010, 7, 143–152. [Google Scholar] [CrossRef]
  16. Berberich, G. Identifikation junger gasführender Störungszonen in der West- und Hocheifel mit Hilfe von Bioindikatoren. Ph.D. Thesis, 2010. [Google Scholar]
  17. Berberich, G.; Klimetzek, D.; Wöhler, C.; Grumpe, A. Statistical Correlation between Red Wood Ant Sites and Neotectonic Strike-Slip Faults. In Proceedings of the EGU General Assembly, Vienna, Austria, 22–27 April 2012.
  18. Jähne, B. Digitale Bildverarbeitung, 6th ed.; Springer: Berlin, Germany, 2005. [Google Scholar]
  19. Geologische Übersichtskarte von Deutschland 1:200 000; Digitale Ausgabe, Blatt Trier CC 6302; Bundesanstalt für Geowissenschaften und Rohstoffe (BGR): Hanover, Germany, 2010.
  20. Erläuterungen zu Blatt C 5502 Aachen; Geologisches Landesamt Nordrhein-Westfalen (GLA-NRW): Krefeld, Germany, 1992; p. 84.
  21. Schreiber, U.; Berberich, G. Red Wood Ant Mounds as Biological Indicators for Earthquake-bearing Fault Systems. In Proceedings of the EGU General Assembly, Vienna, Austria, 4–8 April 2011.
  22. Grant, R.A.; Halliday, T.; Bladerer, W.P.; Leuenberger, F.; Newcomer, M.; Cyr, G.; Freund, F.T. Ground water chemistry changes before major earthquakes and possible effects on animals. Int. J. Environ. Res. Public Health 2011, 8, 1936–1956. [Google Scholar]
  23. Dologlou, E. Recent aspects on possible interrelation between precursory electric signals and anomalous bioeffects. Nat. Hazards Earth Syst. Sci. 2010, 10, 1951–1955. [Google Scholar] [CrossRef]
  24. Weaver, J.C.; Vaughan, T.E.; Astumian, R.D. Biological sensing of small field differences by magnetically sensitive chemical reactions. Nature 2000, 405, 707–709. [Google Scholar]
  25. Weaver, J.C. Understanding conditions for which biological effects of nonionizing electromagnetic fields can be expected. Bioelectrochemistry 2002, 56, 207–209. [Google Scholar] [CrossRef]
  26. Ouzounov, D.; Freund, F.T. Mid-infrared emission prior to strong earthquakes analyzed by remote sensing data. Adv. Pace Res. 2004, 33, 268–273. [Google Scholar]
  27. Freund, F.T. Toward a unified solid state theory for pre-earthquake signals. Acta Geophys. 2010, 58, 719–766. [Google Scholar] [CrossRef]
  28. Carreno, E.; Capote, R.; Yague, A. Observations of Thermal Anomaly Associated to Seismic Activity from Remote Sensing. In Proceedings of General Assembly of European Seismology Commission, Lisbon, Portugal, 10–15 September 2001; pp. 265–269.
  29. Hinzen, K.G. Stress field in the Northern Rhine area, Central Europe, from earthquake fault plane solutions. Tectonophysics 2003, 377, 325–356. [Google Scholar] [CrossRef]
  30. Illies, J.; Baumann, H. Crustal dynamics and morphodynamics of the Western European Rift System. Z. f. Geomorph. 1982, 42 (Suppl.), 135–165. [Google Scholar]
  31. Rummel, F.; Baumgärtner, J. Hydraulic Fracturing in-situ stress and permeability measurements in the research borehole Konzen, Hohes Venn (West Germany). N. Jb. Geol. Paläont. Abh. 1985, 171, 183–195. [Google Scholar]
  32. Schreiber, U.; Rotsch, S. Cenozoic block rotation according to a conjugate shear system in central Europe—Indications form paleomagnetic measurements. Tectonophysics 1998, 299, 111–142. [Google Scholar] [CrossRef]
  33. Walker, K.T.; Bokelmann, G.H.R.; Klemperer, S.L.; Bock, G. Shear-wave splitting around the Eifel hotspot: Evidence for a mantle upwelling. Geophys. J. Int. 2005, 163, 962–980. [Google Scholar] [CrossRef]
  34. Tesauro, M.; Hollenstein, C.; Egli, R.; Geiger, A.; Kahle, H.-G. Analysis of central western Europe deformation using GPS and seismic data. J. Geodynam. 2006, 42, 194–209. [Google Scholar] [CrossRef]
  35. Weisser, D. Tektonik und Barytgänge in der SE-Eifel. Zeitschrift der deutschen geologischen Gesellschaft 1963, 115, 33–68. [Google Scholar]
  36. Redecke, P. Zur Geochemie und Genese variszischer und postvariszischer Buntmetallmineralisation in der Nordeifel und der Niederrheinischen Bucht. Ph.D. Thesis, RWTH Aachen, Aachen, Germany, 1992. [Google Scholar]
  37. Campbell, J.; Kümpel, H.‐J.; Fabian, M.; Fischer, D.; Görres, B.; Keysers, C.J.; Lehmann, K. Recent movement pattern of the Lower Rhine Embayement from tilt, gravity and GPS data. Netherlands J. Geosci. 2002, 81, 223–230. [Google Scholar]
  38. Ahorner, L. Present-day stress field and seismotectonic block movements along major faults in Western Europe. Tectonophysics 1975, 29, 233–249. [Google Scholar] [CrossRef]
  39. Earthquake Data Catalogue. Department of Earthquake Geology of Cologne University. Available online: http://www.seismo.uni-koeln.de/catalog/index.htm (accessed on 9 November 2012).
  40. Earthquake Catalogue. Landesamt für Geologie und Bergbau Rheinland-Pfalz. Available online: http://www.lgb-rlp.de/erdbeben.htm (accessed on 10 September 2009).
  41. Schreiber, U.; Berberich, G. Why does the Size of the Laacher See Magma Chamber and its Caldera Size not go together?—New Findings with regard to Active Tectonics in the East Eifel Volcanic Field. In Proceedings of the EGU General Assembly, Vienna, Austria, 7–12 April 2013.
  42. Earth Tides Vertical Displacements. Astronomisches Institut, University of Bern. Available online: http://www.aiub.unibe.ch/content/services/earth_tides/index_eng.html (accessed on 1 September 2009).
  43. Earth Tides and Volcano Monitoring. Hawaian Volcano Observatory. USGS. Available online: http://hvo.wr.usgs.gov/volcanowatch/98_05_28.html (accessed on 28 November 2012).
  44. Agnew, D.C. Earth Tides. In Treatise on Geophysics; Schubert, G., Herring, T., Eds.; Elsevier: Oxford, UK, 2007; Volume 3, pp. 163–195. [Google Scholar]
  45. Wetter und Klima. Deutscher Wetterdienst. Freie Met. Infos. Available online: http://www.dwd.de/...202758871200642573928&_urlType=action&_pageLabel=_dwdwww_menu2_leistungen_a-z_freiemetinfos (accessed on 1 September 2009).
  46. Hsu, Y.Z.; Nagel, H.H.; Rekers, G. New likelihood test methods for change detection in image sequences. Comput. Vision Graph. Image Process. 1984, 26, 73–106. [Google Scholar] [CrossRef]
  47. Rosin, P.; Ellis, T. Image Difference Threshold Strategies and Shadow Detection. In Proceedings of the Sixth British Machine Vision Conference, Birmingham, UK, 11–14 September 1995; pp. 347–356.
  48. Brown, L.G. A Survey of image registration techniques. ACM Comput. Surv. 1992, 24, 325–376. [Google Scholar] [CrossRef]
  49. Lighton, J.R.B.; Duncan, F.D. Shaken, not stirred: A serendipitous study of ants and earthquakes. J. Exp. Biol. 2005, 208, 3103–3107. [Google Scholar] [CrossRef]
  50. Gößwald, K. Artunterschiede der Waldameisen in Aussehen, Lebensweise, Organisation, Verhalten, Nest- und Straßenbau, Ökologie und Verbreitung. Waldhygiene 1982, 14, 161–192. [Google Scholar]
  51. Goropashnaya, A. Phylogeographic Structure and Genetic Variation in Formica Ants. Ph.D. Thesis, Uppsala University, Uppsala, Sweden, 2003. [Google Scholar]

Share and Cite

MDPI and ACS Style

Berberich, G.; Berberich, M.; Grumpe, A.; Wöhler, C.; Schreiber, U. Early Results of Three-Year Monitoring of Red Wood Ants’ Behavioral Changes and Their Possible Correlation with Earthquake Events. Animals 2013, 3, 63-84. https://doi.org/10.3390/ani3010063

AMA Style

Berberich G, Berberich M, Grumpe A, Wöhler C, Schreiber U. Early Results of Three-Year Monitoring of Red Wood Ants’ Behavioral Changes and Their Possible Correlation with Earthquake Events. Animals. 2013; 3(1):63-84. https://doi.org/10.3390/ani3010063

Chicago/Turabian Style

Berberich, Gabriele, Martin Berberich, Arne Grumpe, Christian Wöhler, and Ulrich Schreiber. 2013. "Early Results of Three-Year Monitoring of Red Wood Ants’ Behavioral Changes and Their Possible Correlation with Earthquake Events" Animals 3, no. 1: 63-84. https://doi.org/10.3390/ani3010063

APA Style

Berberich, G., Berberich, M., Grumpe, A., Wöhler, C., & Schreiber, U. (2013). Early Results of Three-Year Monitoring of Red Wood Ants’ Behavioral Changes and Their Possible Correlation with Earthquake Events. Animals, 3(1), 63-84. https://doi.org/10.3390/ani3010063

Article Metrics

Back to TopTop