FIELD OF THE INVENTION 
       [0001]    The invention relates to the field of weather forecasting and more particularly to the field of short term precipitation forecasting (or precipitation nowcasting). 
       BACKGROUND OF THE INVENTION 
       [0002]    Numerical Weather Predictors (NWPs) are computer systems that predict weather by modeling the atmosphere and using measured weather or atmospheric values as initial input and/or as feedback to predict future weather conditions. NWP&#39;s can be very sophisticated and expensive systems that are commonly used for long term precipitation forecasting (more than 12 hours). NWPs comprise generally a data simulator that takes into account a multiple of variables (such as atmospheric heat, atmospheric pressure, wind motion, etc.) to predict weather conditions. NWPs are often able to predict all weather parameters, including temperature, wind, and precipitation (namely rain, sleet or snow). NWP&#39;s have limitations in predicting weather due to errors in modeling and a lack of precision or detail in the observations of weather parameters used in the model. 
         [0003]    Due to these limitations, NWP&#39;s are not considered efficient for short term precipitation forecasting in comparison to short term precipitation forecasting systems. Known precipitation nowcasting systems are very simple systems that use the precipitation pattern persistence principle in order to predict short term future precipitation patterns. Because they are based on the continuation of what has actually taken place in the very recent past, the immediate short term prediction accuracy is very good. While persistence cannot predict new events, such as when rain will begin to fall, or a change in weather patterns when two fronts collide, absent such abrupt changes, persistence-based nowcasting systems are more accurate than NWP&#39;s for short term weather prediction. 
         [0004]    Doppler weather radar networks can be found today across North America and are able to provide accurate images of spatial distribution of precipitation within large regional areas. The persistence principle involves using observed precipitation patterns from the very recent past (obtained by Doppler radar, satellite images or even weather station observations) to model precipitation pattern motion and project it in the future as a function of wind motion in presuming persistency in time of the observed patterns. Such nowcasting systems are efficient for predicting weather, and most importantly precipitation, for the future few hours with a geographical accuracy corresponding to the observation density. 
         [0005]    Even though such precipitation nowcasting systems are widely used by precipitation forecasters because of their simplicity and low cost, they lack the ability to predict changes to weather that can occur in the short term, such as a sudden development of precipitation by way of thundershower or drop in atmospheric temperature, that only NWPs can attempt to predict. 
       SUMMARY OF THE INVENTION 
       [0006]    It is therefore an object of the present invention to provide a system to enhance accuracy of nowcasting systems, and in some embodiment, precipitation nowcasting systems. 
         [0007]    The present invention takes advantage of the complexity of NWP to enhance accuracy of nowcasting systems. It provides a short term weather forecasting system for nowcasting using a Numerical Weather Predictor (NWP). The present nowcasting system tracks the evolution of differences between NWP and radar based precipitation patterns and adjusts the NWP forecast to account for theses differences. These differences are due to amplitude and phase errors of the NWP, such as model misses, false alarms, intensity errors and position errors. In presuming persistency in time of these errors and in estimating precipitation pattern time evolution due to other weather conditions such as wind motion, the present system corrects the NWP short term predicted patterns to compensate for these errors, thus enhancing nowcasting accuracy. 
         [0008]    According to one aspect of the invention, there is provided a short term weather forecaster for nowcasting comprising a numerical weather predictor receiving measured weather variables and generating, as a function of a first mathematical model and of the measured weather variables, short term and long term predicted weather variables, an error modeling unit receiving the predicted variables and measured weather variable corresponding to the predicted weather variables, and generates a time function error model based on short term past predicted and measured weather variables, and a short term error correction unit, where the unit receives the short term predicted weather variables and the time function error model, corrects the short term weather variables as a function of the time function error model and outputs corrected short term predicted weather variables associated with the short term weather predicted variables. 
         [0009]    In some embodiments, the invention is applied to precipitation, instead of other weather variables such as temperature, pressure, wind, etc. The measured weather variables and the predicted weather variables can be precipitation variables representing precipitation maps, and the error model can correct errors in the precipitation maps. The measured weather variables can in some embodiments include wind motion variables, while the error modeling unit identifies a miss, and the error correction unit generates a prediction of the miss using persistence of the measured wind motion and precipitation variables. Wind motion variables may be estimated using a variational echo tracker. 
         [0010]    In some embodiments, the error correction unit generates a persistence prediction of short term weather based on persistence of the measured weather variables and combines the persistence prediction with the short term weather variables corrected as a function of the time function error model to output the corrected short term predicted weather variables associated with the short term weather predicted variables. 
         [0011]    Short term forecasting is often understood to mean less than twelve hours, and reasonable accuracy of short term forecasting may be limited to about six hours in certain climate conditions. 
         [0012]    In some embodiment, the short term weather forecaster for nowcasting is for use with a numerical weather predictor receiving measured weather variables and generating, as a function of a first mathematical model and of the measured weather variables, short term and long term predicted weather variables. As above, the short term weather forecaster may include an error modeling unit receiving the predicted variables and measured weather variable corresponding to the predicted weather variables, and generates a time function error model based on short term past predicted and measured weather variables, and a short term error correction unit, where the unit receives the short term predicted weather variables and the time function error model, corrects the short term weather variables as a function of the time function error model and outputs corrected short term predicted weather variables associated with the short term weather predicted variables. 
         [0013]    In other embodiments, a numerical weather predictor (NWP) for forecasting comprises a numerical weather prediction module receiving measured weather variables and generating, as a function of a first mathematical model and of the measured weather variables, short term and long term predicted weather variables, an error modeling unit receiving the short term predicted weather variables and corresponding measured weather variables, and generates an error model based on short term past predicted and measured weather variables, and the numerical weather prediction module uses the error model to estimate adjusted initial conditions of the first mathematical model. The numerical prediction module uses the error model to generate short term corrected predicted weather variables, the corrected predicted weather variables being used to estimate the adjusted initial conditions of the first mathematical model. 
         [0014]    In further embodiments, the invention is embodied as a computer program product. This product comprises a data recording medium, such as a CD-ROM, electronic memory device, magnetic storage drive, etc. having recorded thereon executable computer program code that when loaded into a suitable computer and executed provides a short term weather forecaster according to any one of the embodiments described or defined in this specification. 
     
     
       BRIEF DESCRIPTION OF THE DRAWINGS 
         [0015]      FIG. 1  shows a block diagram of a short term weather forecaster for nowcasting using a Numerical Weather Predictor (NWP); 
           [0016]      FIG. 2  illustrates an amplitude error of the NWP due to a model miss; 
           [0017]      FIG. 3  illustrates an amplitude error of the NWP due to a false alarm; 
           [0018]      FIG. 4  illustrates an amplitude error of the NWP due to an intensity error; 
           [0019]      FIG. 5  illustrates a phase error (position error) of the NWP; 
           [0020]      FIG. 6  illustrates the cumulative error effect due to model misses, false alarms, intensity errors and phase errors; 
           [0021]      FIG. 7  illustrates generation of the extrapolation function (tao) using the error distance between observed and measured precipitation patterns; 
           [0022]      FIG. 8  illustrates use of the extrapolation function (tao) to generate a short term forecasting precipitation pattern and compare it to the measured pattern; 
           [0023]      FIG. 9  is a Model-Radar correlation curve illustrating enhancement of accuracy between a traditional nowcasting system and the nowcasting system as provided by the present invention. 
       
    
    
     DETAILED DESCRIPTION OF THE INVENTION 
       [0024]    The functionality of the developed weather nowcasting system could be detailed as follows. A Weather Observation Unit  2 , generally connected to a weather radar, measures weather conditions and transmits Measured Weather Variables  4  to a Numerical Weather Predictor (NWP)  6 . The Measured Weather Variables  4  comprise precipitation variables, wind motion variables and any other measured variables related to weather conditions. The NWP  6  receives the Measured Weather Variables  4 , processes them according to a precipitation-forecasting algorithm and outputs Long Term Predicted Weather Variables  8  and Short Term Predicted Weather Variables  10 . The Short Term Predicted Weather Variables  10  and the corresponding Measured Weather Variables  4  are transmitted to an Error Modeling Unit  12 . In one embodiment of the invention, the Short Term Predicted Weather Variables  10  include the present state of precipitation (at t 0 ) and the state of precipitation just before (at t −1 ) and the Error Modeling Unit  12  measures the error “distance” between predicted precipitation and measured precipitation over the short term past. 
         [0025]    By error distance, it is understood to mean an error function or model describing the spatial or morphological and/or temporal error or difference in the spatial or spatio-temporal map of weather variables. The error model can be a correction in space or distance only (namely a spatial shift of a region within the map corresponding to an event), an error in time only, an account of false event prediction only (such as an incorrect forecast of rain from a group of clouds, while actual measurements do not show rain), an account of a failure or a miss to predict an event only (such as an incorrect forecast of no rain from a group of clouds, while actual measurements show rain), or an error in intensity of a prediction only (such as a level of temperature, wind or precipitation within a local area or event). Of course, the error model can combine one or more of such elements. In the case of a false event prediction, the error function or model strives to remove or greatly attenuate the event from the short term predicted values from the NWP. In the case of a miss, the event is predicted from the persistence of the observed past values without using the NWP&#39;s predicted values, and the NWP&#39;s past values are used merely to identify the miss. 
         [0026]    As shown in  FIG. 2 , there is a model miss when the NWP  6  fails to predict, for a specific spot in space and for a specific time, an existing state of precipitation measured by a Radar. On the contrary, there is a false alarm when the NWP  6  predicts, for a specific spot in space and for a specific time, a non-existing state of precipitation according to the Radar (see  FIG. 3 ). There is an intensity error when the predicted intensity of precipitation is different from the measured one (see  FIG. 4 ). As regards the position error, the latter arises when the NWP  6  predicts a state of precipitation in a certain spot in space different from the real spot where the given state of precipitation arose according to the Radar (see  FIG. 5 ). 
         [0027]    After measuring the error distance, the Error Modeling Unit  12  generates and updates a Time Function Error Model  14  (see  FIG. 6 to 8 ). The Time Function Error Model  14  is then transmitted to a Short Term Error Correction Unit  16  that corrects, according to the Error Model  14 , the Short Term Predicted Variables  10  transmitted by the NWP  6 . 
         [0028]    The Time Function Error Model  14  is preferably an extrapolation function that extrapolates error occurrences in the future time, based on occurred errors and on measured weather variables such as precipitation motion direction and precipitation motion speed. The Short Term Error Correction Unit  16  outputs Corrected Short Term Predicted Variables  18  associated with the received Short Term Predicted Variables  10 .  FIG. 9  illustrates an enhancement of accuracy of precipitation nowcasting (short term prediction) by using the above-illustrated system. 
         [0029]    It is also possible to enhance the accuracy of precipitation forecasting (long term prediction) by using the Corrected Short Term Predicted Variables  18  to set up the initial conditions of the mathematical model of the NWP  6 . 
         [0030]    The present invention is based on correction of errors of NWP (numerical weather prediction) outputs. The correction of model outputs is performed with two components: phase and amplitude (or intensity) errors (τ P  and τ A ). These errors and their tendencies are determined by comparison with current and past observations for each geographical pixel determined by the resolution of the NWP. Then, during a certain forecast time model outputs are modified by correcting the time-dependent errors. For precipitation we apply various corrections: 
       1. Correction of Constant Phase Errors (CPE), τ P (x,t 0 ): 
       [0031]    Establish a homomorphism that is, a pixel-to-pixel correspondence between model precipitation output and radar patterns. For this minimize the following cost function 
         [0000]    
       
         
           
             J 
             = 
             
               
                 J 
                 ψ 
               
               + 
               
                 J 
                 s 
               
             
           
         
       
       
         
           where 
         
       
       
         
           
             
               J 
               ψ 
             
             = 
             
               ∫ 
               
                 
                   ∫ 
                   Ω 
                   
                       
                   
                 
                  
                 
                   
                     
                       a 
                        
                       
                         [ 
                         
                           
                             
                               ψ 
                               R 
                             
                              
                             
                               ( 
                               
                                 
                                   x 
                                   + 
                                   α 
                                 
                                 , 
                                 
                                   y 
                                   + 
                                   β 
                                 
                                 , 
                                 t 
                               
                               ) 
                             
                           
                           - 
                           
                             
                               ψ 
                               M 
                             
                              
                             
                               ( 
                               
                                 x 
                                 , 
                                 y 
                                 , 
                                 t 
                               
                               ) 
                             
                           
                         
                         ] 
                       
                     
                     2 
                   
                    
                   
                       
                   
                    
                   
                      
                     x 
                   
                    
                   
                      
                     y 
                   
                 
               
             
           
         
       
       
         
           
             
               J 
               s 
             
             = 
             
               b 
                
               
                 ∫ 
                 
                   
                     ∫ 
                     Ω 
                     
                         
                     
                   
                    
                   
                     
                       [ 
                       
                         
                           ( 
                           
                             
                               
                                 ∂ 
                                 2 
                               
                                
                               α 
                             
                             
                               ∂ 
                               
                                 x 
                                 2 
                               
                             
                           
                           ) 
                         
                         + 
                         
                           ( 
                           
                             
                               
                                 ∂ 
                                 2 
                               
                                
                               α 
                             
                             
                               ∂ 
                               
                                 y 
                                 2 
                               
                             
                           
                           ) 
                         
                         + 
                         
                           ( 
                           
                             
                               
                                 ∂ 
                                 2 
                               
                                
                               β 
                             
                             
                               ∂ 
                               
                                 x 
                                 2 
                               
                             
                           
                           ) 
                         
                         + 
                         
                           ( 
                           
                             
                               
                                 ∂ 
                                 2 
                               
                                
                               β 
                             
                             
                               ∂ 
                               
                                 y 
                                 2 
                               
                             
                           
                           ) 
                         
                         + 
                         
                           ( 
                           
                             
                               
                                 ∂ 
                                 2 
                               
                                
                               α 
                             
                             
                               
                                 ∂ 
                                 x 
                               
                                
                               
                                 ∂ 
                                 y 
                               
                             
                           
                           ) 
                         
                         + 
                         
                           ( 
                           
                             
                               
                                 ∂ 
                                 2 
                               
                                
                               β 
                             
                             
                               
                                 ∂ 
                                 x 
                               
                                
                               
                                 ∂ 
                                 y 
                               
                             
                           
                           ) 
                         
                       
                        
                       
                           
                       
                       ] 
                     
                      
                     
                        
                       x 
                     
                      
                     
                        
                       y 
                     
                   
                 
               
             
           
         
       
     
         [0032]    Here ψ R  and ψ M  are the precipitation intensity or accumulation values, as a function of space at time t=t 0  of the initiation of the nowcast, measured by radar and given by the NWP model output, respectively; Ω is the domain over which ARMOR is applied. α and β are the control variables of the minimization problem. The solution of the minimization gives, for each pixel, a vector 
         [0000]    
       
         
           
             
               ( 
               
                 
                   
                     α 
                   
                 
                 
                   
                     β 
                   
                 
               
               ) 
             
             . 
           
         
       
     
         [0000]    The ensemble of these vectors gives the x and y components of initial spatial phase errors within a domain Ω, that is, the matrix of errors τ P (x,t 0 ). This matrix represents the full two-dimensional field of vectors necessary to produce the displacements and deformations of the NWP model output to match the observations at t 0 . 
         [0033]    The minimization of the cost function is performed by a variational method using a conjugate-gradient method (although other methods are equally possible). In this manner a field of vectors (τ P ) is determined over the domain, one vector per each resolution pixel of the NWP precipitation output. The parameters a, b are adjustable weights, with a representing the uncertainty in radar measurements and b chosen as an empirical compromise between eliminating noise in the retrieved spatial phase error vectors and the spatial variability in the phase error vectors. To account for time phase errors the cost function is modified to 
         [0000]    
       
         
           
             
               J 
               ψ 
             
             = 
             
               ∫ 
               
                 ∫ 
                 
                   
                     ∫ 
                     Ω 
                     
                         
                     
                   
                    
                   
                     
                       
                         a 
                          
                         
                           [ 
                           
                             
                               
                                 ψ 
                                 R 
                               
                                
                               
                                 ( 
                                 
                                   
                                     x 
                                     + 
                                     α 
                                   
                                   , 
                                   
                                     y 
                                     + 
                                     β 
                                   
                                   , 
                                   
                                     t 
                                     + 
                                     γ 
                                   
                                 
                                 ) 
                               
                             
                             - 
                             
                               
                                 ψ 
                                 M 
                               
                                
                               
                                 ( 
                                 
                                   x 
                                   , 
                                   y 
                                   , 
                                   t 
                                 
                                 ) 
                               
                             
                           
                           ] 
                         
                       
                       2 
                     
                      
                     
                         
                     
                      
                     
                        
                       x 
                     
                      
                     
                        
                       y 
                     
                   
                 
               
             
           
         
       
       
         
           
             
               J 
               s 
             
             = 
             
               
                 b 
                  
                 
                   ∫ 
                   
                     ∫ 
                     
                       
                         ∫ 
                         Ω 
                         
                             
                         
                       
                        
                       
                         
                           [ 
                           
                             
                               ( 
                               
                                 
                                   
                                     ∂ 
                                     2 
                                   
                                    
                                   α 
                                 
                                 
                                   ∂ 
                                   
                                     x 
                                     2 
                                   
                                 
                               
                               ) 
                             
                             + 
                             
                               ( 
                               
                                 
                                   
                                     ∂ 
                                     2 
                                   
                                    
                                   α 
                                 
                                 
                                   ∂ 
                                   
                                     y 
                                     2 
                                   
                                 
                               
                               ) 
                             
                             + 
                             
                               ( 
                               
                                 
                                   
                                     ∂ 
                                     2 
                                   
                                    
                                   β 
                                 
                                 
                                   ∂ 
                                   
                                     x 
                                     2 
                                   
                                 
                               
                               ) 
                             
                             + 
                             
                               ( 
                               
                                 
                                   
                                     ∂ 
                                     2 
                                   
                                    
                                   β 
                                 
                                 
                                   ∂ 
                                   
                                     y 
                                     2 
                                   
                                 
                               
                               ) 
                             
                             + 
                             
                               2 
                                
                               
                                 ( 
                                 
                                   
                                     
                                       ∂ 
                                       2 
                                     
                                      
                                     α 
                                   
                                   
                                     
                                       ∂ 
                                       x 
                                     
                                      
                                     
                                       ∂ 
                                       y 
                                     
                                   
                                 
                                 ) 
                               
                             
                             + 
                             
                               2 
                                
                               
                                 ( 
                                 
                                   
                                     
                                       ∂ 
                                       2 
                                     
                                      
                                     β 
                                   
                                   
                                     
                                       ∂ 
                                       x 
                                     
                                      
                                     
                                       ∂ 
                                       y 
                                     
                                   
                                 
                                 ) 
                               
                             
                           
                            
                           
                               
                           
                           ] 
                         
                          
                         
                             
                         
                          
                         
                            
                           x 
                         
                          
                         
                            
                           y 
                         
                          
                         
                            
                           t 
                         
                       
                     
                   
                 
               
               + 
               
                 b 
                  
                 
                   ∫ 
                   
                     ∫ 
                     
                       
                         ∫ 
                         Ω 
                         
                             
                         
                       
                        
                       
                         
                           [ 
                           
                             
                               ( 
                               
                                 
                                   
                                     ∂ 
                                     2 
                                   
                                    
                                   γ 
                                 
                                 
                                   ∂ 
                                   
                                     t 
                                     2 
                                   
                                 
                               
                               ) 
                             
                             + 
                             
                               ( 
                               
                                 
                                   
                                     ∂ 
                                     2 
                                   
                                    
                                   γ 
                                 
                                 
                                   ∂ 
                                   
                                     x 
                                     2 
                                   
                                 
                               
                               ) 
                             
                              
                             
                                 
                             
                             + 
                             
                               ( 
                               
                                 
                                   
                                     ∂ 
                                     2 
                                   
                                    
                                   γ 
                                 
                                 
                                   ∂ 
                                   
                                     y 
                                     2 
                                   
                                 
                               
                               ) 
                             
                             + 
                             
                               ( 
                               
                                 
                                   
                                     ∂ 
                                     2 
                                   
                                    
                                   α 
                                 
                                 
                                   ∂ 
                                   
                                     t 
                                     2 
                                   
                                 
                               
                               ) 
                             
                             + 
                             
                               ( 
                               
                                 
                                   
                                     ∂ 
                                     2 
                                   
                                    
                                   β 
                                 
                                 
                                   ∂ 
                                   
                                     t 
                                     2 
                                   
                                 
                               
                               ) 
                             
                           
                           ] 
                         
                          
                         
                            
                           x 
                         
                          
                         
                            
                           y 
                         
                          
                         
                            
                           t 
                         
                       
                     
                   
                 
               
               + 
               
                 b 
                  
                 
                   ∫ 
                   
                     ∫ 
                     
                       
                         ∫ 
                         Ω 
                         
                             
                         
                       
                        
                       
                         
                           2 
                            
                           
                             [ 
                             
                               
                                 ( 
                                 
                                   
                                     
                                       ∂ 
                                       2 
                                     
                                      
                                     γ 
                                   
                                   
                                     
                                       ∂ 
                                       x 
                                     
                                      
                                     
                                       ∂ 
                                       y 
                                     
                                   
                                 
                                 ) 
                               
                               + 
                               
                                 ( 
                                 
                                   
                                     
                                       ∂ 
                                       2 
                                     
                                      
                                     γ 
                                   
                                   
                                     
                                       ∂ 
                                       t 
                                     
                                      
                                     
                                       ∂ 
                                       y 
                                     
                                   
                                 
                                 ) 
                               
                               + 
                               
                                 ( 
                                 
                                   
                                     
                                       ∂ 
                                       2 
                                     
                                      
                                     γ 
                                   
                                   
                                     
                                       ∂ 
                                       x 
                                     
                                      
                                     
                                       ∂ 
                                       t 
                                     
                                   
                                 
                                 ) 
                               
                                
                               
                                   
                               
                               + 
                               
                                 ( 
                                 
                                   
                                     
                                       ∂ 
                                       2 
                                     
                                      
                                     α 
                                   
                                   
                                     
                                       ∂ 
                                       t 
                                     
                                      
                                     
                                       ∂ 
                                       y 
                                     
                                   
                                 
                                 ) 
                               
                               + 
                               
                                 ( 
                                 
                                   
                                     
                                       ∂ 
                                       2 
                                     
                                      
                                     α 
                                   
                                   
                                     
                                       ∂ 
                                       t 
                                     
                                      
                                     
                                       ∂ 
                                       y 
                                     
                                   
                                 
                                 ) 
                               
                               + 
                               
                                 ( 
                                 
                                   
                                     
                                       ∂ 
                                       2 
                                     
                                      
                                     β 
                                   
                                   
                                     
                                       ∂ 
                                       t 
                                     
                                      
                                     
                                       ∂ 
                                       y 
                                     
                                   
                                 
                                 ) 
                               
                               + 
                               
                                 ( 
                                 
                                   
                                     
                                       ∂ 
                                       2 
                                     
                                      
                                     β 
                                   
                                   
                                     
                                       ∂ 
                                       t 
                                     
                                      
                                     
                                       ∂ 
                                       y 
                                     
                                   
                                 
                                 ) 
                               
                             
                             ] 
                           
                         
                          
                         
                            
                           x 
                         
                          
                         
                            
                           y 
                         
                          
                         
                            
                           t 
                         
                       
                     
                   
                 
               
             
           
         
       
     
         [0034]    The minimization of the above leads to a field of vectors 
         [0000]    
       
         
           
             ( 
             
               
                 
                   α 
                 
               
               
                 
                   β 
                 
               
               
                 
                   γ 
                 
               
             
             ) 
           
         
       
     
         [0000]    in the three-dimensional space giving the phase error (τ P ) for each pixel are the displacement in space and time necessary to establish the above pixel-to-pixel correspondence. 
         [0035]    Model outputs at a certain forecast time t i  are corrected by backward advection with derived three-dimensional initial phase error τ P (x,t 0 ). A cubic interpolation is done to place a corrected model value at a grid point. Since the correction of model errors is done with τ P (x,t 0 ), the growth and decay predicted by the NWP model is retained. However, τ P (x,t 0 ) does not take into account the time tendency of the phase errors. Thus, the following step is added to the procedure: 
       2. Correction of Lagrangian Time-Dependent Phase Errors (LTPE), τ P (x,t) 
       [0036]    Consider now the time tendency of the phase error along the movement of the precipitation pattern. First, determine the motion pattern, defined by the vector 
         [0000]    
       
         
           
             
               u 
               = 
               
                 ( 
                 
                   
                     
                       u 
                     
                   
                   
                     
                       v 
                     
                   
                 
                 ) 
               
             
             , 
           
         
       
     
         [0000]    of radar derived precipitation by minimizing the difference between radar precipitation patterns at a times t 0  and t 0 −Δt. For this minimize the following cost function: 
         [0000]    
       
         
           
             
               J 
               ψ 
             
             = 
             
               ∫ 
               
                 
                   ∫ 
                   Ω 
                   
                       
                   
                 
                  
                 
                   
                     
                       a 
                        
                       
                         [ 
                         
                           
                             
                               ψ 
                               R 
                             
                              
                             
                               ( 
                               
                                 x 
                                 , 
                                 y 
                                 , 
                                 t 
                               
                               ) 
                             
                           
                           - 
                           
                             
                               ψ 
                               R 
                             
                              
                             
                               ( 
                               
                                 
                                   x 
                                   - 
                                   
                                     u 
                                      
                                     
                                         
                                     
                                      
                                     Δ 
                                      
                                     
                                         
                                     
                                      
                                     t 
                                   
                                 
                                 , 
                                 
                                   y 
                                   - 
                                   
                                     v 
                                      
                                     
                                         
                                     
                                      
                                     Δ 
                                      
                                     
                                         
                                     
                                      
                                     t 
                                   
                                 
                                 , 
                                 
                                   t 
                                   - 
                                   
                                     Δ 
                                      
                                     
                                         
                                     
                                      
                                     t 
                                   
                                 
                               
                               ) 
                             
                           
                         
                         ] 
                       
                     
                     2 
                   
                    
                   
                       
                   
                    
                   
                      
                     x 
                   
                    
                   
                      
                     y 
                   
                 
               
             
           
         
       
       
         
           
             
               J 
               s 
             
             = 
             
               b 
                
               
                 ∫ 
                 
                   
                     ∫ 
                     Ω 
                     
                         
                     
                   
                    
                   
                     
                       [ 
                       
                         
                           ( 
                           
                             
                               
                                 ∂ 
                                 2 
                               
                                
                               u 
                             
                             
                               ∂ 
                               
                                 x 
                                 2 
                               
                             
                           
                           ) 
                         
                          
                         
                             
                         
                         + 
                         
                           ( 
                           
                             
                               
                                 ∂ 
                                 2 
                               
                                
                               u 
                             
                             
                               ∂ 
                               
                                 y 
                                 2 
                               
                             
                           
                           ) 
                         
                         + 
                         
                           ( 
                           
                             
                               
                                 ∂ 
                                 2 
                               
                                
                               v 
                             
                             
                               ∂ 
                               
                                 x 
                                 2 
                               
                             
                           
                           ) 
                         
                         + 
                         
                           ( 
                           
                             
                               
                                 ∂ 
                                 2 
                               
                                
                               v 
                             
                             
                               ∂ 
                               
                                 y 
                                 2 
                               
                             
                           
                           ) 
                         
                         + 
                         
                           ( 
                           
                             
                               
                                 ∂ 
                                 2 
                               
                                
                               u 
                             
                             
                               
                                 ∂ 
                                 x 
                               
                                
                               
                                 ∂ 
                                 y 
                               
                             
                           
                           ) 
                         
                         + 
                         
                           ( 
                           
                             
                               
                                 ∂ 
                                 2 
                               
                                
                               v 
                             
                             
                               
                                 ∂ 
                                 x 
                               
                                
                               
                                 ∂ 
                                 y 
                               
                             
                           
                           ) 
                         
                       
                       ] 
                     
                      
                     
                        
                       x 
                     
                      
                     
                        
                       y 
                     
                   
                 
               
             
           
         
       
     
         [0037]    Here, Δt is the time lag over which pattern motion is determined. 
         [0038]    The total time period over which motion of precipitation from radar and from NWP model outputs are to be compared is nΔt and n patterns of precipitation motion are derived for the period between t=t 0  and t=t 0 −nΔt. A more detailed description of the procedure to derive the motion pattern is given in Germann, U. and I. Zawadzki, 2002: “Scale dependence of predictability of precipitation from continental radar images. Part I: Description of methodology”, Monthly Weather Review, 130, pp 2859-2873. 
         [0039]    The same procedure is then repeated for ψ M  for a period between t=t 0 −nΔt and t=t i  (where t i  is the forecast time) to obtain the motion patterns of the NWP model precipitation outputs for the past nΔt and for the future up to the forecast time t i . 
         [0040]    Then, two corresponding pixels of the radar pattern and NWP output, established by τ P (x,t 0 ), are backward advected from t 0  with their own motion vectors. Now, the time-dependent phase errors ψ P (x,t) are derived by comparing positions of corresponding advected pixels at t=t 0 −Δt, . . . , t 0 −nΔt. At the forecast time, t i , the NWP model forecast fields ψ M  at t=t i  are corrected with ψ P (x,t i ) with each pixel of ψ M (x, y, t i ) tracked from its position at t 0  by following the motion of the NWP model precipitation patterns. 
       3. Correction of Lagrangian Intensity Errors (LIE), τ A (x,t). 
       [0041]    After adjustment of the phase errors, whatever residual differences remain between the intensity of precipitation given by radar and model are classified as NWP false alarms, misses, and amplitude errors. False alarms are pixels for which NWP predicts precipitation but none is observed above the detectable threshold; misses are pixels for which NWP fails to predict the observed precipitation; amplitude errors τ A  are the differences in intensity of precipitation for each pixel for which NWP predicts precipitation and precipitation is observed at that pixel. 
         [0000]      τ A ( x,t )=ψ R   [x+τ   p ( t ), t]−ψ   M   [x,t]   
         [0000]    where t=t 0 , . . . , t 0 −nΔt. 
         [0042]    False alarms are replaced with the value of no rain. Similar to the phase correction, we correct ψ M (x,y,t i ) with the time dependent τ A (x,t). Again, all corrections are applied following the motion of NWP precipitation patterns. 
         [0043]    The corrected ψ M (x,y,t i ) maintains the growth and decay as well as the motion of precipitation patterns at a future time that rely on the performance of the NWP models. 
         [0044]    It will be appreciated that a model of the motion pattern of short-term past observed weather can be used to perform short-term prediction of weather. This technique is presently implemented for precipitation nowcasting at McGill University in a tool called MAPLE. In a further embodiment of the invention, such short-term prediction or nowcasting technique can be combined with the short-term prediction obtained by using the error model and the NWP&#39;s short-term predicted values as described hereinabove to yield a balanced or weighted combination of the two nowcasting techniques that may be more reliable under some conditions.
The present invention provides a short term weather forecaster for nowcasting using a Numerical Weather Predictor (NWP). The system of the present invention tracks the evolution of differences between NWP and radar based precipitation patterns and adjusts the NWP forecast to account for theses differences. These differences are due to amplitude and phase errors of the NWP, such as model misses, false alarms, intensity errors and position errors. In presuming persistency in time of these errors and in estimating precipitation pattern time evolution due to other weather conditions such as wind motion, the present system corrects the NWP short term predicted patterns to compensate for these errors, thus enhancing nowcasting accuracy.