Gary Craig

@gcraig1

DNA testing shows that a male genetic sample on the clothing of Tammy Jo Alexander does not match three men it was checked against.

Livingston County Sheriff's Investigator Brad Schneider said an FBI lab sent the results to the Sheriff's Office last month.

"We do have the DNA results back from the lab and we did not come up with a match," he said.

This was not a complete surprise. Law enforcement asked three men who had known Alexander to provide samples. There was not evidence linking them to the November 1979 homicide, and the three men each agreed to provide a DNA swab for comparison.

Police ID teen found dead in cornfield in '79

► START FROM THE BEGINNING: Finding Tammy Jo podcast

Schneider said tips still come in regularly about the fatal shooting of Alexander, a 16-year-old from Florida found dead in a Caledonia, Livingston County, cornfield on the chilly damp morning of Nov. 10, 1979.

What happened in the decades after is now widely known: The girl, shot twice, would remain unidentified for years, despite the untiring efforts by the Livingston County Sheriff's Office and, particularly, now-retired Sheriff John York, to give her a name and find her killer. An unusual convergence of events — involving Tammy Jo's closest teenage friend, a long-ago boyfriend, her sister, an online community of men and women dedicated to help locate missing people, and law enforcement in two states — led to her identification in January 2015.

Schneider said he regularly talks to investigators with the Hernando County, Florida, Sheriff's Office, which, along with the FBI, also has been part of the investigation.

"As a matter of fact, one (tip) was received this week," he said in a recent interview. "I got one ... last week. The FBI has received a couple as well.

" ... Within the last couple weeks we've received several tips that we're still working on," he said.

All tips are treated as possible leads, and researched until they either are targeted for more investigation or determined to be fruitless or unreliable.

Alexander was from a troubled family in Brooksville, Florida, which is in Hernando County. She occasionally ran away and hitchhiked, and it is possible that's how she traveled to New York.

Sister hoped 'Cali' escaped toxic life

Laurel Nowell, who once hitchhiked with Alexander, began looking for her high school friend several years ago and found no online trail for her. She also communicated with Alexander's high school boyfriend, Kevin Williams, who did not know what had happened to Alexander.

Nowell connected with Alexander's sister, Pamela Dyson, and they then approached the Hernando County Sheriff's Office, which could not locate a missing person report for Alexander, who was last seen in Florida in 1979. Dyson provided a DNA sample, and information about Alexander was posted on the Department of Justice's National Missing and Unidentified Persons System, or NamUs.

NamUs includes information about missing persons or bodies that are found and unidentified. An online community, called Websleuths, has worked for years to try to identify missing persons. When the information about Alexander went online, one of its members, Carl Koppelman of California, believed that Tammy's appearance matched that of the young woman found in the Caledonia cornfield, whose image had been circulated internationally.

He reached out to law enforcement in Florida. A Hernando County Sheriff's Office Detective, George Loydgren, was assigned cold cases and was the contact for the missing person case of Tammy Jo Alexander. Once he saw the photo of the girl found dead in New York, "I thought right away it was her," Loydgren said.

Dyson's DNA sample showed a familial match with the body found in Caledonia.

Dyson still holds out hope that the investigation into her sister's death will end with the identification of a killer.

"Somebody has to know something, in just my personal opinion," she said. "We haven't reached the right person that has the piece of information that will lead to this case being solved."

"I always have hope. And I will never give up till the day I die hoping that they find who killed her."

GCRAIG@Gannett.com

Want to know more?

You can listen and subscribe to the "Finding Tammy Jo" podcast at findingtammyjo.com. The podcast is co-hosted by Gary Craig of the Democrat and Chronicle. Other crime and "cold case" podcasts from the USA TODAY Network can be found at coldcasepodcasts.tumblr.com.


Contents

Spinelli was born in Rome, and joined the Italian Communist Party (PCI) at an early age in order to oppose the regime of Benito Mussolini's National Fascist Party. Following his entry into radical journalism, he was arrested in 1927 and spent ten years in prison and a further six in confinement. During the war he was interned on the island of Ventotene (in the Gulf of Gaeta) along with some eight hundred other political opponents of the regime.[1] During those years, he broke with the Italian Communist Party over Stalin's purges (resulting in him being ostracised by many of the other prisoners), but refused to compromise with the fascist regime, despite offers of early release.[citation needed]

In June 1941, well before the outcome of the war was safely predictable, Spinelli and fellow prisoner Ernesto Rossi completed the Ventotene Manifesto, eventually entitled Per un’Europa libera e unita ("For a Free and United Europe. A Draft Manifesto"), which argued that, if the fight against the fascist powers was successful, it would be in vain if it merely led to the re-establishment of the old European system of sovereign nation-states in shifting alliances. This would inevitably lead to war again. The document called for the establishment of a European federation by the democratic powers after the war. Because of a need for secrecy and a lack of proper materials at the time, the Manifesto was written on cigarette papers, concealed in the false bottom of a tin box and smuggled to the mainland by Ursula Hirschmann. It was then circulated through the Italian Resistance, and was later adopted as the programme of the Movimento Federalista Europeo, which Spinelli, Rossi and some 20 others established, as soon as they were able to leave their internment camp. The founding meeting was held in clandestinity in Milan on the 27/28 August 1943.

The Manifesto was widely circulated in other resistance movements towards the end of the war. Resistance leaders from several countries met clandestinely in Geneva in 1944, a meeting attended by Spinelli.

The Manifesto put forward proposals for creating a European federation of states, the primary aim of which was to tie European countries so closely together that they would no longer be able to go to war with one another. As in many European left-wing political circles, this sort of move towards federalist ideas was argued as a reaction to the destructive excesses of nationalism. The ideological underpinnings for a united Europe can thus be traced to hostility to nationalism: "If a post war order is established in which each State retains its complete national sovereignty, the basis for a Third World War would still exist even after the Nazi attempt to establish the domination of the German race in Europe has been frustrated" (founding meeting of the MFE).

After the war, Spinelli, leading the federalist MFE, played a vanguard role in the early episodes of European integration, criticising the small steps approach and the dominance of intergovernmentalism, feeling even that the chance to unite Europe had been missed as sovereign states were re-established without any common bond other than the functionalist OEEC and the largely symbolic Council of Europe. Even the European Coal and Steel Community (ECSC) was felt to be too sectoral. The MFE believed governments alone would never relinquish their national power without popular pressure. They advocated a European constituent assembly to draft a European Constitution.

This approach eventually had a response from governments when they set up the "ad hoc assembly" of 1952–3. It was Spinelli who persuaded Italian Prime Minister Alcide De Gasperi to insist in the negotiation of the European Defence Community (EDC) treaty on a provision for a parliamentary assembly to draw up plans for placing the EDC, the ECSC and any other development within a global constitutional framework to "replace the present provisional organization" with "a subsequent federal or confederal structure based on the principle of the separation of powers and having, in particular, a two-chamber system of representation". The Assembly was invited to submit its proposals within six months of its constitutive meeting following the entry into force of the EDC treaty. In fact, the Foreign Ministers, meeting three months after the signature of the EDC treaty, invited the ECSC Assembly immediately to draft a "treaty constituting a European Political Authority" without waiting for ratification of the EDC Treaty.

Spinelli played a significant role in advising the drafting of the Assembly's proposal for a European "Statute". However, the failure of France to ratify the EDC treaty meant it was all to no immediate avail. Some of its ideas, however, were taken up in subsequent events.

Following the crisis of the failure of the EDC and the "re-launch" under the Paul-Henri Spaak committee, which led to the 1958 EEC Treaty, Spinelli, recognising that the EEC institutions were the only real existing form of European integration, but still considering that they were insufficient and that they lacked a democratic legitimacy, embarked on a "long march through the institutions". In 1970, he was nominated by the Italian government to be a member of the European Commission[2] from 1970 to 1976, taking responsibility for industrial policy in order to develop European policies in a new field.

Spinelli decided to run in the first direct elections to the European Parliament in 1979. He did so as an independent candidate on the list of the Italian Communist Party, which by then had become a Eurocommunist party and was keen to have prominent independent figures to stand on its list of candidates. He was elected and used the position to urge the first elected parliament to use its democratic legitimacy to propose a radical reform of the European Community, to transform it into a democratic European state.

To this end, he began to gather like-minded Members of the European Parliament around him, taking care to involve Members from different political groups. An initial meeting at the "Crocodile" restaurant in Strasbourg set up the "Crocodile Club", which, once it was of sufficient size, tabled a motion for Parliament to set up a special committee (eventually established in January 1982 as the Committee on Institutional Affairs, with Spinelli as General Rapporteur) to draft a proposal for a new treaty on union.

The idea was that the European Parliament should act as a constituent assembly, although Spinelli was prepared to make compromises on the way to secure broad majorities behind the process. On 14 February 1984, the European Parliament adopted his report and approved the Draft Treaty Establishing the European Union. The decision was taken with 237 votes for and 31 against (43 abstentions).[3]

Spinelli's project was soon buried by the governments of the member states. However, it provided an impetus for the negotiations which led to Single European Act of 1986 and the Maastricht Treaty of 1992. This happened with the help of several National parliaments, which adopted resolutions approving the Draft Treaty, and of French President François Mitterrand who, following a meeting with Spinelli, came to the European Parliament to speak in favour of its approach, thereby reversing France's policy (since Charles De Gaulle) of hostility to anything but an intergovernmental approach to Europe. This momentum was enough to obtain the support of a majority of national governments to trigger the treaty revision procedure.

Although the resultant treaties fell short of what Spinelli would have liked, his efforts did trigger a new momentum in European integration, including a major increase in the powers of the European Parliament within the EU system. In honour of his work, the largest building of the European Parliament complex in Brussels was named after him.

^See Spinelli, Altiero: Come ho tentato di diventare saggio (1987). This autobiographical work was to be in two parts. The first part, which appeared during Spinelli's life-time under the title Io Ulisse (1984), covers the life of the author until the Liberation (1943). Spinelli's writing on the second part (La goccia e la roccia) was interrupted by his death in 1986.

^See Bieber, Jacqué and Weiler: An Ever Closer Union. A critical analysis of the Draft Treaty establishing the European Union. Luxembourg, 1985. This publication also contains the text of the treaty. More recent comment and analysis of the 1984 treaty and Spinelli's role is found in Burgess, Michael: Federalism and European Union: the Building of Europe, 1950–2000. Routledge 2000; Corbett, Richard: The European Parliament's Role in Greater EU Integration.(1998) Palgrave. 424 pages. ISBN 0-333-94938-2


[A case of mesoblastic nephroma in adulthood].
We report a rare adult case of mesoblastic nephroma. A right renal mass was detected in a 51-year-old woman by ultrasonography. The radiographic examination revealed that the mass was the hypovasucular tumor, but there was a possibility of a malignant tumor. We therefore operated on her for right radical nephrectomy. The tumor was 40 X 52 mm in size and had several cysts. Based on the histological analysis, we diagnosed it as congenital mesoblastic nephroma consisting of spindle cells.


Introduction {#Sec1}
============

Air quality is, and will likely continue to be, an important issue over the Eastern Mediterranean and the Middle East (EMME) region. Air pollution influences the quality of life, the number of premature deaths, and the frequency and severity of various respiratory diseases (Giannadaki et al. [@CR21]; Lelieveld et al. [@CR39]; Abdo et al. [@CR2]; Khader et al. [@CR34]; Dayan et al. [@CR10]). The EMME region is a receptor of pollution from multiple sources of anthropogenic and natural origin (Lelieveld et al. [@CR38]; Gerasopoulos et al. [@CR19]; Mihalopoulos et al. [@CR42]; Astitha and Kallos [@CR4]; Im and Kanakidou [@CR29]). The island of Cyprus, located centrally in the EMME region, is ideally located to assess current regulations and evaluate model-benchmarking methodologies. Cyprus is affected by long-range pollution transport from three continents (Europe, Africa, and Asia) and dust aerosols from the two largest desert regions in the world (North Africa and Middle East).

Station observations can only monitor the level of pollutants over specific locations, while models can provide insights into the emission, transport, and transformation of pollution over the whole region of interest. Therefore, models that exhibit satisfactory performance can complement observations or even replace them in cases of disoperation, for policy-related applications and reporting purposes. Thus, it is necessary to harmonize the criteria of model skill and capability to reproduce air quality features over a specific region, for officially reporting national air pollution levels and for examining compliance with regulations.

Model performance criteria are defined to benchmark a model application against agreed and/or regulated quality standards, by comparing statistical indicators against bound values. Specific validation protocols (i.e., Dennis et al. [@CR11] for assessment; Thunis and Clappier [@CR52] for emission scenarios) have been developed to support the use of air quality models when performing various tasks (assessment, forecasting, planning). A methodology for unified model evaluation process has been developed (Thunis et al. [@CR53], [@CR54]) in the framework of the Forum for Air Quality Modelling in Europe (FAIRMODE). The FAIRMODE model evaluation methodology aims at promoting and supporting the harmonized use of models by EU Member States, with emphasis on model application under the European Air Quality Directives. The approach is based on paired modeled and monitored data to offer diagnostics of model performance using various statistical indicators and diagrams. The FAIRMODE IDL-based DELTA software tool (<http://aqm.jrc.ec.europa.eu/DELTA/>) that incorporates this methodology can produce summary reports on performance indicators and guide model development to improve statistical metrics that are not in line with target values.

The FAIRMODE methodology has been used in the past to evaluate chemical transport models over Europe, by Carnevale et al. ([@CR6]) over the Po Valley and Georgieva et al. ([@CR17]) over Bulgaria. Recently, Monteiro et al. ([@CR43]) presented an analysis of the strengths and weaknesses of the FAIRMODE benchmarking approach, focusing on the pollutants regulated by the Air Quality Directive (PM2.5, NO~2~, and O~3~) and based on feedback from different research groups over Europe using continental, regional, and urban scale models. Through a Strengths, Weaknesses, Opportunities and Threats (SWOT) analysis, they identified the main advantages and value of the FAIRMODE approach compared to other methodologies. The main strengths are the successful promotion of harmonized reporting for air quality model applications to meet the relevant directives and the integration of the most essential (statistically and representatively) quality indicators. Weaknesses include among others the lack of a clear definition and use of measurement uncertainty for various pollutants, a significant component of the methodology. Evaluation and SWOT studies are essential and provide motivation and technical basis for a continuous improvement of the methodology. Evaluation and SWOT studies are essential and provide for motivation and technical basis for a continuous improvement of the methodology.

In this work, we apply the FAIRMODE methodology to evaluate the skill of the Weather Research and Forecasting model coupled with Chemistry (WRF-Chem) to simulate air quality for the year 2015. We have selected an integrated modeling system since the impact of aerosols on radiation, clouds, and precipitation remains an open research question, especially in our region of interest. Aerosols greatly influence the atmospheric processes over the EMME due to their abundance in the region from both natural (sea salt, mineral dust) and anthropogenic sources (Real and Sartelet [@CR47]; Gerasopoulos et al. [@CR20]; Kushta et al. [@CR35]; Mailler et al. [@CR41]; Abdelkader et al. [@CR1]; Gkikas et al. [@CR22]). Hence, models that incorporate the links and feedbacks between airborne pollution and atmospheric processes are needed (Grell and Baklanov [@CR23]; Baklanov [@CR5]).

The evaluation of WRF-Chem is based on the comparison of model results with available observations from the national air quality monitoring network. Ambient air quality is regulated in the EU member countries by the Directive 2008/50/EC. We model regulated pollutants, namely ozone (maximum daily 8-h mean O~3~), nitrogen dioxide (hourly NO~2~), and fine particulate matter (daily mean PM2.5 values), and discuss the use of model results in the context of the application of the EU Air Quality Directive for the country of Cyprus. We focus on the island of Cyprus to highlight strengths and weaknesses when applying such a methodology in countries of limited size with no land borders, subject to long-range transport of pollution from multiple sources, and with specific constrains in monitoring sites. Air quality over Cyprus is monitored through a network of 12 stations including five rural stations (four rural background and one industrial stations) and seven urban stations (four traffic and three urban background stations).

The paper is organized as follows: "Model and data" includes a description of the model data and observations used in this study, as well as the statistical background of the FAIRMODE evaluation methodology. In "[Results and discussion](#Sec8){ref-type="sec"}," we present the results of the model and we compare with observation for O~3~, NO~2~, and PM2.5. In "[Concluding remarks](#Sec12){ref-type="sec"}," we give our concluding remarks on the capabilities of the evaluation methodology and its use as a policy tool for air quality and pollutant-level exceedance reporting in European countries.

Model and data {#Sec2}
==============

Model description and input data {#Sec3}
--------------------------------

For this study, we use the coupled meteorological and atmospheric chemistry model WRF-Chem (v3.9.1.1) that simultaneously simulates physical and chemical processes taking into account direct and indirect feedbacks (Grell et al. [@CR25]; Fast et al. [@CR15]). The model is applied for a year-long period (2015) using two domains with respective grid spacing of 50 and 10 km with the nested one, covering the EMME countries, centered over Cyprus (Fig. [1](#Fig1){ref-type="fig"}a). The model physical and chemical parameterizations are summarized in [Table](#MOESM1){ref-type="media"} [1](#Tab1){ref-type="table"}. The parent domain has a large extent to cover important emission sources inside the regional model instead of representing them as lateral forcing from the global modeling system. WRF-Chem has been applied and evaluated for a summer period (July 2014), over Cyprus, in two recent studies by Kushta et al. ([@CR36]) and Georgiou et al. ([@CR18]). The model achieves an overall good representation of the unique geophysical features and atmospheric chemical composition of the region, considering the limitations imposed by uncertainties in the emission inventories and physical and chemical model parameterizations it is based on.Fig. 1**a** The domains covered by this study: the parent domain covers Europe, North Africa, and the Arabian Peninsula at a horizontal resolution of 50 km, and the nested domain of the EMME has a grid spacing of 10 km. **b** The location and type of the air pollution monitoring stations over Cyprus. Rural background stations (3) are shown in purple squares, urban (background and traffic, 7 in total) in red, the Zygi industrial station as black star, and the Troodos rural background station, located above PBL, as blue circleTable 1Physical and chemical configuration of the modeling applicationProcessScheme/parameterizationReferenceMicrophysicsMorrison 2-moment schemeMorrison et al. [@CR44]Land surfaceNOAH land surface modelChen and Dudhia [@CR8]Boundary layerYonsei University (YSU) planetary boundary layerHong et al. [@CR27]CumulusGrell 3D ensemble schemeGrell and Devenyi [@CR24]Surface layerMM5 similarity surface layer schemeZhang and Anthes [@CR60]RadiationRapid Radiative Transfer Model (RRTTM)Iacono et al. [@CR28]Gas-phase chemistryRACM regional atmospheric chemistry mechanismStockwell et al. [@CR51]AerosolsModal Aerosol Dynamics Model for Europe (MADE), Secondary Organic Aerosol Model (SORGAM)Ackermann et al. [@CR3], Schell et al. [@CR48]

The WRF-Chem model meteorology is driven by 6-hourly boundary conditions at 0.5° × 0.5° horizontal resolution from the US National Centers for Environmental Prediction (NCEP) global forecast system (GFS). The modeled chemistry uses boundary conditions from the global Model for Ozone And Related chemical Tracers (MOZART version 4) for the outermost domain (Emmons et al. [@CR13]). The nested domain receives its meteorological and chemical lateral conditions from the coarse domain. The emissions of natural aerosol species (mineral dust, sea salt) are simulated online using the parameterization schemes included in the model. Biogenic emissions of volatile organic species are derived from the Model on Emissions of Gases and Aerosols from Nature (MEGAN), as described in Guenther et al. ([@CR26]).

The anthropogenic emissions are based on the EDGAR-HTAP v2 emission estimates for 2010 at a resolution of 0.1° × 0.1°, with no further modifications to incorporate changes that may have occurred in the recent years in the region (Janssens-Maenhout et al. [@CR30]). The emission inventory is processed for spatial and temporal allocation to the resolution of the model domains and lumped speciation to match the species of the chemical mechanisms used (Kushta et al. [@CR36]).

Observations {#Sec4}
------------

The model performance is evaluated against hourly measurements of the national monitoring network for the year 2015 over Cyprus. The Cypriot Legal Framework on Ambient Air Quality comprises the Quality of Ambient Air Law of 2010 (Ν. 77(Ι)/2010) and two series of Regulations that determine limits for the concentrations of certain pollutants in ambient air, in line with the European regulations on ambient air quality. The implementation and management of the law and regulations are the responsibility of the Ministry of Labour and Social Insurance Department of Labour Inspection (DLI). The DLI air quality monitoring network consists of four rural background stations (one of which located above the atmospheric boundary layer in the Troodos mountain), seven urban stations comprising of three urban background (hereafter referred to as urban) and four urban traffic stations (traffic), and a station located in the vicinity of the coastal, heavy-industry zone of Zygi (industrial). Station classification information, following the Exchange of Information Decision (EOI 97/101/EC [@CR9]) and the Implementation Decision for Reporting (2011/850/EU [@CR12]), and the availability of measurements per station and per species are summarized in Table [2](#Tab2){ref-type="table"}. We focus on the atmospheric pollutants included in the air quality directives and supported by the evaluation software: hourly values regarding NO~2~, mean daily values for PM2.5, and daily maximum 8-h O~3~. The analysis is performed in volume mixing ratio units (ppbV) for gaseous pollutants and in mass concentration units (μg m^−3^) for particulate matter that are directly comparable as they are the original units for both observations and model results.Table 2The air quality monitoring network of Cyprus (station, type, availability of observations for the respective pollutant)StationType of areaType of stationNONO~2~SO~2~O~3~COPM10PM2.5BenzeneAgia MarinaRuralBackgroundXXXXXXXXCavo GrecoRuralBackgroundXXXIniaRuralBackgroundXXXTroodosRuralBackgroundXXXLarnaca ResUrbanBackgroundXXXXXLarnaca TraUrbanTrafficXXXXXXXXLimassol ResUrbanBackgroundXXXXXLimassol TraUrbanTrafficXXXXXXXNicosia ResUrbanBackgroundXXXXXXNicosia TraUrbanTrafficXXXXXXXXPafosUrbanTrafficXXXXXXXZygiRuralIndustrialXXXXXX

As shown in Fig. [1](#Fig1){ref-type="fig"}b, the rural background stations are located along the central west-east axis of the island monitoring air pollution transport over the country under the main atmospheric circulation regimes (westerlies and easterlies). The coastal stations capture the impact of local emissions and/or long-range transport of pollutants depending on the dominating wind direction (Tyrlis and Lelieveld [@CR57]; Tyrlis et al. [@CR58]). The rural background stations are all used in the analysis regarding ozone and nitrogen dioxide. The urban background and traffic stations are located within the major urban areas (Nicosia, Limassol, Larnaca, and Pafos). For our study, we only use the urban background stations for the gaseous pollutants since the traffic stations capture high spikes in local emissions (especially NOx) due to their proximity to main traffic routes that cannot be resolved at the horizontal resolution of a regional-scale air quality model (10--50 km). Regarding fine particulate matter, only one rural background station has available PM2.5 observations (Agia Marina). There are no other stations located in rural areas that can contribute to the present PM2.5 analysis. The urban traffic stations provide additional measurements of particulate matter; thus, we include these stations in our analysis of aerosol pollutants, to investigate the performance criteria stringency, even though the methodology is not designed for application in traffic stations.

Benchmarking evaluation methodology {#Sec5}
-----------------------------------

### Modeling quality indicator {#Sec6}

The overall benchmarking evaluation procedure used in this study is presented in detail in Thunis et al. ([@CR53], [@CR54]). The methodology works on paired series of model and observed values of the respective pollutant for year-long periods. It investigates the model capabilities by introducing an overall indicator, namely, the modeling quality indicator (MQI), taking into account the measurement uncertainty of each pollutant. The measurement uncertainty parameters currently used in the FAIRMODE methodology are discussed in detail in Thunis et al. ([@CR55]) for ozone and in Pernigotti et al. ([@CR46]) for coarse particulate matter and nitrogen dioxide.

MQI defines the deviation between measured (O) and modeled (M) values at a given *i* time (hour or day) as a factor of the measurement uncertainty and a scaling factor that indicates the stringency of the objective to be satisfied (e.g., a scaling factor of 2 means the allowed deviation should be within a factor of 2 of the measurement uncertainty gap):$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym}
                \usepackage{amsfonts}
                \usepackage{amssymb}
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$ MQI=\frac{\mid Oi- Mi\mid }{\beta\ {U}_{95}\ (Oi)},\beta =2 $$\end{document}$$

*O*~i~ and *M*~i~ are the observed and modeled values, respectively, *U*~95~ is the 95th percentile measurement uncertainty of the observed concentration level, and *β* is the coefficient that scales the proportionality of the bias to the measurement uncertainty (Thunis et al. [@CR54]). MQI can be generalized to a yearly time series as follows:$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym}
                \usepackage{amsfonts}
                \usepackage{amssymb}
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$ MQI=\frac{RMSE}{\beta\ RMSu} $$\end{document}$$*where RMSu* is the root mean square of the measurement uncertainty (details in Thunis et al. [@CR55]; Pernigotti et al. [@CR46]) and RMSE is the root mean square error between model and measured concentrations. For annual values, *O*~i~ and *M*~i~ are substituted in the formula by the observed and modeled annual mean, respectively, and the uncertainty is calculated at the 95th percentile upon the annual mean measured concentration. The modeling quality objective (MQO) is the criterion for the value of the MQI to be satisfied for satisfactory model performance in terms of air quality representation for reporting applications (minimum level of quality). MQO is fulfilled if the MQI is less than or equal to unity for at least 90% of the available monitoring data.

The measurement uncertainty parameter is a significant component of the evaluation methodology reflecting our confidence level in the observational data. Initially, the measurement uncertainty of the benchmarking methodology was set to be constant (independent of the concentration level of the respective pollutant) and defined by the data quality objective (DQO) value of the Air Quality Directive to 15, 15, and 25% for O~3~, NO~2~, and PM10, respectively (Thunis et al. [@CR53], [@CR54]). Thunis et al. ([@CR55]) reshaped the measurement uncertainty as a combination of two components, one given as a function of the pollutant concentration and a non-proportional component (independent of measured values) calculated on a concentration level of choice (reference value, RV). The total measurement uncertainty represents the 95th percentile highest value among all uncertainty values calculated for each pollutant. The calculations are performed on data from JRC instrument inter-comparison results (Lagler et al. [@CR37]) for particulate matter, EU AIRBASE stations for series of meteorological years for NO~2~, and analytical relationships for O~3~. Several updates of the uncertainty parameters used to estimate the measurement uncertainty for particulate matter have been made (see "Guidance document on modeling quality objectives and benchmarking" (FAIRMODE [@CR14]) with current values reflecting uncertainties associated to different types of instrumentations (e.g., β-ray measurement technique).

For the visualization of the MQO, a target diagram adapted from Jolliff et al. ([@CR31]) is used (results of current study shown in Fig. [3](#Fig3){ref-type="fig"}, top plots). The horizontal axis represents the central root mean square error (CRMSE) while the vertical axis refers to BIAS, normalized by measurement uncertainty of the respective pollutant levels. For each station, a dot is placed on the target diagram and the distance of the dot from the center of the diagram represents the MQI for that station. Stations with MQI within the green area are identified as stations that fulfill the performance criteria. The dashed line represents the limit outside which (but still within the green area) model results are within the measurement uncertainty range. When the MQI of a specific station is greater than one and falls out of the colored area, there are statistically significant differences between model and observations. Four zones on the plot help the user identify the reasons for model-observation differences in terms of standard deviation (SD), bias, and correlation (*R*) (Thunis et al. [@CR53], [@CR54]; Jolliff et al. [@CR31]). The top (bottom) two zones exhibit positive (negative) bias while the assignment of the station to the left (right) zone of the graph indicates a unsystematic (systematic) RMSE ratio (error dominated by SD or by *R*). Model performance indicators through the "target" approach have been assessed in Pederzoli et al. ([@CR45]).

### Model performance indicators {#Sec7}

The MQI proposed in the framework of the FAIRMODE benchmarking methodology combines the bias error (BIAS), standard deviation (σ), and correlation (R) into a single number, providing a general measure of the model performance associated with RMSE. These three statistical metrics are related as follows:$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym}
                \usepackage{amsfonts}
                \usepackage{amssymb}
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$ {MQI}^2=\frac{RMSE^2}{{\left(\beta\ RMSu\right)}^2}=\frac{BIAS^2}{{\left(\beta\ RMSu\right)}^2}+\frac{{\left({\sigma}_M-{\sigma}_{\mathrm{o}}\right)}^2}{{\left(\beta\ RMSu\right)}^2}+\frac{2{\sigma}_{\mathrm{o}}{\sigma}_{\mathrm{M}}\left(1-R\right)}{{\left(\beta\ RMSu\right)}^2} $$\end{document}$$

Assuming ideal cases where two out of the three statistical indicators perform perfectly, we get model performance indicators (MPI) that depend solely on the third remaining metric. MPIs associated separately with correlation, standard deviation, and bias can then be used to highlight which aspect of the model performance can be improved. Model performance criteria (MPC) are defined for each MPI to be fulfilled during model evaluation.

That is, assuming bias to be zero and the standard deviation of observations equal to that of modeled values leads to the following condition on *R*:$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym}
                \usepackage{amsfonts}
                \usepackage{amssymb}
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$ R>1-0.5\ {\beta}^2\frac{RMSu^2}{\sigma_{\mathrm{o}}{\sigma}_{\mathrm{M}}} $$\end{document}$$where *σ*~o~ and *σ*~M~ are the standard deviation of the observed and modeled concentrations, respectively. Similarly, the MPIs for bias and standard deviation translate in the following conditions:$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym}
                \usepackage{amsfonts}
                \usepackage{amssymb}
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$ \kern2.75em \left| BIAS\right|\le \beta RMSu $$\end{document}$$$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym}
                \usepackage{amsfonts}
                \usepackage{amssymb}
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$ \kern2.25em \left|{\sigma}_{\mathrm{M}}-{\sigma}_{\mathrm{o}}\right|\le \beta RMSu $$\end{document}$$

It is important to note that the model performance criteria for bias, correlation, and standard deviation represent necessary but not sufficient conditions to ensure fulfillment of the MQO.

Results and discussion {#Sec8}
======================

We focus on the pollutant metrics as required by policy applications: the maximum daily 8-hourly mean O~3~, hourly NO~2~, and mean daily PM2.5 concentrations. The tool additionally enables the investigation of other metrics (hourly concentrations for all pollutants). We have not performed an analysis on the coarse particulate matter (PM10) in the current study. The fraction of particles with diameter between 2.5 and 10 μm in the region is dominated by natural aerosols (mineral dust) that can introduce large uncertainties to the evaluation procedure without a separate assessment. As shown in Karanasiou et al. ([@CR33]) and Chalbot et al. ([@CR7]), crustal soil particles typically concentrate into the coarse fraction mode of particulate matter and only a small fraction of them (5--10%) are present in the fine mode (diameter less than 2.5 μm).

All the stations included in the analysis satisfy the 75% data availability threshold, as required by the methodology. The EU Air Quality Policy introduced as Air Quality Standards by the Directive 2008/50/EU for the Member States defines a maximum threshold of 25 (current) and 20 μg m^−3^ (future) for mean annual PM2.5 concentration, 120 μg m^−3^ for maximum daily 8-h mean O~3~ not to be exceeded over 25 days (averaged over a 3-year period), and 40 (200) μg m^−3^ for annual (hourly) NO~2~ concentrations. The target plot analysis and threshold exceedance assessment of the gaseous pollutants are performed separately on the urban and rural background stations as they exhibit dissimilar characteristics. The main statistical indicators (bar plots of mean, standard deviation, bias, etc.) include both categories of stations (rural and urban) for an overview of the model performance.

Ozone {#Sec9}
-----

The bar plots of mean annual values of the daily max 8-h O~3~ concentrations (Online Resource, Fig. [S1](#MOESM1){ref-type="media"}a) indicate that the model is capable of reproducing the distribution of ozone over the island. An overestimation is noticeable in urban areas. This overestimation may be related to the lateral forcing from the global model and low local emission fluxes of ozone precursors such as nitrogen oxides, as shown in Kushta et al. ([@CR36]). There is a difference of daily max 8-h O~3~ values between rural and urban background stations that leads to lower annual mean max 8-h concentrations in the urban sites as a result of ozone titration by nitric oxide from traffic and residential activities. This outcome highlights a local-scale intense conversion effect (of O~3~ by NO into NO~2~) near emission sources where higher levels of NOx emissions are present relative to rural areas. The standard deviation ranges from 6.5 to 15.2 ppbV, with areas under "clean air" influence, like the Inia station located in the westernmost part of the island (affected by the prevailing north--northwest wind flow) and the Troodos station (mountainous site above boundary layer) exhibiting the largest variations in max 8-h O~3~ levels (Online Resource, Fig. [S1](#MOESM1){ref-type="media"}b). The model shows a smoother variation in terms of standard deviation of about 9 ± 2 ppbV. The mean bias and correlation coefficients are larger over the urban stations (Online Resource, Figs. [S1](#MOESM1){ref-type="media"}c and [S1](#MOESM1){ref-type="media"}d). While the model cannot reproduce the magnitude of ozone in urban background sites (large mean bias), the variation correlates better due to the more pronounced photochemical processes. In rural background stations, the mean bias is lower (modeled values are closer to the observed ones) but the correlations lack accuracy mostly due to the influence of long-range transported pollution with no distinguishable temporal variation in the production and depletion life cycle. The fact that the Troodos station (not influenced by the boundary layer) exhibits smaller mean bias despite its large observed standard deviation underlines the significance of accurately representing the boundary layer in the simulation of pollutant advection and diffusion.

We next examine the compliance with the model performance objectives as defined by the evaluation methodology. Figure [2](#Fig2){ref-type="fig"} (top plots) presents the assessment target plot of the MQI for rural background (Fig. [2](#Fig2){ref-type="fig"}a), urban background stations (Fig. [2](#Fig2){ref-type="fig"}b), and all stations combined (Fig. [2](#Fig2){ref-type="fig"}c). The model quality objective requires a MQI of less than unity for at least 90% of the stations included in the analysis. As depicted on the top panel plots, the MQI for both group of stations (four background and three residential stations) is satisfied with MQI~rural~ = 0.387 and MQI~urban~ = 0.391 for hourly values and MQI~rural~ = 0.231 and MQI~urban~ = 0.519 for annual means under a measurement uncertainty of 18%. For all stations combined, the MQI~all~ is 0.4 for the hourly assessment and 0.5 for the annual mean values, again satisfying the MQO of less than one (Fig. [2](#Fig2){ref-type="fig"}c). A small positive bias of the model is evident due to the location of the station dots in the upper part of the target diagram. The statistical metric of correlation is the one with less satisfactory performance.Fig. 2Assessment target plots (top row) for **a** rural background, **b** urban background, and **c** all stations. **d** Model performance criteria (bottom row) associated with correlation coefficient and **e** scatter plot of mean annual values of 8-h max ozone

The model performance criteria related to the MPI for correlation is given in Fig. [2](#Fig2){ref-type="fig"}d for all types of stations (rural and urban background sites). All stations satisfy the MPI for correlation (as well as the MPI for standard deviation and mean bias, not shown). There is a clear distinction between observed mean annual values for averaged max 8-h O~3~ concentrations in rural versus urban background stations (Fig. [2](#Fig2){ref-type="fig"}e), but they remain close in their respective modeled values (\~ 50 ppbV).

Another component of the benchmarking methodology is the assessment of threshold exceedances and the provision of a summary report that includes information on the model performance (temporal and spatial statistical indicators) and exceedances for all stations in one diagram. As seen in Fig. [3](#Fig3){ref-type="fig"}, for Cyprus in 2015, the mean max 8-h O~3~ over the rural background stations is approximately 50--53 ppbV while in the urban background stations, it is almost constant in all cities at about 45 ppbV (row 1 of summary statistics plot in Fig. [3](#Fig3){ref-type="fig"}). Row 2 shows the number of days with exceedances of the threshold set by the user. We performed a threshold analysis for two values, 40 ppbV and 60 ppbV, to account for the impact of ozone on both vegetation and human health, respectively. The value of 40 ppbV has been proposed in the second edition of the World Health Organization Air Quality Guidelines for Europe (WHO [@CR59]). The assessment of threshold exceedances over urban areas is important for human health; the analysis of the respective situation in vegetated areas is of equal importance for ecosystems.Fig. 3Summary report of threshold exceedances (with threshold concentrations set to 40 ppbV for top row and 60 ppbV for bottom row) and model performance statistical metrics of 8-h max ozone for the rural background stations on the left column (Cavo Greco, Inia, Agia Marina, and Troodos) and urban background stations on the right column (Larnaca, Limassol, and Nicosia)

The rural background concentrations of ozone over Cyprus are consistently 5--10 ppbV higher than urban levels. From the summary plot (Fig. [3](#Fig3){ref-type="fig"}, second row), we can see that in all rural background stations (important for vegetation and agricultural crops), the threshold of 40 ppbV is exceeded at all stations for more than 100 days. In urban background sites, where population density is larger, there occur 10--15 days when the 60-ppbV threshold is exceeded that are within the 25 exceedance days limit set by the EU. We note that the EU directive refers to the average---over 3 years---value; therefore, a continuous longer-term assessment should be performed for regulatory purposes.

Regarding model performance, rows 3 to 6 of the summary plot provide an overview of the temporal statistics for bias, correlation, standard deviation and a metric of the ability of the model to capture the highest range of concentration values (the latter still under development). Rows 7 and 8, in turn, provide metrics of the spatial performance of the model in terms of correlation and standard deviation. Average values for each station (over the selected time period) are calculated, and subsequently, the spatial correlation and standard deviation are computed. All indicators are normalized by the measurement uncertainty. In each statistical metric (row), the fulfillment of the performance criteria is shown with the green-shaded area. The orange area depicts regions where the criteria are fulfilled but model performance is dominated by the respective error. As shown in Fig. [3](#Fig3){ref-type="fig"}, the temporal and spatial statistical indicators are satisfied for both rural and urban background stations, for all metrics (bias, correlation, and standard deviation in time; correlation and standard deviation in space). The red dot in the temporal correlation MPI indicates that the time correlation--associated error is the one that dominates the model performance quality.

Nitrogen dioxide {#Sec10}
----------------

A similar analysis is performed for hourly nitrogen dioxide concentrations showing an overall fulfillment of the objectives set by the methodology. The mean, standard deviation, mean bias, and correlation reveal that there is a significant underestimation of the mean concentrations and variation of NO~2~ over the urban background stations with respective modeled values being 50% less than observed (Online Resource, Fig. [S2](#MOESM1){ref-type="media"}). The mean bias is low in the rural background stations; in urban areas, where the horizontal resolution of the regional model (10 km) cannot accurately reproduce the magnitude and variation of the local emissions, the bias is mostly negative between − 5 and − 6 ppbV. The correlation of the observed and modeled hourly values varies from 0.2 in rural background stations to 0.6 in urban background stations. The low correlation in rural stations is an outcome of the low and relatively constant NO~2~ concentrations in these areas, with no distinct diurnal cycle and background intensity. In urban areas, the correlations are higher but still dominated by possible discrepancies in the spatial and temporal evolution of emission fluxes related to residential activities. Overall, the NO~2~ levels over Cyprus are relatively low and well within the regulation limits.

The MQI are satisfied in both rural and urban stations as shown in Fig. [4](#Fig4){ref-type="fig"}. At the rural background stations (Fig. [4](#Fig4){ref-type="fig"}a), the objectives are met more closely with both annual and hourly MQI less than 0.1. At the urban stations (Fig. [4](#Fig4){ref-type="fig"}b), the MQI is again below unity with the hourly MQI = 0.5 and the yearly MQI = 0.7. The overall MQIs of all stations combined are closer to the values of the urban stations (0.5 and 0.66 for hourly and annual data, respectively). The two coastal stations (Larnaca and Limassol) are located on the left side of the plot indicating that correlation is the lowest performing metric, while Nicosia is on the right part indicating that standard deviation is captured poorly. Despite the low correlation, the model performance criteria for all stations are within the green area in the MPC plots (Fig. [4](#Fig4){ref-type="fig"}d) with smaller RMSu to *σ*~o~ ratio at urban stations as expected, mostly due to the low NO~2~ values over the island. There are no hourly or annual threshold exceedances at any of the stations for 2015. Figure [4](#Fig4){ref-type="fig"}e shows the large variation in observed nitrogen dioxide levels varying between \~ 1 ppbV for the rural background stations (well reproduced by the model) and \~ 8--10 ppbV over urban areas (strongly underestimated by the model with a value of \~ 2--4 ppbV). Meeting the target set by the methodology, despite the evident discrepancies between modeled and observed values in the urban stations, is a result of the low NO~2~ levels at these stations. Low concentrations are indeed associated to large measurement uncertainties. In this methodology, the performance criteria are strongly linked to the uncertainty gap, meaning that low concentrations with large measurement uncertainty lead to easily satisfied criteria. For policy applications, however, there is no need for a more rigorous assessment since annual and hourly concentrations are below the threshold set by the air quality directive.Fig. 4As in Fig. [2](#Fig2){ref-type="fig"} but for hourly NO~2~

The underestimation in NO~2~ levels may result from limitations in the resolution of national emissions in these areas as represented in the emission database (\~ 10-km resolution, same as the nested grid spacing). Several other studies have highlighted the discrepancies between different emission inventories, obtained with similar methodology and resolution (top-down regional or global inventories) or with different scales and approaches. Lopez et al. ([@CR40]) compared fine scale bottom-up urban emission inventories with regional top-down datasets and found large discrepancies in NOx and fine and coarse particulate matter, both in totals and sectorial emissions, with the regional emission inventories underestimating by 20--80% the NOx and 50--90% the PM10 emission fluxes. The authors indicated that the activity data based on fuel sales and population, used in the compilation of the regional inventories, versus actual traffic volume used in the bottom-up inventories, is the dominating reason for the underestimation of NO~2~ emissions. Moreover, even regional emission inventories based on top-down approaches exhibit substantial differences over urban areas, especially for NOx and VOCs, in terms of total, sectorial emissions, and spatial distribution, possibly due to downscaling approaches and choice of spatial proxies (Trombetti et al. [@CR56]). Thus, the accurate modeling of NOx fluxes at the urban scale requires higher resolution in both emission data and simulations and preferably bottom-up emission inventories.

Specifically for the country of Cyprus, studies indicate that the global top-down anthropogenic inventories do not adequately represent the gradient of emission fluxes between rural and populated areas, leading to an underestimation in modeled nitrogen oxide levels and overestimation of ozone levels in urban zones (Kushta et al. [@CR36]; Georgiou et al. [@CR18]). It is anticipated that the use of an up-to-date high-resolution (\< 2 km) and temporally resolved (diurnal to seasonal) national emission inventory encompassing ground-based and satellite information and an accompanying model configuration with similar grid spacing can help capture the magnitude of emission fluxes and ozone levels in urban areas.

Fine particulate matter {#Sec11}
-----------------------

The fine particulate matter (PM2.5) analysis is performed using daily mean values. Cyprus is frequently affected by dust storms from both North Africa and the Middle East (Kallos et al. [@CR32]; Gkikas et al. [@CR22]; Solomos et al. [@CR49], [@CR50]), and these events can contribute to the overall exceedances of the health safety thresholds set by the Air Quality Directive. The WRF-Chem model includes an online dust mobilization and transport module that simulates the sources and transport of mineral dust particles in the atmosphere. We note that an evaluation and fine tuning of its performance for the region of this study have not been performed in an integrated way, other than event-based analyses. The situation in the Eastern Mediterranean is complex due the proximity (and influence) of two large dust sources, North Africa and the Middle East.

During September 2015, a large dust event over the EMME region originated in Syria through the inert forcing of a thermal low developed in the area and strengthened by the convective activity over northern Iraq (haboob forming mechanism) that merged with and enhanced the dust cloud mobilization and transport (Gasch et al. [@CR16]; Solomos et al. [@CR49]). The study of Solomos et al. ([@CR49]) showed that the default land use--type dataset used by the model (USGS Global Land Cover Characteristics Data Base Version 2) does not include very active dust sources of uncultivated fields (especially around the Euphrates River) that are a result of a disruption in summer agriculture activities either due to climatic forcing in the region and/or, mainly, military conflict. These changes in land use due to war in the areas of northern Iraq and Syria significantly affect model performance. Additionally, relatively higher resolution (finer than 4 km over dust sources) is needed to accurately simulate convective processes and vertical motion for dust uplifting. Over Cyprus, the study also showed that the model cannot reproduce the magnitude of the observed fine and coarse particulate matter concentrations measured at the coastal and inland stations, as well as ground-based LIDAR vertical distributions. The authors attributed the discrepancies to several dust and atmospheric processes such as an underestimation of the intense downward mixing and emission fluxes at the source area, model limitations due to the size distribution and deposition rates. Thus, to avoid misleading results regarding the model performance over the rest of the year, we exclude the 2 days of the dust episode from our analysis (by omitting only PM2.5 concentrations above 100 μg m^−3^ on 8 and 9 September). This is the only exclusion performed on the observational data on a year during which other dust episodes also occurred.

For the PM2.5 assessment, we include traffic stations in the analysis since aerosol measurements are not available at urban background stations. We also include measurements performed at the Zygi station despite its classification as an industrial site. These decisions are driven by the scarcity of measurements. We however need to remember these choices when analyzing the results. From the rural background stations, only Agia Marina has a set of complete observations for the year 2015.

The model captures the magnitude and standard deviation of the fine particulate matter in the area with a slight overestimation over Agia Marina (rural background) and an underestimation of urban levels by 2--3 μg m^−3^ (Online Resource, Fig. [S3](#MOESM1){ref-type="media"}). Modeled yearly averaged PM2.5 concentrations are above the 10-μg m^−3^ threshold of the World Health Organization (WHO) for all four stations, but within the annual threshold of 25 (current) and 20 μg m^−3^ (future) as defined in the Air Quality Standards by the Directive 2008/50/EU. The observed annual mean at Agia Marina is below this value indicating that the model shows a non-measured annual exceedance. Exceedances at the other stations are adequately simulated by the model. The largest biases (negative in sign) and the largest correlations are over the urban stations of Larnaca and Nicosia (Online Resource, Fig. [S3](#MOESM1){ref-type="media"}c, d).

The model fulfills the performance criteria for normalized standard deviation and correlation coefficient (Fig. [5](#Fig5){ref-type="fig"}a). The error related to the correlation coefficient dominates the model performance for both background and industrial stations. This is also shown in the MPC plots (Fig. [5](#Fig5){ref-type="fig"}, bottom plots) with the MPC for correlation in the orange area (Agia Marina and Zygi). The scatter plot of the mean annual values reveals a smaller range of model variation of PM2.5 distribution from site to site than in the observations (Fig. [5](#Fig5){ref-type="fig"}b). All model results fall within the green target area and significantly close to the 1:1 line. We note that even though the model evaluation benchmarking methodology is not fit for traffic and industrial sites, the model objectives are also satisfied over these sites. As seen in Online Resource, Fig. [S3](#MOESM1){ref-type="media"}, there is no common pattern in the PM2.5 pollution characteristics of the traffic, industrial, and rural background stations. Therefore, the objective fulfillment highlights the need to scrutinize whether the MQO are stringent enough regarding daily fine particulate matter.Fig. 5**a** Target plot and **b** scatter plot of mean annual PM2.5 concentrations and model performance criteria for **c** correlation and **d** standard deviation calculated on mean daily PM2.5 values for two urban background stations (Nicosia and Larnaca), one rural background station (Agia Marina), and one industrial station (Zygi)

We also present the summary report and threshold exceedance analysis for PM2.5 (Fig. [6](#Fig6){ref-type="fig"}). There is no threshold value in the Air Quality Directive regarding mean daily PM2.5 exceedance days (row 2). We use the values of 20, 25, and 30 μg m^−3^ to assess the consistency of the methodology (the summary plot shown in Fig. [6](#Fig6){ref-type="fig"}a refers to 25 μg m^−3^). The thresholds set by the National Ambient Air Quality Standards (NAAQS) for the USA (65 μg m^−3^) and SEPA, Vehicle Emission Control Center for China in residential areas (50 μg m^−3^), are only exceeded for 1 and 2 days, respectively, during the dust event described above. As shown in the threshold exceedance plots (Fig. [6](#Fig6){ref-type="fig"}b for 20, Fig. [6](#Fig6){ref-type="fig"}c for 25, and Fig. [6](#Fig6){ref-type="fig"}d for 30 μg m^−3^), the model misses a number of exceedance days in the urban sites (Nicosia and Larnaca), overestimates the threshold exceedances in the rural background station of Agia Marina and the industrial site of Zygi for the thresholds below 25 μg m^−3^, and overestimates exceedance days in all stations for the upper threshold of 30 μg m^−3^. These results indicate a need for a more detailed investigation of the sources of pollution during days with high PM2.5 concentrations and the influence of natural versus anthropogenic pollution. It is anticipated that local emission information at higher resolution as well as a dedicated analysis of the natural aerosol sources in the region and their representation in the model will contribute to the improvement of the model performance regarding particulate matter. This is a necessary step towards creating methodologies to assess the contribution of such sources for policy-related applications.Fig. 6**a**--**d** Summary report (top panel) and threshold exceedance plot for 20 μg m^−3^ (left bottom plot), 25 μg m^−3^ (middle bottom plot), and 30 μg m^−3^ (right bottom plot) for mean daily PM2.5 concentrations for two urban background stations (Nicosia and Larnaca), one rural background station (Agia Marina), and one industrial station (Zygi)

Concluding remarks {#Sec12}
==================

We performed a statistical analysis of the regional atmospheric and chemistry model WRF-Chem using the benchmarking evaluation methodology developed in the framework of FAIRMODE, driven by the necessity of harmonizing model evaluation criteria for regulatory purposes. We applied the evaluation procedure over Cyprus, an EU member country located centrally in the Eastern Mediterranean and Middle East region, being an ideal case study due to its relatively isolated position at the crossroads of pollution from multiple sources.

The target plot of the model performance indicators delivers a combined matrix of the bias, standard deviation, and correlation of the model results with observations. The separate model performance criteria for each of these statistical metrics provide an in-depth analysis, highlighting the components that require improvement. Based on the case for Cyprus, we showed that the model is capable of reproducing the spatiotemporal distribution of ozone measured in the monitoring network that includes rural and urban areas, fulfilling the objectives set by the FAIRMODE evaluation methodology. Regarding nitrogen dioxide, there is a need for an improved, high-resolution representation of local emission fluxes to highlight the large urban-rural gradient in the observations. Particulate matter levels over the region are adequately simulated. However, the correlation between mean daily modeled and observed values needs further improvement. Overall, the model exhibits less variability in the pollutant concentration fields than the observations. The stringency of MQO related to PM2.5 must be further assessed.

The atmospheric aerosol distribution in the region is influenced by frequent mineral dust episodes, as well as other natural contributors such as sea salt and, especially during summer, black and organic carbon from accidental forest fires mostly upwind of the region (e.g., in southeastern Europe). For policy-related applications, the contribution of natural sources to mean daily and annual levels of pollutants must be well defined and the option to exclude the periods of episodes from the analysis must be provided, especially with regard to threshold exceedances. If this methodology is to be applied in forecasting applications, it must provide for solutions to separately assess events of episodic nature and distinguish between threshold exceedances from natural and anthropogenic sources.

Overall, the benchmarking evaluation methodology can provide comprehensive and detailed insights into the model performance related to regulated pollutants for scientific and policy applications. Since a large number of the objectives (and indicators) depend on the measurement uncertainty of the different pollutants, it is necessary to perform an assessment of this component for the region, due to its unique features regarding pollution, such as multiple source regions, chemical aging of atmospheric components during long-range transport, and regional climatic parameters.

Electronic supplementary material
=================================

 {#Sec13}

ESM 1(PDF 248 kb)

This work was performed within FAIRMODE (<http://fairmode.ew.eea.europa.eu>/); the community members of which are acknowledged for their indirect contribution through technical meetings. The observations over Cyprus were kindly provided by the Department of Labour Inspection of Cyprus, and the simulations were performed with computational resources from Cy-TERA under the European Union Horizon 2020 research and innovation programme VI-SEEM project (grant agreement 675121). The authors acknowledge the support of the InDUST Cost Action (action number CA16202).

This study received funding from the Government of Cyprus through the ESA Contract under the PECS (Plan for European Cooperating States) no. 4000124009/18NL/SC.

Disclaimer {#FPar1}
==========

The view expressed herein can in no way be taken to reflect the official opinion of the European Space Agency.


Comparison of anterior segment parameter values obtained with Scheimpflug-Placido topographer, optical low coherence reflectometry and noncontact specular microscopy in morbid obesity.
To investigate the measurement of anterior segment parameters using Sirius Scheimpflug-Placido topographer, Lenstar optical low coherence reflectometry (OLCR), and noncontact specular microscopy (SM) in morbidly obese and nonobese subjects. Twenty-eight morbidly obese subjects (BMI ≥ 40; Group 1) and 28 age- sex-matched healthy nonobese subjects (BMI 18.50-24.99; Group 2) were included in this study. Anterior segment parameters were measured by Scheimpflug-Placido topographer and OLCR. Corneal endothelial cell parameters were measured by non-contact SM. The group data were analyzed using the Mann-Whitney U test and Student's t-test. Bland-Altman plots were used to assess agreement among the instruments, and 95% limits of agreement (LoA) for each comparison were calculated. In group 1, the mean CCT by Scheimpflug-Placido topographer, OLCR, and noncontact SM were 549.44±30.10 µm, 544.15±31.48 µm, and 541.59±29.87 µm respectively. In group 2, the mean CCT by Scheimpflug-Placido topographer, OLCR, and noncontact SM were 531.0±22.09 µm, 523.15±21.39 µm, and 521.12±21.70 µm respectively. Mean CCT values obtained by the three methods were significantly higher in the morbidly obese than the nonobese subjects. In both groups, mean CCT was significantly higher when measured by Scheimpflug-Placido topographer than by OLCR and noncontact SM, and mean AD and ACD were significantly higher when measured by Scheimpflug-Placido topographer than OLCR. No significant differences were found between mean corneal curvature and corneal astigmatism when measured by Scheimpflug-Placido topographer and OLCR. The mean CCT of the morbidly obese subjects were significantly higher than the nonobese subjects when measured by all three methods. The CCT values obtained by Scheimpflug-Placido topographer were significantly higher than those by OLCR and SM.


Battleworld

Marvel has a pretty strange approach to ending its long standing universe and launching a new one simultaneously. It’s been a while since the main Secret Wars tale unfolded, with their focus shifting to all the new #1’s that have been pouring out of Marvel each week. There have been a few patches of confusion, […]

Jonathan Hickman is and has always been a strong writer in the world of comic books, and his contributions to Marvel since 2009 have been among his best. Heavy on science fiction theory and with deep character motivation and development, his incarnations of long standing heroes rank high above most. Secret Wars, so far, is […]

Post-apocalyptic stories where all but a handful of your favorite heroes are dead is a setting that has been around forever in comics. It’s a setting that, when done properly, can easily become a reader’s favorite. Marvel’s previous visit with the gray haired and broken Wolverine was a welcomed one and Old Man Logan #1 […]

With all the excitement and drama going on with Secret Wars, Marvel is still wrapping up a few characters in their final 616-universe moments. Of those characters, Magneto still has a lot going on. Recently the “Master of Magnetism” has gone pretty dark, almost as dark as his current costume. Now that the world knows […]

Some may think that Secret Wars is really starting to feel like a thinly veiled excuse for writers to create lengthy What If? plots. In the first issue of Amazing Spider-Man: Renew Your Vows, Dan Slott gets to dive into a reality where Peter and Mary Jane didn’t separate and we never had to live through One More […]

Welcome to Kabooooom! — where we review, rant, and rave about comics, movies, TV and all they inspire.


Polygenic risk for schizophrenia, transition and cortical gyrification: a high-risk study.
Schizophrenia is a highly heritable disorder, linked to several structural abnormalities of the brain. More specifically, previous findings have suggested that increased gyrification in frontal and temporal regions are implicated in the pathogenesis of schizophrenia. The current study included participants at high familial risk of schizophrenia who remained well (n = 31), who developed sub-diagnostic symptoms (n = 28) and who developed schizophrenia (n = 9) as well as healthy controls (HC) (n = 16). We first tested whether individuals at high familial risk of schizophrenia carried an increased burden of trait-associated alleles using polygenic risk score analysis. We then assessed the extent to which polygenic risk was associated with gyral folding in the frontal and temporal lobes. We found that individuals at high familial risk of schizophrenia who developed schizophrenia carried a significantly greater burden of risk-conferring variants for the disorder compared to those at high risk (HR) who developed sub-diagnostic symptoms or remained well and HC. Furthermore, within the HR cohort, there was a significant and positive association between schizophrenia polygenic risk score and bilateral frontal gyrification. These results suggest that polygenic risk for schizophrenia impacts upon early neurodevelopment to confer greater gyral folding in adulthood and an increased risk of developing the disorder.


The answer to that question is, apparently, a lot more:

Alec Baldwin apologizes and speaks out in support of the gay community http://t.co/dkdWjYLO4R — GLAAD (@glaad) June 28, 2013

Wow, that was easy.

Even TMZ wasn’t buying GLAAD’s quick acceptance of Baldwin’s “apology”:

GLAAD left many others bewildered as well:

@JimmyLaSalvia so @glaad plays the FOG (Friend of Gays) card and thinks that entitles Baldwin to a pass? #hypocrisywatch — Stacey-SisterToldjah (@sistertoldjah) June 28, 2013

Would love to see look on @rolandsmartin & @BrettRatner's face when they learn @GLAAD gave Alec Baldwin a TOTAL pass for FAR worse. — John Nolte (@NolteNC) June 28, 2013

So it's all good now? Alec Baldwin pardoned without so much as a bad boy slap?// http://t.co/0hFKyqDV84 — PJ Swan (@pjswan) June 28, 2013

What I learned from Alec Baldwin and @glaad today. It's no big if straight men refer to gay men as queens. Thanks GLAAD and Alec! — Rob Tam (@robtr2) June 28, 2013

It doesn't matter how you actually treat people, if you pretend to be ideologically compliant. http://t.co/fcjGsbxZ52 — Blasphemous Hat (@Popehat) June 29, 2013

When will he begin making the penitent PSA's? RT @glaad: Alec Baldwin apologizes and speaks out in support of the gay community — Kelly (@flyoverangel) June 28, 2013

In between shoots for Capital One commercials, perhaps?

Update:

GLAAD insists it’s not a free pass. It’s that Baldwin is well on his way to getting a pass.

Baldwin's apology is a step, not the end of the dialogue. He should now work to fight anti-gay violence http://t.co/dkdWjYLO4R — GLAAD (@glaad) June 29, 2013

Related:

‘Almost all of my best friends are gay’; Alec Baldwin’s wife defends ‘queen’ slur

Hilary Rosen knows Alec Baldwin’s heart (and mouth): ‘He’s not a homophobe’


Click to email this to a friend (Opens in new window)

Hillary came over. Ask not Hillary Who. Secret Service types stayed outside. We were alone. Just we two. An off-the-record old-time longtime friendship. No interview. No Q&A. No notes. No quotes. I’m just remembering. Should my memory prove faulty, you can cancel my order for a new Studebaker.

She favors Cuomo. Appreciates his stance on gun control.

About our next mayor. Ain’t her. She thinks somebody with a name and portfolio — such as Bloomberg, or someone like him.

About the conga line for president. No good words about no good Bernie Sanders. Stuff like: Anyone overtaking him in a district considered his, he’ll burn the place down.

She acknowledges Kamala Harris’ ability, calls non-winner Elizabeth Warren “smart as a whip,” likes Montana’s newest newie, Gov. Steve Bullock, and says each dwarf (my word, not hers) asked her advice.

Quote: “I remember when politics used to be fun.” And: “Russia will always meddle in our political system.”

She looked marvelous. Hair perfect. Makeup perfect. Outfit perfect. Carried only her phone. Three-piece, navy Akris pantsuit, matching rhinestone-studded Jimmy Choo pumps. Gold bracelet. No earrings.

Slim. “I’m dieting. I’m trying. I take walks. I’m working out. I have a beautiful backyard with surrounding trees, and flowers and greens, and I enjoy sitting outside reading.”

Coffee, fruit, a bagel, many laughs, and we told each other stories. Over a cream cheese schmear, she recalled her 1974 lawyer days when she was working on a Nixon memo describing exactly “what’s an impeachable event.”

A book on the horizon

October, Simon & Schuster will publish Hillary and daughter Chelsea’s “Book of Gutsy Women.” A compendium of Earth’s big-time females starting from B.C. — before computers — such as France’s Marie Antoinette, radium’s Marie Curie, from back when only boobs flashed their boobs and nobody could spell Kardashian. The book favors such first ladies as Eleanor Roosevelt, Betty Ford and Barbara Bush, yet somehow neglects to include Nancy Reagan. Discussing today’s world, she recalled how after 9/11 she and Chuck Schumer were the ones who pressed then President Bush’s $20 billion allocation to NYC. And she’ll be a grandma again. Her daughter’s having baby No. 3. It’s a busy household. Chelsea even has a Yorkie, age 12.

That’s not right

A publicist’s release says Chris Mann’s new single debuts “May 31th.” Congratulationth . . . Another says I’m in John King’s new paperback, and I don’t know who/what that is . . . Another: “ti action star” Paul Mormando’s in “Silent Partners” with “Django Unchained’s” James Russo . . . P.R. stuff’s as exciting as seeing Anthony Weiner in the shower.

Congrats, Isabella

Isabella Rossellini says she owes her new Master’s degree from Hunter College in Animal Behavior to her father, director Roberto Rossellini, and Oscar-winner mom, Ingrid Bergman . . . High-class, elegant, supremely VIP diner Dan Cathy was at the 21 Club. He owns Chick-fil-A. It’s at airports, restaurants, hospitals, etc. So this supremely elegant gent, what was he eating? Chicken hash.

Why do people condemn Russia? It’s really growing quite democratic. Pravda even began a letters-to-the-editor section this year. Not exactly a New York Post version. Theirs publishes all beefs, criticisms and complaints — but you have to give your name, address and next of kin.

Said Only in DC, kids, Only in DC.


/*
 * Minecraft Forge
 * Copyright (c) 2016-2020.
 *
 * This library is free software; you can redistribute it and/or
 * modify it under the terms of the GNU Lesser General Public
 * License as published by the Free Software Foundation version 2.1
 * of the License.
 *
 * This library is distributed in the hope that it will be useful,
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
 * Lesser General Public License for more details.
 *
 * You should have received a copy of the GNU Lesser General Public
 * License along with this library; if not, write to the Free Software
 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA  02110-1301  USA
 */

package net.minecraftforge.fml.common.gameevent;

import javax.annotation.Nonnull;
import net.minecraft.entity.item.EntityItem;
import net.minecraft.entity.player.EntityPlayer;
import net.minecraft.inventory.IInventory;
import net.minecraft.item.ItemStack;
import net.minecraftforge.fml.common.eventhandler.Event;

public class PlayerEvent extends Event {
    public final EntityPlayer player;
    private PlayerEvent(EntityPlayer player)
    {
        this.player = player;
    }

    public static class ItemPickupEvent extends PlayerEvent {
        @Deprecated
        public final EntityItem pickedUp;
        /**
        * Original EntityItem with current remaining stack size
        */
        private final EntityItem originalEntity;
        /**
         * Clone item stack, containing the item and amount picked up
         */
        private final ItemStack stack;
        public ItemPickupEvent(EntityPlayer player, EntityItem entPickedUp, ItemStack stack)
        {
            super(player);
            this.originalEntity = entPickedUp;
            this.pickedUp = entPickedUp;
            this.stack = stack;
        }

        public ItemStack getStack() {
            return stack;
        }

        public EntityItem getOriginalEntity() {
            return originalEntity;
        }
    }

    public static class ItemCraftedEvent extends PlayerEvent {
        @Nonnull
        public final ItemStack crafting;
        public final IInventory craftMatrix;
        public ItemCraftedEvent(EntityPlayer player, @Nonnull ItemStack crafting, IInventory craftMatrix)
        {
            super(player);
            this.crafting = crafting;
            this.craftMatrix = craftMatrix;
        }
    }
    public static class ItemSmeltedEvent extends PlayerEvent {
        @Nonnull
        public final ItemStack smelting;
        public ItemSmeltedEvent(EntityPlayer player, @Nonnull ItemStack crafting)
        {
            super(player);
            this.smelting = crafting;
        }
    }

    public static class PlayerLoggedInEvent extends PlayerEvent {
        public PlayerLoggedInEvent(EntityPlayer player)
        {
            super(player);
        }
    }

    public static class PlayerLoggedOutEvent extends PlayerEvent {
        public PlayerLoggedOutEvent(EntityPlayer player)
        {
            super(player);
        }
    }

    public static class PlayerRespawnEvent extends PlayerEvent {
        private final boolean endConquered;

        public PlayerRespawnEvent(EntityPlayer player, boolean endConquered)
        {
            super(player);
            this.endConquered = endConquered;
        }

        /**
         * Did this respawn event come from the player conquering the end?
         * @return if this respawn was because the player conquered the end
         */
        public boolean isEndConquered()
        {
            return this.endConquered;
        }

    }

    public static class PlayerChangedDimensionEvent extends PlayerEvent {
        public final int fromDim;
        public final int toDim;
        public PlayerChangedDimensionEvent(EntityPlayer player, int fromDim, int toDim)
        {
            super(player);
            this.fromDim = fromDim;
            this.toDim = toDim;
        }
    }
}


FW Trend Watch: Chokers

If this frenzied collage of AW 2014 runway head shots hasn’t convinced you of how necessary this look is, you should have seen both aerielist and lolitawoolf reppin solid black and diamante chokers at the White Material party this weekend. Click through for more proof of how omnipresent chokers, collars, and neck braces were at every fashion week this season…


For more info Contact

Details

Climb this distinctive peak on the Divide from Dry Gulch. Follow use trails and elk trails up steep slopes to the southeast ridge. Traverse west to climb scree/snow couloir and scramble north to summit. Helmet and ice axe required. Climb is intended to allow Alpine Scrambling Course students to practice field skills.


Effects of lidocaine on random skin flap survival in rats.
Use of a random skin flap is common for repairing wounds and for reconstruction. Lidocaine is a traditional local anesthetic that blocks sodium channels and has positive effects on ischemia-reperfusion injury. To investigate the effects of lidocaine on random skin flap survival in rats. McFarlane flaps were established in 20 rats divided into 2 groups. Lidocaine was injected in the lidocaine group, and the same concentration of saline was injected in the control group. The survival area of the flaps was measured on Day 7. Levels of inflammation were evaluated by hematoxylin and eosin (H&E)-stained slices, and superoxide dismutase and malonyldialdehyde contents were examined. The mean survival area of the flaps in the lidocaine group was significantly larger than that in the control group. Superoxide dismutase activity increased significantly in the lidocaine group compared with that in the control group. Malonyldialdehyde level in the lidocaine group was significantly lower than that in the control group. The H&E-stained slices showed that inflammation was clearly inhibited in the lidocaine group. Lidocaine improved the survival of random skin flaps.