STAM-LSGRU: a spatiotemporal radar echo extrapolation algorithm with edge computing for short-term forecasting

With the advent of Mobile Edge Computing (MEC), shifting data processing from cloud centers to the network edge presents an advanced computational paradigm for addressing latency-sensitive applications. Specifically, in radar systems, the real-time processing and prediction of radar echo data pose significant challenges in dynamic and resource-constrained environments. MEC, by processing data near its source, not only significantly reduces communication latency and enhances bandwidth utilization but also diminishes the necessity of transmitting large volumes of data to the cloud, which is crucial for improving the timeliness and efficiency of radar data processing. To meet this demand, this paper proposes a model that integrates a spatiotemporal Attention Module (STAM) with a Long Short-Term Memory Gated Recurrent Unit (ST-ConvLSGRU) to enhance the accuracy of radar echo prediction while leveraging the advantages of MEC. STAM, by extending the spatiotemporal receptive field of the prediction units, effectively captures key inter-frame motion information, while optimizations to the convolutional structure and loss function further boost the model’s predictive performance. Experimental results demonstrate that our approach significantly improves the accuracy of short-term weather forecasting in a mobile edge computing environment, showcasing an efficient and practical solution for processing radar echo data under dynamic, resource-limited conditions.


Introduction
The precision and timeliness of weather forecasting are crucial for addressing extreme weather events, agricultural production, aviation safety, and many other domains.Radar echo extrapolation, as a vital weather prediction technique, provides essential information on short-term weather changes.However, the effectiveness of this method greatly depends on the accurate capture and analysis of spatiotemporal features in radar data [1,2].
Traditional radar extrapolation methods primarily rely on linear or simple mathematical models to predict weather patterns, which often perform poorly in handling complex weather systems [3].With the advancement of MEC and Artificial Intelligence (AI) technologies, new solutions have emerged for radar echo extrapolation.The low-latency characteristic of MEC allows for rapid processing of data near its point of origin, while AI, particularly deep learning technologies, demonstrate immense potential in analyzing large-scale, complex datasets [4,5].
Firstly, MEC plays a pivotal role in processing radar data.Traditionally, radar data required transmission to remote servers for processing and analysis, which was not only time-consuming but could also lead to data delays [6].MEC significantly reduces data transmission time by providing computational resources near the data source, thus accelerating data processing [7].This type of near-source processing is particularly well-suited for weather forecasting, as it necessitates rapid response and real-time analysis [8].Secondly, AI technologies, especially machine learning and deep learning, have proven highly effective in interpreting radar data and enhancing forecast accuracy.Deep learning models can learn from historical weather data and predict future changes in weather patterns [9].These models are particularly adept at handling large volumes of radar data and extracting meaningful insights, aiding meteorologists in making more accurate predictions [10].
In recent years, there has been growing interest in developing algorithms to infer radar echoes beyond the instrument range and to forecast the evolution of echoes over time [11][12][13][14].Currently, weather forecasting methods can be broadly categorized into two main approaches: numerical weather prediction (NWP) methods and radar echo extrapolation.NWP methods utilize fluid dynamics and thermodynamic laws to simulate the physical processes of the lower atmosphere, providing predictions based on complex physical state equations and supercomputers [15].While NWP methods offer valuable insights, they face challenges such as prediction delays, low resolution, and limitations in forecasting sudden severe weather events [16,17].On the other hand, radar echo extrapolation methods, such as artificial neural networks [18], support vector machines [19], and decision trees [20], leverage radar data to understand the relationship between radar echoes and other variables, enabling the prediction of future weather conditions.These data-driven methods have gained attention in recent years due to the availability of large amounts of historical data and have shown superior performance in various fields [21,22].
The purpose of radar echo extrapolation is to predict future radar echo maps for a specific area based on previously observed radar echoes.This prediction task poses significant challenges as it requires spatiotemporal modeling of radar data to accommodate high resolution, thereby rendering it a spatiotemporal forecasting problem.Convolutional Neural Networks and Recurrent Neural Networks have been extensively utilized in such spatiotemporal forecasting tasks [23,24].However, existing models still confront challenges in handling high spatiotemporal resolution and complex non-stationary information, particularly in the context of convection formation and dissipation.Moreover, theoretical models for radar echo extrapolation are capable of generating prediction sequences of any length.Yet, as the prediction length increases, error accumulation can lead to image blurriness and loss of details.
In this context, the integration of MEC and AI offers a new perspective for radar echo extrapolation.The core advantage of MEC lies in its low-latency characteristics, enabling rapid processing near the data generation point, which is particularly crucial for real-time radar data analysis.This capability for rapid response, coupled with advanced abilities in handling and analyzing large-scale complex datasets, provides robust support for enhancing the efficiency and accuracy of radar data processing.Therefore, this paper proposes a STAM-LSGRU network to address key challenges in radar echo extrapolation, including error accumulation and effective extraction of high-order non-stationary information.The main contributions of this paper can be summarized as follows: and Heidke Skill Score (HSS) evaluation metrics, improvements are made to the loss function, reducing the ambiguity and distortion of the prediction results, and enhancing predictive performance in heavy rainfall regions.
The remainder of this article is organized as follows: Related Work, Methodology, Experiments, and Conclusion.The Related Work section provides an overview of previous studies, highlighting existing methods and findings in the field.The Methodology section details the theoretical frameworks and techniques used in the research, outlining the design of the proposed model.The Experiments section introduces the experimental setup, dataset description, and obtained results, offering empirical validation and comparison.Finally, the Conclusion section summarizes the contributions and proposes directions for future research.

Related work
The MEC with AI technologies demonstrates significant potential in fields such as radar echo extrapolation and weather forecasting.As an emerging computational paradigm, MEC shifts computational tasks from the cloud to the network edge, achieving low-latency and highefficiency data processing.For instance, zhou [25] noted that 'Edge Intelligence' is a product of the convergence of MEC and AI.This concept aims to provide superior solutions for key issues in edge computing.It also explores how to establish AI models on edge devices, including model training and inference processes.This approach exemplifies the innovative strides being made in combining AI with edge computing, to optimize computational efficiency and enhance the capabilities of edge devices in processing complex tasks.
MEC technology enhances data processing efficiency by providing computational resources at the network edge, thereby significantly reducing latency and bringing processing closer to the data source.The integration of AI technologies further elevates MEC's data processing capabilities.AI algorithms, particularly deep learning models, have demonstrated exceptional performance in areas like image recognition, pattern detection, and predictive analytics.Al-Habob and Dobre [26] explored the symbiotic relationship between MEC and AI, highlighting AI's critical role in the MEC offloading process, such as resource management and scheduling.Huang [27] proposed an infrastructure for executing machine learning tasks on MEC servers, assisted by Reconfigurable Intelligent Surfaces.Deng [28] discussed the role of AI with software orchestration and hardware acceleration in reducing edge computing latency.Yazid [29] provided a comprehensive review of Unmanned Aerial Vehicles (UAVs) in the application of MEC and AI, exploring their role in enhancing the efficiency of IoT applications.Dahmane [30] introduced a blockchain-based AI paradigm for secure implementation of UAVs in MEC.Wang [31] surveyed the convergence of MEC, Metaverse, 6G wireless communications, AI, and blockchain, and their impact on modern applications.Chakraborty and Sukapuram [32] examined the application of MEC in urban informatics, emphasizing its contribution to the development of smart cities.
Recently, deep learning-based radar echo extrapolation models have been proposed that are more accurate than traditional methods.In 2015, Shi [21] introduced the Convolutional Long-Short Term Memory (ConvL-STM) model for precipitation nowcasting.This model is designed to handle time series data with spatial structures and uses convolutional operations instead of the Hadamard product in the FC-LSTM [33].In 2016, Shi [34] improved the ConvGRU model and proposed the TrajGRU model, which can dynamically learn the network recursive structure.Wang [35,36] proposed Pre-dRNN and PredRNN++ models based on ConvLSTM.The LSTM unit was rebuilt by the team, and the Spatiotemporal LSTM (ST-LSTM) unit was created, allowing memory state to propagate in both vertical and horizontal directions and no longer being restricted to each individual LSTM unit.They developed the Gradient Highway Units (GHU), which were added between the first and second layers of the model at each time step, and the Causal LSTM unit in subsequent research.This greatly shortened the gradient propagation path and solved the issue of information loss in long-term predictions.Due to the relatively simple state transition function and the ineffective differential signal processing in the majority of RNNs used for spatiotemporal prediction, it is challenging for the model to learn complicated spatiotemporal changes.They proposed the MIM structure [37], which is used to extract stationary features and non-stationary features (MIM-S layer and MIM-N layer, respectively).The model achieved better performance on radar datasets.Lin [38] proposed the Sa-ConvLSTM, which added a self-attention mechanism at the output end of ConvL-STM.By using an additional memory unit M and a selfattentive feature aggregation mechanism, it computed pairwise similarity scores to fuse the previous features that contain the global spatial receptive field.Wu [39] has made further advances in the utilization of spatiotemporal information by proposing the MotionRNN architecture and designing the MotionGRU unit.This unit can model transient changes and motion trends in a unified manner, and a new motion highway has been introduced, which significantly enhances the ability to predict variable motion and avoids the issue of vanishing motion when stacking multiple prediction models.Chang [40] proposed a spatial-temporal residual prediction model applicable for high-resolution video prediction.The model employs a spatial-temporal encoding-decoding scheme to capture complex motion information in highresolution videos.Jin [41] proposed a novel spatiotemporal graph neural network model called BGGRU, which integrates spatial and temporal information to explore the temporal patterns and spatial propagation effects of time series, aiming to enhance prediction accuracy.However, in the above methods, the global spatiotemporal dependencies of radar echoes have not been fully explored.This paper analyzes existing spatiotemporal prediction models and proposes an STAM module that addresses the issue of error accumulation.Additionally, the convolutional structure and loss function of the basic unit are improved, resulting in more accurate predictions of high-echo regions at different scales.

Methodology
The task of radar echo extrapolation aims to learn the mapping from input sequences to a latent space.To achieve this objective, we constructed a convolutional recurrent neural network, STAM-LSGRU.As shown in Fig. 1, its recurrent connections bestow it with memory functionality, enabling the capture and storage of previously inputted information.This memory capability allows the network to consider the information of the entire sequence comprehensively, rather than being limited to the inputs at the current timestep.Within STAM-LSGRU, this memory mechanism is crucial, allowing the network to effectively learn patterns and regularities within sequence data, thereby achieving the mapping from input sequences to latent space.By stacking three ST-ConvLSGRUs and one STAM-LSGRU, an encoder-decoder network is formed, which is currently the mainstream method for spatiotemporal sequence prediction.The ST-ConvLSGRU integrates the temporal information flow processing capabilities of the conventional GRU with a newly added spatial memory flow propagation mechanism, achieving simultaneous capture of temporal and spatial dependencies.The STAM-LSGRU predicts the next time step's radar echo image, not solely relying on the output of the previous time step, allowing for better handling of long sequence inputs and outputs.At a single time step, the vertical arrows represent the direction of memory and state updates along the spatial dimension, while the horizontal arrows represent the direction of updates along the time dimension.The spatiotemporal memory M is transferred from the lowest layer of the recurrent unit to the highest layer within a single time step, and is then transferred to the lowest layer of the following time step in a "Z"-shaped direction, first along the spatial dimension and then along the time dimension.The input radar data is downsampled by a factor of four and goes through three layers of ST-ConvLSGRU for information extraction and transformation before being input into the STAM-LSGRU.This enables the model to focus on past input states, avoiding error accumulation.The output is then upsampled to obtain the final prediction results.In order to improve the prediction of high echo areas, this paper proposes an enhanced loss function.In addition, the gating mechanisms has been optimized with an inception module for all the basic units of RNNs.Experimental results demonstrate that employing the STAM-LSGRU network leads to a significant improvement in prediction accuracy, thereby enabling more precise forecasting of future echo image sequences at different time points.

ST-ConvLSGRU
Inspired by ST-LSTM [35], this paper introduces the concept of LSTM memory units into the ConvGRU model, resulting in the ST-ConvLSGRU model depicted in Fig. 2 This model serves as the foundation for subsequent improvements.The original ST-LSTM model utilized a dual LSTM structure to process images with high spatiotemporal resolution effectively, storing and transmitting spatiotemporal information within and outside the memory cells.However, this structure tends to increase the model's complexity and the number of parameters, leading to overfitting issues.To address this challenge, this work integrates LSTM memory units into ConvGRU, aiming to reduce the model's parameter count while maintaining effective spatiotemporal information processing.This enhanced ST-ConvLSGRU model, compared to ST-LSTM, not only reduces the model's complexity but also enhances its capability to handle spatiotemporal information, avoiding overfitting and thus enabling more accurate predictions.As shown in Fig. 2, Z t , R t and X t are update gates, reset gates and input states, respectively.ht , i and f are new information, input gates, and forget gates, respectively.g is used as a temporary variable to update M. t denotes the t th time step and l denotes that the loop cell is located at the l th level of the stacked structure.
For a single ST-ConvLSGRU unit at time t, if the unit is located in the first layer (i.e., when l=1), the input state X l t is a tensor converted from the radar echo map input at the current time.If the unit is not in the first layer (i.e., l >1), the hidden state H l−1 t output at time t is used as the input state X l t for the unit.The ST-ConvLSGRU unit first passes the input state X l t and the hidden state H l t−1 output by the unit in the same layer at time t-1 through a gating structure.Two different convolution filters are applied to obtain the reset gate R t , the update gate Z t , and the new information vector ht .The calculation method for R t and Z t is consistent with that of ConvGRU, and is shown as follows, '*' denotes the convolution operation: Similar to LSTM, the input state X l t and spatiotemporal memory M l t−1 are fed into a gated structure.Three different convolution filters are applied to obtain the forget gate f t , input gate i t , and input modulation gate g t .The forget gate f t and element-wise Hadamard product are used to forget unimportant features of past time steps in the temporal memory M. The input gate and input modulation gate are used to update the features in memory through element multiplication, resulting in the updated spatiotemporal memory M l t .The ' • ' represents the Hadamard product.This process can be represented as follows: Next, the input state X l t , hidden state H l−1 t , and updated spatiotemporal memory M l t are convolved using convolution filters to obtain new information ht .The hid- den state H l t is updated using the reset gate R t and update gate Z t , which extract abundant and rich spatiotemporal features through two gating mechanisms.As a result, the extrapolation network is able to accurately model the motion of radar echoes and make precise predictions on whether they will continue to expand or dissipate in the future.This process can be represented as follows:

Spatiotemporal attention memory module
The extrapolation of radar echoes can also be classified as a regression problem.The extrapolation model theoretically has the ability to generate prediction sequences (1) of arbitrary lengths.However, as the prediction length increases, there is a strong interdependence between adjacent frames, resulting in cumulative errors, which leads to blurred and distorted extrapolated images with missing details.To address this issue and enable the model to review the historical input sequence at each predicted time step, a STAM is constructed to utilize the input H l t of the l th layer, to recall the historical input X l h .The model can adaptively learn the mapping from X 0:n to X n+1:T based on a rich history of data.
The specific implementation is illustrated in Fig. 3, and its design is inspired by the dual-attention mechanism, which is named as STAM.STAM receives three inputs, including the predicted results of past time series H l t−τ :t−1 , the multi-layer results H l−τ :l−1 t of the current time series, and the low-level results H l t of the current time series, where τ represents the step size.STAM consists of two modules, the attention module and the fusion module.The convolutional layer passes the current hidden state H l t ∈ R C×H ×W , and generates the query Q ∈ R C×H ×W through a 1x1 convolution operation, where C, H, and W represent the number of channels, the length and width of the input data, respectively.Similarly, the key value K t ∈ R τ ×C×H ×W and value V t ∈ R τ ×C×H ×W can be obtained from the predicted results of the past time series H l t−τ :t−1 ∈ R τ ×C×H ×W through two independent convolutions.The weight matrix A t is obtained by mul- tiplying Q and K t , and then applying sum and softmax operations: Subsequently, the new temporal state T l t can be calculated according to the formula of temporal attention: Finally, the reshaped T l t is resized to the same size as the original hidden state, and is used as the input of the fusion module.
The above approach enables adaptive learning of the historical input X l t in the temporal dimension.In order to address the problem of information loss during the propagation process from the low-level to high-level layers, the output of each convolutional neural network layer is kept in the multi-layer state.Then, the top-level hidden state H l t is used to recall H l−τ :l−1 t and generate a new spatial hidden state S l t as the input of the fusion module. ( The computation process is akin to temporal attention, whereby the current hidden state H l t ∈ R C×H ×W can be transformed into a query Q ∈ R C×H ×W by reshap- ing the convolutional layer.Subsequently, the keys K s ∈ R τ ×C×H ×W and values V s ∈ R τ ×C×H ×W can be generated via two independent 1x1 convolutions using H 0:l−1 t ∈ R τ ×C×H ×W .The weight matrix A s is obtained by multiplying Q and K s , followed by summation and softmax operations.
Subsequently, the new spatial state S l t can be calculated according to the formula for spatial attention: Finally, the reshaped S l t is resized to the same size as the original hidden state and serves as the input to the fusion module.
The fusion module aggregates the temporal state T l t and the spatial state S l t , and uses gating mechanisms to control the output of the current time sequence.First, T l t and S l t are concatenated along the channel dimension, and the number of channels is adjusted through convolutional operations to obtain the fusion features: Subsequently, to effectively control the fusion of historical attention information and the current hidden state H l t , two gating mechanisms are used: (7) The fusion features G are used to generate the input gate e i and the forget gate e f through convolution, and the output of the STAM, represented as H l t , is given by: In the STAM-LSTM unit, the STAM will be embedded into the RNN unit shown in Fig. 4, forming the STAM-LSGRU.

Convolutional Inception optimization
The radar echo image exhibits varying strengths and sizes of echoes in different regions, and the commonly used 5x5 convolution is insufficient for capturing the multiscale radar echoes and high-order non-stationary (10) Fig. 3 The STAM consists of two components: the attention module and the fusion module This capacity for information aggregation enables the network to better integrate abstract features at different levels, further improving its understanding of image content.
As shown in Fig. 5, the enhanced convolutional network structure comprises three branches, each undergoing 1x1 convolution, 3x3 convolution, and two consecutive 3x3 convolution operations, respectively.By concatenating two consecutive 3x3 convolutions, our structure not only achieves the same receptive field as a single 5x5 convolution but also significantly reduces the parameter count compared to the latter, offering a more efficient computational approach.Furthermore, the improved Inception convolutional structure effectively captures features at different scales by synthesizing convolution kernels of different sizes, thereby demonstrating superior performance in capturing image features compared to a single-size 5x5 convolution.This design optimizes parameter usage, reduces computational burden, enhances the network's ability to capture multi-scale information in images, and improves the model's expressive power in handling complex image tasks.

Loss function optimization
Currently, most of the existing deep learning radar echo algorithms use mean squared error (MSE) loss as a loss function, which is a common loss function used to evaluate the difference between the predicted value and the real value of the model, and is suitable for regression problems.The smaller the value, the smaller the difference between the predicted and true results.The calculation of the MSE is as follows: In this equation, H and W are the length and width of the radar image, and MSE loss is the sum of the squared errors for each pixel in every extrapolated radar image ŷ and its corresponding true image y.However, since MSE is sensitive to outliers and penalizes large prediction errors, using MSE as the loss function on datasets with outliers can be influenced by those outliers.In practical images, noise or other interference may exist, which could greatly affect MSE and thereby impact the model's predictive ability.MSE only considers the difference between the predicted and actual values of each pixel, without considering the correlation between pixels.In image prediction, there is usually some correlation (12 5 The structure of the proposed convolutional inception optimization between pixels, and ignoring this correlation may lead to a decrease in prediction performance.Therefore, to improve upon MSE, commonly used meteorological indicators such as the CSI and the HSS can be incorporated into the evaluation metrics.The improved loss function is as follows: Considering that the CSI and the HSS may not be differentiable, they can be made differentiable by applying the sigmoid function, and incorporated into the final differentiable loss function.

Dataset
This paper utilizes the data from the China Central Meteorological Observatory's radar network in the eastern region of China from 2020 to 2022, with a spatial resolution of 0.01 • and a temporal resolution of 6 minutes.The radar data is sliced at the central point to achieve a spatial resolution of 400x400.dBZ represents the radar echo value, with larger dBZ values indicating a higher likelihood and intensity of severe convective weather.Atmospheric motion exhibits periodicity, particularly in the same region, with many similar samples in both real-time observations and model forecasts, leading to overfitting.Additionally, severe convective weather occurs only a few days out of the year, thus filtering out some relatively mediocre data is necessary.Ultimately, 10,000 sequential samples are selected, with 6,000 sequences serving as the training set, 2,000 sequences as the validation set, and 2,000 sequences as the test set.The radar echo values range from 0 to 70, and normalization is employed to distribute the values between 0 and 1, facilitating better model convergence.

Implementation
For all experiments, a Nvidia GeForce RTX 3090 GPU was used for training.The default hyperparameters and experimental setup were as follows: the model was trained using a batch size of four image sequences, the Adam optimizer with an initial learning rate of 0.001, and momentum decay set to 0.90.The four-layer model was configured with 64 channels, and a total of 70,000 training steps were performed.At every 5000 training steps, evaluation metrics were recorded for both the training and validation sets.During training, the model was designed to predict the next 10 time sequences, with the attention mechanism of the STAM module focusing on a time step of 5.All experiments were conducted using the same set of hyperparameter values to ensure consistency and comparability.To enhance training performance (13) and explore different strategies, various techniques were employed.These included the use of teacher forcing, which involved providing the correct sequence as input during training, and bidirectional training, where the model was trained in both forward and backward directions.To prevent overfitting, early stopping techniques were utilized.Specifically, training was terminated if the validation loss did not decrease for 10,000 consecutive steps, indicating that the model's performance had plateaued.This ensured that the model was not trained excessively and retained its generalization ability.

Evaluation indicators
This paper evaluates the effectiveness of various models based on MSE, CSI, HSS, and Structural Similarity (SSIM).MSE measures the average difference between model predictions and true values, reflecting the model's precision in predicting future radar echoes.CSI assesses the model's detection capability for precipitation events, taking into account false alarms and missed detections.HSS compares the correctness of the model's predictions to random forecasting, used to measure the predictive capability of the model.The SSIM is used to evaluate the degree of similarity between the model-predicted radar echo images and the actual radar echo images, considering aspects such as luminance, contrast, and structure.Considering the correlation between radar echo values and actual weather, three thresholds of 20dBZ, 35dBZ, and 45dBZ were chosen to evaluate the efficacy of the algorithm for radar extrapolation.Under these thresholds, radar echo images were binarized by designating a value of 1 if the echo value exceeded the threshold and a value of 0 otherwise.TP represents a predicted event that actually occurred, FP represents a predicted event that did not occur, FN represents an unpredicted event that did occur, and TN represents a predicted event that did not occur.Following are the formulas used to calculate the evaluation metrics in this paper: where µ x and µ y are the means of x and y, respectively.σ 2 x and σ 2 y denote the variances.σ xy is the covariance of x and y. c1 and c2 are constants.

Results and analysis
Currently, radar echo extrapolation models are primarily based on stacking multiple layers of basic Convolutional Recurrent Units.There is no fixed standard for the number of layers in the stacked radar echo extrapolation (14) model, as it is contingent on variables such as data volume, data complexity, network structure, and hardware.
Increasing the number of layers can improve the network's expressive power and the model's capacity to represent and abstract data, but it may also increase the network's computational and storage burden, resulting in overfitting and other issues.The number of layers is typically determined by the extent and complexity of the dataset as well as the training effect of the network.For smaller datasets and relatively straightforward problems, a shallower network structure may be optimal, whereas a deeper network structure may be preferable for larger datasets and more complex problems.Using experimental methods, this paper determines the optimal number of stacked layers to avoid overfitting and underfitting issues.Table 1 presents the MSE scores of various models at different numbers of layers.It is observed that as the number of layers increases, the MSE value of each model typically decreases to reach a minimum value before starting to rise again.For most models, the MSE value reaches its minimum when the number of layers is four, indicating that the models perform best at this depth.As the depth of the model increases, it is able to learn more complex features and deeper data representations.There comes a point where this learning capability is optimized, and further increases in the number of layers make the model more complex, increasing the number of parameters.This complexity can lead to gradients gradually vanishing or exploding during the backpropagation process, making it difficult to train the model.The STAM-LSGRU model achieves its best performance at four layers; when the number of layers is less than or greater than four, the model's performance metrics decrease.Therefore, this paper sets the number of layers of the STAM-LSGRU model to four.
Through ablation experiments and comparative experiments, we validated the best performance of radar echo extrapolation at different thresholds.The results are shown in Tables 2 and 3, where CSI and HSS scores were calculated at the thresholds of 20 dBZ, 35 dBZ, and 45 dBZ, while MSE and SSIM were averaged across all thresholds.The symbol ↑ indicates that higher values indicate better performance in radar echo extrapolation, while ↓ indicates the opposite.
The ablation experiment results are presented in Table 2, where evaluation metrics of the original ST-Con-vLSGRU model and various improvement methods are compared.Across all metrics and thresholds, ST-Con-vLSGRU-1, STAM-LSGRU-0, and STAM-LSGRU-1 outperform ST-ConvLSGRU.Specifically, STAM-LSGRU* demonstrates improvements of 6.87%, 6.45%, 5.8%, and 7.7% in CSI, HSS, SSIM, and MSE, respectively, compared to the ST-ConvLSGRU network.Figure 6 illustrates the radar echo extrapolation results in order to visually compare the experimental outcomes of various enhancement techniques.Indicating that these modules can enhance the performance of intensity prediction, the prediction results of STAM-LSGRU* are clearer, have more distinct edges, and pay closer attention to high echo regions.As shown in Table 3, different evaluation metrics have been improved.Compared to the state-of-the-art MotionRNN model, the CSI score has increased by an average of 1.6%, the HSS score has increased by 1.1%, the SSIM score has increased by 2.7%, and the mean squared error has increased by 3.2%, under different thresholds.Figure 7 displays the hourly scores of each model in the prediction of the next hour.The models exhibit similar performances, but as time progresses, the scores of other models sharply decline, whereas STAM-LSGRU shows a more gradual decrease in scores.
For a more intuitive observation, a visual example is also shown in Fig. 8. From Fig. 8, it can be observed that the STAM-LSGRU model designed in this research outperformed the other five models.The ConvLSTM model produced smoother results than other methods and suffered from severe detail loss and prediction errors in the high reflectivity regions of radar images.TrajGRU and PredRNN had poor performance in predicting the central echo region and also suffered from distortion to some extent.Since radar image changes are a high-order nonstationary process, these methods were able to effectively predict the radar motion trend.MIM and STAM-LSGRU better captured the overall changing trend characteristics of the echo region, but STAM-LSGRU captured the high and low echo characteristics better than MIM.The prediction results of MotionRNN and STAM-LSGRU were similar, but the predicted image blocks of STAM-LSGRU were closer to the actual observation results.From the edge of the echo main body, its predicted image edges were more in line with the actual ones, and its blurring degree was lower, showing greater consistency with the real image.

Conclusion
This research introduces a neural network-based radar echo extrapolation algorithm named STAM-LSGRU.By deploying the STAM-LSGRU model in an edge computing environment, we not only achieve enhanced real-time data processing capabilities but also significantly reduce data transmission delays.Compared with traditional radar echo extrapolation algorithms and other deep learning-based algorithms, STAM-LSGRU exhibits markedly improved predictive performance in complex environments, particularly in heavy rain areas.This paper designs STAM to capture reliable inter-frame motion information by expanding the temporal and spatial receptive fields of the prediction units.The convolutional structure and loss function of the basic unit have been improved to  3 The predictive performance of different models on radar echo data was compared.All models were trained using 5 time series and evaluated based on their ability to predict the next 10 time series.The evaluation metric was the average of the 10 predicted time series enhance the robustness of model predictions.Compared to the MotionRNN model, the CSI score has increased by an average of 1.6%, the HSS score by 1.1%, and the SSIM score by 2.7%.In the future, we plan to further advance meteorological forecasting by integrating more observational data and model outputs, aiming to improve the accuracy and timeliness of weather predictions.With the continuous advancements in MEC and AI technologies, along with the increasing abundance of meteorological observation data, we anticipate that the STAM-LSGRU model will demonstrate higher predictive capabilities in a wider range of meteorological scenarios, bringing new breakthroughs to the field of weather forecasting.

Fig. 4
Fig.4 The STAM-LSGRU unit incorporates the STAM module and optimized convolutional patterns on the basis of the ST-ConvLSGRU

Fig. 6
Fig. 6 Comparison of results of different improved versions based on ST-ConvLSGRU

Fig. 7 Fig. 8
Fig. 7 Extrapolated objective indicators of 10 time series of different models

Table 1
MSE score of 10 frame advance prediction for different depth models

Table 2
The improved components were evaluated through a series of ablation experiments.The ST-ConvLSGRU served as the baseline model, while the ST-ConvLSGRU-1 incorporated an inception optimization.The STAM-LSGRU introduced an STAM module, and the loss function of STAM-LSGRU-1 was further optimized.Finally, the STAM-LSGRU* model integrated all the improvement strategies