The impact of biofuels for road traffic on air quality : a modelling approach

The selection of the energy source to power the transport sector is one of the main current concerns, not only relative with the energy paradigm but also due to the strong influence of road traffic in urban areas, which highly affects human exposure to air pollutants and human health and quality of life. Due to current important technical limitations of advanced energy sources for transportation purposes, biofuels are seen as an alternative way to power the world’s motor vehicles in a near-future, helping to reduce GHG emissions while at the same time stimulating rural development. Motivated by European strategies, Portugal, has been betting on biofuels to meet the Directive 2009/28/CE goals for road transports using biofuels, especially biodiesel, even though, there is unawareness regarding its impacts on air quality. In this sense, this work intends to clarify this issue by trying to answer the following question: can biodiesel use contribute to a better air quality over Portugal, particularly over urban areas? The first step of this work consisted on the characterization of the national biodiesel supply chain, which allows verifying that the biodiesel chain has problems of sustainability as it depends on raw materials importation, therefore not contributing to reduce the external energy dependence. Next, atmospheric pollutant emissions and air quality impacts associated to the biodiesel use on road transports were assessed, over Portugal and in particular over the Porto urban area, making use of the WRF-EURAD mesoscale numerical modelling system. For that, two emission scenarios were defined: a reference situation without biodiesel use and a scenario reflecting the use of a B20 fuel. Through the comparison of both scenarios, it was verified that the use of B20 fuels helps in controlling air pollution, promoting reductions on PM10, PM2.5, CO and total NMVOC concentrations. It was also verified that NO2 concentrations decrease over the mainland Portugal, but increase in the Porto urban area, as well as formaldehyde, acetaldehyde and acrolein emissions in the both case studies. However, the use of pure diesel is more injurious for human health due to its dominant VOC which have higher chronic hazard quotients and hazard indices when compared to B20.

Regional and urban air quality modelling under a climate change scenario

The better understanding of the interactions between climate change and air quality is an emerging priority for research and policy. Climate change will bring changes in the climate system, which will affect the concentration and dispersion of air pollutants. The main objective of the current study is to assess the impacts of climate change on air quality in 2050 over Portugal and Porto urban area. First, an evaluation and characterization of the air quality over mainland Portugal was performed for the period between 2002 and 2012. The results show that NO2, PM10 and O3 are the critical pollutants in Portugal. Also, the influence of meteorology on O3, NO2 and PM10 levels was investigate in the national main urban areas (Porto and Lisboa) and was verified that O3 has a statistically significant relationship with temperature in most of the components. The results also indicate that emission control strategies are primary regulators for NO2 and PM10 levels. After, understanding the national air quality problems and the influence that meteorology had in the historical air quality levels, the air quality modelling system WRF-CAMx was tested and the required inputs for the simulations were prepared to fulfil the main goal of this work. For the required air quality modelling inputs, an Emission Projections under RCP scenarios (EmiPro-RCP) model was developed to assist the estimation of future emission inventories for GHG and common air pollutants. Also, the current emissions were estimated for Portugal with a higher detailed disaggregation to improve the performance of the air quality simulations. The air quality modelling system WRF/CAMx was tested and evaluated over Portugal and Porto urban area and the results point out that is an adequate tool for the analysis of air quality under climate change. For this purpose, regional simulations of air quality during historical period and future (2045-2050) were conducted with CAMx version 6.0 to evaluate the impacts of simulated future climate and anthropogenic emission projections on air quality over the study area. The climate and the emission projections were produced under the RCP8.5 scenario. The results from the simulations point out, that if the anthropogenic emissions keep the same in 2050, the concentrations of NO2, PM10 and O3 will increase in Portugal. When, besides the climate change effects, is consider the projected anthropogenic emissions the annual mean concentrations of NO2 decrease significantly in Portugal and Porto urban area, and on the contrary the annual mean PM10 concentrations increases in Portugal and decrease in Porto urban area. The O3 results are mainly caused by the reduction of ozone precursors, getting the higher reductions in urban areas and increases in the surrounding areas. All the analysis performed for both simulations for Porto urban area support that, for PM10 and O3, there will be an increase in the occurrence of extreme values, surpassing the annual legislated parameters and having more daily exceedances. This study constitutes an innovative scientific tool to help in future air quality management in order to mitigate future climate change impacts on air quality.

Influência de plumas de água doce na hidrodinâmica do estuário da Ria de Arousa (NW Espanha), uma abordagem de modelação matemática

Localizada na parte sul da costa da Galiza, a ria de Arousa, tem um interesse particular devido à sua importante indústria de exploração de marisco, especialmente culturas de bivalves. O crescimento, qualidade e mortalidade destes bivalves não só são afetados por possíveis incidentes de poluição, mas também por episódios de grande run-off, que podem produzir decréscimos sazonais significativos de salinidade. Neste contexto, a presente dissertação pretende contribuir para uma melhor compreensão de como as plumas de água doce provenientes dos rios vão influenciar os gradientes de salinidade e temperatura no estuário da ria de Arousa. O objetivo é proceder a uma melhoria da previsão hidrodinâmica do modelo operacional oceanográfico desta ria, que é executado diariamente pela MeteoGalicia. Para atingir esse objectivo executaram-se simulações numéricas para diferentes cenários de forçamento. Foram escolhidos dois períodos com diferentes condições meteorológicas: um período de inverno, coincidindo com chuvas fortes e um período de verão, coincidindo com a ausência de precipitação. Foi aplicado o modelo hidrodinâmico MOHID (Hydrodynamic Model) para simular os processos físicos da ria de Arousa, utilizando a metodologia de modelos encaixados. O forçamento para a superfície foi imposto através do modelo atmosférico WRF (Weather Regional Forecast), que fornece as variáveis meteorológicas necessárias. Para fronteira terra, uma ferramenta hidrológica foi implementada, o modelo SWAT (Soil water Assessment Tool), este providencia as descargas de água doce para os principais rios da região. É feita neste trabalho uma descrição da implementação deste sistema, assim como a caracterização da área de estudo. Os resultados do modelo foram comparados com medidas provenientes de bóias oceanográficas e perfis de CTDs (Conductivity-Temperature-Depth instrument) em diferentes localizações.

Dynamics and variability of the alongshore flows on the northwestern iberian margin

Tese de doutoramento, Ciências Geofísicas e da Geoinformação (Oceanografia), Universidade de Lisboa, Faculdade de Ciências, 2014

Aerosol-Radiation Feedback and PM10 Air Concentrations Over Poland

We have implemented the WRF-Chem model version 3.5 over Poland to quantify the direct and indirect feedback effects of aerosols on simulated meteorology and aerosol concentrations. Observations were compared with results from three simulations at high spatial resolutions of 5 × 5 km: (1) BASE—without any aerosol feedback effects; (2) DIR—with direct aerosol-radiative effects (3) INDIR—with direct and indirect aerosol-radiative effects. We study the overall effect during January 2011 as well as selected episodes of the highest differences in PM10 concentrations between the three simulations. For the DIR simulation, the decrease in monthly mean incoming solar radiation (SWDOWN) appears for the entire study area. It changes geographically, from about −8.0 to −2.0 W m−2, respectively for the southern and northern parts of the country. The highest changes do not correspond to the highest PM10 concentration. Due to the solar radiation changes, the surface mean monthly temperature (T2) decreases for 96 % of the area of Poland, but not more than 1.0 °C. Monthly mean PBLH changes by more than ±5 m for 53 % of the domain. Locally the differences in PBLH between the DIR and BASE are higher than ± 20 m. Due to the direct effect, for 84 % of the domain, the mean monthly PM10 concentrations increase by up to 1.9 µg m−3. For the INDIR simulation the spatial distribution of changes in incoming solar radiation as well as air temperature is similar to the DIR simulation. The decrease of SWDOWN is noticed for the entire domain and for 23 % of the domain is higher than −5.0 W m−2. The absolute differences of PBLH are slightly higher for INDIR than DIR but similarly distributed spatially. For daily episodes, the differences between the simulations are higher, both for meteorology and PM10 concentrations, and the pattern of changes is usually more complex. The results indicate the potential importance of the aerosol feedback effects on modelled meteorology and PM10 concentrations.

Mudanças climáticas nas regiões insulares

Tese de Doutoramento, Física, 27 de Novembro de 2013, Universidade dos Açores.

Climatology of the Iberia coastal low-level wind jet: weather research forecasting model high-resolution results

Coastal low-level jets (CLLJ) are a low-tropospheric wind feature driven by the pressure gradient produced by a sharp contrast between high temperatures over land and lower temperatures over the sea. This contrast between the cold ocean and the warm land in the summer is intensified by the impact of the coastal parallel winds on the ocean generating upwelling currents, sharpening the temperature gradient close to the coast and giving rise to strong baroclinic structures at the coast. During summertime, the Iberian Peninsula is often under the effect of the Azores High and of a thermal low pressure system inland, leading to a seasonal wind, in the west coast, called the Nortada (northerly wind). This study presents a regional climatology of the CLLJ off the west coast of the Iberian Peninsula, based on a 9km resolution downscaling dataset, produced using the Weather Research and Forecasting (WRF) mesoscale model, forced by 19 years of ERA-Interim reanalysis (1989-2007). The simulation results show that the jet hourly frequency of occurrence in the summer is above 30% and decreases to about 10% during spring and autumn. The monthly frequencies of occurrence can reach higher values, around 40% in summer months, and reveal large inter-annual variability in all three seasons. In the summer, at a daily base, the CLLJ is present in almost 70% of the days. The CLLJ wind direction is mostly from north-northeasterly and occurs more persistently in three areas where the interaction of the jet flow with local capes and headlands is more pronounced. The coastal jets in this area occur at heights between 300 and 400 m, and its speed has a mean around 15 m/s, reaching maximum speeds of 25 m/s.

Wind resource modelling in complex terrain using different mesoscale–microscale coupling techniques

Wind resource evaluation in two sites located in Portugal was performed using the mesoscale modelling system Weather Research and Forecasting (WRF) and the wind resource analysis tool commonly used within the wind power industry, the Wind Atlas Analysis and Application Program (WAsP) microscale model. Wind measurement campaigns were conducted in the selected sites, allowing for a comparison between in situ measurements and simulated wind, in terms of flow characteristics and energy yields estimates. Three different methodologies were tested, aiming to provide an overview of the benefits and limitations of these methodologies for wind resource estimation. In the first methodology the mesoscale model acts like “virtual” wind measuring stations, where wind data was computed by WRF for both sites and inserted directly as input in WAsP. In the second approach, the same procedure was followed but here the terrain influences induced by the mesoscale model low resolution terrain data were removed from the simulated wind data. In the third methodology, the simulated wind data is extracted at the top of the planetary boundary layer height for both sites, aiming to assess if the use of geostrophic winds (which, by definition, are not influenced by the local terrain) can bring any improvement in the models performance. The obtained results for the abovementioned methodologies were compared with those resulting from in situ measurements, in terms of mean wind speed, Weibull probability density function parameters and production estimates, considering the installation of one wind turbine in each site. Results showed that the second tested approach is the one that produces values closest to the measured ones, and fairly acceptable deviations were found using this coupling technique in terms of estimated annual production. However, mesoscale output should not be used directly in wind farm sitting projects, mainly due to the mesoscale model terrain data poor resolution. Instead, the use of mesoscale output in microscale models should be seen as a valid alternative to in situ data mainly for preliminary wind resource assessments, although the application of mesoscale and microscale coupling in areas with complex topography should be done with extreme caution.

Computational models for the prediction of severe thunderstorms over east indian region

Thunderstorm is one of the most spectacular weather phenomena in the atmosphere. Many parts over the Indian region experience thunderstorms at higher frequency during pre-monsoon months (March- May), when the atmosphere is highly unstable because of high temperatures prevailing at lower levels. Most dominant feature of the weather during the pre-monsoon season over the eastern Indo-Gangetic plain and northeast India is the outburst of severe local convective storms, commonly known as ‘Nor’wester’ or ‘Kalbaishakhi’. The severe thunderstorms associated with thunder, squall line, lightning and hail cause extensive losses in agriculture, damage to structure and also loss of life. The casualty due to lightning associated with thunderstorms in this region is the highest in the world. The highest numbers of aviation hazards are reported during occurrence of these thunderstorms. In India, 72% of tornadoes are associated with this thunderstorm.

Ice formation and development in aged, wintertime cumulus over the UK: observations and modelling

In situ high resolution aircraft measurements of cloud microphysical properties were made in coordination with ground based remote sensing observations of a line of small cumulus clouds, using Radar and Lidar, as part of the Aerosol Properties, PRocesses And InfluenceS on the Earth's climate (APPRAISE) project. A narrow but extensive line (~100 km long) of shallow convective clouds over the southern UK was studied. Cloud top temperatures were observed to be higher than −8 °C, but the clouds were seen to consist of supercooled droplets and varying concentrations of ice particles. No ice particles were observed to be falling into the cloud tops from above. Current parameterisations of ice nuclei (IN) numbers predict too few particles will be active as ice nuclei to account for ice particle concentrations at the observed, near cloud top, temperatures (−7.5 °C). The role of mineral dust particles, consistent with concentrations observed near the surface, acting as high temperature IN is considered important in this case. It was found that very high concentrations of ice particles (up to 100 L−1) could be produced by secondary ice particle production providing the observed small amount of primary ice (about 0.01 L−1) was present to initiate it. This emphasises the need to understand primary ice formation in slightly supercooled clouds. It is shown using simple calculations that the Hallett-Mossop process (HM) is the likely source of the secondary ice. Model simulations of the case study were performed with the Aerosol Cloud and Precipitation Interactions Model (ACPIM). These parcel model investigations confirmed the HM process to be a very important mechanism for producing the observed high ice concentrations. A key step in generating the high concentrations was the process of collision and coalescence of rain drops, which once formed fell rapidly through the cloud, collecting ice particles which caused them to freeze and form instant large riming particles. The broadening of the droplet size-distribution by collision-coalescence was, therefore, a vital step in this process as this was required to generate the large number of ice crystals observed in the time available. Simulations were also performed with the WRF (Weather, Research and Forecasting) model. The results showed that while HM does act to increase the mass and number concentration of ice particles in these model simulations it was not found to be critical for the formation of precipitation. However, the WRF simulations produced a cloud top that was too cold and this, combined with the assumption of continual replenishing of ice nuclei removed by ice crystal formation, resulted in too many ice crystals forming by primary nucleation compared to the observations and parcel modelling.

Effects of spatial resolution in the simulation of daily and subdaily precipitation in the southwestern US

We evaluate the effects of spatial resolution on the ability of a regional climate model to reproduce observed extreme precipitation for a region in the Southwestern United States. A total of 73 National Climate Data Center observational sites spread throughout Arizona and New Mexico are compared with regional climate simulations at the spatial resolutions of 50 km and 10 km for a 31 year period from 1980 to 2010. We analyze mean, 3-hourly and 24-hourly extreme precipitation events using WRF regional model simulations driven by NCEP-2 reanalysis. The mean climatological spatial structure of precipitation in the Southwest is well represented by the 10 km resolution but missing in the coarse (50 km resolution) simulation. However, the fine grid has a larger positive bias in mean summer precipitation than the coarse-resolution grid. The large overestimation in the simulation is in part due to scale-dependent deficiencies in the Kain-Fritsch convective parameterization scheme that generate excessive precipitation and induce a slow eastward propagation of the moist convective summer systems in the high-resolution simulation. Despite this overestimation in the mean, the 10 km simulation captures individual extreme summer precipitation events better than the 50 km simulation. In winter, however, the two simulations appear to perform equally in simulating extremes.

High resolution simulation of the variability of surface energy balance fluxes across central London with urban zones for energy partitioning

The parameterization of surface heat-flux variability in urban areas relies on adequate representation of surface characteristics. Given the horizontal resolutions (e.g. ≈0.1–1km) currently used in numerical weather prediction (NWP) models, properties of the urban surface (e.g. vegetated/built surfaces, street-canyon geometries) often have large spatial variability. Here, a new approach based on Urban Zones to characterize Energy partitioning (UZE) is tested within a NWP model (Weather Research and Forecasting model;WRF v3.2.1) for Greater London. The urban land-surface scheme is the Noah/Single-Layer Urban Canopy Model (SLUCM). Detailed surface information (horizontal resolution 1 km)in central London shows that the UZE offers better characterization of surface properties and their variability compared to default WRF-SLUCM input parameters. In situ observations of the surface energy fluxes and near-surface meteorological variables are used to select the radiation and turbulence parameterization schemes and to evaluate the land-surface scheme

Multi-site evaluation of an urban land-surface model: intra-urban heterogeneity, seasonality and parameter complexity requirements

An extensive off-line evaluation of the Noah/Single Layer Urban Canopy Model (Noah/SLUCM) urban land-surface model is presented using data from 15 sites to assess (1) the ability of the scheme to reproduce the surface energy balance observed in a range of urban environments, including seasonal changes, and (2) the impact of increasing complexity of input parameter information. Model performance is found to be most dependent on representation of vegetated surface area cover; refinement of other parameter values leads to smaller improvements. Model biases in net all-wave radiation and trade-offs between turbulent heat fluxes are highlighted using an optimization algorithm. Here we use the Urban Zones to characterize Energy partitioning (UZE) as the basis to assign default SLUCM parameter values. A methodology (FRAISE) to assign sites (or areas) to one of these categories based on surface characteristics is evaluated. Using three urban sites from the Basel Urban Boundary Layer Experiment (BUBBLE) dataset, an independent evaluation of the model performance with the parameter values representative of each class is performed. The scheme copes well with both seasonal changes in the surface characteristics and intra-urban heterogeneities in energy flux partitioning, with RMSE performance comparable to similar state-of-the-art models for all fluxes, sites and seasons. The potential of the methodology for high-resolution atmospheric modelling application using the Weather Research and Forecasting (WRF) model is highlighted. This analysis supports the recommendations that (1) three classes are appropriate to characterize the urban environment, and (2) that the parameter values identified should be adopted as default values in WRF.

Contribution of land use changes to near-surface air temperatures during recent summer extreme heat events in the Phoenix Metropolitan Area

The impact of 1973–2005 land use–land cover (LULC) changes on near-surface air temperatures during four recent summer extreme heat events (EHEs) are investigated for the arid Phoenix, Arizona, metropolitan area using the Weather Research and Forecasting Model (WRF) in conjunction with the Noah Urban Canopy Model. WRF simulations were carried out for each EHE using LULC for the years 1973, 1985, 1998, and 2005. Comparison of measured near-surface air temperatures and wind speeds for 18 surface stations in the region show a good agreement between observed and simulated data for all simulation periods. The results indicate consistent significant contributions of urban development and accompanying LULC changes to extreme temperatures for the four EHEs. Simulations suggest new urban developments caused an intensification and expansion of the area experiencing extreme temperatures but mainly influenced nighttime temperatures with an increase of up to 10 K. Nighttime temperatures in the existing urban core showed changes of up to 2 K with the ongoing LULC changes. Daytime temperatures were not significantly affected where urban development replaced desert land (increase by 1 K); however, maximum temperatures increased by 2–4 K when irrigated agricultural land was converted to suburban development. According to the model simulations, urban landscaping irrigation contributed to cooling by 0.5–1 K in maximum daytime as well as minimum nighttime 2-m air temperatures in most parts of the urban region. Furthermore, urban development led to a reduction of the already relatively weak nighttime winds and therefore a reduction in advection of cooler air into the city.

Toward seamless prediction: calibration of climate change projections using seasonal forecasts

The CWRF is developed as a climate extension of the Weather Research and Forecasting model (WRF) by incorporating numerous improvements in the representation of physical processes and integration of external (top, surface, lateral) forcings that are crucial to climate scales, including interactions between land, atmosphere, and ocean; convection and microphysics; and cloud, aerosol, and radiation; and system consistency throughout all process modules. This extension inherits all WRF functionalities for numerical weather prediction while enhancing the capability for climate modeling. As such, CWRF can be applied seamlessly to weather forecast and climate prediction. The CWRF is built with a comprehensive ensemble of alternative parameterization schemes for each of the key physical processes, including surface (land, ocean), planetary boundary layer, cumulus (deep, shallow), microphysics, cloud, aerosol, and radiation, and their interactions. This facilitates the use of an optimized physics ensemble approach to improve weather or climate prediction along with a reliable uncertainty estimate. The CWRF also emphasizes the societal service capability to provide impactrelevant information by coupling with detailed models of terrestrial hydrology, coastal ocean, crop growth, air quality, and a recently expanded interactive water quality and ecosystem model. This study provides a general CWRF description and basic skill evaluation based on a continuous integration for the period 1979– 2009 as compared with that of WRF, using a 30-km grid spacing over a domain that includes the contiguous United States plus southern Canada and northern Mexico. In addition to advantages of greater application capability, CWRF improves performance in radiation and terrestrial hydrology over WRF and other regional models. Precipitation simulation, however, remains a challenge for all of the tested models.