2.1.1. General Consideration

A GSM (Global System for Mobile communications) relay site is an infrastructure comprising a base station (BTS) [20]-[23] with its antennas, a tower, a technical shelter for equipment, an air conditioning system, power systems (batteries, generator) and transmission links to the core network [24][25], managing the coverage of a radio cell for mobile communication, with an architecture that also includes BSCs (base station controllers) and the MSC (switching center) [26][27].

2.1.2. Physical and Technical Elements, Operation, and Network Architecture

According to information from the literature [28]-[30], the physical and technical elements of a GSM relay site are as follows:

Pylon/Mast: Metal structure (guyed or self-supporting) supporting the antennas at an optimal height for radio coverage.Antennas: Transmit and receive radio signals (frequencies such as 900 MHz, 1800 MHz) to and from mobile devices.BTS (Base Transceiver Station): The electronic heart of the site, managing radio communication with mobile devices, composed of several TRX (transceivers).Technical shelter (Shelter): Container housing electronic equipment (BTS, power supplies, batteries, air conditioning).Air conditioning system: Maintains a stable temperature for the proper functioning of sensitive equipment.Power supply: Ensures continuous operation with backup batteries and often a generator.Grounding & protection: Overvoltage and discharge protection systems (lightning rods).Transmission links: Connection of the site to the central network (fiber optic, microwave link).

According to the work proposed by [31]-[34], such a system has the following operations and architecture:

Cell: Geographic area covered by a base station (BTS).BSC (Base Station Controller): Controls several BTS, manages radio resources and call transfers (handovers) between cells.MSC (Mobile Switching Center): A switching center that manages calls, authenticates users and connects them to the fixed network or other networks.Radio dialogue: The mobile connects to the BTS, which transmits the information to the BSC, then to the MSC to route the call, with appropriate frequencies and power.Authentication: The MSC queries databases (HLR/VLR) to verify the subscriber via the SIM card.To achieve this [34]-[36], the key technical indicators of a GSM relay antenna are:Frequencies: GSM bands 900 MHz, 1800 MHz, etc.Initial speed: 9.6 kbps (GSM 2G), evolving towards higher speeds with GPRS, EDGE.Frame structure: 8 time slots per channel to multiplex users.

A telecommunications site, also known as a GSM site, is a location where a telecommunications operator has installed equipment to form part of its network, thereby providing a range of services to populations living within a certain radius of the site. This location may be shared with other operators [37]. It is characterised by its configuration, which is linked to its urban environment, and by the infrastructure put in place by one or more mobile operators [38]. Open, elevated locations are preferentially chosen to allow for optimal signal propagation and an optimised network [39]. Depending on the situation, this is achieved either by using a metal support structure (pylon) or by utilising an existing building in the surrounding area [40].

Besides the supporting structure, a telecommunications site primarily comprises the following elements, schematically represented in Figures 1-2[40].

Figure 1. Composition of a GSM site [40].

Figure 2. Detailed description of the elements present in a GSM site [41].

2.1.3. GSM Antennas and Essential Characteristics

An antenna is a system that converts electrical energy into electromagnetic energy for reception and transmission. Several criteria are used to describe the characteristics and performance of antennas, such as input impedance, reflection coefficient, directivity, gain, and radiation patterns.

The antenna has multiple roles, the main ones being:

To allow for a proper adaptation between the radio system and the propagation medium.To ensure the propagation or reception of energy in preferred directions.To transmit information in the most perfectly feasible way. Furthermore, to present the performance of antennas, several criteria are used [20][23].

These criteria are classified into two groups. The first group characterizes the antenna as an electrical circuit component with an input impedance and a reflection coefficient (Z_in and S₁), and the other group shows great interest in its radiation characteristics, such as the radiation pattern, directivity, and gain.

Finally, it should be noted that the concept of power (absorbed or radiated) is essential to the study of antennas. An antenna is characterized by various factors that can be ordered either into electrical characteristics or into technical radiation specifications.

2.1.4. Principle of Electromagnetic Radiation from GSM Antennas

According to wave theory, all electromagnetic radiation possesses fundamental properties and behaves predictably [42]. Electromagnetic radiation is composed of an electric field (E) and a magnetic field (H). This radiation can be intentional, as is the case for antennas that emit EM energy to establish wireless communication with radio receivers, or unintentional, as with conductors carrying electrical energy. The electric field varies in magnitude and is oriented perpendicular to the direction of radiation propagation (Figure 3).

The magnetic field is oriented perpendicular to the electric field. Both fields travel at the speed of light (c). The radiation is characterised by several parameters, including the radiation pattern, radiated power, directivity, gain, radiation resistance, polarisation, bandwidth, and quality factor [43].

These parameters are respectively modeled by expressions 1 - 10.

Figure 3. Plane electromagnetic wave [44].

The radiated field at long distances is a function of θ (Site angle: vertical plane) and ϕ (Azimuth angle: horizontal plane). It can therefore be written, up to a factor, in the form:

Or F ( θ , ϕ ) is called the characteristic function of radiation.

The radiation pattern represents the radiation intensity K(θ, θ) as defined later, based on the deflection angles (θ, θ) in space. This representation provides us with the antenna’s most efficient radiation directions.

Radiated power is the power that passes through a sphere of infinite radius. It is determined by integrating the Poynting vector over a spherical surface. The Poynting vector in the radiation zone is defined by:

The power radiated through a sphere of infinite radius is given by:

With 𝑛 being a normal vector to any point on the surface of the sphere. The radiation intensity is given by the expression:

With dΩ, the unit of solid angle. The total radiated power is defined by:

The directivity of an antenna characterizes how that antenna concentrates its radiation in certain directions in space. Directivity is the quotient of the radiation intensity in a direction Δ ( θ , ϕ ) by the average value of this radiation intensity for all directions in space.

An isotropic antenna radiates the same power density uniformly regardless of direction. Gain is a quantity that describes the performance of an antenna. The gain of an isotropic antenna is taken as a unit reference (0 dB).

The gain of an antenna in a given direction is the ratio of the radiated intensity to that of an isotropic antenna.

The direction of maximum radiation is often taken as the direction for deducing the power gain. If ηis the radiation efficiency of an antenna, then P_IN = μP_rad, where P_rad is the total radiated power. The gain is thus written as:

The gain is associated with an equivalent radiation area S_r in the direction u defined by relation (9).

These rules are equally valid for Wi-Fi and ISM.

Let P_r be the active power radiated by an antenna. If it is possible to know the current I_Q at a point Q of this antenna, we define the radiation resistance at this point with respect to:

Electromagnetic compatibility (EMC) is defined as the discipline that deals with the problems of electromagnetic coexistence of electrical equipment located in the same environment [42]. This is to ensure the proper functioning of this equipment without producing disturbances that could disrupt the normal operation of other nearby electrical devices. In EMC terminology, a system is said to be compatible if, on the one hand, it does not generate excessive disturbances in its environment, and on the other hand, if it is able to function correctly in the presence of nearby disturbances.

This definition leads to two EMC concepts, immunity and emissivity (Figure 4), which are associated with conducted and radiated interference. The distinction between these two modes is made according to the transmission medium, which is either the electrical conductor or the ambient air [43].

Figure 4. Concepts of EMC [41]-[43].

According to existing literature [43]-[46], the concept of Electromagnetic Compatibility (EMC) in telecommunications, especially with GSM antennas, concerns the ability of a device to function well in its electromagnetic environment without disturbing others, and vice versa, via radio frequency (RF) waves, governing good mobile-antenna link, the management of disturbances (noise, interference) and health risks, notably through the control of the transmission power of antennas and phones to minimize exposure to non-ionizing electromagnetic fields (EMC).

EMC in the GSM context is crucial to guarantee quality of service (no interruptions, good throughput) while managing exposure to waves, using engineering techniques (antennas, power) and respecting regulatory thresholds to minimize nuisances on the network and on health [47].

Established in 1870, Maxwell’s equations are fundamental laws of physics. They constitute the basic theories of electromagnetism. These equations provide relationships between the variations of the four characteristic vectors of the electromagnetic field: the electric field E(V/m), the magnetic field H(A/m), the electric flux density D (C/m²), and the magnetic flux density B(T) at any point in space.

Time-dependent variations are expressed by the partial derivative with respect to time, and spatial variations by the differential operators: curl and divergence.

At any point in space, which is not located on a surface separating two media, that is to say, in a linear, homogeneous, and isotopic (LHI) medium, Maxwell’s general equations specify that [48]-[51]:

The basic variables of these equations are [49]:

B : Magnetic induction (Tesla, T)

H : Magnetic field strength (Ampere/ meter², Am⁻²)

D : Electric flux density (colomb/meter², cm⁻²)

E : Electric field density (volts/meter, Vm⁻¹)

J : Electric current density (ampere/ meter², Am⁻²)

ρ: Electric charge density (coulomb/meter³, cm⁻³)

With [50]:

And [51]:

µ: Magnetic permeability;

ε: Electrical permeability;

σ: Electrical conductivity.

In this form, called local or differential, Maxwell’s equations express relationships between spatial variations of some fields and temporal variations of other fields.

Into the void:

Then Maxwell’s equations become:

The differential operator nabla ∇ is used to express the curl operation ∇ × = rot and the divergence operation ∇ ⋅ = div .

Maxwell’s equations can also be expressed in “global form” as follows [51][52]:

Using the divergence theorem after integrating Equations (23) and (24), Equations (13) and (14) then become:

The boundary conditions for fields [52] are:

With:

n : is the normal to the separation surface, going from midpoint 2 to midpoint 1;

j S : is the surface current density;

ρ S : is the surface charge density.

Starting from the wave equations below obtained from Maxwell’s equations, for the electric and magnetic fields of interest, the wave equations, at a point r and at time t, are given respectively by [47]-[51]:

where E is the electric field H is the magnetic field and μ₀ and ε₀ are respectively the magnetic permeability and electric permittivity of air (vacuum). The wave equations are written as follows [51]:

where 𝜀 and 𝜇 are characteristic of the medium considered and represent, respectively, the permittivity and permeability of the medium. These quantities allow us to define the propagation speed vand the characteristic impedance of the medium 𝜂 by the following relations:

3.2.1. FDTD Digital Model

However, depending on the type of problem one wishes to study, there are several ways to rewrite Maxwell’s equations [52]-[55]. Today, with the increasing power of computing tools, several numerical tools and methods exist for translating the behavior of these equations. In this study, we implement the FDTD method [53] because it can simulate the behavior of an electromagnetic wave in any type of medium (dielectric, vacuum, metal, plasma, etc.), while taking into account the most complex geometric shapes of the objects that may constitute the system [54].

It does not involve any matrix inversion. Its extremely simple theoretical formulation provides highly accurate predictions for a wide range of problems in the electromagnetic domain [52]-[55]. It is characterized by a broad bandwidth; an impulse excitation in the time domain is sufficient to give the response of a system over a wide frequency range via a Fourier transform [55]. In this method, the unknowns are the electric and magnetic B fields E .

The principle consists of approaching the spatial and temporal derivatives with finite differences using an explicit scheme: this means that at each time step, it is possible to calculate all the derivatives without having to invert matrices [52]-[55].

The points where E and H are calculated are offset by half a step, in space and in time. For each half-step of time, the values of the E and H fields are updated alternately. Using a uniform grid in Cartesian coordinates, the following Figure 6: Δx= Δy= Δz=Δ, the indices in the x, y, z directionsare i, j, k; the time step is Δt and therefore, the time is: n.Δt.

Figure 6. Yee cell.

The centered difference approximation is applied subsequently [52]-[55]:

The dielectric parameters ε, μand σ aretime-independent, so the two previous equations can be written as follows in the system (x, y, z) [52]-[55]:

E_x,E_y, E_z and H_x, H_y,H_z are the components of E and H respectively. Subsequently, we calculate the three components of E and H:

Or:

3.2.2. Hertzian Dipole Model

The use of dipole theory not only simplifies the general expressions of electromagnetic fields radiated by GSM relay antennas, but also allows modeling of EM radiation emitted by a radiating structure from knowledge of the currents involved [56]-[59].

The EM fields radiated by a GSM relay antenna can be estimated by considering an infinitesimal portion of length dl of the line carrying a constant current as a radio dipole. The EM fields created by this antenna are then considered to correspond to the superposition of the fields produced by all the dipoles constituting the wire structure [56]. It is assumed that the current distribution is variable along the antenna, but that the current is constant for each elementary dipole (see Figure 7) [57]. The accuracy of the method is then governed by the number of dipoles required to adequately represent this current distribution [56]-[59]. The greater the number of dipoles, the better the accuracy [56]-[59].

Figure 7. Geometry of a GSM relay antenna [56]-[59].

Consider a wire antenna of length, centered at the origin of the coordinate system and parallel to the axis (𝑧). The antenna in Figure 7 carries a sinusoidal current (I), ω = 2πf where f represents the frequency [56]-[59].

This antenna is considered as a superposition of N elementary dipoles of identical lengths dz such that: L = N∙dz. The observation point P is located at any position in space [58]. The vector potential created by this structure at a given point is simply the contribution of the elementary vector potentials constituting the wire antenna, i.e. [59]:

Similarly, the expressions for electromagnetic fields obtained for an infinitesimal dipole can be generalized to the case of an antenna satisfying the thin-wire assumptions [56]-[59]. These assumptions assume that the current does not undergo azimuthal variation and is uniformly distributed over the surface of the cable. This is valid if the radius is small compared to the wavelength and the distance from the interface. Thus, the electromagnetic fields radiated by this antenna are simply deduced by summing the elementary fields, hence [56]-[59]:

The expansion of the expressions for the elementary fields of Equation (30), obtained in Figure 7, the non-zero components of the electromagnetic fields in a cylindrical coordinate system ( e r , e φ , e z ) , are written for a wire antenna:

3.2.3. Empirical Evaluation Methods

The evaluation of electromagnetic fields (EMF) from GSM antennas using empirical models relies on propagation formulas (free space, Okumura-Hata) that take into account the EIRP (Equivalent Isotropically Radiated Power) [60], antenna gain, distance, and environment to estimate the power density and electric field (V/m), and then deduce the magnetic field. These models predict far-field signal attenuation [61]. These methods are generally used for pre-installation estimations or theoretical checks, with actual on-site measurements required for full accuracy [62].

Empirical models for evaluating the electromagnetic fields (EMF) of GSM antennas, such as propagation models (Hata, COST-231), estimate exposure based on real measurements rather than complex physical simulations. They use simplified formulas incorporating power, antenna gain, frequency, and distance [60]-[62].

Advantages [60]-[62]:

Speed: Fast calculations with low computing resource consumption.Simplicity: Usable for general coverage assessments without complex 3D modeling.Adaptability: Suitable for GSM frequencies (900/1800 MHz) in urban or rural environments.

Weaknesses [63]-[65]:

Limited accuracy: Less accurate than deterministic models (ray tracing) in complex environments.Data dependency: Requires precise technical information (antenna diagrams, tilt) which is sometimes difficult to obtain.Specific cases: Less effective at predicting the far field or the impact of unique obstacles.These models are useful for an initial estimation of overall exposure, but often require validation through field measurements.

The power density vector, Poynting vector S, of an electromagnetic field is given by the cross product of the electric component E and the magnetic component H of the field [63]-[65]:

Ideal conditions [63]-[65], that is, when it is important that the ground or other obstacles have no influence, this equation can be simplified because the electric field, the magnetic field, and the direction of propagation are perpendicular to each other. Furthermore, the ratio of the amplitudes of the electric field, E, and the magnetic field, H, is a constant Z₀_,called the characteristic impedance of free space, equal to approximately 377 Ω (or 120π Ω).

Thus, in the far-field region, the free-space power density, S, is given by the following non-vector equation [63]-[65]:

In general, the characteristic impedance of a medium is given by the formula: Z = μ ε where µ is the magnetic permeability (µ = 1.2566 × 10⁶ F/m in free space) and ε is the permittivity (= 8.85418 × 10¹² H/m in free space) [63]-[65].

The power density regardless of distance and direction can be calculated in the far-field region using the following equation:

where [63]-[65]:

S: power density (W/m²) in a given direction

P: power (W) supplied to the radiation source, assuming a lossless system

G_i: gain factor of the radiation source in the direction considered, relative to an isotropic radiating element

r: distance (m) from the radiation source.

P∙G_i product In Equation (44), EIRP is called the EIRP and represents the power that a fictitious isotropic radiating element would have to emit to produce the same field strength at the receiving point. The antenna pattern must be taken into account for power densities in other directions [66]-[68].

To use Equation (44) with an antenna whose gain G_a is related to a reference antenna of isotropic gain G_r, for example a half-wave dipole or a short unipolar antenna, the gain factor G_i must be replaced by the product G_r·G_a as in Equation (44). The relevant factor G_r is given in [60]; it is 1.64 for half-wave dipole applications (television, broadcasting on metric waves and sometimes on decametric waves) and 3.0 for short unipolar antennas (broadcasting on kilometric, hectometric, and sometimes decametric waves) [66]-[69]:

Thus, when the antenna gain G_d (G_a = G_d) is referred to that of a half-wave dipole [63]-[65]:

where: G_d: antenna gain relative to a half-wave dipole.

Similarly, when the antenna gain G_a = G_m is related to that of a short unipolar antenna [63]-[65]:

where: G_m: antenna gain relative to a short unipolar antenna.

Equations (42) to (46) assume that the conditions are those specific to the far-field region where the field is in the form of a plane wave; they do not concern calculations for the near-field region. If Equation (42) is inserted into Equation (44) to eliminate S and a factor C is introduced to account for the direction characteristic of the radiation source, Equation (47) is obtained, which allows the electric field of a radiation source to be calculated in the far-field region [63]-[65]:

Or:

E: electric field (V/m);

Z₀ = 377 Ω, characteristic impedance of free space;

P: power supplied to the radiation source (W), assuming a lossless system;

C: factor (0 ≤ C ≤ 1) which takes into account the direction characteristic of the radiation source (in the principal direction of the radiation, C = 1).

If the antenna gain is referred to that of a half-wave dipole or a short unipolar antenna and not to that of an isotropic radiating element, it is appropriate to use instead of G_i the factors G_d and G_m respectively in Equations (49) and (50) [63]-[65].

To calculate the magnetic field, in the far-field region, of a radiation source, Equation (51) [63]-[65] is used:

Or:

E: electric field (V/m);

H: magnetic field (A/m);

Z₀ = 377 Ω (120π), characteristic impedance of free space.

The equipment used for geolocating measurement points and quantifying electromagnetic fields, taken from existing literature [65]-[68], is:

Location systems: GPS to geolocate measurement points [65]. For confirmation of the choice of measurement areas, the HF DIGIMETER Endotronic D8826 can be used, which is an electromagnetic radiation meter focused on measuring noise voltage.Electric (E) and magnetic (B) field probes: Specific to telecommunications frequencies (from a few MHz to several GHz), often with triaxial probes to cover all directions [66].Spectrum analyzers: To identify the different sources (operators, technologies) and their respective frequencies [67].Software: For data acquisition, processing and results mapping [67].

According to the ANFR/COMSIS framework [67], a good RF/GSM electromagnetic measurement protocol should include the following aspects:

Definition of zones: Residential premises, public places (parks, shops), places accessible to the public.Antenna location: Use of Cartoradio (ANFR website) to identify sources.On-site measurements [65]-[67]:Measurement points: At height (1.5 m) and in the corners of rooms for the interior; at various locations outside.Scenarios: Measurements in real-world conditions (low/high number of users) and theoretical extrapolation (all antennas at maximum power).Frequencies: Taking into account the different bands (GSM, 3G, 4G, 5G).Calculation of limit values: The results are compared to the limit values set by decree (e.g., 5 V/m at 900 MHz for the public).Analysis of results [65]-[67]:Comparison with thresholds: Verification that exposure remains below regulatory limit values.Mapping: Creation of exposure level maps by area.Reports: Preparation of detailed reports by accredited bodies (such as Apave or others) to assess compliance with standards.

The metric method describes Data Analysis (or Statistical Analysis), which uses systematic techniques (statistical and logical) to describe, structure, condense, visualize (images, tables, graphs) and evaluate biases in data, in order to draw meaningful conclusions, often combining descriptive and inferential statistics to understand phenomena and make predictions [69].

We reiterate that one of the objectives of this article is to contribute to the definition of an exposure indicator by jointly exploiting aspects related to measurements and simulations. One of the means used to achieve this objective is to focus on indicators for evaluating errors between measurements and simulations [69]. Given the number of exposure situations that may need to be studied, finding a relevant comparison indicator suitable for all exposure cases is not a trivial matter. The question of the criterion for evaluating prediction error has therefore become a crucial problem in the search for an exposure indicator to enable [69]:

to conduct a mutual validation of the methods and protocols used;analyze the uncertainties related to the calculation methods;to best characterize all exposure configurations;develop lower cost monitoring systems.

In this section, we will present a state-of-the-art review of reliability indicators for calculation methods commonly used in radio wave simulation. From this, we will identify a relevant criterion for evaluating prediction errors to meet all our requirements.

The criteria for comparing prediction methods with measurements generally used in radio wave simulation can be classified according to 3 types of criteria:

The “error” type criterion;The “link” type criterion;The criterion of type “characteristics of the data distribution”.

Equation (52) presents the equations modeling the criteria described above. For further details, we recommend that researchers consult the bibliographic references [69]. These indicators allow for the search for a linear relationship between measurements and simulations. However, depending on the variability of the data distribution, the results may lead to erroneous conclusions. They are better suited to normally distributed data without exceptional values.

Therefore, in the context of evaluating methods for simulating exposure, due to the high variability of exposure, their use is not appropriate, especially in exposure cases.

This is because a very high value and a very low value can be found at two relatively close points, as everything depends on the point’s position relative to the antenna and the medium through which the wave propagates.

5.1.1. GSM Sites

See Table 1 below:

Table 1. GSM site characteristics.

5.1.2. Field Measurements

Beyond the responses, most of which were subjective, obtained in the user survey, a more objective quantification of the exposure level seemed necessary. We therefore conducted field measurements in targeted areas, without which comparison to relevant standards would be impossible.

5.1.3. Selection of Measurement Zones

The measurements were carried out in the Gombe district, identified not only as the city center but also as the area most exposed to non-ionizing radiation (NIR) due to the presence of numerous telecommunications antennas. This location is of particular interest because it encompasses a high concentration of socio-economic activity. Consequently, the number of potential subscribers exceeds that of other districts and even the entire Lukunga district. Numerous GSM sites are installed there, justifying the importance of this study area.

The choice of measurement zones is based on two fundamental criteria: a geometric criterion and a demographic or urban planning criterion. On the one hand, the cellular subdivision inherent to GSM technology allows for the concentration of non-invasive radiofrequency (NIR) in certain specific areas of the municipality. On the other hand, areas with high population density, as well as those exhibiting high sensitivity characteristics, such as hospitals, were also selected due to their increased vulnerability to the effects of NIR.

To delineate these zones, a satellite image of the Gombe district was obtained using Google Earth (version 5.2.1.1588). Two orthogonal axes were drawn, one parallel to Tshatshi Avenue, and two diagonals connecting the four corners of the Golf neighborhood. These axes defined the areas of interest for the measurements.

On each axis and diagonal, measurement points were selected, with between two and four points per axis, excluding the central point (zone 9). Thus, zones 1, 5, 10 and 12 are located on diagonal A; zones 3, 7, 11 and 13 on diagonal B; zones 2, 6 and 14 on the minor axis; and zones 4 and 8 on the major axis.

The geodetic coordinates of the measurement points were extracted directly from Google Earth. A field validation campaign was then carried out using a Garmin eTrex H GPS receiver, allowing for confirmation or adjustment of the coordinates, particularly based on the accessibility of the selected sites. Figure 8 shows the cadastral map detailing the measurement areas.

Radiation measurements require substantial equipment. In addition to the GPS mentioned earlier for confirming the choice of measurement areas, we used the HF DIGIMETEREndotronicD8826. This is an electromagnetic radiation meter focused on measuring noise voltage.

The technical specifications of this equipment are presented below:

Receiver principle: Wideband linear receiver measuring the sum of peak values of resonant amplitude modulations with 3D sound performance.Frequency range, measured electromagnetic noise voltage and linearity: 5 KHz - 8 GHz distributed across 3 inputs:Middle BNC input: 5 kHz to 8 GHz (linear 4 GHz) if their intensity is greater than 10 mV;Left BNC input: Same frequency range with a voltage greater than 30 mV;On the right BNC input: very low noise (minimum 40 μV) can still be measured from 10 MHz (linear), and from 3 MHz (non-linear).RF signal sensitivity: ± 3 μV, linear for electromagnetic fields;Acoustic output with 50 amplifier;Antennas: the device is equipped with 5 antennas including a 15 × 1 mm, a 90 × 20 mm, a 45 × 1 mm, a 980 × 10 mm and the AS2002 antenna;Battery life: 9 V, 12 mA alkaline battery with approximately 2 hours of performance;Power consumption: between 40 and 80 mA, depending on the sound level;Weight: ± 200 g;Dimensions: 145 × 76 × 30 mm (L × W × H);Uncertainty;Calibration: The device was supplied to us already calibrated. A conversion chart between the displayed values (VA) for the three inputs and the corresponding electromagnetic noise voltage values (in volts) was included in the user manual. Since no standard is expressed in terms of noise voltage, we sought the assistance of overseas experts to obtain conversion curves between the displayed values and the power density and electric field values obtained by comparison with other field strength meters and spectrum analyzers.

This study considers 50 measurement points, comprising 29 outdoor and 21 indoor points, taken from 30 selected sites. At each point identified by its geographic coordinates, we delimited an area of 8 m²; the GPS used for location had a resolution of ± 4 m. Once the perimeter was defined, we first traversed the corresponding area, device in hand, in a vertical position at a distance of 1 to 2 m from our support (ground or concrete, depending on the location), guided by sound, in search of the point with the highest level of electromagnetic noise. In certain accessible areas, we took several measurements, both indoors and outdoors, but sometimes also at different heights, particularly on buildings.

We then measured the highest value displayed on the screen using the 45 × 1 mm antenna (vertically polarized), first placing it at the right input and then at the center input if the value exceeded the limit measured at the first input. Finally, we used calibration curves to determine the electric field strength and power density values corresponding to each displayed value.

Electromagnetic noise levels, and therefore electric field and power density levels, were measured in all selected areas between July 6th and 14th, 2022, from 9:00 AM to 5:00 PM. Due to the lack of meteorological equipment, we relied on the websites of TV5 and MSN, whose weekly weather forecasts (parameters relevant to this study) for the period in question were described as follows:

Temperature: 21˚C - 33˚C;Relative humidity: 70% - 80% from 9 a.m. to 5 p.m.;Sunrise: 06:06;Sunset: 17:59.

See Figures 8-16 below:

Figure 8. 2D simulation results of the spatial distribution of the electric field on 21 points of the external sample.

Figure 9. 2D simulation results of the spatial distribution of the electric field on 29 points of the internal sample.

Figure 10. Estimated 2D electric field magnitude as a function of simulation positions [Achraf Liakouti, 2018] [67].

Figure 11. 3D simulation result of the electric field using Hertzian dipole theory.

Figure 12. 3D electric field magnitude estimated as a function of measurement position using dipole theory [Achraf Liakouti, 2018] [67].

Considering the characteristics of the GSM relay antennas used, those of the measuring instruments used and the climatic conditions of the sample in the chosen geographical area, we obtained the results presented in Figures 13-23, Tables 1-3 below:

Figure 13. Calibration curve of the right-hand input.

Figure 14. Calibration curve of the center input.

Figure 15. Modeled mapping for assessing the distance between the site and sensitive areas of use.

Figure 16. Modeled map of the distribution of measurement points.

Table 2.Statistics of the 50 power density and electric field measurements corresponding.

Figure 17. Statistics of the 50 measurements of power density and corresponding electric field.

Figure 18. 2D electric field residual graph (external measurement).

Figure 19. 2D spatial distribution of the electric field (external measurement).

Figure 20. 2D electric field probability distribution (external measurement).

Figure 21. 2D residual graph of the electric field (Indoor Measurement).

Figure 22. 2D spatial distribution of the electric field (internal measurement).

Figure 23. 2D electric field probability distribution: (internal measurement).

Table 3. Joint comparison of 2D simulation and external measurement E (V/m).

Table 4. Joint comparison of 2D simulation and interior measurement E (V/m).