Search

The roadmap to enable operational soil moisture (SM) monitoring for meteorologic and hydrological early warning is challenged by the uncertainty within the available remote sensing and modelling products. This study addressed two relevant uncertainties: the residual trends in the series and the spatial consistency. While the latter has been often revisited to validate remote sensing and modelling products against in-situ data, the former is often disregarded in studies addressing SM changes. This study evaluated three SM products: (1) the Satellite Application Facility on Support to Operational Hydrology and Water Management (H SAF) active Advanced SCATterometer (ASCAT)-derived dataset, (2) the passive subset of the European Space Agency (ESA) - Climate Change Initiative (CCIp), and (3) the modelled dataset from the European Drought Observatory (EDO). The analysis was carried out over Europe in the period 2007–2022 at 10-day temporal scales. We obtained that even these popular datasets are subject to patches of spatial inconsistency and residual trends when compared to the in-situ data from the International Soil Moisture Network (ISMN). In view of the great complementarity shown by the active and passive remote sensing and the modelled SM estimates, two merged products are proposed and tested against in-situ data. Results indicate that combining H SAF ASCAT, CCIp and EDO equals or surpasses the spatial and temporal consistency of the individual SM products alone, even when only the near-real-time products of H SAF ASCAT and EDO are combined. Thus, merging remote sensing and modelled SM products is advantageous for enhanced spatial and temporal operational monitoring of SM at the European scale. [ABSTRACT FROM AUTHOR]

The impact of assimilating GNSS-ZTD (global navigation satellite system–zenith total delay) on the precipitable water vapor and precipitation forecast over Italy is studied for the month of October 2019, which was characterized by several moderate to intense precipitation events, especially over northwestern Italy. The WRF (Weather Research and Forecasting) model, version 4.1.3, is used with its 3D-Var data assimilation system to assimilate ZTD observations from 388 GNSS receivers distributed over the country. The dataset was built collecting data from all the major national and regional GNSS permanent networks, achieving dense coverage over the whole area. The water vapor forecast is verified for the forecast hours of 1–6 h after the last data assimilation time. Results show that WRF underestimates the atmospheric water vapor content for the period, and GNSS-ZTD data assimilation improves this underestimation. The precipitation forecast is verified in the phases of 0–3 and 3–6 h after the last data assimilation time using more than 3000 rain gauges spread over Italy. The application of GNSS-ZTD data assimilation to a case study improved the precipitation forecast by increasing the rainfall maximum and by better focusing the precipitation pattern over northeastern Italy, with the main drawback being the prediction of false alarms. Considering the study over the whole period, GNSS-ZTD data assimilation had a positive impact on rainfall forecast, with an improvement in the performance up to 6 h and with statistically significant results for moderate to intense rainfall thresholds (25–30 mm (3 h) -1). [ABSTRACT FROM AUTHOR]

The use of SAR data for flood mapping is well established. However, the problem of the missed detection of rapidly evolving floods remains. Recently, the Italian Space Agency deployed the COSMO-SkyMed Second Generation (CSG) constellation, with an on-demand capability that makes it possible to reduce the number of missed floods. However, for on-demand SAR acquisitions, pre-flood images are generally not available, and change-detection methods, commonly adopted for flood mapping using SAR, cannot be applied. This study focused on the high-resolution CSG images of a flood that occurred in central Italy. An accurate analysis of the radar responses of the different targets included in the scenes observed by GSG was performed. Then, a methodology to detect floods using high-resolution single SAR images was developed. The methodology was based on image segmentation and fuzzy logic. Image segmentation allowed us to analyze homogeneous areas in the CSG images. Fuzzy logic was used to integrate the SAR data with ancillary information that was useful when change-detection methods could not be applied. A comparison with the maps produced by the Copernicus Emergency Service, using high-resolution optical images, demonstrated the reliability of the methodology. [ABSTRACT FROM AUTHOR]

The impact of assimilating GNSS-ZTD (Global Navigation Satellite Systems - Zenith Total Delay) on the precipitable water vapor and precipitation forecast over Italy is studied for the month of October 2019, characterized by several moderate to intense precipitation events. The WRF (Weather Research and Forecasting) model 4.1.3 is used with its 3DVar data assimilation system to assimilate ZTD observations from 388 GNSS receivers distributed over the country. The dataset was built collecting data from all the major national and regional GNSS permanent networks, achieving dense coverage over the whole area. Results show that WRF underestimates the atmospheric water vapor content for the period, and GNSS-ZTD data assimilation improves this underestimation. The water vapor forecast is verified for the forecast hours 1-6h after the last data assimilation time, while the precipitation forecast is verified in the phases 0-3h and 3-6h after the last data assimilation time. More than 3000 rain gauges spread over Italy were used to verify the precipitation forecast. The application of GNSS-ZTD data assimilation to a case study showed an improvement of the precipitation forecast in different ways, the main drawback being the prediction of false alarms. Considering the study over the whole period, GNSS-ZTD data assimilation had a positive impact on the precipitable water vapor and rainfall forecast, with an improvement of the performance up to 6 hours, and with statistically significant results for moderate to intense rainfall thresholds (25-30 mm/3h). [ABSTRACT FROM AUTHOR]

We present IT-SNOW, a serially complete and multi-year snow reanalysis for Italy (∼ 301 × 10 3 km2) – a transitional continental-to-Mediterranean region where snow plays an important but still poorly constrained societal and ecological role. IT-SNOW provides ∼ 500 m daily maps of snow water equivalent (SWE), snow depth, bulk snow density, and liquid water content for the initial period 1 September 2010–31 August 2021, with future updates envisaged on a regular basis. As the output of an operational chain employed in real-world civil protection applications (S3M Italy), IT-SNOW ingests input data from thousands of automatic weather stations, snow-covered-area maps from Sentinel-2, MODIS (Moderate Resolution Imaging Spectroradiometer), and H SAF products, as well as maps of snow depth from the spatialization of over 350 on-the-ground snow depth sensors. Validation using Sentinel-1-based maps of snow depth and a variety of independent, in situ snow data from three focus regions (Aosta Valley, Lombardy, and Molise) show little to no mean bias compared to the former, and root mean square errors are of the typical order of 30–60 cm and 90–300 mm for in situ, measured snow depth and snow water equivalent, respectively. Estimates of peak SWE by IT-SNOW are also well correlated with annual streamflow at the closure section of 102 basins across Italy (0.87), with ratios between peak water volume in snow and annual streamflow that are in line with expectations for this mixed rain–snow region (22 % on average and 12 % median). Examples of use allowed us to estimate 13.70 ± 4.9 Gm3 of water volume stored in snow across the Italian landscape at peak accumulation, which on average occurs on 4 March ± 10 d. Nearly 52 % of the mean seasonal SWE is accumulated across the Po river basin, followed by the Adige river (23 %), and central Apennines (5 %). IT-SNOW is freely available at 10.5281/zenodo.7034956 and can contribute to better constraining the role of snow for seasonal to annual water resources – a crucial endeavor in a warming and drier climate. [ABSTRACT FROM AUTHOR]

A comprehensive analysis of the July 2021 event that occurred on Lake Como (Italy), during which heavy hailstorms and floods affected the surroundings of Lake, is presented. The study provides a detailed analysis of the event using different observation sources currently available. The employed techniques include both conventional (rain gauges, radar, atmospheric sounding) and non-conventional (satellite-based Earth observation products, GNSS, and lightning detection network) observations for hydro-meteorological analysis. The study is split in three main topics: event description by satellite-based observations; long-term analysis by the ERA5 model and ASCAT soil water index; and short-term analysis by lightning data, GNSS delays and radar-VIL. The added value of the work is the near-real-time analysis of some of the datasets used, which opens up the potential for use in alerting systems, showing considerable application possibilities in NWP modeling, where it can also be useful for the implementation of early warning systems. The results highlight the validity of the different techniques and the consistency among the observations. This result, therefore, leads to the conclusion that a joint use of the innovative techniques with the operational ones can bring reliability in the description of events. [ABSTRACT FROM AUTHOR]

We present IT-SNOW, a serially complete and multi-year snow reanalysis for Italy (300k+ km²) covering a transitional continental-to-Mediterranean region where snow plays an important, but still poorly constrained societal and ecological role. IT-SNOW provides ∼500-m, daily maps of Snow Water Equivalent (SWE), snow depth, bulk-snow density, and liquid water content for the 5 period 01/09/2010 - 31/08/2021, with future updates envisaged on a regular basis. As the output of an operational chain employed in real-world civil-protection applications (S3M Italy), IT-SNOW ingests input data from thousands of automatic weather stations, snow-covered-area maps from Sentinel 2, MODIS, and H-SAF products, and maps of snow depth from the spazialization of 350+ on-the-ground snow-depth sensors. Validation using Sentinel-1-based maps of snow depth and a variety of independent, in-situ snow data from three focus regions (Aosta Valley, Lombardia, and Molise) shows little to none mean bias compared to the former, and Root Mean Square Errors on the order of 30 to 60 cm and 90 to 300 mm for in-situ, measured snow depth and Snow Water Equivalent, respectively. Estimates of peak SWE by IT-SNOW are also well correlated with annual streamflow at the closure section of 102 basins across Italy (0.87), with ratios between peak SWE and annual streamflow that are in line with expectations for this mixed rain-snow region (22% on average). Examples of use allowed us to estimate 13.70 ± 4.9 Gm3 of SWE across the Italian landscape at peak accumulation, which on average occurs on the 4th of March. Nearly 52% of mean seasonal SWE is accumulated across the Po river basin, followed by the Adige river (23%), and central Apennines (5%). IT-SNOW is freely available with the following DOI: https://doi.org/10.5281/zenodo.7034956 (Avanzi et al., 2022b) and can contribute to better constraining the role of snow for seasonal to annual water resources - a crucial endevor in a warming and drier climate. [ABSTRACT FROM AUTHOR]

Lightning is an important threat to life and properties and its forecast is important for several applications. In this paper, we show the performance of the "dynamic lightning scheme" for next-day total strokes forecast. The predictions were compared against strokes recorded by a ground observational network for a forecast period spanning one year. Specifically, a total of 162 case studies were selected between 1 March 2020 and 28 February 2021, characterized by at least 3000 observed strokes over Italy. The events span a broad range of lightning intensity from about 3000 to 600,000 strokes in one day: 69 cases occurred in summer, 46 in fall, 18 in winter, and 29 in spring. The meteorological driver was the Weather Research and Forecasting (WRF) model (version 4.1) and we focused on the next-day forecast. Strokes were simulated by adding three extra variables to WRF, namely, the potential energies for positive and negative cloud to ground flashes and intracloud strokes. Each potential energy is advected by WRF, it is built by the electrification processes occurring into the cloud, and it is dissipated by lightning. Observed strokes were remapped onto the WRF model grid with a 3 km horizontal resolution for comparison with the strokes forecast. Results are discussed for the whole year and for different seasons. Moreover, statistics are presented for the land and the sea. In general, the results of this study show that lightning forecast with the dynamic lightning scheme and WRF model was successful for Italy; nevertheless, a careful inspection of forecast performance is necessary for tuning the scheme. This tuning is dependent on the season. A numerical experiment changing the microphysics scheme used in WRF shows the sensitivity of the results according to the choice of the microphysics scheme. [ABSTRACT FROM AUTHOR]

Because of the ongoing climate change, the frequency of extreme rainfall events at the global scale is expected to increase, resulting in higher social and economic impacts. Thus, improving the forecast accuracy and the risk communication is a fundamental goal to limit social and economic damages. Both Numerical Weather Prediction (NWP) and radar-based nowcasting systems still have open issues, mainly in terms of precipitation correct time/space localization predictability and rapid forecast accuracy decay, respectively. Trying to overcome these issues, this work aims to present a nowcasting system combining an NWP model (WRF), using a 3 h rapid update cycling 3DVAR assimilation of radar reflectivity data, with the radar-based nowcasting system PhaSt through a blending technique. Moreover, an innovative post-processing algorithm named SWING (Score-Weighted Improved NowcastinG) has been developed in order to take into account the timely and spatial uncertainty in the convective field simulation. The overarching goal is to pave the way for an easy and automatic communication of the heavy rainfall warning derived by the nowcasting procedure. The results obtained applying the SWING algorithm over a case study of 22 days in the fall 2019 season suggest that the algorithm could improve the predictive capability of a traditional deterministic nowcasting forecast system, keeping a useful forecast timing and thus integrating the current forecast procedures. Eventually, the main advantage of the SWING algorithm is also its very high versatility, since it could be used with any meteorological model also in a multi-model forecast approach. [ABSTRACT FROM AUTHOR]

A Multivariate Neural Network (MNN) algorithm is proposed for the reconstruction of significant wave height time series, without any increase of the error of the MNN output with the number of modelled data. The algorithm uses a weighted error function during the learning phase, to improve the modelling of the higher significant wave height. The ability of the MNN to reconstruct sea storms is tested by applying the equivalent triangular storm model. Finally an application to the NOAA buoys moored off California shows a good performance of the MNN algorithm, both during sea storms and calm time periods.

Heavy and localized summer events are very hard to predict and, at the same time, potentially dangerous for people and properties. This paper focuses on an event occurred on 15 July 2020 in Palermo, the largest city of Sicily, causing about 120 mm of rainfall in 3 h. The aim is to investigate the event predictability and a potential way to improve the precipitation forecast. To reach this aim, lightning (LDA) and radar reflectivity data assimilation (RDA) was applied. LDA was able to trigger deep convection over Palermo, with high precision, whereas the RDA had a key role in the prediction of the amount of rainfall. The simultaneous assimilation of both data sources gave the best results. An alert for a moderate–intense forecast could have been issued one hour and a half before the storm developed over the city, even if predicting only half of the total rainfall. A satisfactory prediction of the amount of rainfall could have been issued at 14:30 UTC, when precipitation was already affecting the city. Although the study is centered on a single event, it highlights the need for rapidly updated forecast cycles with data assimilation at the local scale, for a better prediction of similar events. [ABSTRACT FROM AUTHOR]

Many severe precipitation events occurred in Europe during the last decade and caused casualties and damages to the historical heritage and natural environment. Protection of civilians and safeguard of the territory are the kay mandates of a civil Protection (CP) agency. This is achieved through activities that mitigate hydro-meteorological risks, such as flooding and droughts. European CP agencies, supported by meteorological and hydrological institutes, aim to assess risk scenarios, to monitor and supervise events and risk levels, providing early warning to National and local authorities.Hydro meteorological risk management consists of three phases: pre-event (forecast), event (early warning, and monitoring) and post-event (recovery and survey).Near real time accurate estimations of hydrological variables such as precipitation and soil moisture are invaluable to the CP agencies, enable them to issue early warnings and plan for disaster relief at the local level.Besides measurements of key hydrological variables by ground-based instruments, often affected by a limited spatial coverage, advanced satellite-based precipitation and soil moisture products developed within different international programs, are becoming available and accessible to users in near-real time. The assessment of the accuracy and reliability of such products is necessary in order to be able to optimally exploit them for hydro-meteorological applications.The comparison of satellite-derived rainfall and soil moisture estimates with respect to ground-based measurements is a challenging task, both because of the temporal and spatial variability of the fields and the inherently different measurement types. This activity is particularly complex in Europe due to the heterogeneity in the ground-based networks of the different countries, due to sampling strategy and processing methodologies. For the validation and assessment of satellite-derived soil moisture products, these issues are exacerbated by the limited availability of ground observation networks that can be used as benchmark. Moreover, the measuring stations cover only small areas over Europe. As an example, in Central Italy only 14 probes are distributed over an area of about 8000 km2.In this study, a comparison between satellite derived products and ground data is presented, according to the EUMETSAT H SAF (Satellite Application Facility on Support to Operational Hydrology and Water Management) Validation protocols. Since 2005, eleven European countries joined in the framework of the H SAF to perform the satellite products monitoring and ground verification in their own regions of interest. Moreover, recent severe rainfall events are selected in order to understand how the main satellite product characteristics, i.e. accuracy, spatial pattern and resolution, update frequency and latency, impact the efficiency of a hydro-meteorological early warning system at a local level in an operational framework. State-of-the-art satellite rainfall and soil moisture products, obtained through international space agencies and programmes (EUMETSAT H SAF, ESA, NASA, and JAXA), are used over different European countries as forcing input for an early warning system for meteorological/hydrological events. [ABSTRACT FROM AUTHOR]

An automated tool for pre-operational mapping of floods and inland waters using Sentinel-1 data is presented. The acronym AUTOWADE (AUTOmatic Water Areas DEtector) is used to denote it. The tool provides the end user (Italian Department of Civil Protection) with a continuous, near real-time (NRT) monitoring of the extent of inland water surfaces (floodwater and permanent water). It implements the following operations: downloading of Sentinel-1 products; preprocessing of the products and storage of the resulting geocoded and calibrated data; generation of the intermediate products, such as the exclusion mask; application of a floodwater/permanent water mapping algorithm; generation of the output layer, i.e., a map of floodwater/permanent water; delivery of the output layer to the end user. The open floodwater/permanent water mapping algorithm implemented in AUTOWADE is based on a new approach, denoted as buffer-from-edge (BFE), which combines different techniques, such as clustering, edge filtering, automatic thresholding and region growing. AUTOWADE copes also with the typical presence of gaps in the flood maps caused by undetected flooded vegetation. An attempt to partially fill these gaps by analyzing vegetated areas adjacent to open water is performed by another algorithm implemented in the tool, based on the fuzzy logic. The BFE approach has been validated offline using maps produced by the Copernicus Emergency Management Service. Validation has given good results with a F1-score larger than 0.87 and a kappa coefficient larger than 0.80. The algorithm to detect flooded vegetation has been visually compared with optical data and aerial photos; its capability to fill some of the gaps present in flood maps has been confirmed. [ABSTRACT FROM AUTHOR]

In this paper, precipitation estimates derived from the Italian ground radar network (IT GR) are used in conjunction with Spinning Enhanced Visible and InfraRed Imager (SEVIRI) measurements to develop an operational oriented algorithm (RAdar INfrared Blending algorithm for Operational Weather monitoring (RAINBOW)) able to provide precipitation pattern and intensity. The algorithm evaluates surface precipitation over five geographical boxes (in which the study area is divided). It is composed of two main modules that exploit a second-degree polynomial relationship between the SEVIRI brightness temperature at 10.8 µm TB10.8 and the precipitation rate estimates from IT GR. These relationships are applied to each acquisition of SEVIRI in order to provide a surface precipitation map. The results, based on a number of case studies, show good performance of RAINBOW when it is compared with ground reference (precipitation rate map from interpolated rain gauge measurements), with high Probability of Detection (POD) and low False Alarm Ratio (FAR) values, especially for light to moderate precipitation range. At the same time, the mean error (ME) values are about 0 mmh−1, while root mean square error (RMSE) is about 2 mmh−1, highlighting a limited variability of the RAINBOW estimations. The precipitation retrievals from RAINBOW have been also compared with the European Organization for the Exploitation of Meteorological Satellites (EUMETSAT) Satellite Application Facility on Support to Operational Hydrology and Water Management (H SAF) official microwave (MW)/infrared (IR) combined product (P-IN-SEVIRI). RAINBOW shows better performances than P-IN-SEVIRI, in terms of both detection and estimates of precipitation fields when they are compared to the ground reference. RAINBOW has been designed as an operational product, to provide complementary information to that of the national radar network where the IT GR coverage is absent, or the quality (expressed in terms of Quality Index (QI)) of the RAINBOW estimates is low. The aim of RAINBOW is to complement the radar and rain gauge network supporting the operational precipitation monitoring. [ABSTRACT FROM AUTHOR]

On 10 October 2018 an intense storm, characterized by heavy rainfall, hit the Sardinia island, reaching a peak of 452 mm of rain measured in 24 h. Among others, two particularly intense phases were registered between 3 and 6 UTC (Universal Coordinated Time), and between 18 and 24 UTC. The forecast of this case study is challenging because the precipitation was heavy and localized. In particular, the meteorological model used in this paper, provides a good prediction only for the second period over the eastern part of the Sardinia island. In this work, we study the impact of lightning data assimilation and horizontal grid resolution on the Very Short-term Forecast (VSF, 3 and 1 h) for this challenging case, using the RAMS@ISAC meteorological model. The comparison between the 3 h VSF control run and the simulations with lightning data assimilation shows the considerable improvement given by lightning data assimilation, especially for the precipitation that occurred in the eastern part of the island. Reducing the VSF range to 1 h, resulted in higher model performance with a good precipitation prediction over eastern and south-central Sardinia. In addition, the comparison between simulated and observed reflectivity shows an important improvement of simulations with lightning data assimilation compared to the control forecast. However, simulations assimilating lightning overestimated the precipitation in the last part of the day. The increasing of the horizontal resolution to 2 km grid spacing reduces the false alarms and improves the model performance. [ABSTRACT FROM AUTHOR]

A fully automated processing chain for near real-time mapping of burned forest areas using Sentinel-2 multispectral data is presented. The acronym AUTOBAM (AUTOmatic Burned Areas Mapper) is used to denote it. AUTOBAM is conceived to work daily at a national scale for the Italian territory to support the Italian Civil Protection Department in the management of one of the major natural hazards, which affects the territory. The processing chain includes a Sentinel-2 data procurement component, an image processing algorithm, and the delivery of the map to the end-user. The data procurement component searches every day for the most updated products into different archives. The image processing part represents the core of AUTOBAM and implements an algorithm for burned forest areas mapping that uses, as fundamental parameters, the relativized form of the delta normalized burn ratio and the normalized difference vegetation index. The minimum mapping unit is 1 ha. The algorithm implemented in the image processing block is validated off-line using maps of burned areas produced by the Copernicus Emergency Management Service. The results of the validation shows an overall accuracy (considering the classes of burned and unburned areas) larger than 95% and a kappa coefficient larger than 80%. For what concerns the class of burned areas, the commission error is around 1%−3%, except for one case where it reaches 25%, while the omission error ranges between 6% and 25%. [ABSTRACT FROM AUTHOR]

Information on snow properties is of critical relevance for a wide range of scientific studies and operational applications, mainly for hydrological purposes. However, the ground-based monitoring of snow dynamics is a challenging task, especially over complex topography and under harsh environmental conditions. Remote sensing is a powerful resource providing snow observations at a large scale. This study addresses the potential of using Sentinel-2 high-resolution imagery to assess moderate-resolution snow products, namely H10—Snow detection (SN-OBS-1) and H12—Effective snow cover (SN-OBS-3) supplied by the Satellite Application Facility on Support to Operational Hydrology and Water Management (H-SAF) project of the European Organisation for the Exploitation of Meteorological Satellites (EUMETSAT). With the aim of investigating the reliability of reference data, the consistency of Sentinel-2 observations is evaluated against both in-situ snow measurements and webcam digital imagery. The study area encompasses three different regions, located in Finland, the Italian Alps and Turkey, to comprehensively analyze the selected satellite products over both mountainous and flat areas having different snow seasonality. The results over the winter seasons 2016/17 and 2017/18 show a satisfying agreement between Sentinel-2 data and ground-based observations, both in terms of snow extent and fractional snow cover. H-SAF products prove to be consistent with the high-resolution imagery, especially over flat areas. Indeed, while vegetation only slightly affects the detection of snow cover, the complex topography more strongly impacts product performances. [ABSTRACT FROM AUTHOR]

The uncertainties associated with rainfall estimates comprise various measurement scales: from rain gauges and ground-based radars to the satellite rainfall retrievals. The quality of satellite rainfall products has improved significantly in recent decades; however, such algorithms require validation studies using observational rainfall data. For this reason, this study aims to apply the H-SAF consolidated radar data processing to the X-band radar used in the CHUVA campaigns and apply the well established H-SAF validation procedure to these data and verify the quality of EUMETSAT H-SAF operational passive microwave precipitation products in two regions of Brazil (Vale do Paraíba and Manaus). These products are based on two rainfall retrieval algorithms: the physically based Bayesian Cloud Dynamics and Radiation Database (CDRD algorithm) for SSMI/S sensors and the Passive microwave Neural network Precipitation Retrieval algorithm (PNPR) for cross-track scanning radiometers (AMSU-A/AMSU-B/MHS sensors) and for the ATMS sensor. These algorithms, optimized for Europe, Africa and the Southern Atlantic region, provide estimates for the MSG full disk area. Firstly, the radar data was treated with an overall quality index which includes corrections for different error sources like ground clutter, range distance, rain-induced attenuation, among others. Different polarimetric and non-polarimetric QPE algorithms have been tested and the Vulpiani algorithm (hereafter, R q 2 V u 15 ) presents the best precipitation retrievals when compared with independent rain gauges. Regarding the results from satellite-based algorithms, generally, all rainfall retrievals tend to detect a larger precipitation area than the ground-based radar and overestimate intense rain rates for the Manaus region. Such behavior is related to the fact that the environmental and meteorological conditions of the Amazon region are not well represented in the algorithms. Differently, for the Vale do Paraíba region, the precipitation patterns were well detected and the estimates are in accordance with the reference as indicated by the low mean bias values. [ABSTRACT FROM AUTHOR]