Analysis of the Arbovirosis Potential Occurrence in Dobrogea, Romania

Abstract

Climate change creates new challenges for preventing and protecting human health against different diseases that could appear and propagate. The Aedes albopictus mosquito species is an important vector for different diseases like dengue fever or zika. Although this species is not “indigenous” in Europe, its presence is noticed in many countries on the continent. The Ae. albopictus establishment is conditioned by the species’ characteristics and environmental factors. To assess the possible spread of Ae. albopictus in the Dobrogea region (situated in the Southeast of Romania), we conducted the following analysis: (1) Investigation of the current distribution and climatic factors favoring Ae. albopictus’ establishment in Europe; (2) Analysis of climate dynamics in Dobrogea in terms of the parameters identified at stage (1); (3) Testing the hypothesis that the climate from Dobrogea favors Ae. albopictus’ establishment in the region; (4) Building a Geographic Information System (GIS)-based model of the potential geographic distribution of Ae. albopictus in Dobrogea. Results show that the climate of Dobrogea favors the apparition of the investigated species and its proliferation.

1. Introduction

During the last years, weather and climate changes have a notable impact on human health [1,2,3,4]. According to the CDC data recently published, climate change has led to an augmentation of the incidence of emerging and re-emerging diseases, out of which vector-borne, food- and water-borne disease are the most important [5].

According to the WHO [6], vector-borne diseases (whose most common vectors are mosquitos) account for more than 17% of all infectious illnesses, causing more than 700,000 deaths annually. Globally, 146 (58.4%) countries/territories reported at least one arboviral disease, while 123 (49.2%) reported more than one arboviral disease, in many cases, local outbreaks [7,8,9].

For instance, dengue fever (caused by mosquito bites, like Ae. albopictus (Skuse, 1895)) is common in many countries from tropical and subtropical areas. The Invasive Species Specialist Group considers that Ae. albopictus is one of the worst 100 invasive species due to its adaptability [10,11]. Thus it comes as no surprise that autochthonous dengue cases were documented in southern France in 2010, 2013, 2014, and 2015 [12,13,14]. In Spain, the first local outbreak of Dengue was registered in 2018. The second one was reported in September 2019, when the local health authorities announced a laboratory-confirmed autochthonous dengue case in Barcelona. The presence of Ae. albopictus in the same country has been reported since 2004 [15]. Autochthonous dengue cases have been reported in Europe (Croatia, France, Madeira Islands) and the United States (Hawaii, Florida, and Texas) [12,13], as well. The number of imported cases has also increased in Germany, Belgium, Italy, Spain, and the United Kingdom, due to the tourism intensification in endemic countries [14,15,16,17,18]. Imported dengue fever cases have been notified in Romania since 2008 [16], their number significantly increasing in 2019 (Figure 1).

Following the outbreaks of dengue and chikungunya virus in Europe, European Centre for Disease Prevention and Control (ECDC) attempted to achieve a current distribution map and provide a historical spread of these mosquito species in Europe. According to the report published in 2009 [19], the North Mediterranean coast was already infested. At that time, there was no information concerning the Ae. albopictus spreading in Romania.

Scientists proposed different approaches to investigate the Ae. albopictus spatial distribution in Europe. Fisher et al. [20] published a review containing the principal models used to predict this species’ future distribution and classified them into two categories—mechanist and correlative. Since none of them is better than the other one, ECDC utilized the Random Forest model (correlative approach) and Multi-Criteria Decision Analyses (MDCA—mechanist approach) in its technical report [16]. It was shown that the Random Forest is suitable to describe mosquitos’ spread in the Mediterranean area or similar zones. According to both the maximum impact short- and long-term scenarios, entire Europe is suitable for Ae. albopictus establishment. The MDCA result (based on Intergovernmental Panel on Climate Change—IPPC’s long-term climate change scenario for 2030 with minimal impact) shows that Romania becomes unsuitable for Ae. albopictus spreading [19].

Recently, Kraemer et al. [21,22] and Kamal et al. [23] mapped the spreading zones of the Ae. albopictus vectors responsible for transmitting the major human arboviral diseases. According to [23], the Southern part of Romania is 50% suitable for the existence of this species, but the situation from Dobrogea (southeastern part of Romania) has yet to be reported. The reconstructed global distribution of Ae. albopictus [22] provides 50% suitability for the establishment of this mosquito species in Romania.

Prioteasa et al. [24] reported for the first time the existence of Ae. albopictus in Romania, in Bucharest, in 2012. Since then, its presence has been noticed every year. Recently, Fălcuţă et al. [25] investigated the existence of Ae. albopictus in 53 localities from Romania for the period 2017–2018. Their results confirm the presence of this invasive Aedes mosquito species, especially in the Bucharest metropolitan zone. ECDC 2020 map [26] shows the current distribution of Ae. albopictus, which takes into account the results from [25]. According to this document, Ae. albopictus is already established in Constanta county (situated in the south of Dobrogea) but was not recorded in Tulcea county (situated in the northern part of Dobrogea). That means that the vector has been confined to the country, and there is a risk of a local outbreak of Dengue fever soon, in the absence of drastic control measures.

In this context, this paper aims to (i) update the understanding of the current distribution and environmental factors favoring the Ae. albopictus’ establishment in different European regions by compiling the data from the literature and (ii) investigate the possibility of the future establishment of this mosquito species in Dobrogea, Romania, given (that rural and urban zones represent 30% in Constanta county and 60% in Tulcea county, that belong to the study region) and the climate could foster this species survival [27,28]. A Geographic Information System (GIS) model based on five climate criteria is proposed to identify the favorable areas for Ae. albopictus establishment.

2. Materials and Methods

2.1. Study Region

Dobrogea region (Figure 2) is situated in the Southeast of Romania between the Black Sea (East), lower Danube (West), and Danube Delta (North). From an administrative viewpoint, this territory, with 16,501 km², is divided into two counties, Constanța and Tulcea. From a geomorphological point of view, the relief of Dobrogea contains the Dobrogea Plateau, the Danube floodplain, and the Danube Delta. The average altitude is about 200–300 m, the highest point being Tutuianu Peak (467 m).

From a hydrological viewpoint, Dobrogea is divided into two basins, one tributary to the Danube River and the other tributary to the Black Sea basin. According to the Management Plan of Dobrogea-Littoral Watershed, the total length of rivers in Dobrogea is 1623 km of rivers; there are 18 lakes and four reservoirs. Wetland areas represent 23% of the region.

The climate of Dobrogea is temperate continental, with two zones. The first one contains the Danube Delta, the two lagoons (Razelm and Sinoe), and the Black Sea Littoral. This unit covers the territory between 20 and 50 km from the Romanian Littoral (depending on the warm or cold season) and has a continental climate. The Black Sea influences are felt during the cold season, the temperatures remaining positive up to an altitude of 100 m. During the warm season, the climate is affected by sea breezes. The second unit covers the western part of Dobogea, where the moderate continental belt influences are felt [29].

The mean annual temperature varies between 9.9 °C at Jurilovca village and 11.7 °C at Constanta. The mean multiannual precipitation varies from 261 mm at Sulina station to 488 mm at Cernavoda.

Different climatic parameters are investigated in the present study, covering the period after 1965. Data series was recorded at the ten principal meteorological stations in the region.

2.2. Materials and Methodology

As presented in the Introduction, Ae. albopictus is an important vector for different diseases. This mosquito species, native to the Southeast Asia region, appeared in Europe in the last decades of the XX century. As Lockwood et al. [30] mentioned, four stages are necessary to declare a species as “invasive”: introduction, establishment, expansion, and impact on the environment. A species can be transported from its original area to a new ecosystem by natural phenomena or humans. This article will address only the second stage (establishment), which is conditioned by environmental and species characteristics. To assess the possible spread of Ae. albopictus in Romania, especially in the Dobrogea region, we conducted the following steps.

• Analysis of the current distribution and climate factors favoring the establishment of Ae. albopictus in Europe.

To achieve this objective and evaluate the habitat of the Ae. albopictus in Europe, we investigated 67 locations where the presence of this species was reported. The selection of location is based on the scientific reports and literature [11,19,25,26,31,32,33,34,35]. The following climatic parameters are analyzed: mean multiannual temperature (°C), mean maximum multiannual temperature in the warmest month (°C), mean multiannual minimum temperature (°C) of the coldest month (generally January), mean multiannual precipitation (mm), and mean multiannual relative humidity (%). The climatic parameters investigated were obtained via [36]. For each of them, we computed the basic statistics and built spatial distribution maps in ArcGIS 10.7 (ESRI).

2.

Understanding the physiological plasticity of Ae. albopictus.

The objective of this direction is to identify the principal climate parameters favoring the establishment of Ae. albopictus. A review of the scientific literature concerning the physiological plasticity of Ae. albopictus [37,38,39,40,41,42] demonstrates this species ability to adapt to different climate conditions.

3.

Analysis of climate dynamics in Dobrogea in terms of the parameters identified in the first two steps. The aim is to determine the suitability of the climate for Ae. albopictus establishment.

Studies regarding the temperature and precipitation dynamics in the Dobrogea region [43,44,45,46,47,48] demonstrate that after 1997–1998 and 1995, respectively, increasing temperature (+0.8 °C on average) and rainfall (88–98 mm) trends are noticed as a result of climate change. Consequently, in this paper, the data records after 1997 (for temperature) and 1995 (for precipitation) are used. Data series provided by the Romanian National Meteorological Agency, registered at ten meteorological stations (Adamclisi, Cernavodă, Constanța, Corugea, Hârșova, Jurilovca, Mangalia, Medgidia, Sulina, and Tulcea) are utilized in this study. The geographical positions of the study locations are presented in Figure 2.

Firstly, we tested the hypothesis that the mean annual temperature series are normally distributed by using the Kolmogorov-Smirnov test. Secondly, we performed statistical tests to check different hypotheses related to the favorable climate conditions for Ae. albopictus establishment.

4.

Building a GIS-based model for analyzing the spatial correlation between the climatic parameters identified.

Figure 3 contains a scheme of the proposed model. Discrete climatic data processing was conducted via IDW (Inverse Distance Weight) interpolation, with the default parameter (beta = 2) in order to obtain the spatial maps. The analysis is based on the Weighted Overlay, a component of spatial modeling that uses the Multicriteria Evaluation (MCE) method integrated with GIS. To create an integrated analysis, a common measurement scale must first be established. Each raster file was reclassified into suitability units (for example, 1 to 6). The value 1 was assigned to the most suitable location and the highest value to the most disadvantageous one. After data homogenization, the Weighted Overlay operations were used. Weights were assigned to each class of each data set. Finally, the spatial distribution of the suitable zones for Ae. albopictus establishment was obtained. In this study, we built the model assuming that all datasets were equally important for identifying the most suitable area.

3. Results and Discussion

3.1. Analysis of the Current Distribution and Climate Factors Favoring the Establishment of Ae. albopictus in Europe

To evaluate the habitat availability for Ae. albopictus proliferation in Europe, we investigated 67 locations where this species’ presence was reported.

Table 1. Geographical distribution and the climatic parameters for locations where the presence of Ae. albopictus was reported. Mean multiannual temperature (MAT) (°C), mean multiannual maximum temperature (MAMaxT), mean multiannual minimum temperature (MAMimT) (°C), mean multiannual precipitation (MAP) (mm), mean multiannual humidity RH (%). The minimum and maximum values of the investigated parameters are highlighted in blue and yellow respectively.

No.	Location	Latitude, Longitude	MAT (°C)	MAMaxT (°C)	MAMinT (°C)	MAP (mm)	RH (%)
1	Durrës	41.3246° N, 19.4565° E	15.9	28.8	4.8	1064	67
2	Genoa	44°24′40″ N 8°55′58″ E	14.7	27.3	3.5	1086	66
3	Veneto	45.4415° N, 12.3153° E	13.2	27.4	0.2	830	76
4	Lombardia	45.4791° N, 9.8452° E	13.1	29.3	−0.8	1013	75
5	Emilia-Romagna	44.5968° N, 11.2186° E	14	30.3	−0.6	774	70
6	Toscana	43.7711° N, 11.2486° E	14.5	30.9	1.9	864	70
7	Piemonte	45.0522° N, 7.5154° E	12.6	29.3	−1.7	846	75
8	Sardegna	40.1209° N, 9.0129° E	16.2	29.2	6.1	419	73
9	Roma	41.9028° N, 12.4964° E	15.7	30.6	3.8	798	75
10	Nice	43.2154 N, 6. 5059 E	10.4	23.3	0.3	721	76
11	Corsica	41.9192° N, 8.7386° E	15.2	26.9	5.8	638	79
12	Orne	48.6389° N, 0.0848° E	11.9	27.3	−0.5	781	79
13	Vienne	45.5256° N, 4.8743° E	9.9	26.1	−3.7	623	76
14	Paris	48.8566° N, 2.3522° E	11.3	24.8	0.7	637	78
15	Vrasene	51.2194° N, 4.1945° E	10.2	21.9	0.1	770	83
16	Oost-Vlaanderen	51.0362° N, 3.7373° E	10.2	21.5	0.4	754
17	Podgorica	42.4304° N, 19.2594° E	15.4	32.6	1.5	1631	64
18	Ticino	46.3317° N, 8.8005° E	11.4	25.4	−0.6	1364	67
19	Chiasso	45.8367° N, 9.0246° E	11.9	26.8	−1.1	1267	67
20	Sant Cugat del Vallès	41°28′ N, 2°4′ E	16.1	27.9	6.0	596	72
21	Zagreb	45.8150° N, 15.9819° E	11	27.3	−3.2	930	75
22	Podobuče	42.9471° N, 17.2865° E	16.3	30.5	4.9	1073
23	Orebić	42.9758° N, 17.1779° E	16	30.1	4.7	1035
24	Korčula	42.9297° N, 16.8886° E	16.3	30.3	5.1	1004
25	Pelješac	42.8653° N, 17.5505° E	14.2	27.9	2.9	1290	64
26	Corfu	39.6243° N, 19.9217° E	16.9	31.2	5.7	1146	70
27	Igoumenidsa	39.5061° N, 20.2655° E	16.7	31.9	4.7	1108	70
28	Haarlemmermeer	52.3004° N, 4.6744° E	9.2	20.6	−0.1	805	84
29	South-Holland	52.0208° N, 4.4938° E	9.6	20.9	0.1	682
30	Noord-Holland	52.5206° N, 4.7885° E	9.2	20.6	−0.1	805	84
31	Utrecht	52.0907° N, 5.1214° E	9.3	21.3	−0.5	804	82
32	Heijningen	51.6559° N, 4.4125° E	9.8	21.3	0.1	776	81
33	Oosterhout	51.6410° N, 4.8617° E	9.6	21.5	−0.2	791	81
34	Montfoort	52.0362° N, 4.9519° E	9.4	21.3	−0.3	802	82
35	Weert	51.2439° N, 5.7142° E	9.7	21.9	−0.4	775	82
36	Banja Luka	44.7722° N, 17.1910° E	11	27.9	−3.2	996	75
37	Lanžhot	48.7244° N, 16.9670° E	9.5	25.6	−4.5	620	75
38	Ladná	48.8054° N, 16.8723°E	9.5	25.5	−4.4	591	75
39	Mikulov	48.8053° N, 16.6377° E	9.2	25.2	−4.6	586	75
40	Baden-Württemberg	47.9958° N, 7.8522° E	9.3	23.4	−2.7	674	77
41	Monaco	43.7384° N, 7.4246° E	14.8	26.6	4.4	811	76
42	San Marino	43.9424° N, 12.4578° E	11.8	26	0.1	805	78
43	Bucuresti	44.4268° N, 26.1025° E	10.8	28.6	−5.5	598	70
44	Sochi	43.6028° N, 39.7342° E	14.5	27.1	3.2	1514	75
45	Wurzburg	49° 47′ 0″ N, 9° 56′ 0″ E	9.5	24.2	−2.6	603	77
46	Freudenstadt	48°27′48″ N, 8°24′40″ E	7.8	22	−4.4	1024	77
47	Essen	51°27′ 3″ N, 7°0′47″ E	9.9	22.7	−0.6	843	78
48	Jena-Lobeda	50°55′38″ N, 11°35′10″ E	8.6	22.8	−2.5	565	78
49	Heidelberg-West	49°25′0″ N, 8°43′0″ E	10.2	24.6	−1.7	666	77
50	Freiburg-East	47°59′0″ N, 7°51′0″ E	10.4	25	−1.1	887	77
51	Penafiel	41°12′0″ N, 8°17′0″ W	12.1	28.5	−0.1	434	77
52	Guilhufe	41°12′14″ N, 8°16′38″ W	12.1	28.5	−0.1	434	77
53	Urrô	40°55′34″ N, 8°17′35″ W	14.5	26.4	4.5	1162	77
54	Gibraltar	36°8′0″ N, 5°21′0″ W	17.8	28.3	9	729	72
55	Batumi	41°38′45″ N, 41°38′30″ E	14.2	26.2	2.8	2393	74.5
56	Bagratashen	41°14′45″ N, 44°49′16″ E	13.3	31.4	−3.1	444	64
57	Tel Aviv	32°4′0″ N, 34°47′0″ E	20.2	32	8.7	562	71
58	Jerusalem	31°47′0″ N, 35°13′0″ E	17.2	30.5	4.9	474	53
59	Haifa	32°49′0″ N, 34°59′0″ E	21	32.2	9.3	525	63
60	Artvin	41°11′0″ N, 41°49′5″ E	11.8	25.7	−1.8	1168
61	Rize	41°1′29″ N, 40°31′20″ E	14.3	25.7	3.4	1860
62	Trabzon	41°0′ 8″ N, 39°43′21″ E	14.4	25.7	3.9	891	72
63	Oradea	47°4′20″ N, 21°55′16″ E	10.6	27.4	−4.9	600	77
64	Ploiesti	44°56′0″ N, 26°2′0″ E	10.3	27.6	−6.2	588	77
65	Dubova	44°37′0″ N, 22°16′0″ E	11.1	27.9	−3.3	619	75
66	Constanta	44°10′0″ N, 28°38′0″ E	11.6	25.8	−2.3	423	81
67	Petru Rares	44°05′50″ N, 25°47′29″ E	10.9	28.6	−5.3	595	77
	Average		12.56	26.59	0.58	850.5	74.81
	Max		21	32.6	9.3	2393	84
	Min		7.8	20.6	−6.2	419	53

presents the geographical distribution and climatic parameters for the territories where the presence of Ae. albopictus was noticed.

Our study shows that Ae. albopictus survives in Europe in the following conditions:

➢

According to Figure 4, the geographical areal boundary is between the latitudes of 52.3004° N (Netherlands—Haarlemmermeer region) and 36°8′0″ N (Gibraltar), and the longitudes of 8°17′0″ W (Penafiel area—Portugal) [32] and 44°49′16″ E (Bagratashen—Armenia) [34]. The limits are debatable, especially in the North and South.

●

In The Netherlands, Ae. albopictus was discovered in greenhouses at 52.3004° N, even if some adults were found outdoors [49]. In Belgium, Shaffner et al. [50] consider that the area surrounding the observation site (Oost-Vlaanderen province) favors the mosquitos’ spread because they found immature stages of this species during the study period. The authors specified that the larvas and pupas came from The Netherlands and “have been on-site for 4–5 months” [50]. Considering the map provided by Kraemer [22], we should consider Jena (50°55′38″ N latitude) or Essen (51°27′3″ N), Germany as the northern limit of Ae. albopictus areal.

●

Information about the presence of the study species South of Gibraltar is not available. However, Leshem et al. [33] reported in 2012 suitable conditions for the autochthonous transmission of dengue in Israel, meaning that the Ae. albopictus was already installed.

➢

According to [26], most areas are located in the northern part of the Mediterranean Sea coast and on the Adriatic Sea Littoral (Figure 4).

➢

The mean annual temperature in the zones where the Ae. albopictus appeared in Europe is 12.56 °C. The territories where the Ae. albopictus were detected are located between the 7.8 °C and 17.8 °C isotherms (Freudenstadt-Germany and Gibraltar, respectively). In most of them, the temperatures vary between 9 °C and 13.8 °C (Figure 4).

●

The average multiannual maximum temperature is 26.59 °C, with variations between 32.6 °C (Podgorica—Montenegro) and 20.6 °C.

●

The average multiannual minimum temperature is 0.58 °C, with variations between –6.2 °C (Ploiesti, Romania) and 7 °C (Gibraltar) or 9.3 °C (when considering the data series from Israel, as well).

Notice the existence of Ae. albopictus in Bucharest, in a temperate continental climate. Prioteasa et al. [21] show that the eggs survived two winters consecutively during the observation period (2013–2015).

➢

The average multiannual precipitation is 850.5 mm, with variations from 419 mm (Sardinia—Italy) to the 2393 mm (Batumi—Georgia) isohyets (Figure 5).

➢

The average multiannual humidity is about 75%, varying between 64% (Croatia and Bagratashen—Armenia) and 84% (when considering The Netherlands as well) or 81% (Constanta—Romania, without The Netherlands data). Taking into account Israel, the minimum is registered at Jerusalem—53%.

3.2. Understanding the Physiological Plasticity of Ae. albopictus

Table 2. Climate condition for development stages of Ae. albopictus.

Phase	Temperature (°C) (4)				Total Rainfall (2) (mm)	Day Length (3) (h)
	No	Min	Optimal	No
Eggs hatching	<10	10–11	20–25	>25	>500	11–11.5
Larval development	<5–10		25–30	>40		13–14
Eggs survival (1)		(−2)–0
Adults survival			25	>35		9
Reproduction			25–30	>35

⁽¹⁾ Thomas et al. [42] proved that the disposed eggs resist 24 h at a temperature of −10 °C in laboratory conditions. ⁽²⁾ Mitchell [40] considers that the Ae. albopictus could not survive in a climate with less than 300 mm annual precipitation. ⁽³⁾ According to [40,51], the optimal photoperiod (day length) is 13–14 h but eggs could hatch at 11–11.5 h daily length. The adults could survive at 9 h day length, at the temperature of 10 °C [51]. ⁽⁴⁾ According to [37,38,39,40,49].

summarizes the climate conditions for the different development stages of Ae. albopictus, retrieved from the above-mentioned literature.

3.3. Analysis of Climate Dynamics in Dobrogea

Figure 4 and Figure 5 show that Romania is located in the region where Ae. albopictus can proliferate. In the context of the analysis presented below, firstly, we investigated the mean multiannual temperature (tm), the mean multiannual minimum temperature in January (TJan), the mean multiannual maximum temperature in July (TJul), the mean multiannual precipitation (Rm), and the mean multiannual of relative humidity (RH).

Figure 6 presents the dynamics of temperature and precipitation in Dobrogea. The mean multiannual temperature (tm) varied between 10 °C (at Corugea station) and 12.3 °C (Constanta), with a regional average of 11.6 °C. TJan varied between −4.1 °C (Corugea station) and −1 °C (Mangalia). TJul took values between 30.2 °C (Cernavoda) and 26.6 °C (Sulina). The warmest month was July and the coldest one was January.

The multiannual precipitation (Rm) varied between 205 mm (Sulina) and 568 mm (Adamclisi). The average relative humidity (RH) in the region was 78.4% (77.8–84.2%). The day length varied between 8.9 h in December and 15 h in July.

To deeply investigate the suitability of climate conditions for the Ae. albopictus existence in the Dobrogea region, we performed some statistical tests on data series.

The Kolmogorov–Smirnov normality test on the annual precipitations series could not reject the hypothesis that the data series are normally distributed because all the p-values (p-val) were greater than 0.05 (which was chosen as significance level).

We used the t-test for checking the null hypothesis H₀: The annual precipitation is equal to 500, against its alternative, H₁: The annual precipitation is greater than 500. The value of 500 (mm) is the minimum value of annual rainfall at which the eggs survive (Table 2).

The null hypothesis was rejected for Jurilovca (t(24) = −2.41, p-val = 0.076) and Sulina series (t(24) = −11.89, p-val <0.001), whereas there was not enough evidence to reject the hypothesis that the annual precipitation is higher than 500 mm for the other eight series (all the p-val being greater than 0.05). This means that the precipitation conditions favor the installment of this species in eight out of ten of the studied locations.

We fitted the best probability density distributions for the mean annual temperature series, and, based on them, we computed the probabilities that the mean annual temperatures fall in the interval 10–25 °C (minimum and maximum optimal temperatures for eggs hatching). After fitting probability density distributions for the mean temperatures in July, we computed the probabilities that the temperatures fall between 20 °C and 25 °C (optimal temperatures for the eggs hatching).

Table 3. The probability that the mean annual temperature (T-mean annual) is between 10 °C and 25 °C, and the probability that the mean temperatures in July are in the interval 20 °C and 25 °C.

	P(10 < T-Mean Annual < 25)	P(20 < T-Mean July < 25)
Adamclisi	0.97633	0.979
Cernavoda	0.98311	0.4053
Constanta	0.992	0.8268
Corugea	0.8714	0.8361
Harsova	0.9892	0.9781
Jurilovca	0.9176	0.8855
Mangalia	0.9934	1
Medgidia	0.9668	0.9908
Sulina	0.995	0.8551
Tulcea	0.9829	0.8435

contains the results of the tests.

Since the probabilities that the mean annual temperatures belong to the interval 10–25 °C have high values, it results in that they should be suitable for the eggs hatching. The probabilities that the mean temperatures in July are above 0.8268 (so, very high) at nine out of ten study locations indicating that the temperatures are optimal for the eggs hatching. P(20 < T-mean July < 25) is 0.4053 at Cernavoda, so the mean temperature in July is less suitable for the eggs hatching. Taking into account that the maximum temperatures in July are greater than 25 °C (as mentioned above), so optimal for reproduction, it results that there are climatic conditions for the Ae. albopictus establishment in the Dobrogea region.

3.4. Building a GIS-Based Model for Analyzing the Spatial Correlation between the Climatic Parameters Identified

Based on the proposed methodology (Figure 3), a spatial map (raster file) has been built for each climatic data set presented in

Table 4. Scores assigned to the climatic parameters.

No.	Criteria	Range	Score	Weight
1	tm (°C)	10.5 to 10.8	6	20%
		10.8 to 11.1	5
		11.1 to 11.4	4
		11.4 to 11.7	3
		11.7 to 12.0	2
		12.0 to 2.3	1
2	TJan (°C)	−4.1 to −3.6	6	20%
		−3.6 to −3.1	5
		−3.1 to −2.5	4
		−2.5 to 2.0	3
		−2.0 to −1.5	2
		−1.5 to −1.0	1
3	TJul (°C)	26.6 to 27.2	1	20%
		27.2 to 27.8	2
		27.8 to 28.4	3
		28.4 to 29.0	4
		29.0 to 29.6	5
		29.6 to 30.2	6
4	Rm (mm)	519 to 568	1	20%
		478 to 519	2
		427 to 478	3
		361 to 427	4
		283 to 361	5
		205 to 283	6
5	RH (%)	72.7 to 74.6	6	20%
		74.6 to 76.5	5
		76.5 to 78.4	4
		78.4 to 80.4	3
		80.4 to 82.3	2
		82.3 to 84.2	1

. An integrated analysis was carried out to create a suitable GIS model for identifying the potential locations for Ae. albopictus establishment. In this respect, all spatial maps were combined. Since the five criteria have different measurement units for each spatial data set, a common scale was defined, ranging from 1 to 6. The proposed scale is presented in Table 4, where the score “1” is assigned to the most favorable class, whereas “6” is assigned to the most disadvantageous one. For example, the values of Rm (mm) vary between 205 and 568. This interval was divided into six subintervals with the same lengths. According to the values presented in Table 2, Table 4, and Figure 4, a score of “1” was attached to the subinterval 519–568 (Table 4), whereas and score of “6” was attached to the more unfavorable class (205–283).

Figure 7 shows that most favorable classes for precipitation cover all Dobrogea, but the Danube Delta (north–east of the region). This is sustained by the rejection of the hypothesis that the annual mean precipitation is less than 500 mm at Sulina and Jurilovca.

After been reclassified to a common measurement scale, the datasets were finally combined using the Weighted Overlay tool. Each input raster was weighted by its respective importance. All datasets were considered equally important, so they have an equal influence on the determination of favorable Ae. albopictus establishment (Table 4).

Figure 8 shows that the most suitable area for the Ae. albopictus establishment is on the Black Sea Littoral and the Danube Delta (the classes 1 and 2). The most unfavorable location is near Corugea, where the average temperature in January is −3.5 °C.

We should remark that the condition that the annual precipitation is greater than 500 is not fulfilled for Sulina, which belongs to the most suitable location based on the map’s interpretation. This is only an apparent discordance since the weight associated with the precipitation is only 20% in the map creation.

According to Skolka [52], Ae. albopictus has been identified in Constanta city (since 2016–2017), in Murfatlar, Cumpăna (approximately 10 km Constanta), the lower Danube, and Danube Delta. In 2016, some dengue fever cases were reported in Mangalia, but there is no indication if those cases were autochthonous or imported [52].

4. Conclusions

Ae. albopictus mosquitos play an important role in speeding infectious diseases. This species was introduced in Europe, and it has a strong adaptive capacity. In this context, the study’s main objective was to develop a tool that could map the suitable locations for Ae. albopictus establishment in Dobrogea, Romania, using the available datasets of the current climate with respect to the Multicriteria Evaluation method. We performed statistical analyses to test the climate variables’ suitability for the establishment of this species in Dobrogea. The statistical tests and the GIS map should be interpreted together.

The previous studies carried out so far in Romania [24,25,52] had the goal of proving the mosquitoes’ presence in the country and not investigating the favorable conditions for their establishment and proliferation. For this reason, we consider that the assumption that all climatic parameters have an equal influence on identifying the most suitable area is reasonable.

Our approach differs from those found in the literature. It is easy to use, can be integrated into open source GIS (such as QGIS), and utilized in monitoring the mosquitoes’ habitat for placing the surveillance tools. This approach will give better results when a higher number of stations and a higher number of record data will be used.

A sensitive point of this study refers to the climatic datasets utilized to evaluate the Ae. albopictus spread in Europe. The direct and immediate interest in this study was to observe if Romania is situated in this area, from the climatic point of view. The objective was achieved by investigating the scientific literature related to (i) European area were Ae. albopictus was discovered and (ii) its climatic parameters (Table 2). However, utilizing a set of gridded data layers provided by different services like Worldclim (https://www.worldclim.org/), E-OBS (https://www.ecad.eu/download/ensembles/download.php), or CHELSA (https://chelsa-climate.org/) can significantly improve our findings (Figure 4 and Figure 5). The results obtained offer the possibility to choose the adequate climatic factors for the investigation performed for the Dobrogea region. Furthermore, other parameters could be introduced to improve the study accuracy. These will be done in the next step of our research.

To develop the model, we must better understand the correlation between the climatic factors, vegetation, soil type, water abundance, land use, etc., and the processes of Ae. albopictus establishment and development. More investigations are necessary to incorporate all information and reach a stage when the proposed approach’s findings could be compared with those of other used models. Still, we obtained encouraging results, and we consider this model a useful tool for a monitoring strategy.

To predict the Ae. albopictus establishment a combined approach of statistical, GIS analysis, and remote sensing imagery represents a task for all stakeholders involved.

According to climate change impact analyses [53], a continuous surveillance program is recommended, especially in areas with the highest suitability predicted by the model developed herein. Since the Ae. albopictus’ presence in the Dobrogea region was noticed, the development of a Protocol for Surveillance and Control of arboviruses transmitted by mosquitos in Romania is necessary.