Ensemble prediction for nowcasting with a convection-permitting model—I: description of the system and the impact of radar-derived surface precipitation rates

A key strategy to improve the skill of quantitative predictions of precipitation, as well as hazardous weather such as severe thunderstorms and flash floods is to exploit the use of observations of convective activity (e.g. from radar). In this paper, a convection-permitting ensemble prediction system (EPS) aimed at addressing the problems of forecasting localized weather events with relatively short predictability time scale and based on a 1.5 km grid-length version of the Met Office Unified Model is presented. Particular attention is given to the impact of using predicted observations of radar-derived precipitation intensity in the ensemble transform Kalman filter (ETKF) used within the EPS. Our initial results based on the use of a 24-member ensemble of forecasts for two summer case studies show that the convective-scale EPS produces fairly reliable forecasts of temperature, horizontal winds and relative humidity at 1 h lead time, as evident from the inspection of rank histograms. On the other hand, the rank histograms seem also to show that the EPS generates too much spread for forecasts of (i) surface pressure and (ii) surface precipitation intensity. These may indicate that for (i) the value of surface pressure observation error standard deviation used to generate surface pressure rank histograms is too large and for (ii) may be the result of non-Gaussian precipitation observation errors. However, further investigations are needed to better understand these findings. Finally, the inclusion of predicted observations of precipitation from radar in the 24-member EPS considered in this paper does not seem to improve the 1-h lead time forecast skill.

Suspended and bed load sediment transport dynamics in two lowland UK streams—storm integrated monitoring

A great deal of work recently has focused on suspended and bedload sediment transport, driven primarily by interest in contaminant transfer. However, uncertainties regarding the role of storm events, macrophyte beds and interactions between the two phases of sediment still exist. This paper compares two study sites within the same catchment whose geology varies significantly. The differences in hydrology, suspended sediment (SS) transport and bed load transport that this causes are examined. In addition, a method to predict the mobilization of different size fractions of sediment during given flows is investigated using critical entrainment thresholds.

Emission: a simple new technique to correct rainfall estimates from attenuation due to both the radome and heavy rainfall

We present a new technique for correcting errors in radar estimates of rainfall due to attenuation which is based on the fact that any attenuating target will itself emit, and that this emission can be detected by the increased noise level in the radar receiver. The technique is being installed on the UK operational network, and for the first time, allows radome attenuation to be monitored using the increased noise at the higher beam elevations. This attenuation has a large azimuthal dependence but for an old radome can be up to 4 dB for rainfall rates of just 2–4 mm/h. This effect has been neglected in the past, but may be responsible for significant errors in rainfall estimates and in radar calibrations using gauges. The extra noise at low radar elevations provides an estimate of the total path integrated attenuation of nearby storms; this total attenuation can then be used as a constraint for gate-by-gate or polarimetric correction algorithms.

Temporal and spatial variability of auroral forms in the 10–14 MLT sector: Relationship to plasma convection and solar wind-magnetosphere coupling

Ground-based observations of dayside auroral forms and magnetic perturbations in the arctic sectors of Svalbard and Greenland, in combination with the high-resolution measurements of ionospheric ion drift and temperature by the EISCAT radar, are used to study temporal/spatial structures of cusp-type auroral forms in relation to convection. Large-scale patterns of equivalent convection in the dayside polar ionosphere are derived from the magnetic observations in Greenland and Svalbard. This information is used to estimate the ionospheric convection pattern in the vicinity of the cusp/cleft aurora. The reported observations, covering the period 0700-1130 UT, on January 11, 1993, are separated into four intervals according to the observed characteristics of the aurora and ionospheric convection. The morphology and intensity of the aurora are very different in quiet and disturbed intervals. A latitudinally narrow zone of intense and dynamical 630.0 nm emission equatorward of 75 degrees MLAT, was observed during periods of enhanced antisunward convection in the cusp region. This (type 1 cusp aurora) is considered to be the signature of plasma entry via magnetopause reconnection at low magnetopause latitudes, i.e. the low-latitude boundary layer (LLB I,). Another zone of weak 630.0 nm emission (type 2 cusp aurora) was observed to extend up to high latitudes (similar to 79 degrees MLAT) during relatively quiet magnetic conditions, when indications of reverse (sunward) convection was observed in the dayside polar cap. This is postulated to be a signature of merging between a northward directed IMF (B-z > 0) and the geomagnetic field poleward of the cusp. The coexistence of type 1 and 2 auroras was observed under intermediate circumstances. The optical observations from Svalbard and Greenland were also used to determine the temporal and spatial evolution of type 1 auroral forms, i.e. poleward-moving auroral events occurring in the vicinity of a rotational convection reversal in the early post-noon sector. Each event appeared as a local brightening at the equatorward boundary of the pre-existing type 1 cusp aurora, followed by poleward and eastward expansions of luminosity. The auroral events were associated with poleward-moving surges of enhanced ionospheric convection and F-layer ion temperature as observed by the EISCAT radar in Tromso. The EISCAT ion flow data in combination with the auroral observations show strong evidence for plasma flow across the open/closed field line boundary.

Causes of plasma flow bursts and dayside auroral transients: An evaluation of two models Invoking reconnection pulses and changes in the Y component of the magnetosheath field

Longitudinal flow bursts observed by the European Incoherent Scatter (EISCAT) radar, in association with dayside auroral transients observed from Svalbard, have been interpreted as resulting from pulses of enhanced reconnection at the dayside magnetopause. However, an alternative model has recently been proposed for a steady rate of magnetopause reconnection, in which the bursts of longitudinal flow are due to increases in the field line curvature force, associated with the By component of the magnetosheath field. We here evaluate these two models, using observations on January 20, 1990, by EISCAT and a 630-nm all-sky camera at Ny Ålesund. For both models, we predict the behavior of both the dayside flows and the 630-nm emissions on newly opened field lines. It is shown that the signatures of steady reconnection and magnetosheath By changes could possibly resemble the observed 630-nm auroral events, but only for certain locations of the observing site, relative to the ionospheric projection of the reconnection X line: however, in such cases, the flow bursts would be seen between the 630-nm transients and not within them. On the other hand, the model of reconnection rate pulses predicts that the flows will be enhanced within each 630-nm transient auroral event. The observations on January 20, 1990, are shown to be consistent with the model of enhanced reconnection rate pulses over a background level and inconsistent with the effects of periodic enhancements of the magnitude of the magnetosheath By component. We estimate that the reconnection rate within the pulses would have to be at least an order of magnitude larger than the background level between the pulses.

EISCAT observations of bursts of rapid flow in the high latitude dayside ionosphere

Observations are presented of short-lived, highly structured bursts of rapid plasma flow observed with the EISCAT radar in the high latitude dayside ionosphere. It is shown that the properties of the bursts are consistent with ionospheric perturbations caused by impulsive, localized reconnection at the Earth's magnetopause, i.e. by flux transfer events.

The DYMECS project: a statistical approach for the evaluation of convective storms in high-resolution NWP models

A new frontier in weather forecasting is emerging by operational forecast models now being run at convection-permitting resolutions at many national weather services. However, this is not a panacea; significant systematic errors remain in the character of convective storms and rainfall distributions. The DYMECS project (Dynamical and Microphysical Evolution of Convective Storms) is taking a fundamentally new approach to evaluate and improve such models: rather than relying on a limited number of cases, which may not be representative, we have gathered a large database of 3D storm structures on 40 convective days using the Chilbolton radar in southern England. We have related these structures to storm life-cycles derived by tracking features in the rainfall from the UK radar network, and compared them statistically to storm structures in the Met Office model, which we ran at horizontal grid length between 1.5 km and 100 m, including simulations with different subgrid mixing length. We also evaluated the scale and intensity of convective updrafts using a new radar technique. We find that the horizontal size of simulated convective storms and the updrafts within them is much too large at 1.5-km resolution, such that the convective mass flux of individual updrafts can be too large by an order of magnitude. The scale of precipitation cores and updrafts decreases steadily with decreasing grid lengths, as does the typical storm lifetime. The 200-m grid-length simulation with standard mixing length performs best over all diagnostics, although a greater mixing length improves the representation of deep convective storms.

NDVI With Artificial Neural Networks For SRTM Elevation Model Improvement – Hydrological Model Application

Digital elevation model (DEM) plays a substantial role in hydrological study, from understanding the catchment characteristics, setting up a hydrological model to mapping the flood risk in the region. Depending on the nature of study and its objectives, high resolution and reliable DEM is often desired to set up a sound hydrological model. However, such source of good DEM is not always available and it is generally high-priced. Obtained through radar based remote sensing, Shuttle Radar Topography Mission (SRTM) is a publicly available DEM with resolution of 92m outside US. It is a great source of DEM where no surveyed DEM is available. However, apart from the coarse resolution, SRTM suffers from inaccuracy especially on area with dense vegetation coverage due to the limitation of radar signals not penetrating through canopy. This will lead to the improper setup of the model as well as the erroneous mapping of flood risk. This paper attempts on improving SRTM dataset, using Normalised Difference Vegetation Index (NDVI), derived from Visible Red and Near Infra-Red band obtained from Landsat with resolution of 30m, and Artificial Neural Networks (ANN). The assessment of the improvement and the applicability of this method in hydrology would be highlighted and discussed.

Performance Evaluation Of SMOS Soil Moisture Retrieval Parameters For Hydrological Application

Microwave remote sensing has high potential for soil moisture retrieval. However, the efficient retrieval of soil moisture depends on optimally choosing the soil moisture retrieval parameters. In this study first the initial evaluation of SMOS L2 product is performed and then four approaches regarding soil moisture retrieval from SMOS brightness temperature are reported. The radiative transfer equation based tau-omega rationale is used in this study for the soil moisture retrievals. The single channel algorithms (SCA) using H polarisation is implemented with modifications, which includes the effective temperatures simulated from ECMWF (downscaled using WRF-NOAH Land Surface Model (LSM)) and MODIS. The retrieved soil moisture is then utilized for soil moisture deficit (SMD) estimation using empirical relationships with Probability Distributed Model based SMD as a benchmark. The square of correlation during the calibration indicates a value of R2 =0.359 for approach 4 (WRF-NOAH LSM based LST with optimized roughness parameters) followed by the approach 2 (optimized roughness parameters and MODIS based LST) (R2 =0.293), approach 3 (WRF-NOAH LSM based LST with no optimization) (R2 =0.267) and approach 1(MODIS based LST with no optimization) (R2 =0.163). Similarly, during the validation a highest performance is reported by approach 4. The other approaches are also following a similar trend as calibration. All the performances are depicted through Taylor diagram which indicates that the H polarisation using ECMWF based LST is giving a better performance for SMD estimation than the original SMOS L2 products at a catchment scale.

An intelligent system to real time rainfall prediction using radar data

This work presents a new approach for rainfall measurements making use of weather radar data for real time application to the radar systems operated by institute of Meteorological Research (IPMET) - UNESP - Bauru - SP-Brazil. Several real time adjustment techniques has been presented being most of them based on surface rain-gauge network. However, some of these methods do not regard the effect of the integration area, time integration and distance rainfall-radar. In this paper, artificial neural networks have been applied for generate a radar reflectivity-rain relationships which regard all effects described above. To evaluate prediction procedure, cross validation was performed using data from IPMET weather Doppler radar and rain-gauge network under the radar umbrella. The preliminary results were acceptable for rainfalls prediction. The small errors observed result from the spatial density and the time resolution of the rain-gauges networks used to calibrate the radar.

Geoestatística aplicada à acumulação da precipitação pluviométrica com radar meteorológico

Pós-graduação em Agronomia (Energia na Agricultura) - FCA

Analisi di dati telerilevati ottici e radar per la gestione dei disastri: le alluvioni nel Bangladesh.

Abstract L’utilizzo dei dati satellitari per la gestione dei disastri naturali è fondamentale nei paesi in via di sviluppo, dove raramente esiste un censimento ed è difficile per i governi aggiornare le proprie banche dati con le tecniche di rilevamento e metodi di mappatura tradizionali che sono entrambe lunghe e onerose. A supporto dell’importanza dell’impiego del telerilevamento e per favorirne l’uso nel caso di catastrofi, vi è l’operato di diverse organizzazioni internazionali promosse da enti di ricerca, da agenzie governative o da organismi sopranazionali, le quali svolgono un lavoro di cruciale valore, fornendo sostegno tecnico a chi si occupa di far giungere alle popolazioni colpite gli aiuti umanitari e i soccorsi nel più breve tempo possibile. L’attività di tesi è nata proprio dalla collaborazione con una di esse, ITHACA (Information Technology for Humanitarian Assistance, Cooperation and Action), organizzazione no-profit, fondata dal Politecnico di Torino e SiTI (Istituto Superiore sui Sistemi Territoriali per l’Innovazione), la quale a sua volta collabora con il WFP (World Food Programme) delle Nazioni Unite, realizzando cartografie speditive necessarie per la valutazione delle conseguenze di un evento catastrofico, attraverso l’impiego di dati acquisiti da satellite. Su questo tema si è inserito il presente lavoro che ha come obiettivo quello di dimostrare la valenza dei dati telerilevati, siano essi di tipo ottico o Radar, nel caso di alcuni dei disastri naturali più catastrofici, le alluvioni. In particolare è stata studiata la vulnerabilità del Bangladesh, il quale annualmente si trova ad affrontare eventi alluvionali, spesso di grave intensità. Preliminarmente allo studio, è stata condotta una ricerca bibliografica al fine di avere una buona conoscenza dell’area sia in termini geografici e fisici che di sviluppo e tipologia di urbanizzazione. E’stata indagata in particolare l’alluvione che ha colpito il paese nel Luglio del 2004, attraverso delle immagini satellitari multispettrali, in particolare Landsat 7, per un inquadramento pre-evento, ed ASTER per studiare la situazione a distanza di tre mesi dall’accaduto (immagine rilevata il 20 Ottobre 2004). Su tali immagini sono state condotte delle classificazioni supervisionate con il metodo della massima verosimiglianza che hanno portato la suddivisione del territorio in quattro classi di destinazione d’uso del suolo: urbano (Build-up), campi e vegetazione (Crops&Vegetation), sabbia e scavi (Sand&Excavation), idrografia e zone alluvionate (Water). Dalla sperimentazione è emerso come tali immagini multispettrali si prestino molto bene per l’analisi delle differenti caratteristiche del territorio, difatti la validazione condotta sulla mappa tematica derivata dall’immagine Landsat 7 ha portato ad un’accuratezza del 93% circa, mentre la validazione dell’immagine ASTER è stata solo di tipo qualitativo, in quanto, considerata l’entità della situazione rilevata, non è stato possibile avere un confronto con dei punti da assumere come verità a terra. Un’interpretazione della mappa tematica derivante dalla classificazione dell’immagine ASTER è stata elaborata incrociandola in ambiente GIS con dati forniti dal CEGIS (Center for Environmental and Geographic Information Services) riguardanti il landuse della zona in esame; da ciò è emerso che le zone destinate alla coltivazione del riso sono più vulnerabili alle inondazioni ed in particolare nell’Ottobre 2004 il 95% delle aree esondate ha interessato tali colture. Le immagini ottiche presentano un grosso limite nel caso delle alluvioni: la rilevante copertura nuvolosa che spesso accompagna siffatti eventi impedisce ai sensori satellitari operanti nel campo dell’ottico di rilevare il territorio, e per questo di frequente essi non si prestano ad essere impiegati per un’indagine nella fase di prima emergenza. In questa circostanza, un valido aiuto giunge dall’impiego di immagini Radar, le quali permettono osservazioni ad ogni ora del giorno e della notte, anche in presenza di nuvole, rendendole di fondamentale importanza nelle situazioni descritte. Per dimostrare la validità di questi sensori si sono analizzati due subset derivanti da un mosaico di immagini della nuova costellazione italiana ad alta risoluzione CosmoSkymed: il primo va dalla città di Dhaka al Golfo del Bengala ed il secondo copre la zona più a Nord nel distretto di Sylhet. Dalla sperimentazione condotta su tali immagini radar, che ha comportato come ovvio problematiche del tutto differenti rispetto alle elaborazioni tradizionalmente condotte su immagini nel campo dell’ottico, si è potuto verificare come l’estrazione dei corpi d’acqua e più in generale dell’idrografia risulti valida e di veloce computazione. Sono emersi tuttavia dei problemi, per esempio per quanto riguarda la classificazione dell’acqua in presenza di rilievi montuosi; tali complicazioni sono dovute alla presenza di zone d’ombra che risultano erroneamente assegnate alla classe water, ma è stato possibile correggere tali errori di attribuzione mascherando i rilievi con l’ausilio di una mappa delle pendenze ricavata da modelli di elevazione SRTM (Shuttle Radar Topographic Mission). La validazione dei risultati della classificazione, condotta con un grande numero di check points, ha fornito risultati molto incoraggianti (ca. 90%). Nonostante le problematiche riscontrate, il Radar, in sé o in accoppiamento con altri dati di diversa origine, si presta dunque a fornire in breve tempo informazioni sull’estensione dell’esondazione, sul grado di devastazione, sulle caratteristiche delle aree esondate, sulle vie di fuga più adatte, diventando un’importante risorsa per chi si occupa di gestire l’emergenza in caso di eventi calamitosi. L’integrazione con i dati di tipo ottico è inoltre essenziale per pervenire ad una migliore caratterizzazione del fenomeno, sia in termini di change detection che di monitoraggio post-evento.

EXPLORING CHANGES IN DETRITAL FLOCCULENT LAYER DYNAMICS DUE TO SHIFTS IN MACROPHYTE COMMUNITIES IN THE NORTHERN EVERGLADES

A shift in plant communities of the Water Conservation Areas (WCAs) within the Everglades has been linked to changes in hydrology and high levels of nutrient loading from surrounding agicultural areas. This has resulted in the encroachment of dense cattail stands (Typha domingensis) into areas that had previously been a ridge and slough landscape populated primarily by native sawgrass (Cladium jamaicense). In order to study ecological management solutions in this area, WCA-2A was broken into study plots; several of which became open water areas through the application of herbicide and burning regimens. The open water areas allowed for Chara spp (a submersed algal species) to replace Typha domingensis as the dominant macrophyte. This study investigated the polymer and ionic profiles of Chara spp, Typha domingensis and Cladium jamaicense and their contributions to detrital flocculent (floc) in the study plots where they are the dominant macrophytes. Floc is not only an important food source for aquatic species; it also supports many algal, fungal and bacterial communities. Data gathered in this study indicated that the floc sample from a phosphorus enriched open water study plot (EO1) where Chara spp was the dominant macrophyte may contain cell wall polymers from sources other than Chara spp (most likely Typha domingensis), while the chemical and polymeric profile of the floc of the study plot where Typha domingensis is the dominant macrophyte (EC1) suggests that the floc layer has contributions from algal sources as well as Typha domingensis. Additionally, monoclonal antibodies to Arabinoglalactan protein (AGP) and (1,4)-β-D galactan were identified as possible biomarkers for distinguishing algal dominated floc layers from layers dominated by emergent vegetation. Calcium labeling could be a useful tool for this as well because of the high amount of Ca2+ associated with Chara spp cell walls. When looking into the soluble phosphorus content of the macrophytes and paired floc samples of WCA-2A, it was found that Chara spp may be contributing a greater amount of Ca-bound phosphorus to floc layers where it is the dominant macrophyte when compared to floc layers from study plots dominated by emergent macrophytes. Floc layers also appear to be acting as a nutrient sink for soluble phosphorus. The findings of this study support the overall hypothesis that the shift from native emergent macrophyte communities to submersed macrophyte communities in study sites of the northern Everglades is affecting the polymeric/chemical profile and ionic content of detrital floc layers. The effects of this shift may contribute to changes in complex flocculent community dynamics.

Origin and evolution of reactive and noble gases dissolved in matrix pore water

Reactive and noble gases dissolved in matrix pore water of low permeable crystalline bedrock were successfully extracted and characterized for the fist time based on drillcore samples from the Olkiluoto investigation site (SW Finland). Interaction between matrix pore water and fracture groundwater occurs predominately by diffusion. Changes in the chemical and isotopic composition of gases dissolved in fracture groundwater are transmitted and preserved in the pore water. Absolute concentrations, their ratios and the stable carbon isotope signature of hydrocarbon gases dissolved in pore water give valuable indications about the evolution of these gases in the nearby-flowing fracture groundwaters. Inert noble gases dissolved in matrix pore water and their isotopes combined with their in-situ production and accumulation rates deliver information about the residence time of pore water.

Matrix pore water in low-permeable crystalline bedrock: An archive for the palaeohydrogeological evolution of the Olkiluoto Investigation Site

Matrix pore water in the connected inter- and intragranular pore space of low-permeable crystalline bedrock interacts with flowing fracture groundwater predominately by diffusion. Based on the slow exchange between the two water reservoirs, matrix pore water acts as an archive of past changes in fracture groundwater compositions and thus of the palaeohydrological history of a site. Matrix pore water of crystalline bedrock from the Olkiluoto investigation site (SW Finland) was characterised using the stable water isotopes (δ18O, δ2H), combined with the concentrations of dissolved chloride and bromide as natural tracers. The comparison of tracer concentrations in pore water and present-day fracture groundwater suggest for the pore water the presence of old, dilute meteoric water components that infiltrated into the fractures during various warm climate stages. These different meteoric components can be discerned based on the diffusion distance between the two reservoirs and be brought into context with the palaeohydrological evolution of the site.