The GELATICA project (GEorgian Large-area Angle and TIme Coincidence Array) consists of several cosmic ray detector arrays located in two cities in Georgia. Designed to detect Extensive Air Showers (EAS) with high-precision timing, these observatories aim to reconstruct the arrival directions of cosmic ray showers. GELATICA is the part of the international CREDO collaboration (Cosmic Ray Extremely Distributed Observatory), with a primary goal of searching for Cosmic Ray Ensembles (CRE)—correlated cosmic ray events potentially originating from a common astrophysical or exotic source. This paper presents an overview of GELATICA's instrumentation and methodology, along with an analysis of multi-year data. We evaluate detector performance, characterize the distribution of recorded signals, and assess the potential for detecting CREs via coincident signals across spatially separated arrays. We report the detection of two pairs of events coinciding in time with 0.5 millisecond accuracy. However, the arrival directions of these coincident events do not suggest a common origin. Nevertheless, further data collection and an increased number of observatories could enhance the likelihood of identifying true CRE candidates.

We compared two methods for determining the detector performance and the cosmic ray small-surface array aperture. The comparison was performed using the GELATICA network station at Telavi Iakob Gogebashvili State University (hereafter TEL) as an example. The first method is a standard analytical method. It is based on mean values of variables and averaged distributions. This analytical approach to data analysis was the focus of the research carried out within the GELATICA project. The project is a member of the international CREDO (Cosmic-Ray Extremely Distributed Observatory) Collaboration, whose main goal is the detection and global analysis of cosmic ray ensembles. In contrast to the traditional approach, which focuses on the detection of individual cosmic ray extensive air showers, CREDO aims to connect existing and yet-to-be built cosmic ray arrays into a worldwide network, thus creating a scientific tool for the global analysis of ultra-high energy cosmic rays. The present work was the first to compare the determination of the parameters of the extensive air showers recorded by the TEL array with simulation results using the CORSIKA code. The interpretation of the results of the analytical method for the evaluation of detector registration thresholds is generally found to be inconclusive. In order to obtain definitive results, we propose additional measurements and a new method of array detector performance. We show that the energy spectrum obtained analytically is nearly in agreement with that obtained from simulations. Differences are apparent for the primary particle energy threshold. The difference in the overall counting rate for the TEL array is of the order of 4%.