In many countries the number of wind turbines is growing rapidly as a response to the increasing demandfor renewable energy.Modern wind turbines are large structures, many reach more than 150 meters above theground. Clusters of densely spaced wind turbines, so called wind farms, are being built both on- and offshore. Wind farm installations relatively near to radar systems generate clutter returns that usually affect the normal operation of these radars. Interference caused by wind turbines is more severe for many radar systems than interference caused by stationary objects such as masts or towers. This is due to the rotating blades of the wind turbines. Many Doppler radars use a filter that removes echoes originating from objects with no or little radial velocity. However, these filters do not work for rotating objects such as the rotating blades of wind turbines. Wind turbines located around the line of sight of Doppler radars can cause clutter, blockage, and erroneous velocity measurements, affecting the performance of both military and civilian radar systems. As a result, the unwanted radar return from wind farms, known as Wind Turbine Clutter (WTC), is considered to be dynamic clutter due to the nonzero Doppler return created by rotating wind turbine blades. Nowadays numerous radar systems are developed in order to exploit the diverse information obtained through transmission of waves with different polarizations. This technique is widely known as polarimetry. Many targets of interest exhibit Radar Cross Sections which vary with different transmitted and received polarizations. Wind Turbines also experiences this variability. In this thesis we propose a method to optimal detect the presence of WTC with the use of radar polarimetry. Since the crucial part of this interference comes from the blades rotation, we initially propose a method to estimate the angular velocity of these blades. The estimation of this parameter is derived with the use of proper combination of maximum likelihood estimation theory and radar polarimetry. As there is absence of Micro-Doppler when the radar beam axis and rotation coincide, a separate estimator for this case is pro-posed. In the final part of this thesis, we present a detection approach based on the same signal model used for angular velocity estimation. Again we define a detection rule for the case when radar beam axis and rotation axis coincide and one when they do not. Although at some extent the used model for the second case is valid for low frequencies (f<1 GHz), both estimator and detector derivations can be further applied for higher frequencies signal models. All these mathematical derivations are accompanied with proper simulations.

During a solar eclipse, sunlight incident on the Earth is reduced due to the (partial) shadow of the Moon. Atmospheric trace gas concentrations which are influenced by the amount of available sunlight, such as nitrogen dioxide (NO2), may be affected due to the disrupted photolysis processes. Large-scale observations of the increased NO2 concentrations caused by the solar eclipse would improve our understanding of the sensitivity of NO2 in the atmosphere to short-term variations in sunlight. Spaceborne measurements can provide valuable information about the large-scale spatial distribution of NO2, which is provided daily by the TROPOMI instrument aboard the Sentinel-5 Precursor satellite by measuring and retrieving locally reflected sunlight. However, the TROPOMI NO2 retrieval is unable to derive reliable concentrations during a solar eclipse, as solar eclipses are not taken into account in its retrieval algorithm. In this research, we have adjusted the NO2 retrieval of TROPOMI such that it can handle solar eclipses and study the large-scale response of NO2 during two solar eclipses over Europe in 2021 and 2022. We found a large-scale increase of NO2 in the adjusted measurements, which linearly correlated with the degree of obscuration. We compared the measured NO2 increase with the values from the atmospheric chemistry model TM5 including an applied eclipse implementation and we found a close agreement in most areas that are not highly polluted. Our measurements and model predict a NO2 increase of 60%±12% and 70%±7% for an obscuration fraction of 1, respectively. More advanced chemistry modelling work is needed to explain the measurements in highly populated areas. We conclude that our results demonstrate that the TROPOMI algorithm is capable of correctly measuring NO2 after an adjustment of the NO2 retrieval. We have shown that it is possible to adjust an atmospheric trace gas retrieval for the influence of a solar eclipse. Moreover, we are the first to provide evidence for an increase in NO2 during a solar eclipse using space-based measurement techniques and to quantify this increase on a large scale with the same instrument. Our measurements can be used to test atmospheric chemistry models, possibly improving their sensitivity to solar eclipses but also artificial shadows on the Earth induced by sunlight-intercepting geoengineering approaches.

Lightning is a natural phenomena that can be dangerous to humans. It is however challenging to study thunderstorm clouds using direct observations since it can be dangerous to fly into thunderstorm clouds. In this study, cloud radar with millimeter wavelength is used to study the properties and dynamics of thunderstorm clouds. It is based on a case of thunderstorm on 2021-06-18 from 16:10 to 17:45 UTC near Cabauw. Polarimetric radar variables are used to investigate possible hydrometeors in the clouds and look for vertical alignment of ice crystals that is expected due to electric torque. The technique of Doppler spectra analysis, which has not been used in previous studies about thunderstorms so far, is used to help understand the behaviours of different types of particles within a radar resolution volume. Due to challenges posed by Mie scattering, scattering simulations are carried out to aid the interpretation of spectral polarimetric variables. From the results, there is a high chance that supercooled liquid water and conical graupel are present in thunderstorm clouds. There is also a possibility of ice crystals arranged in chains at the cloud top. Ice crystals become vertically aligned a few seconds before lightning and return to their usual horizontal alignment afterwards. However, this phenomenon has been witnessed in only a few cases, specifically when the lightning strike is in close proximity to the radar's line of sight or when the lightning is exceptionally strong. Doppler analyses show that updrafts are found near the core of the thunderstorm cloud, while downdrafts are observed at the edges. Strong turbulence is also observed as reflected by the large Doppler spectrum width.

Severe wind gusts associated with mid-latitude convective storms contribute to an increasing amount of natural hazard related losses in Europe. Modifications of the atmosphere associated with anthropogenic climate change are projected to increase the frequency of favorable conditions for convective storms, and it is absolutely critical to understand the implications of these changes to prepare for a resilient future. To this end, the fate of convective gusts in Europe in a future warmer climate is addressed through the output of two high resolution regional climate models (RCMs). One RCM includes convection permitting (CP) physics, and the other assumes hydrostatic conditions in which deep convection is parameterized. It is found that the magnitude and characteristics of extreme straight line gusts from mesoscale convective systems are well resolved in the CP-RCM but not in the hydrostatic RCM. The RCMs are forced with a high carbon emission scenario for the end of the 21st century, and the CP-RCM shows an increase in the frequency and magnitude of extreme convective gusts over mainland Europe. These changes are likely related to an increase in conditions where strong wind shear (≥15 m/s) simultaneously occurs with unstable environments (lifted index ≤-2). However, the inherent low frequency of extreme convective storms requires continued investigation to draw more robust conclusions. In addition, the environmental conditions in which the severe gusts take place indicate the importance of addressing gusts separately from the more commonly studied effects of climate change on extreme convective precipitation. This thesis provides an important bridge to understand the fate of extreme gusts from convective storms in a future warmer climate and the high potential of CP-RCMs in such studies.

Oceans cover a significant part of the Earth's surface. The coupling between the upper ocean and the atmosphere is very complicated with defied theoretical understanding, while it is essential for climate studies, weather prediction, and marine ecosystems. With the advent of spaceborne Synthetic Aperture Radar (SAR) systems, surface signatures of ocean and atmospheric processes have been revealed. As winds blowing over the ocean excite the wind waves, all undulations of the ocean surface are assumed as waves in this study. The primary sources for ocean surface signatures in SAR images are waves that are created by the exertion of the local wind stress. Wind waves cause changes in the backscattered power due to three mechanisms: specular reflections, Bragg scattering, and a contribution from wave breaking. A statistical multi-static normalized radar cross-section (NRCS) background model in terms of the directional wave spectrum is developed, considering both Bragg and non-Bragg mechanisms for various polarization states. As the qualitative comparison between optical and SAR data reveals a significant correlation in sea surface signatures, a synthetic attempt is made to estimate the SAR signals from optical signatures. This is realized with the transformation of the wave spectrum in a nonuniform medium, as a consequence of surface currents, and varying near-surface wind fields. A comparison between modeled NRCS and observations is presented. This modulated NRCS model advances the quantitative interpretation of the upper ocean dynamics from satellite measurements.

Stratiform precipitation is one of the most important precipitation systems in the mid-latitudes. To gain understanding of the melting layer, which is part of a stratiform precipitating system, an algorithm which is able to select melting layer data from large datasets is required. The goal of this thesis work is to create and deliver a program which is able to performthis work. For development the available data was separated into three groups: no melting layer, only melting layer and the rest. Based on the second group, typical melting layer signatures are analysed. In particular the different signatures (reflectivity, polarimetric and Doppler) occur at different heights. Based on the knowledge gained from the analysis and focussing on the reflectivity and polarimetric signatures, three different approaches were taken to detect and characterise the melting layer fromthe data. The first approach is based on an existing method in literature, this approach acts as the reference method. The second approach is based on image processing techniques, while the third approach is based on machine learning techniques. The second approach is later abandoned because of limitations of the techniques investigated. The third approach, machine learning, is the main contribution of this thesis work. For analysis of the performance of the different algorithms an annotated test dataset is created which represents the entire dataset. The performance is determined by the melting layer detection probability using a confusion matrix and by the determination of errors of the melting layer boundaries. The reference method proved to show an extremely low false positive rate (0%) on the test dataset. This means that if the method detected a melting layer, it is almost certain that there is one. The overall detection probability was 80%. The method fails in detecting the melting layer when the peak reflectivity is below the threshold (30dBZ) used in the method. The detected upper melting layer boundary of the reference method is on average 365 m lower compared to the ground truth. The lower boundary is on average 149 m higher than the ground truth. This means that the reference method only selects a part of the entire thickness of the ML. The correlation between the lower boundary and the ground truth is higher than the upper boundary (0.987 vs. 0.965). The proposed machine learning method has a detection performance of almost 94%, which is higher than the reference method, but some false positives occur (3%). The upper boundary is on average 20 m above the ground truth, while the lower boundary is 66 m lower than the ground truth. The machine learned method thus overestimates the thickness of the ML. The correlation between the lower boundary and the true boundary is 0.976 and the correlation between the upper boundary and the true boundary is 0.966. The upper boundaries of both methods have a very similar correlation, while the lower boundary of the machine learning method has a lower correlation. However, despite the lower correlation of the lower boundary and the introduction of false positives, the machine learning method is an improvement over the reference method since it has a significantly higher detection probability and it captures the entire thickness of theML much better. It is therefore suited forML analysis.

The knowledge of the raindrop size distribution is key for characterizing precipitation. It is however still a challenge to retrieve it with radars. Several polarimetric and spectral techniques are proposed for cm-wavelength radars (weather radars). What about the mm-wavelength radars (cloud radars), which have a better spatial and time resolution and can still measure light and moderate rain? Knowing that 90% of the rain volume in Europe comes from rainfall rates between 0.1 mm/h and 10 mm/h, this is worthwhile to investigate. The goal of this thesis is to retrieve 1 of the 3 parameters of the modelled gamma raindrop size distribution, the median volume diameter (D0), during stratiform rainfall events using a slantwise profiling dual-frequency polarimetric cloud radar. Focus is given to phase measurements, which are not affected by attenuation. Simulations show that the differential backscatter phase (δco) strongly depends on D0. At mm-wavelength, backscattering and propagation effects need to be disentangled first. To achieve this, an algorithm to detect and characterize Rayleigh plateaus is proposed and implemented. After the application of this algorithm, a methodology to estimate the differential backscatter phase and its error is given. The 95% confidence interval of δco is estimated with the re-sampling method bootstrapping.  Using simulation results, an attempt is made to find combinations of D0 and the raindrop size distribution shape parameter μ that match with the confidence interval of δco. The confidence interval of δco restricts D0, but not μ in most cases. This proposed technique is applied for both the 35 and 94 GHz frequency band of the new cloud radar at Cabauw (Ruisdael Observatory site near Utrecht). The resulting 95% confidence intervals of D0 with 35 and 94 GHz and their overlap are compared with in-situ disdrometer measurements of the mass-weighted mean diameter (Dm) which is closely related to D0. The median volume diameter retrieved with the 35 and 94 GHz frequency bands both shows a normalized cross correlation coefficient of 0.845 with the measured Dm of the disdrometer. Therefore, the cloud radar seems to have the capability to provide the detailed variations of the raindrops mean/median diameter like a local disdrometer, but at different heights. Nonetheless, the values differ. The disdrometer provides higher values than the cloud radar. One possible explanation is the inability of the disdrometer to measure raindrops smaller than 0.25 mm and the expected underestimation of the number of raindrops with sizes between 0.25 and 0.375 mm. However, because D0 values retrieved from 35 GHz data are also higher than the ones at 94 GHz, further research, which can use all the methodologies proposed in this master thesis work, is needed to examine the quantitative values of the median volume diameter retrieval.  These techniques can be implemented for all the single-frequency cloud radars (94 GHz) of the national Ruisdael Observatory (cloud and precipitation profiling mobile station and Lutjewad site above Groningen).

Wind is an important indicator of circulation related processes in the atmosphere. Accurate wind in- formation is used as an input for weather and circulation models. Wind data with a high temporal and spatial resolution are useful for research on the microphysics of the atmosphere, which is the area of application for the three-beam Transportable Atmospheric Radar (TARA) located at Cabauw Ex- perimental Site for Atmospheric Research (CESAR) in the Netherlands. Although some comparisons were performed using radiosondes, no thorough quality assessment of the wind estimation results has been done so far. During this research, the quality assessment is performed for the estimation of the wind speed and direction during the Analysis of the Composition of Clouds with Extended Polarization Techniques (ACCEPT) campaign, which took place in October and November 2014. This research shows that the wind estimation is working well during precipitation and a TARA-based criteria for the data selection to guarantee the quality of the wind retrievals is identified in the form of the coefficient of variation for the mean Doppler velocities. The quality assessment of the results of the TARA wind estimation is done by comparing the wind retrievals to the measurements of the meteorological tower at a height of 200 m above the surface at minute resolution. For the assessment to be performed, an effective way of removing clear air measurements is needed, which is done for the whole ACCEPT-campaign using a rough selection of time steps including rain. This approach showed that the wind estimation is working well during precipitation and shows that an improved data selection is needed. Several approaches towards improving the data selection for the wind estimation of TARA are per- formed. The most successful is the use of the coefficient of variation, which is defined as the standard deviation divided by the mean. This coefficient is calculated for the mean Doppler velocities, which are the input for the wind estimation algorithm. Comparing the results of the coefficient of variation method to the one based on rain selection shows that both return good results. This leads to the conclusion that the coefficient of variation is useful to improve the data selection.

Mixed-phase clouds, which have a significant impact on the global climate, are complex systems where liquid water and various types of ice particles coexist at temperatures below the freezing point. A key process in mixed-phase clouds is riming which alters microphysical and scattering properties of ice particles. Cloud radar is a powerful instrument for observing and understanding the processes that occur within mixed-phase clouds. Observations from multi-frequency radars and simulation results were combined in recent research to retrieve microphysical properties of ice particles in snowfall and ice clouds. This report presents an ambitious attempt to retrieve all common microphysical properties of ice particles, such as maximum dimension, density, aspect ratio and number concentration in slight rime condition using Doppler spectra. Two mixed-phase cloud events with low liquid water path are studied for such purpose. Spectral dual-wavelength ratio is introduced to retrieve maximum dimension of particles. An iteration process is developed in order to retrieve aspect ratio and density of ice particles from observation of spectral differential reflectivity. The number concentration of particles is retrieved with additional spectral reflectivity. With all the retrieved microphysical properties, ice water content and particle size distribution can be further derived. Ice water content is compared with results from an empirical model. The retrieved properties obtained from using three distinct mass-size relations are compared. Also the bulk and spectral retrieved profiles are compared. The retrieval process can provide consistent microphysical properties of ice particles. It is found that the retrieved ice water content is generally smaller than that from empirical model. Besides, the mass-size relation has significant impact on all retrieved microphysical properties except maximum dimension. The resulting profiles from bulk retrieval are smoother, while spectral retrieval can provide values in regions where the former cannot. The possible error from different sources are discussed or estimated, including the effect on dual-wavelength ratio from the elevation angle of radar, the neglect of differential attenuation caused by liquid and the usage of soft spheroid model. Recommendations are discussed, which include the usage of the latest microphysical models for ice aggregates and Discrete Dipole Approximation for electromagnetic wave scattering simulation.