We present a method to estimate spatially and temporally variable uncertainty of areal precipitation data. The aim of the method is to merge measurements from different sources, remote sensing and in ...situ, into a combined precipitation product and to provide an associated dynamic uncertainty estimate. This estimate should provide an accurate representation of uncertainty both in time and space, an adjustment to additional observations merged into the product through data assimilation, and flow dependency. Such a detailed uncertainty description is important for example to generate precipitation ensembles for probabilistic hydrological modelling or to specify accurate error covariances when using precipitation observations for data assimilation into numerical weather prediction models. The presented method uses the Local Ensemble Transform Kalman Filter and an ensemble nowcasting model. The model provides information about the precipitation displacement over time and is continuously updated by assimilation of observations. In this way, the precipitation product and its uncertainty estimate provided by the nowcasting ensemble evolve consistently in time and become flow-dependent. The method is evaluated in a proof of concept study focusing on weather radar data of four precipitation events. The study demonstrates that the dynamic areal uncertainty estimate outperforms a constant benchmark uncertainty value in all cases for one of the evaluated scores, and in half the number of cases for the other score. Thus, the flow dependency introduced by the coupling of data assimilation and nowcasting enables a more accurate spatial and temporal distribution of uncertainty. The mixed results achieved in the second score point out the importance of a good probabilistic nowcasting scheme for the performance of the method.
This publication intends to prove that a network of low-cost local area weather radars (LAWR) is a reliable and scientifically valuable complement to nationwide radar networks. A network of four ...LAWRs has been installed in northern Germany within the framework of the Precipitation and Attenuation Estimates from a High-Resolution Weather Radar Network (PATTERN) project observing precipitation with a temporal resolution of 30 s, a range resolution of 60 m and a sampling resolution of 1° in the azimuthal direction. The network covers an area of 60 km × 80 km. In this paper, algorithms used to obtain undisturbed precipitation fields from raw reflectivity data are described, and their performance is analysed. In order to correct operationally for background noise in reflectivity measurements, noise level estimates from the measured reflectivity field are combined with noise levels from the last 10 time steps. For detection of non-meteorological echoes, two different kinds of clutter algorithms are applied: single-radar algorithms and network-based algorithms. Besides well-established algorithms based on the texture of the logarithmic reflectivity field (TDBZ) or sign changes in the reflectivity gradient (SPIN), the advantage of the unique features of the high temporal and spatial resolution of the network is used for clutter detection. Overall, the network-based clutter algorithm works best with a detection rate of up to 70%, followed by the classic TDBZ filter using the texture of the logarithmic reflectivity field. A comparison of a reflectivity field from the PATTERN network with the product from a C-band radar operated by the German Meteorological Service indicates high spatial accordance of both systems in the geographical position of the rain event as well as reflectivity maxima. Long-term statistics from May to September 2013 prove very good accordance of the X-band radar of the network with C-band radar, but, especially at the border of precipitation events, higher-resolved X-band radar measurements provide more detailed information on precipitation structure because the 1 km range gate of C-band radars is only partially covered with rain. The standard deviation within a range gate of the C-band radar with a range resolution of 1 km is up to 3 dBZ at the borders of rain events. The probability of detection is at least 90%, the false alarm ratio less than 10% for both systems. Therefore, a network of high-resolution low-cost LAWRs can give valuable information on the small-scale structure of rain events in areas of special interest, e.g. urban regions, in addition to the nationwide radar networks.
An advanced package of microwave remote sensing instrumentation has been developed for the operation on the new German High Altitude LOng range research aircraft (HALO). The HALO Microwave Package, ...HAMP, consists of two nadir-looking instruments: a cloud radar at 36 GHz and a suite of passive microwave radiometers with 26 frequencies in different bands between 22.24 and 183.31 ± 12.5 GHz. We present a description of HAMP's instrumentation together with an illustration of its potential. To demonstrate this potential, synthetic measurements for the implemented passive microwave frequencies and the cloud radar based on cloud-resolving and radiative transfer model calculations were performed. These illustrate the advantage of HAMP's chosen frequency coverage, which allows for improved detection of hydrometeors both via the emission and scattering of radiation. Regression algorithms compare HAMP retrieval with standard satellite instruments from polar orbiters and show its advantages particularly for the lower atmosphere with a root-mean-square error reduced by 5 and 15% for temperature and humidity, respectively. HAMP's main advantage is the high spatial resolution of about 1 km, which is illustrated by first measurements from test flights. Together these qualities make it an exciting tool for gaining a better understanding of cloud processes, testing retrieval algorithms, defining future satellite instrument specifications, and validating platforms after they have been placed in orbit.
The estimation of CO2 exchange between the ocean and the atmosphere is essential to understand the global carbon cycle. The eddy-covariance technique offers a very direct approach to observe these ...fluxes. The turbulent CO2 flux is measured, as well as the sensible and latent heat flux and the momentum flux, a few meters above the ocean in the atmosphere. Assuming a constant-flux layer in the near-surface part of the atmospheric boundary layer, this flux equals the exchange flux between ocean and atmosphere. The purpose of this paper is the comparison of long-term flux measurements at two different heights above the Baltic Sea to investigate this assumption. The results are based on a 1.5-year record of quality-controlled eddy-covariance measurements. Concerning the flux of momentum and of sensible and latent heat, the constant-flux layer theory can be confirmed because flux differences between the two heights are insignificantly small more than 95 % of the time. In contrast, significant differences, which are larger than the measurement error, occur in the CO2 flux about 35 % of the time. Data used for this paper are published at http://doi.pangaea.de/10.1594/PANGAEA.808714 .
The theoretical framework of a novel approach for absolute radar calibration is presented and its potential analysed by means of synthetic data to lay out a solid basis for future practical ...application. The method presents the advantage of an absolute calibration with respect to the directly measured reflectivity, without needing a previously calibrated reference device. It requires a setup comprising three radars: two devices oriented towards each other, measuring reflectivity along the same horizontal beam and operating within a strongly attenuated frequency range (e.g. K or X band), and one vertical reflectivity and drop size distribution (DSD) profiler below this connecting line, which is to be calibrated. The absolute determination of the calibration factor is based on attenuation estimates. Using synthetic, smooth and geometrically idealised data, calibration is found to perform best using homogeneous precipitation events with rain rates high enough to ensure a distinct attenuation signal (reflectivity above ca. 30 dBZ). Furthermore, the choice of the interval width (in measuring range gates) around the vertically pointing radar, needed for attenuation estimation, is found to have an impact on the calibration results. Further analysis is done by means of synthetic data with realistic, inhomogeneous precipitation fields taken from measurements. A calibration factor is calculated for each considered case using the presented method. Based on the distribution of the calculated calibration factors, the most probable value is determined by estimating the mode of a fitted shifted logarithmic normal distribution function. After filtering the data set with respect to rain rate and inhomogeneity and choosing an appropriate length of the considered attenuation path, the estimated uncertainty of the calibration factor is of the order of 1 to 11 %, depending on the chosen interval width. Considering stability and accuracy of the method, an interval of eight range gates on both sides of the vertically pointing radar is most appropriate for calibration in the presented setup.
The influence of spatial surface temperature changes over the Arctic Ocean on the 2-m air temperature variability is estimated using backward trajectories based on ERA-Interim and JRA25 wind fields. ...They are initiated at Alert, Barrow and at the Tara drifting station. Three different methods are used. The first one compares mean ice surface temperatures along the trajectories to the observed 2-m air temperatures at the stations. The second one correlates the observed temperatures to air temperatures obtained using a simple Lagrangian box model that only includes the effect of sensible heat fluxes. For the third method, mean sensible heat fluxes from the model are correlated with the difference of the air temperatures at the model starting point and the observed temperatures at the stations. The calculations are based on MODIS ice surface temperatures and four different sets of ice concentration derived from SSM/I (Special Sensor Microwave Imager) and AMSR-E (Advanced Microwave Scanning Radiometer for EOS) data. Under nearly cloud-free conditions, up to 90% of the 2-m air temperature variance can be explained for Alert, and 70% for Barrow, using these methods. The differences are attributed to the different ice conditions, which are characterized by high ice concentration around Alert and lower ice concentration near Barrow. These results are robust for the different sets of reanalyses and ice concentration data. Trajectories based on 10-m wind fields from both reanalyses show large spatial differences in the Central Arctic, leading to differences in the correlations between modeled and observed 2-m air temperatures. They are most pronounced at Tara, where explained variances amount to 70% using JRA and 80% using ERA. The results also suggest that near-surface temperatures at a given site are influenced by the variability of surface temperatures in a domain of about 200 km radius around the site.
The earth’s surface is characterized by small-scale heterogeneity attributable to variability in land cover, soil characteristics and orography. In atmospheric models, this small-scale variability ...can be partially accounted for by the so-called mosaic approach, i.e., by computing the land-surface processes on a grid with an explicit higher horizontal resolution than the atmosphere. The mosaic approach does, however, not account for the subgrid-scale variability in the screen-level atmospheric parameters, part of which might be related to land-surface heterogeneity itself. In this study, simulations with the numerical weather prediction model COSMO are shown, employing the mosaic approach together with a spatial disaggregation of the atmospheric forcing by the screen-level variables to the subgrid-scale. The atmospheric model is run with a 2.8 km horizontal grid resolution while the land surface processes are computed on a 400-m horizontal grid. The disaggregation of the driving atmospheric variables at screen-level is achieved by a three-step statistical downscaling with rules learnt from high-resolution fully coupled COSMO simulations, where both, atmosphere and surface, were simulated on a 400-m grid. The steps encompass spline interpolation of the grid scale variables, conditional regression based on the high-resolution runs, and an optional stochastic noise generator which restores the variability of the downscaled variables. Simulations for a number of case studies have been carried out, with or without mosaic surface representation and with or without atmospheric disaggregation, and evaluated with respect to the surface state variables and the turbulent surface exchange fluxes of sensible and latent heat. The results are compared with the high-resolution fully coupled COSMO simulations. The results clearly demonstrate the high importance of accounting for subgrid-scale surface heterogeneity. It is shown that the atmospheric disaggregation leads to clear additional improvements in the structures of the two-dimensional surface state variable fields, but to only marginally impacts on the simulation of the turbulent surface exchange fluxes. A detailed analysis of these results identifies strongly correlated errors in atmospheric and surface variables in the mosaic approach as the main reason for the latter. The effects of these errors largely cancel out in the flux parameterization, and thus explain the comparably good results for the fluxes in the mosaic approach without atmospheric disaggregation despite inferior performance for the surface state variables themselves. Inserting noise in the disaggregation scheme leads to a deterioration of the results.
A stochastic version of the Iterative Amplitude Adjusted Fourier Transform (IAAFT) algorithm is presented. This algorithm is able to generate so-called surrogate time series, which have the amplitude ...distribution and the power spectrum of measured time series or fields. The key difference between the new algorithm and the original IAAFT method is the treatment of the amplitude adjustment: it is not performed for all values in each iterative step, but only for a fraction of the values. This new algorithm achieves a better accuracy, i.e. the power spectra of the measurement and its surrogate are more similar. We demonstrate the improvement by applying the IAAFT algorithm and the new one to 13 different test signals ranging from rain time series and 3-dimensional clouds to fractal time series and theoretical input. The improved accuracy can be important for generating high-quality geophysical time series and fields. The traditional application of the IAAFT algorithm is statistical nonlinearity testing. Reassuringly, we found that in most cases the accuracy of the original IAAFT algorithm is sufficient for this application.
For driving soil-vegetation-transfer models or hydrological models, high-resolution atmospheric forcing data is needed. For most applications the resolution of atmospheric model output is too coarse. ...To avoid biases due to the non-linear processes, a downscaling system should predict the unresolved variability of the atmospheric forcing. For this purpose we derived a disaggregation system consisting of three steps: (1) a bi-quadratic spline-interpolation of the low-resolution data, (2) a so-called 'deterministic' part, based on statistical rules between high-resolution surface variables and the desired atmospheric near-surface variables and (3) an autoregressive noise-generation step. The disaggregation system has been developed and tested based on high-resolution model output (400mhorizontal grid spacing).Anovel automatic search-algorithm has been developed for deriving the deterministic downscaling rules of step 2. When applied to the atmospheric variables of the lowest layer of the atmospheric COSMO-model, the disaggregation is able to adequately reconstruct the reference fields. Applying downscaling step 1 and 2, root mean square errors are decreased. Step 3 finally leads to a close match of the subgrid variability and temporal autocorrelation with the reference fields. The scheme can be applied to the output of atmospheric models, both for stand-alone offline simulations, and a fully coupled model system.