Global transport demand is increasing and the associated emissions of GHGs are becoming a key factor in the environmental planning of networks and the actual operational management stage. We describe the development of a novel methodology for the calculation of the carbon footprint of a major motorway (Egnatia Odos) in Northern Greece. Total GHG emissions were calculated for 2014 as the base year, from all corporate activities, including stationary and mobile combustion, fugitive gases, purchased electricity, personnel travel, waste disposal and contractor (third party) emissions. End user vehicle emissions constitute a significant part of the carbon footprint and are estimated from annual traffic data supplied by Egnatia Odos S.A., combined with emission factors extracted from the TREMOVE model, for different vehicle categories and fuel types. The emissions model developed uses additional input data, concerning fuel and electricity consumption, fugitive gases from a/c equipment, personnel commuting and travel, waste disposal, and contractor data, also provided by Egnatia Odos S.A., for the calculation of direct and indirect GHG emissions for all three Scopes of the GHG Protocol. The final product was developed in the form of a specialized electronic platform, which is user-friendly and can be updated with the most recent input data available, in order to be used by highways’ personnel for the calculation of their carbon footprint in future years. The overall carbon footprint of Egnatia Odos S.A. in 2014 has been estimated at 770 ktCO₂e, with 91% of emissions sourced by the end user category (i.e. vehicle emissions).

Mediterranean tropical-like cyclones, or Medicanes, present unique challenges for precipitation estimations due to their rapid development and localized impacts. This study evaluates the performance of satellite precipitation products in capturing the precipitation associated with Medicane Daniel that struck Greece in early September 2023. Utilizing a combination of ground-based observations, reanalysis, and satellite-derived precipitation data, we assess the accuracy and spatial distribution of the satellite precipitation products GPM IMERG, GSMaP, and CMOPRH during the cyclone event, which formed in the Eastern Mediterranean from 4 to 7 September 2023, hitting with unprecedented, enormous amounts of rainfall, especially in the region of Thessaly in central Greece. The results indicate that, while satellite precipitation products demonstrate overall skill in capturing the broad-scale precipitation patterns associated with Medicane Daniel, discrepancies exist in estimating localized intense rainfall rates, particularly in convective cells within the cyclone’s core. Indeed, most of the satellite precipitation products studied in this work showed a misplacement of the highest amounts of associated rainfall, a significant underestimation of the event, and large unbiased root mean square error in the areas of heavy precipitation. The total precipitation field from IMERG Late Run and CMORPH showed the smallest bias (but significant) and good temporal correlation against rain gauges and ERA5-Land reanalysis data as a reference, while IMERG Final Run and GSMaP showed the largest underestimation and overestimation, respectively. Further investigation is needed to improve the representation of extreme precipitation events associated with tropical-like cyclones in satellite precipitation products.

In recent decades, large-scale urbanisation has developed rapidly, resulting in significant changes in the local and regional environment and climate. Large metropolitan areas worldwide induce significant changes in local atmospheric circulation and boundary layer meteorology by modifying the underlying surface characteristics and through the emission of anthropogenic heat and pollutants into the atmosphere. We investigate the urban heat island (UHI) characteristics in the city of Athens, Greece, which is one of Europe’s largest metropolitan complexes with a population of approximately 3.7 million inhabitants. The UHI effect is intense due to the city’s size, dense construction, high incident solar radiation, and almost complete lack of natural vegetation, with previous studies suggesting a temperature rise of 4 °C on average in the city centre compared to summer background conditions. We used high-resolution WRF simulations (1-km horizontal grid) driven with ERA5 reanalysis data to produce surface temperature maps in the city of Athens and the surrounding areas (Region of Attiki) during the summer period of 1 July–20 August 2021. Different model parameterizations were tested, both with respect to urban characteristics and physical parameters. The daily minimum and maximum temperatures (T_min and T_max) derived from the model were validated against observational data from a dense network of weather stations covering metropolitan Athens and surrounding locations. We further investigate the influence of different meteorological conditions on the UHI gradients as produced by the model and the observational dataset, including the extreme heat wave of 28 July–5 August 2021, during which persistent maximum temperatures of >40 °C were recorded for nine consecutive days. The results indicate a strong correlation between WRF output and recorded minimum and maximum temperatures throughout the test period (R ranges from 0.80 to 0.93), with an average mean absolute bias (MAB) of 1.5 °C, and reveal the intensity and spatiotemporal variability of the UHI phenomenon in the city of Athens, with UHI magnitude reaching 8–9 °C at times. Our work aims to maximise the potential of using high-resolution WRF modelling for simulating extreme heat events and mapping the UHI effect in large metropolitan complexes.

Synthetic Aperture Radar Interferometry (InSAR) is a space geodetic technique used for mapping deformations of the Earth’s surface. It has been developed and used increasingly during the last thirty years to measure displacements produced by earthquakes, volcanic activity and other crustal deformations. A limiting factor to this technique is the effect of the troposphere, as spatial and temporal variations in temperature, pressure, and relative humidity introduce significant phase delays in the microwave imagery, thus “masking” surface displacements due to tectonic or other geophysical processes. The use of Numerical Weather Prediction (NWP) models as a tropospheric correction method in InSAR can tackle several of the problems faced with other correction techniques (such as timing, spatial coverage and data availability issues). High-resolution tropospheric modelling is particularly useful in the case of single interferograms, where the removal of the atmospheric phase screen (and especially the highly variable turbulent component) can reveal large-amplitude deformation signals (as in the case of an earthquake). In the western Gulf of Corinth, prominent topography makes the removal of both the stratified and turbulent atmospheric phase screens a challenging task. Here, we investigate the extent to which a high-resolution WRF 1-km re-analysis can produce detailed tropospheric delay maps of the required accuracy by coupling its output (in terms of Zenith Total Delay or ZTD) with the vertical delay component in GNSS measurements. The model is operated with varying physical parameterization in order to identify the best configuration, and validated with GNSS zenithal tropospheric delays, providing a benchmark of real atmospheric conditions. We correct sixteen Sentinel-1A interferograms with differential delay maps at the line-of-sight (LOS) produced by WRF re-analysis. In most cases, corrections lead to a decrease in the phase gradient, with average root-mean-square (RMS) and standard deviation (SD) reductions in the wrapped phase of 6.0% and 19.3%, respectively. Results suggest a high potential of the model to reproduce both the long-wavelength stratified atmospheric signal and the short-wave turbulent atmospheric component which are evident in the interferograms.

In the last thirty years, Synthetic Aperture Radar interferometry (InSAR) and the Global Navigation Satellite System (GNSS) have become fundamental space geodetic techniques for mapping surface deformations due to tectonic movements. One major limiting factor to those techniques is the effect of the troposphere, as surface velocities are of the order of a few mm yr⁻¹, and high accuracy (to mm level) is required. The troposphere introduces a path delay in the microwave signal, which, in the case of GNSS Precise Point Positioning (PPP), can nowadays be partly removed with the use of specialized mapping functions. Moreover, tropospheric stratification and short wavelength spatial turbulences produce an additive noise to the low amplitude ground deformations calculated by the (multitemporal) InSAR methodology. InSAR atmospheric phase delay corrections are much more challenging, as opposed to GNSS PPP, due to the single pass geometry and the gridded nature of the acquired data. Thus, the precise knowledge of the tropospheric parameters along the propagation medium is extremely useful for the estimation and correction of the atmospheric phase delay. In this context, the PaTrop experiment aims to maximize the potential of using a high-resolution Limited-Area Model for the calculation and removal of the tropospheric noise from InSAR data, by following a synergistic approach and integrating all the latest advances in the fields of remote sensing meteorology (GNSS and InSAR) and Numerical Weather Forecasting (WRF). In the first phase of the experiment, presented in the current paper, we investigate the extent to which a high-resolution 1 km WRF weather re-analysis can produce detailed tropospheric delay maps of the required accuracy, by coupling its output (in terms of Zenith Total Delay or ZTD) with the vertical delay component in GNSS measurements. The model is initially operated with varying parameterization, with GNSS measurements providing a benchmark of real atmospheric conditions. Subsequently, the final WRF daily re-analysis run covers an extended period of one year, based on the optimum model parameterization scheme demonstrated by the parametric analysis. The two datasets (predicted and observed) are compared and statistically evaluated, in order to investigate the extent to which meteorological parameters that affect ZTD can be simulated accurately by the model under different weather conditions. Results demonstrate a strong correlation between predicted and observed ZTDs at the 19 GNSS stations throughout the year (R ranges from 0.91 to 0.93), with an average mean bias (MB) of –19.2 mm, indicating that the model tends to slightly underestimate the tropospheric ZTD as compared to the GNSS derived values. With respect to the seasonal component, model performance is better during the autumn period (October–December), followed by spring (April–June). Setting the acceptable bias range at ±23 mm (equal to the amplitude of one Sentinel-1 C-band phase cycle when projected to the zenithal distance), it is demonstrated that the model produces satisfactory results, with a percentage of ZTD values within the bias margin ranging from 57% in summer to 63% in autumn.