Empirical dynamic modelling and enhanced causal analysis of short-length Culex abundance timeseries with vector correlation metrics.

Employing Empirical Dynamic Modelling we investigate whether model free methods could be applied in the study of Culex mosquitoes in Northern Greece. Applying Simplex Projection and S-Map algorithms on yearly timeseries of maximum abundances from 2011 to 2020 we successfully predict the decreasing trend in the maximum number of mosquitoes which was observed in the rural area of Thessaloniki during 2021. Leveraging the use of vector correlation metrics we were able to deduce the main environmental factors driving mosquito abundance such as temperature, rain and wind during 2012 and study the causal interaction between neighbouring populations in the industrial area of Thessaloniki between 2019 and 2020. In all three cases a chaotic and non-linear behaviour of the underlying system was observed. Given the health risk associated with the presence of mosquitoes as vectors of viral diseases these results hint to the usefulness of EDM methods in entomological studies as guides for the construction of more accurate and realistic mechanistic models which are indispensable to public health authorities for the design of targeted control strategies and health prevention measures.

Effects of climatic and environmental factors on mosquito population inferred from West Nile virus surveillance in Greece.

Mosquito-borne diseases' impact on human health is among the most prominent of all communicable diseases. With limited pool of tools to contrast these diseases, public health focus remains preventing mosquito-human contacts. Applying a hierarchical spatio-temporal Bayesian model on West Nile virus (WNV) surveillance data from Greece, we aimed to investigate the impact of climatic and environmental factors on Culex mosquitoes' population. Our spatio-temporal analysis confirmed climatic factors as major drivers of WNV-transmitting-Culex mosquitoes population dynamics, with temperature and long periods of moderate-to-warm climate having the strongest positive effect on mosquito abundance. Conversely, rainfall, high humidity, and wind showed a negative impact. The results suggest the presence of statistically significant differences in the effect of regional and seasonal characteristics, highlighting the complex interplay between climatic, geographical and environmental factors in the dynamics of mosquito populations. This study may represent a relevant tool to inform public health policymakers in planning preventive measures.

An epidemiological model for mosquito host selection and temperature-dependent transmission of West Nile virus.

We extend a previously developed epidemiological model for West Nile virus (WNV) infection in humans in Greece, employing laboratory-confirmed WNV cases and mosquito-specific characteristics of transmission, such as host selection and temperature-dependent transmission of the virus. Host selection was defined by bird host selection and human host selection, the latter accounting only for the fraction of humans that develop symptoms after the virus is acquired. To model the role of temperature on virus transmission, we considered five temperature intervals (≤ 19.25 °C; > 19.25 and < 21.75 °C; ≥ 21.75 and < 24.25 °C; ≥ 24.25 and < 26.75 °C; and > 26.75 °C). The capacity of the new model to fit human cases and the week of first case occurrence was compared with the original model and showed improved performance. The model was also used to infer further quantities of interest, such as the force of infection for different temperatures as well as mosquito and bird abundances. Our results indicate that the inclusion of mosquito-specific characteristics in epidemiological models of mosquito-borne diseases leads to improved modelling capacity.

Technical note: AQMEII4 Activity 1: evaluation of wet and dry deposition schemes as an integral part of regional-scale air quality models.

We present in this technical note the research protocol for phase 4 of the Air Quality Model Evaluation International Initiative (AQMEII4). This research initiative is divided into two activities, collectively having three goals: (i) to define the current state of the science with respect to representations of wet and especially dry deposition in regional models, (ii) to quantify the extent to which different dry deposition parameterizations influence retrospective air pollutant concentration and flux predictions, and (iii) to identify, through the use of a common set of detailed diagnostics, sensitivity simulations, model evaluation, and reduction of input uncertainty, the specific causes for the current range of these predictions. Activity 1 is dedicated to the diagnostic evaluation of wet and dry deposition processes in regional air quality models (described in this paper), and Activity 2 to the evaluation of dry deposition point models against ozone flux measurements at multiple towers with multiyear observations (to be described in future submissions as part of the special issue on AQMEII4). The scope of this paper is to present the scientific protocols for Activity 1, as well as to summarize the technical information associated with the different dry deposition approaches used by the participating research groups of AQMEII4. In addition to describing all common aspects and data used for this multi-model evaluation activity, most importantly, we present the strategy devised to allow a common process-level comparison of dry deposition obtained from models using sometimes very different dry deposition schemes. The strategy is based on adding detailed diagnostics to the algorithms used in the dry deposition modules of existing regional air quality models, in particular archiving diagnostics specific to land use-land cover (LULC) and creating standardized LULC categories to facilitate cross-comparison of LULC-specific dry deposition parameters and processes, as well as archiving effective conductance and effective flux as means for comparing the relative influence of different pathways towards the net or total dry deposition. This new approach, along with an analysis of precipitation and wet deposition fields, will provide an unprecedented process-oriented comparison of deposition in regional air quality models. Examples of how specific dry deposition schemes used in participating models have been reduced to the common set of comparable diagnostics defined for AQMEII4 are also presented.

A climate-dependent spatial epidemiological model for the transmission risk of West Nile virus at local scale.

In this study, initial elements of a modelling framework aimed to become a spatial forecasting model for the transmission risk of West Nile virus (WNV) are presented. The model describes the dynamics of a WNV epidemic in population health states of mosquitoes, birds and humans and was applied to the case of Greece for the period 2010-2019. Calibration was performed with the available epidemiological data from the Hellenic Centre for Disease Control and Prevention and the environmental data from the European Union's earth observation program, Copernicus. Numerical results of the model for each municipality were evaluated against observations. Specifically, the occurrence of WNV, the number of infected humans and the week of incidence predicted from the model were compared to the corresponding numbers from observations. The results suggest that dynamic downscaling of a climate-dependent epidemiological model is feasible down-to roughly 80km2. This below nomenclature of territorial units for statistics (NUTS) 3 level represents the municipalities being the lowest level of administrative units, able to cope with WNV and take actions. The average detection probability in hindcast mode was 72%, improving strongly as the number of infected humans increased. Using the developed model, we were also able to show the fundamental importance of the May temperatures in shaping the WNV dynamics. The modeling framework couples epidemiological and environmental dynamical variables with surveillance data producing risk maps downscaled at a local level. The approach can be expanded to provide targeted early warning probabilistic forecasts that can be used to inform public health policy decision making.

On the Transmission Dynamics of SARS-CoV-2 in a Temperate Climate.

An epidemiological model, which describes the transmission dynamics of SARS-CoV-2 under specific consideration of the incubation period including the population with subclinical infections and being infective is presented. The COVID-19 epidemic in Greece was explored through a Monte Carlo uncertainty analysis framework, and the optimal values for the parameters that determined the transmission dynamics could be obtained before, during, and after the interventions to control the epidemic. The dynamic change of the fraction of asymptomatic individuals was shown. The analysis of the modelling results at the intra-annual climatic scale allowed for in depth investigation of the transmission dynamics of SARS-CoV-2 and the significance and relative importance of the model parameters. Moreover, the analysis at this scale incorporated the exploration of the forecast horizon and its variability. Three discrete peaks were found in the transmission rates throughout the investigated period (15 February-15 December 2020). Two of them corresponded to the timing of the spring and autumn epidemic waves while the third one occurred in mid-summer, implying that relaxation of social distancing and increased mobility may have a strong effect on rekindling the epidemic dynamics offsetting positive effects from factors such as decreased household crowding and increased environmental ultraviolet radiation. In addition, the epidemiological state was found to constitute a significant indicator of the forecast reliability horizon, spanning from as low as few days to more than four weeks. Embedding the model in an ensemble framework may extend the predictability horizon. Therefore, it may contribute to the accuracy of health risk assessment and inform public health decision making of more efficient control measures.

Assessment of West nile virus transmission risk from a weather-dependent epidemiological model and a global sensitivity analysis framework.

West Nile virus (WNV) transmission risk is strongly related to weather conditions due to the sensitivity of the mosquitoes to climatic factors. We assess the WNV transmission risk of humans to seasonal weather conditions and the relative effects of parameters affecting the transmission dynamics. The assessment involves a known epidemiological model we extend to account for temperature and precipitation and a global uncertainty and sensitivity analysis framework. We focus on three relevant quantities, the basic reproduction number (R0), the minimum infection rate (MIR), and the number of infected individuals. The highest-priority weather-related WNV transmission risks can be attributed to the birth and death rate of mosquitoes, the biting rate of mosquitoes to birds, and the probability of transmission from birds to mosquitoes. Global sensitivity analysis indicates that these parameters make up a big part of the explained variance in R0 and MIR. The analysis allows for a dynamic assessment over time capturing the period parameters are more relevant than others. Global uncertainty and sensitivity analysis of WNV transmission risk to humans enable insights into the relative importance of individual parameters of the transmission cycle of the virus facilitating the understanding of the dynamics and the implementation of tailored control strategies.

Assessment and economic valuation of air pollution impacts on human health over Europe and the United States as calculated by a multi-model ensemble in the framework of AQMEII3.

The impact of air pollution on human health and the associated external costs in Europe and the United States (US) for the year 2010 are modeled by a multi-model ensemble of regional models in the frame of the third phase of the Air Quality Modelling Evaluation International Initiative (AQMEII3). The modeled surface concentrations of O3, CO, SO2 and PM2.5 are used as input to the Economic Valuation of Air Pollution (EVA) system to calculate the resulting health impacts and the associated external costs from each individual model. Along with a base case simulation, additional runs were performed introducing 20 % anthropogenic emission reductions both globally and regionally in Europe, North America and east Asia, as defined by the second phase of the Task Force on Hemispheric Transport of Air Pollution (TF-HTAP2). Health impacts estimated by using concentration inputs from different chemistry-transport models (CTMs) to the EVA system can vary up to a factor of 3 in Europe (12 models) and the United States (3 models). In Europe, the multi-model mean total number of premature deaths (acute and chronic) is calculated to be 414 000, while in the US, it is estimated to be 160 000, in agreement with previous global and regional studies. The economic valuation of these health impacts is calculated to be EUR 300 billion and 145 billion in Europe and the US, respectively. A subset of models that produce the smallest error compared to the surface observations at each time step against an all-model mean ensemble results in increase of health impacts by up to 30 % in Europe, while in the US, the optimal ensemble mean led to a decrease in the calculated health impacts by ~ 11 %. A total of 54 000 and 27 500 premature deaths can be avoided by a 20 % reduction of global anthropogenic emissions in Europe and the US, respectively. A 20 % reduction of North American anthropogenic emissions avoids a total of ~ 1000 premature deaths in Europe and 25 000 total premature deaths in the US. A 20 % decrease of anthropogenic emissions within the European source region avoids a total of 47 000 premature deaths in Europe. Reducing the east Asian anthropogenic emissions by 20 % avoids ~ 2000 total premature deaths in the US. These results show that the domestic anthropogenic emissions make the largest impacts on premature deaths on a continental scale, while foreign sources make a minor contribution to adverse impacts of air pollution.

Two-scale multi-model ensemble: is a hybrid ensemble of opportunity telling us more?

In this study we introduce a hybrid ensemble consisting of air quality models operating at both the global and regional scale. The work is motivated by the fact that these different types of models treat specific portions of the atmospheric spectrum with different levels of detail, and it is hypothesized that their combination can generate an ensemble that performs better than mono-scale ensembles. A detailed analysis of the hybrid ensemble is carried out in the attempt to investigate this hypothesis and determine the real benefit it produces compared to ensembles constructed from only global-scale or only regional-scale models. The study utilizes 13 regional and 7 global models participating in the Hemispheric Transport of Air Pollutants phase 2 (HTAP2)-Air Quality Model Evaluation International Initiative phase 3 (AQMEII3) activity and focuses on surface ozone concentrations over Europe for the year 2010. Observations from 405 monitoring rural stations are used for the evaluation of the ensemble performance. The analysis first compares the modelled and measured power spectra of all models and then assesses the properties of the mono-scale ensembles, particularly their level of redundancy, in order to inform the process of constructing the hybrid ensemble. This study has been conducted in the attempt to identify that the improvements obtained by the hybrid ensemble relative to the mono-scale ensembles can be attributed to its hybrid nature. The improvements are visible in a slight increase of the diversity (4 % for the hourly time series, 10 % for the daily maximum time series) and a smaller improvement of the accuracy compared to diversity. Root mean square error (RMSE) improved by 13-16 % compared to G and by 2-3 % compared to R. Probability of detection (POD) and false-alarm rate (FAR) show a remarkable improvement, with a steep increase in the largest POD values and smallest values of FAR across the concentration ranges. The results show that the optimal set is constructed from an equal number of global and regional models at only 15 % of the stations. This implies that for the majority of the cases the regional-scale set of models governs the ensemble. However given the high degree of redundancy that characterizes the regional-scale models, no further improvement could be expected in the ensemble performance by adding yet more regional models to it. Therefore the improvement obtained with the hybrid set can confidently be attributed to the different nature of the global models. The study strongly reaffirms the importance of an in-depth inspection of any ensemble of opportunity in order to extract the maximum amount of information and to have full control over the data used in the construction of the ensemble.

Seasonal ozone vertical profiles over North America using the AQMEII3 group of air quality models: model inter-comparison and stratospheric intrusions.

This study evaluates simulated vertical ozone profiles produced in the framework of the third phase of the Air Quality Model Evaluation International Initiative (AQMEII3) against ozonesonde observations in North America for the year 2010. Four research groups from the United States (US) and Europe have provided modeled ozone vertical profiles to conduct this analysis. Because some of the modeling systems differ in their meteorological drivers, wind speed and temperature are also included in the analysis. In addition to the seasonal ozone profile evaluation for 2010, we also analyze chemically inert tracers designed to track the influence of lateral boundary conditions on simulated ozone profiles within the modeling domain. Finally, cases of stratospheric ozone intrusions during May-June 2010 are investigated by analyzing ozonesonde measurements and the corresponding model simulations at Intercontinental Chemical Transport Experiment Ozonesonde Network Study (IONS) experiment sites in the western United States. The evaluation of the seasonal ozone profiles reveals that, at a majority of the stations, ozone mixing ratios are underestimated in the 1-6 km range. The seasonal change noted in the errors follows the one seen in the variance of ozone mixing ratios, with the majority of the models exhibiting less variability than the observations. The analysis of chemically inert tracers highlights the importance of lateral boundary conditions up to 250 hPa for the lower-tropospheric ozone mixing ratios (0-2 km). Finally, for the stratospheric intrusions, the models are generally able to reproduce the location and timing of most intrusions but underestimate the magnitude of the maximum mixing ratios in the 2-6 km range and overestimate ozone up to the first kilometer possibly due to marine air influences that are not accurately described by the models. The choice of meteorological driver appears to be a greater predictor of model skill in this altitude range than the choice of air quality model.