删除或更新信息,请邮件至freekaoyan#163.com(#换成@)

Characteristics of Fengyun-4A Satellite Atmospheric Motion Vectors and Their Impacts on Data Assimil

本站小编 Free考研考试/2022-01-02

Yaodeng CHEN1,,,
Jie SHEN1,
Shuiyong FAN2,,,
Deming MENG1,
Cheng WANG1

Corresponding author: Yaodeng CHEN,keyu@nuist.edu.cn;
Corresponding author: Shuiyong FAN,syfan@ium.cn;
1.Key Laboratory of Meteorological Disaster of Ministry of Education/Collaborative Innovation Center on Forecast and Evaluation of Meteorological Disasters, Nanjing University of Information Science and Technology, Nanjing 210044, China
2.Institute of Urban Meteorology, China Meteorological Administration, Beijing 100089, China
Manuscript received: 2020-04-19
Manuscript revised: 2020-07-10
Manuscript accepted: 2020-07-27
Abstract:The high observation efficiency, scanning speed and observation frequency of the Fengyun-4A (FY-4A) satellite indicates the progress of Chinese geostationary meteorological satellites. The characteristics of FY-4A atmospheric motion vectors (AMVs) derived from the high-level water vapor (WV-High) channel, mid-level water vapor (WV-Mid) channel, and infrared (IR) channel of FY-4A are analyzed, and their corresponding observation errors estimated. Then, the impacts of single-channel and multi-channel FY-4A AMVs on RMAPS-ST (the Rapid-refresh Multi-scale Analysis and Prediction System—Short Term) are evaluated based on one-month data assimilation cycling and forecasting experiments. Results show that the observation errors of FY-4A AMVs from the three channels have an explicit vertical structure. Results from the cycling experiments indicate that the assimilation of AMVs from WV-High produces more apparent improvement of the wind in the upper layer, while a more positive effect in the lower layer is achieved by the assimilation of AMVs from IR. Furthermore, the assimilation of AMVs from IR is more skillful for medium and moderate precipitation than from other channels owing to the good quality of data in the lower layer in the AMVs from IR. Assimilation of FY-4A AMVs from the three channels could combine the advantages of assimilation from each individual channel to improve the wind in the upper, middle and lower layers simultaneously.
Keywords: data assimilation,
FY-4A satellite,
atmospheric motion vector,
observation error
摘要:本研究选取FY-4A卫星的大气运动矢量产品(AMVs)作为研究对象,在对各通道产品数据进行特征分析的基础上,计算得到了各通道产品的观测误差,并基于睿图业务模式(RMAPS-ST)开展了为期一个月的同化及预报试验。特征分析表明:FY-4A中各通道AMV产品的观测误差垂直结构特征显著。各通道的单通道同化试验表明,水汽通道的AMVs对高层风场的改善相对明显;红外通道由于存在部分低层风信息,在滚动预报一段时间后逐渐改善低层风场和水汽的预报,从而使得降水中心的降雨强度预报更接近观测。AMVs多通道同时同化试验看出,多通道联合同化,可以综合三个通道的各自特点,改进对大多数高度层的风场预报。
关键词:资料同化,
FY-4A卫星,
大气运动矢量,
观测误差





--> --> -->
1. Introduction
Atmospheric motion vectors (AMVs) are wind observations derived from consecutive satellite images by tracking features such as cirrus edges, gradients in water vapor, or small cumulus clusters (Nieman et al., 1993; Velden et al., 1997). They mainly describe the wind field in the middle and upper troposphere, which could effectively make up for the deficiency of the upper observation data. AMVs have been investigated since the 1970s and applied in operational numerical weather prediction (NWP) systems since the 1990s (Schmetz et al., 1993; Menzel, 2001). Various studies have pointed out that the assimilation of AMVs can improve the atmospheric dynamics condition, which in turn increases the accuracy of NWP (Wang et al., 2006; Goerss, 2009; Berger et al., 2011; Joo et al., 2013; Wu et al., 2014).
Compared to other wind measurements, AMVs from geostationary satellites can provide superior temporal and spatial coverage in the middle and upper troposphere (Kaur et al., 2015), due to the capability of geostationary satellites in achieving frequent and continuous observations over a specific region (Honda et al., 2018). Various studies on the assimilation of AMVs derived from different satellites have been carried out to take advantage of such exceptional temporal and spatial overage, and positive results obtained. For example, AMVs from the Geostationary Operational Environment Satellite (GOES) have been assimilated in hurricane research, and the forecasts of hurricane tracks and intensity improved (Soden et al., 2001; Velden et al., 2017; Li et al., 2020). The assimilation of Kalpana-1 AMVs also reduced the initial position errors and improved track forecasts of tropical cyclones (Deb et al., 2011). Besides, studies on the assimilation of AMVs from MTSAT (Multi-function Transmission Satellite) and Himawari-8 showed positive impacts on the timing and intensity of rainfall forecasts, and even flood risk prediction (Yamashita, 2012; Otsuka et al., 2015; Kunii et al., 2016).
In addition to the satellites mentioned above, Fengyun (FY) geostationary satellites from China have also provided AMV products for operational NWP systems and the research community. Since the launch of the first Chinese geostationary meteorological satellite, FY-2A, in 1997, studies and applications of AMVs from Chinese geostationary meteorological satellites have been gradually increasing (Feng et al., 2008; Ren et al., 2014; Wan et al., 2018). Compared to the FY-2 series, the spatial and temporal resolutions have been dramatically improved in the FY-4A satellite, and two visible channels and seven infrared (IR) channels included. Moreover, it can produce full-disk imagery every 15 minutes, and the temporal resolution can be increased to up to 5 minutes for the coverage of China (Dong, 2016; Zhang et al., 2017; Wang and Shen, 2018; Zhao et al., 2019). The capability of obtaining AMVs has been greatly increased. However, studies on the AMV products of FY-4A are relatively rare (Wan et al., 2019). Therefore, in order to further understand the characteristics of FY-4A AMVs, and to promote the application of FY-4A AMVs in NWP, the characteristics of AMV products are analyzed in this study, and their impacts on data assimilation discussed.
The rest of this paper is organized as follows: Section 2 analyzes the characteristics of FY-4A AMVs, and the observation error of AMV products is estimated in section 3. Section 4 describes the framework of the assimilation and forecasting experiments, and results are discussed in section 5. Section 6 provides a summary of the study’s key findings.

2. FY-4A AMVs and their characteristics
2
2.1. AMV data from FY-4A
--> The monitoring and prediction of severe weather systems is one of the most important and complicated parts of operational weather forecasting. The Advanced Geosynchronous Radiation Imager (AGRI) onboard the FY-4A geostationary meteorological satellite plays a critical role in the monitoring of severe weather systems in China and its surrounding areas. There are 14 channels in AGRI, and AMV products used in this study are derived from channel 09 (λ = 6.5 μm), channel 10 (λ = 7.2 μm) and channel 12 (λ = 11.0 μm); namely, the high-level water vapor (WV-High) channel, the mid-level water vapor (WV-Mid) channel, and the IR channel, respectively. According to the documentation for FY-4A AMV products, the spatial resolution corresponding to the sub-satellite point is 64 km, and products with 3-h resolution can be provided. The horizontal distribution of FY-4A AMVs from WV-High, WV-Mid and IR at 0000 UTC 1 November 2018 are displayed in Figs. 1a-c. It is shown that the AMVs cover the West Pacific and Indian oceans, where conventional observations are sparse. Meanwhile, the horizontal distribution of AMVs from WV-High and WV-Mid are similar but different from that of IR. The IR channel supplements the AMV data of some regions where two water vapor channels are not observed. Figure 1d shows the calculation of the vertical profile of the total number of FY-4A AMVs for the three channels from 1?30 November 2018. The vertical distribution of AMVs derived from the two water vapor channels is obviously different from that from the IR channel, which is due to the difference in absorption characteristics of each channel. More specifically, WV-High and WV-Mid are in the water vapor absorption bands (Velden et al., 1997), and their different weighting function peak heights lead to the discrepancy in heights at which the channel AMV count is the most. In contrast, the IR channel is in the IR window band (Yang et al., 2014, Lu et al., 2017), so it can obtain AMVs of various heights, even data on the middle and lower troposphere.
Figure1. Horizontal distributions of AMVs derived from (a) WV-High, (b) WV-Mid, and (c) IR at 0000 UTC 1 November 2018, and (d) vertical profiles of the total number of single AMV data from the three channels in each layer for the period from 1?30 November 2018.



2
2.2. Quality indicator for AMVs
--> The quality indicator (QI) is a weighted average of several normalized QI components, which is calculated to represent the quality of each FY-4A AMV (Holmlund, 1998). The QI value of FY-4A AMV is provided by the National Satellite Meteorological Center of China of the China Meteorological Administration (CMA). It is calculated as an empirical function of a spatial and temporal consistency check of wind speed, wind direction, vector and horizontal wind components. The value decreases as the agreement of a given vector with its surroundings decreases. As a result, a larger QI value means a higher credibility (Velden et al., 1998; Deb et al., 2015).
Figure 2 shows the vertical distributions of FY-4A AMVs derived from the three channels from 1?30 November 2018 with different QI values. It can be seen that AMVs with different QI thresholds have similar vertical structure but the proportion of total amount is different. For example, AMVs with QI greater than 70 account for more than 80%, AMVs with QI greater than 80 make up about 60%, and AMVs with QI greater than 90 share only approximately 30% of the total data quantity.
Figure2. Vertical distributions of AMVs derived from (a) WV-High, (b) WV-Mid, and (c) IR for the period from 1?30 November 2018 with different QI values (the values in brackets represent the percentage of the total amount of data under this threshold).


Next, the wind speed bias (Fig. 3) and root-mean-square error (RMSE) (Fig. 4) of FY-4A AMVs in November 2018 are calculated against the Final Operational Global Analysis (FNL) for FY-4A AMVs with different QI thresholds. As displayed in Fig. 3 and Fig. 4, smaller speed bias and RMSE of AMVs are obtained with a larger QI value, which is consistent with the conclusion obtained by Kim and Kim (2018). However, as the primary source of FY-4A AMVs below 500 hPa, data from the IR channel with QI greater than 80 (red line in Fig. 3c) have a slightly smaller wind speed bias than that with other thresholds, which is related to the AMV products themselves.
Figure3. Average wind speed bias against FNL AMVs derived from (a) WV-High, (b) WV-Mid, and (c) IR in November 2018 (the black solid line denotes that the bias of AMV data is zero).


Figure4. Average wind speed RMSE against FNL AMVs derived from (a) WV-High, (b) WV-Mid, and (c) IR in November 2018.


Besides, the probability density function (PDF) of deviation for AMV data with the three QI values against FNL is shown in Fig. 5. In order to analyze FY-4A AMVs conveniently, AMV data are divided into three bins according to height; namely, the upper layer (400?50 hPa), middle layer (700?400 hPa) and lower layer (1000?700 hPa). It can be seen that fewer outliers with negative deviation are found in a larger QI threshold, making the deviation distribution more symmetric. At the same time, FY-4A AMVs from the two water vapor channels have a poor PDF distribution in the lower layer owing to its small number of samples.
Figure5. PDF of deviation of AMVs derived from (a?c) WV-High, (d?f) WV-Mid and (g?i) IR in the upper layer (above 400 hPa), middle layer (400?700 hPa) and lower layer (below 700 hPa).


Taking the data amount and quality (bias, RMSE and PDF) of the three channels into consideration, AMVs with QI greater than 80 will be extracted from products for subsequent research, and AMVs below 700 hPa in the two water vapor channels will not be assimilated.

3. Observation error estimation
Both observation error and background error determine the effect of the observation on the background in the assimilation. It is of great importance to estimate the reasonable observation error before assimilating a new kind of observation (Errico et al., 2001). Otsuka et al. (2018) analyzed the error characteristics of Rapid Scan AMVs from Himawari-8, and found that AMVs derived from different channels have different characteristics, suggesting the need to calculate the observation error of each channel.
Conventional observation data or reanalysis data are often used as reference fields in the observation error statistics. The RMSE is usually regarded as the observation error on the premise that the observation is not related to the reference field (Benjamin et al., 1998; Gao et al., 2012). Most previous studies estimated AMV observation error through direct comparison between AMVs with radiosonde data (Hollingsworth and L?nnberg, 1986; Bormann et al., 2003; Bedka et al., 2008). As radiosonde data are too sparse in terms of spatial and temporal resolution compared to FY-4A AMVs, if selected as reference fields, fewer samples of data pairs could be selected and therefore representative errors would be caused. Therefore, the FNL reanalysis is used as the reference field in this study.
Referring to the layered criterion in the default observation error file from the Weather Research and Forecasting Model Data Assimilation (WRFDA) system and the statistical method used in Cordoba et al. (2017) and Desroziers et al. (2005), the model layer is vertically divided into n layers at an interval of 50 hPa. The number of AMV data at each layer is recorded as $ {m}_{j}\left(j=\mathrm{1,2},3,\cdots,n\right) $. An assumption is made that the error statistics at each layer are homogeneous. Hence, the corresponding observation error at the $ j\mathrm{t}\mathrm{h} $ layer $ \left(j=\mathrm{1,2},3,\cdots,n\right) $ is recorded as
where $ {y} $, $ {{x}}_{\mathrm{b}} $ and $ H $ denote the FY-4A AMV observations, reference background fields, and nonlinear observation operator, respectively. It should be noted that the observation passing quality control is used here for error statistics.
Figure 6 shows the AMV observation error profiles for the three channels averaged from 1?30 November 2018. Compared to default observation errors, the new observation errors vary with height obviously. The error profiles of the three channels have similar structure and comparable values above 400 hPa. However, between 700 hPa and 400 hPa, the observation errors of WV-High ($ {R}_{\mathrm{w}\mathrm{v}\mathrm{h}} $) increase with height, while the observation errors of WV-Mid ($ {R}_{\mathrm{w}\mathrm{v}\mathrm{m}} $) and IR ($ {R}_{\mathrm{i}\mathrm{r}} $) decrease with height.
Figure6. Vertical profiles of default observation error (black dotted line) and new observation error (red line) of AMVs for (a) WV-High, (b) WV-Mid and (c) IR (units: m s?1).


In terms of the specific values, $ {R}_{\mathrm{w}\mathrm{v}\mathrm{h}} $ is greater than the default error above 400 hPa, but less than the default error below 400 hPa. $ {R}_{\mathrm{w}\mathrm{v}\mathrm{m}} $ is greater than the default error in the upper and middle layers, which will weaken the weight of observation. As for $ {R}_{\mathrm{i}\mathrm{r}} $, it has smaller values below 800 hPa and larger values above 800 hPa. In addition, it is shown that the minimum values of the new observation errors exist at 200 hPa in $ {R}_{\mathrm{w}\mathrm{v}\mathrm{h}} $, $ {R}_{\mathrm{w}\mathrm{v}\mathrm{m}} $ and $ {R}_{\mathrm{i}\mathrm{r}} $, and all of them are smaller than the default error in WRFDA (4.5 m s?1).
A more accurate observation error value is conducive to a more reasonable estimation of the weight of AMV observations in the process of assimilation, leading to a better 3D wind field analysis. Thus, new observation errors of the three channels will be applied to the subsequent assimilation experiments.

4. Experimental framework
2
4.1. Model configuration
--> This study is based on the Rapid-refresh Multi-scale Analysis and Prediction System—Short Term (RMAPS-ST), which runs operationally at the Beijing Meteorological Bureau (Fan et al., 2013; Chen et al., 2014; He et al., 2019; Xie et al., 2019). In RMAPS-ST, the Advanced Research version of the WRF model, version 3.9, and the WRFDA system are employed. Assimilation is performed on a large domain with a grid spacing of 9 km and $ 649\times 500 $ grid points, using 51 vertical levels laid out with an upper limit of 50 hPa.
The experimental domain is shown in Fig. 7 with the distribution of the main conventional observations at 0000 UTC 1 November 2018. AMVs derived from the WV-High, WV-Mid and IR channels are assimilated in this study by merging with conventional observations. Both the initial and boundary conditions are obtained from the Global Forecast System produced by the National Centers for Environmental Prediction. The Thompson microphysics scheme (Thompson et al., 2008) and RRTMG long- and shortwave radiation scheme (Iacono et al., 2008) are used. In addition, the Kain?Fritsch cumulus parameterization scheme (Kain, 2004) is also adopted in this domain.
Figure7. Domain of the model and the locations of the (a) main types of conventional observation data, (b) WV-High AMVs, (c) WV-Mid AMVs and (d) IR AMVs (SYNOP denotes surface station observations; SHIP denotes ship observations; BUOY denotes buoy observations, and SOUND denotes radiosonde observations).



2
4.2. Experimental setup
--> Capability to assimilate AMVs derived from the three channels with new observation errors has been implemented within RMAPS-ST in this study. Five experiments (Table 1) are configured to investigate the added value from the assimilation of FY-4A AMVs on analyses and forecasts with new observation error. Firstly, the control experiment only assimilates the conventional observations (hereafter CTRL). Then, AMVs from WV-High (CONWVh), WV-Mid (CONWVm) and IR (CONIR) are assimilated together with conventional data, respectively. Moreover, an experiment with AMVs from the three channels assimilated simultaneously (CONALL) is added as well. Additionally, to better display the impact of the new observation error on analyses and forecasts, an experiment assimilating three-channel AMVs with default observation error in WRFDA is necessary, and is named CONALL-DEF.
ExperimentObservation errorAssimilation data
CTRLDefault value defined in WRFDAConventional observations
CONWVhNew observation error RwvhAMVs from WV-High Conventional observations
CONWVmNew observation error RwvmAMVs from WV-Mid Conventional observations
CONIRNew observation error RirAMVs from IR Conventional observations
CONALLNew observation error Rwvh, Rwvm and Rir for each channelConventional observations AMVs from three channels
CONALL-DEFDefault observation error in WRFDAConventional observations AMVs from three channels


Table1. Experiment design schemes.


In the six experiments, partial cycle operation (Hsiao et al., 2012) is adopted at 0000 UTC every day, which means that the initial field of 0000 UTC comes from a 6-h forecast beginning at 1800 UTC of the previous day. Then, the assimilation runs eight times at 0000, 0300, 0600, 0900, 1200, 1500, 1800 and 2100 UTC per day with 3-h cycling. At each analysis time, data within a 3-h window centered on the synoptic time will be ingested and a 24-h forecast will be obtained.

5. Results and discussion
2
5.1. Data absorption ratio
--> Table 2 calculates the average proportion of the AMVs absorbed by the assimilation system to the total AMVs used in each assimilation process from experiments with new observation errors first. On average, more than 90% of the data are absorbed in each assimilation cycle, indicating the utilization rate of FY-4A AMVs from the three channels with QI greater than 80 is relatively high. Since observations from radiosondes are obtained twice a day at 0000 UTC and 1200 UTC, the time series of AMV observation minus background (OMB) and AMV observation minus analysis (OMA) of horizontal wind components at 1200 UTC are statistically analyzed in Fig. 8. It can be seen from the time series of OMB (Figs. 8a and c), the U component in FY-4A AMVs is smaller than the background, while the V component is larger than the background. However, the assimilation clearly reduces the deviation between the analysis and observation (Figs. 8b and d), with a more stable time series and a weakened amplitude.
ExperimentLayers (hPa)
900?1000800?900800?700700?600600?400400?300300?200200?100100?50
CONWVh???99.699.9799.9499.9599.9999.99
CONWVm???99.9999.9899.9699.9799.9999.99
CONIR97.598.799.8499.9699.999.9899.9899.9999.99
CONALL97.598.799.7999.8999.9399.9699.9699.9799.99


Table2. Average absorption ratios (%) of AMVs in different layers (hPa) from experiments assimilating with new observation errors.


Figure8. Time series of the (a, c) OMB and (b, d) OMA statistics of the U and V components for assimilating FY-4A AMVs at 1200 UTC from the experiment assimilating FY-4A AMVs with new observation errors.


Moreover, the OMB sequence obtained by simultaneous assimilation of AMVs from the three channels is stable with no obvious outliers under the mutual restriction of data from the three channels. Also, the OMA is closer to 0 than that of the single-channel assimilation experiments, which means that the wind analysis obtained from the multi-channel assimilation experiment is better than that from the single-channel assimilation experiments.

2
5.2. Verification for single-channel assimilation
--> The RMSEs for the analysis and forecast against sounding observations are calculated by the Model Evaluation Tools (MET) verification package (Kirby, 2004), which is designed to be a highly configurable, state-of-the-art suite of verification tools. MET provides a variety of verification techniques, including standard verification scores comparing gridded model data obtained from the WRF modeling system with point-based observations.

3
5.2.1. Wind fields
--> Since the AMV products mainly contain wind field information, only horizontal wind components (U and V) are discussed. Figure 9 shows the vertical profiles of the mean analysis RMSEs for horizontal winds (U, V) (solid lines) and the forecast RMSEs at t+6 h (dashed lines) and t+24 h (dotted lines) of CTRL and the experiments assimilating the single-channel AMVs. The experiments assimilating FY-4A AMVs achieve comparable analysis with CTRL, and even slightly inferior to CTRL, especially for the U component at 300 hPa, which may be related to the sparse resolution of the reference radiosonde observations.
Figure9. Averaged vertical RMSE pro?les for analyses and 6-h and 24-h forecasts in comparison with radiosonde data for CTRL and single-channel assimilation experiments.


After 6-h forecasts, the experiments assimilating FY-4A AMVs reduce the RMSEs in wind fields compared to CTRL. Positive effects mainly occur in the upper and middle layers, which is in accordance with the layer where FY-4A AMVs are distributed (Fig. 1d). Besides, among the three channels, the assimilation of WV-High AMVs produces more obvious improvement, while AMVs from the other two channels have a relatively weak effect on the 6-h forecast. After 24-h forecasts, the RMSEs of the horizontal wind components become larger, and the advantage of assimilating FY-4A AMVs gradually diminishes, but still appears in the U component below 500 hPa and V component above 500 hPa. Moreover, AMVs from the IR channel have a clearer positive impact on the 24-h forecast than that from the two water vapor channels.
At the same time, the mean RMSEs in the three layers for different experiments are calculated for various lead times. The results (Fig. 10) show that the assimilation of FY-4A AMVs reduces the RMSEs of the U and V components in the upper layer, and the improvement can even be seen in the 24-h forecast. Focusing in the middle and lower layers, AMVs from the two water vapor channels produce little positive effect because of no observations below 700 hPa. However, the assimilation of AMVs from IR improves the forecast of the V component in the middle layer, and improvement of the horizontal components in the lower layers in CONIR is gradually shown after 18 h of forecasting.
Figure10. Average RMSEs at different lead times for the (a, b) upper layer (above 400 hPa), (c, d) middle layer (400?700 hPa) and (e, f) lower layer (below 700 hPa) for the CTRL and single-channel assimilation experiments.


In terms of the vertical distribution of FY-4A AMVs, a large number of AMVs from the three channels concentrate at 300 hPa, while only AMVs from IR are available below 700 hPa. Figure 11 displays the trend of RMSEs at 300 hPa and 700 hPa with lead time. It is clear that AMVs derived from the two water vapor channels have an advantage in improving the forecast in the upper layer, especially the data from WV-High, while AMVs from IR can produce noticeable improvement in the lower layer after 12 h of forecasting.
Figure11. RMSEs at 300 hPa and 700 hPa at different lead times for the CTRL and single-channel assimilation experiments.



3
5.2.2. Precipitation forecast
--> Several precipitation processes occur during the research period, so the impact of FY-4A AMVs on the precipitation forecasting skill is evaluated. The bias score (BS), threat score (TS), probability of detection (POD), and false alarm ratio (FAR) (formulae shown in Table 3) are used. The closer the BS, TS and POD are to 1, the higher the prediction accuracy of precipitation; and the smaller the FAR, the fewer false alarms in the precipitation forecast. The precipitation data for comparison are the hourly observations from surface stations provided by the CMA.
IndexFormula
Bias Score (BS)(hits + false alarms)/(hits + misses)
Threat Score (TS)hits/(hits + false alarms + misses)
Probability of Detection (POD)hits/(hits + misses)
False Alarm Ratio (FAR)false alarms/(hits + false alarms)


Table3. Four statistical indices for rainfall forecast skill.


The average TS, BS, POD and FAR for the 6-h accumulated precipitation with thresholds of 0.1 mm, 1.0 mm, 5.0 mm 10.0 mm and 25.0 mm are shown in Fig. 12. Compared to CTRL, the assimilation of the FY-4A AMVs results in a decrease in BS for 6-h accumulated precipitation with intensity between 1.0 mm and 10.0 mm, and CONIR has the smallest values. Besides, the average TS and POD in CONWVm and CONIR show a slight improvement in precipitation with thresholds between 1.0 mm and 10.0 mm. In addition, the FAR of the experiments assimilating FY-4A AMVs decrease when compared to CTRL, and CONIR produces the lowest FAR value. However, owing to the actual intensity of the 6-h accumulated precipitation from these rainfall events being lower than 25.0 mm, the added AMVs lead to an under-prediction with the threshold of 25.0 mm, overall.
Figure12. Average BS (a), TS (b), POD (c) and FAR (d) for the 6-h accumulated precipitation with different thresholds for the CTRL and single-channel assimilation experiments.


Generally, the results suggest that AMVs derived from the IR channel produce the largest improvement among the three channels. The reason may be that the precipitation intensity of these processes during the experiment is relatively weak and the vertical development relatively low. As a result, the wind of the lower layer in the IR channel has a positive impact on the precipitation.

2
5.3. Verification for multi-channel assimilation
--> Through the above discussion of the impact of assimilating single-channel AMVs on the analysis and forecasts, it can be concluded that the FY-4A AMVs from the water vapor channels mainly improve the wind field in the upper layer. Meanwhile, assimilating AMVs from IR has a slightly positive impact in the middle and lower layers. Next, to combine the advantages of the three channels, FY-4A AMVs from the three channels are simultaneously assimilated, and the impact on the forecast further analyzed.
Figure 13 displays the average RMSEs of the analysis and forecasts for the horizontal winds (U, V) of CTRL, CONALL and CONALL-DEF. Similar to Fig. 9, the assimilation of FY-4A AMVs derived from the three channels produces analyses a little inferior to that of CTRL, especially for the upper layer where a large number of AMVs are distributed, which is due to the rough spatial resolution of the radiosonde data. Also, the improvement from assimilating FY-4A AMVs will be seen after integrating for several hours. For example, the differences between the RMSEs of the experiments assimilating AMVs and that of CTRL decrease after six hours of model integration, and their RMSEs are relatively smaller than those of CTRL for the V component under 400 hPa. For the 24-h forecast, the RMSEs of U and V in CONALL and CONALL-DEF are smaller than those in CTRL at almost all levels. The main advantage of assimilating FY-4A AMVs in the U component gradually shifts from the upper layer to the middle and lower layer, while the improvement of the V component is still shown in the upper layer.
Figure13. Averaged vertical RMSE pro?les for analyses and 6-h and 24-h forecasts in comparison with radiosonde data for CTRL, CONALL and CONALL-DEF.


In addition, Fig. 13 also shows the impact of different observation errors on assimilation and prediction results. As can be seen from the RMSE profiles, when compared to CONALL-DEF (gray lines), assimilating FY-4A AMVs with new observation errors (purple lines) can clearly reduce the RMSEs of the horizontal wind components in almost all layers both in the analyses and 6-h and 24-h forecasts, which can serve as a reference for the significance of estimated observation error in new observation datasets.
Besides, according to the precipitation forecasting skill scores displayed in Fig. 14, except for the threshold of 0.1 mm, CONALL results in a decrease of BS and FAR when compared to CTRL. Moreover, the average TS and POD in CONALL show a slight improvement for precipitation with thresholds between 1.0 mm and 10.0 mm, which is similar to the result obtained in the single-channel assimilation.
Figure14. Average BS (a), TS (b), POD (c) and FAR (d) for the 6-h accumulated precipitation with different thresholds for CTRL and CONALL.



6. Conclusions
In order to further promote the practical application of FY-4A satellite data in NWP, the features of FY-4A AMV products derived from the WV-High, WV-Mid and IR channels are analyzed in this paper and their observation errors estimated. Then, one-month single-channel and multi-channel data assimilation and forecasting experiments are carried out based on the operational RMAPS-ST system. The main results are as follows:
(1) FY-4A AMVs derived from the two water vapor channels are concentrated in the upper and middle layers, especially between 400 hPa and 200 hPa, while AMVs from IR have a maximum total number in the upper and lower layers. In addition, AMVs derived from the three channels show relatively normal PDF distributions of deviation from the reference field in the upper and middle layers, as do FY-4A AMVs in the lower layer derived from the IR channel.
(2) Observation errors of FY-4A AMV products have obvious vertical structures. The profiles of the three channels above 400 hPa are similar in structure and are quantitatively comparable. Between 400 hPa and 700 hPa, the observation errors of WV-High increase with height, while those of WV-Mid and IR decrease with height.
(3) The one-month single-channel data assimilation cycling and forecasting experiments demonstrate that assimilating FY-4A AMVs improves wind forecasts. Moreover, AMVs from WV-High produce more apparent improvement in the wind field in the upper layer, and AMVs from IR can achieve more obvious positive effects in the lower layer. Furthermore, assimilating AMVs from IR produces higher skill scores for medium and moderate precipitation than the other channels owing to its good quality of data in the lower layer.
(4) Simultaneous assimilation of FY-4A AMVs from the three channels combines the advantages of assimilating the three individual channels to improve the forecasts in the upper, middle and lower layers at the same time. Moreover, simultaneous assimilation has an obvious positive impact on the precipitation forecast with thresholds between 1.0 mm and 10.0 mm.
Considering that the AMV products derived from the FY-4A satellite are still being studied, more comprehensive and detailed comparisons and evaluations should be conducted in the future by combining the AMV products derived from other satellites, such as those from Himawari-8. In addition, stricter quality-control procedures for FY-4A AMVs should be implemented to reduce the correlation between AMVs. Moreover, the experimental scheme of simultaneous assimilation of FY-4A AMVs of the three channels still needs to be improved, and duplication of the AMVs between different channels for assimilation should be further considered.
Acknowledgements. This research was funded by the National Key Research and Development Plan (Grant No. 2018YFC1507105). We wish to thank Dr. Xiaohu ZHANG from the National Satellite Meteorological Center of China for providing the AMV data and other relevant help. We also thank the Editor and reviewers for their comments and suggestions for improving the manuscript.

闂備浇顕х€涒晠顢欓弽顓炵獥闁圭儤顨呯壕濠氭煙閸撗呭笡闁抽攱鍨堕妵鍕箛閸撲焦鍋у銈忕到瀵墎鎹㈠┑瀣劦妞ゆ帒瀚粻銉︺亜閺冣偓閸庡啿鐣烽妷銉㈡斀妞ゆ梻鐡斿▓妯汇亜閿斿墽銆楅梻鍌氬€烽悞锔锯偓绗涘懐鐭欓柟鐑橆殕閸庢銇勯弬鍨挃闁告宀搁幃妤€鈽夊▎瀣窗闂佸磭绮Λ鍐蓟濞戙垹唯闁靛繆鍓濋悵锕€鈹戦悙鎻掓倯鐎光偓閹间礁钃熼柕鍫濐槸缁狙囨煕椤垵鏋熸繛鍫濈埣濮婅櫣绱掑Ο鐓庘拡闂佺ǹ绨洪崐妤冨垝濞嗘劕绶為柟閭﹀墰椤旀帡姊洪懖鈹炬嫛闁告挻鐩獮濠傗攽鐎n偆鍘介梺缁樏鑸靛緞閸曨厸鍋撶憴鍕闁挎洏鍨介弫鎰版倷绾版ê浜鹃柨婵嗛娴滃墽绱掗幉瀣М闁哄瞼鍠愬ḿ蹇斻偅閸愨晩鈧秴顪冮妶鍡樺暗闁硅櫕鎹囨俊鎾箳閹搭厽鍍甸梺鎯ф禋閸嬪棝濡存繝鍌楁斀闁绘﹩鍋勬禍楣冩⒑閸涘﹦鈽夐柣掳鍔戝鏌ヮ敆閸曨剙鈧爼鏌i幇顓炵祷闁抽攱姊规穱濠勪沪閸屾艾顫囬梺鍝勭焿缁绘繂鐣锋總鍓叉晝闁靛繒濮寸粻鎵磽閸屾瑧鍔嶆い銊ユ嚇钘濋梻鍫熺〒閺嗭附鎱ㄥ鍡楀伎缂佹唻绲介湁闁挎繂瀚鐔访归悡搴㈠枠婵﹥妞藉畷銊︾節閸屾粎鎳嗘繝纰樻閸嬪懐鎹㈤崼銉у祦闁圭増婢橀柋鍥煛閸モ晛鏋戞繛鍛墬缁绘稒娼忛崜褏袣濠电偛鎷戠紞浣哥暦閳ユ枼鍫柛顐ゅ枎閳ь剛鏁婚弻锝夊閻樺吀绮甸梺鎼炲€愰崑鎾绘⒒娴h鍋犻柛鏂跨Ф閹广垽宕熼鐐茬亰闁荤姴娲ゅ顒勫籍閸屾浜鹃柨婵嗛婢т即鏌ㄥ☉铏
547闂傚倸鍊风粈浣革耿闁秲鈧倹绂掔€n亞锛涢梺鐟板⒔缁垶鍩涢幒妤佺厱闁哄洦顨嗗▍鍛存煟韫囷絼閭柡宀嬬到铻栭柍褜鍓熼弻濠囨晲婢跺﹦鐤囬梺瑙勫礃椤曆呯矆閸愨斂浜滈煫鍥ㄦ尰椤ユ粓鏌i幘瀵告创婵﹥妞藉畷妤呭礂閼测晛鍤掗梻浣告憸閸c儵宕戦崟顖涘仼闁绘垼濮ら弲鎼佹煥閻曞倹瀚�1130缂傚倸鍊搁崐椋庣矆娓氣偓钘濇い鏍ㄧ矌閻捇鏌i姀銏╂毌闁稿鎹囧畷褰掝敃閿濆洦顓婚梻渚€鈧偛鑻晶鍙夈亜椤愩埄妲洪柛鐘诧工椤撳吋寰勬繝鍐发闂備線鈧偛鑻晶鎾煛鐏炵ǹ澧茬€垫澘瀚埀顒婄秵娴滄粓锝炴惔鈾€鏀介柣鎰皺婢ф盯鏌涢妸銉т虎闁伙絿鍏橀弫鎾绘偐閼碱剦妲伴柣鐔哥矊缁绘帞鍒掑▎鎾崇倞闁冲搫鍟伴敍婊勭節閵忥絾纭鹃柨鏇檮閺呫儵姊绘担鍛婂暈閼垦囨煕婵犲啰澧辩紓鍌涙崌閹稿﹥绔熼埡鍌涱棃鐎规洦浜濋幏鍛喆閸曨厽袙闂傚倷绀佸﹢閬嶅储瑜旈幃娲Ω閳轰胶顔囬梺鍓插亝濞叉牠鎮為崹顐犱簻闁瑰鍋涢婊呯磼娴e搫顣奸柕鍥у椤㈡洟濮€閳轰礁澹夐梻浣瑰缁嬫帡宕濆▎鎾跺祦閻庯綆鍣弫宥嗙節婵犲倹濯兼俊鏌ユ敱缁绘繈鎮介棃娑楁勃闂佹悶鍔岀紞濠傤嚕椤愩倗鐭欐俊顐ゆ暩閸庛倗鎹㈠┑瀣倞闁靛ǹ鍎查娲⒑閸濆嫷妲搁柣妤佺矒瀹曟垿骞掑Δ鈧崥褰掓煛閸愩劎澧涢柣鎾崇箻閺屾盯鍩勯崘鈺冾槷闂佺ǹ绻愰惉鑲╂閹烘鏁婇柛蹇擃槸娴滈箖鎮洪幒宥嗙グ妞ゎ偄顦甸獮鍡涘棘鎼存挻顫嶉梺鍦劋閸ㄨ櫕绔熼弴鐔虹閺夊牆澧介幃濂告煛閸滀礁浜滈崡閬嶆煕椤愮姴鍔滈柣鎾存礋閺岀喖鎮滃Ο璇查瀺闂佺ǹ顑囨繛鈧柡灞界Х椤т線鏌涜箛鏃傘€掓繛鍡愬灩椤繄鎹勫ú顏咃紬濠德板€х徊浠嬪疮椤栫偛绠犳慨妯诲閸嬫捇鐛崹顔煎闂佸摜鍠撴繛鈧€规洘鍨块獮妯尖偓娑櫭鎾绘⒑缂佹ê鐏﹂柨姘舵煃瑜滈崜娑㈠磻閻旂厧鐒垫い鎺戝枤濞兼劙鏌熼鑲╁煟鐎规洘娲熼、娑㈡倷閼碱剙濮︽俊鐐€栭崹濂稿磿婵犳艾绀傞悘鐐板嫎娴滄粓鏌熼幆褏鎽犻柛濠冨姍閺岀喖顢欓懖鈺佺厽閻庤娲熸禍鍫曠嵁閸ヮ剙惟闁靛鍊曢ˉ姘攽閻愯埖褰х紒鑼亾缁傚秹鎳為妷褜娴勯梺鎸庢磵閸嬫捇妫佹径鎰厱闊洦娲栫敮鍫曟煙閸忓吋鍊愰柡灞界Х椤т線鏌涜箛鏃傗槈闂囧鏌涢妷顔煎缁炬儳銈搁弻鏇熺節韫囨洜鏆犲銈嗘礉妞存悂骞堥妸銉庣喓绮欓崹顔碱潕闂備胶枪椤戝棝骞愭繝姘闁告侗鍨遍崰鍡涙煕閺囥劌澧伴柡鍡樺哺濮婄粯鎷呴崫鍕粯閻庢鍠楅崕濂稿Φ閹版澘绀冩い蹇撴婢跺嫭淇婇妶蹇曞埌闁哥噥鍨跺畷鎴︽偐缂佹ḿ鍘遍梺鏂ユ櫅閸燁偅鎱ㄩ埀顒勬⒑缁嬫鍎愰柟绋款煼楠炲繘宕ㄧ€涙ê浠惧銈嗙墬閼瑰墽绮婚幋锔解拻濞达絽鎲¢幆鍫ユ煕婵犲媱鍦弲闂侀潧绻堥崐鏇㈡倿閸偁浜滈柟瀵稿仜椤曟粎绱掓担瑙勭凡妞ゎ叀娉曢幉鎾礋椤掆偓绾炬娊鎮楀▓鍨灓闁轰礁顭烽妴浣糕槈濡粎鍠庨悾鈩冿紣娴e壊妫滄繝鐢靛仩閹活亞绱為埀顒併亜椤愩埄妯€闁诡噯绻濋崺鈧い鎺戝閳锋帡鏌涚仦鎹愬闁逞屽厸缁瑩銆佸▎鎰瘈闁告洦鍓﹀ḿ鐔兼⒑閸撴彃浜濇繛鍙夌墱缁崵绱掑Ο闀愮盎闂佸湱鍋撳ḿ娆撍夊鑸电厱閻庯綆鍋呯亸浼存煙瀹勭増鍤囩€规洜鍏橀、妯衡攦閹傚婵犵數濮村ú锕傚磹閸偆绠鹃柟瀛樼箘閺嬪啰绱掗埀顒傗偓锝庡枟閻撳啰鎲稿⿰鍫濈婵炴垶姘ㄩ惌鍡椕归敐鍫熴€冮柣鎺戯躬閺岀喎鈻撻崹顔界亐闂佺ǹ顑嗛幐鎼佸煡婢跺备鍋撻崷顓炐ユい锔垮嵆濮婃椽鏌呴悙鑼跺闁告ê鎲$换娑㈠级閹寸偛鏋犲銈冨灪閿曘垹鐣烽锕€唯鐟滃瞼绮径濞炬斀閹烘娊宕愰幘缁樺€块柨鏃€宕樻慨铏叏濡灝鐓愰柣鎾崇箻閺屾盯鍩勯崘鈺冾槶濡炪倧璁f俊鍥焵椤掆偓濠€杈ㄧ仚濡炪値鍘奸崲鏌ユ偩閻戣姤鏅查柛銉檮閸曞啴姊虹粙鎸庢拱妞ゃ劌鎳忕粋宥夊箚椤€崇秺閺佹劙宕ㄩ鍏兼畼闂備礁鎲″ú鐔虹不閺嶎厼钃熼柨婵嗩槹閸婄兘鏌涘▎蹇f▓婵☆偅鍨垮娲焻閻愯尪瀚伴柛妯绘倐閺岋綁骞掗悙鐢垫殼閻庢鍣崜鐔风暦閹烘埈娼╅柨婵嗘閸欏啯绻濆▓鍨灓闁硅櫕鎸哥叅闁绘ǹ顕х壕濠氭煙閹呮憼濠殿垱鎸抽弻娑樷槈濮楀牊鏁鹃梺闈╃稻閹倿寮婚敐鍡樺劅闁靛繒濮撮弸娆撴煟閻樺啿濮夐柛鐘崇墵瀹曟椽鍩€椤掍降浜滈柟瀵稿仜椤曟粍銇勯敃鍌ゆ缂佽鲸鎸搁濂稿椽娴gǹ澹庨梻浣侯攰濞呮洟鎮ч悩璇茬疇闁绘ɑ妞块弫鍡涙煕閹邦喖浜鹃柦鎴濐槺缁辨捇宕掑▎鎴濆闂佸憡鍔曢…鐑界嵁韫囨稒鏅搁柨鐕傛嫹28缂傚倸鍊搁崐椋庢閿熺姴纾诲鑸靛姦閺佸鎲搁弮鍫涒偓浣肝旈崨顔间簻闂佹儳绻愬﹢閬嶆偂鐎n喗鈷戦悷娆忓閸庢鏌涢妸銉э紞濠㈣娲熼弫鎾绘晸閿燂拷
相关话题/Characteristics Fengyun Satellite

闂傚倸鍊烽懗鍓佸垝椤栫偛绀夋俊銈呮噹缁犵娀鏌熼幑鎰靛殭闁告俺顫夐妵鍕棘閸喒鎸冮梺鍛婎殕瀹€鎼佸箖濡ゅ懏鏅查幖瀛樼箘閻╁海绱撴担椋庤窗闁革綇缍佸濠氭偄閻撳海顦ч梺鍏肩ゴ閺呮繈鎮¢崒鐐粹拺缂佸娉曢悞鍧楁煙閸戙倖瀚�2濠电姷鏁搁崑鐐哄垂閸洖绠扮紒瀣紩濞差亜惟闁冲搫顑囩粙蹇涙⒑閸︻厼鍔嬫い銊ユ瀹曠敻鍩€椤掑嫭鈷戦柛娑橈工婵箑霉濠婂懎浠辩€规洘妞介弫鎾绘偐瀹曞洤骞楅梻渚€娼х换鍫ュ磹閺嵮€妲堢憸鏃堝蓟閿濆鐒洪柛鎰典簼閸d即姊虹拠鈥虫殭闁搞儜鍥ф暪闂備焦瀵х换鍌毭洪姀銈呯;闁圭儤顨嗛埛鎴︽煕濠靛棗顏╅柍褜鍓欑紞濠囧箖闁垮缍囬柍鍝勫亞濞肩喖姊虹捄銊ユ珢闁瑰嚖鎷�
濠电姷鏁告慨浼村垂瑜版帗鍋夐柕蹇嬪€曢悙濠勬喐瀹ュ棙鍙忛柕鍫濐槹閳锋垿鏌涘☉姗堝伐缂佹甯楁穱濠囶敃閿濆洦鍒涘銈冨灪濡啯鎱ㄩ埀顒勬煏閸繃锛嶆俊顐㈠閺岋絾鎯旈婊呅i梺鍝ュУ閻楃姴鐣烽姀銈呯妞ゆ梻鏅崢鎼佹⒑閸涘﹥绀嬫繛浣冲洦鍊堕柨婵嗘娴滄粓鏌熺€涙ḿ绠栧璺哄缁辨帞鈧綆浜跺Ο鈧梺绯曟杹閸嬫挸顪冮妶鍡楀潑闁稿鎹囬弻宥囨嫚閺屻儱寮板┑锟犵畺娴滃爼寮诲鍫闂佸憡鎸婚悷鈺佺暦椤栨稑顕遍悗娑櫭禍顖氣攽閻愬弶鈻曞ù婊勭箞瀵煡顢楅崟顒€鈧爼鏌i幇顔芥毄闁硅棄鍊块弻娑㈠Χ閸ヮ灝銏ゆ婢跺绡€濠电姴鍊搁弳锝嗐亜鎼淬埄娈曢柕鍥у閸╃偤顢橀悙宸痪婵犳鍨遍幐鎶藉蓟閿熺姴鐐婇柍杞扮劍閻忎線姊洪崨濠勬喛闁稿鎹囧缁樻媴閸濄儳楔濠电偘鍖犻崱鎰睏闂佺粯鍔楅弫鍝ョ不閺冨牊鐓欓柟顖嗗苯娈堕梺宕囩帛濮婂綊骞堥妸銉庣喓鎷犻幓鎺濇浇闂備焦鎮堕崐褏绮婚幘璇茶摕闁绘棁娅i惌娆撴煙缁嬪灝顒㈤柟顔界懇濮婄儤瀵煎▎鎴犘滅紓浣哄У閻楁洟顢氶敐澶樻晝闁冲灈鏅滈悗濠氭⒑瑜版帒浜伴柛妯哄⒔缁瑩宕熼娑掓嫼闂佸湱枪濞寸兘鍩ユ径鎰厸闁割偒鍋勬晶瀵糕偓瑙勬礀缂嶅﹥淇婂宀婃Ъ婵犳鍨伴妶鎼佸蓟濞戞ǚ妲堟慨妤€鐗婇弫鍓х磽娴e搫校閻㈩垳鍋ら崺鈧い鎺嗗亾闁诲繑鑹鹃…鍨潨閳ь剟骞冭瀹曞崬霉閺夋寧鍠樼€规洜枪铻栧ù锝夋櫜閻ヮ亪姊绘担渚敯闁规椿浜浠嬪礋椤栨稒娅栭梺鍝勭▉閸樹粙鎮¤箛娑欑厱闁斥晛鍟粈鈧銈忕岛閺嗘竼e濠电姷鏁告慨浼村垂閸︻厾绀婂┑鐘叉搐閻掑灚銇勯幒宥堝厡闁愁垱娲熼弻鏇㈠幢濡も偓閺嗭綁鏌$仦鍓ф创妤犵偞甯¢獮瀣倻閸℃﹩妫у┑锛勫亼閸婃牜鏁悙鍝勭獥闁归偊鍠氶惌娆忊攽閻樺弶澶勯柛瀣姍閺岋綁濮€閵忊剝姣勯柡浣哥墦濮婄粯鎷呯粙鎸庡€┑鐘灪閿曘垹鐣烽娑橆嚤閻庢稒锚娴滎垶姊洪崨濠勭畵濠殿垵椴搁幆鏃堝閿涘嫮肖婵$偑鍊栭崝妤呭窗鎼淬垻顩插Δ锝呭暞閻撴盯鏌涢妷锝呭闁汇劍鍨块弻锝夋偄閸欏鐝旈梺瀹犳椤︾敻鐛Ο鑲╃闁绘ê宕銏′繆閻愵亜鈧牕煤濠靛棌鏋嶉柡鍥╁亶缂傛岸鐓崶銊р槈鐎瑰憡绻冮妵鍕箻濡も偓閸燁垶顢欓敓锟�20婵犲痉鏉库偓妤佹叏閻戣棄纾婚柣妯款嚙缁犲灚銇勮箛鎾搭棡妞ゎ偅娲樼换婵嬫濞戝崬鍓扮紓浣哄У閸ㄥ潡寮婚妶鍡樺弿闁归偊鍏橀崑鎾澄旈埀顒勫煝閺冨牆顫呴柣妯烘閹虫捇銈导鏉戠闁冲搫锕ラ敍鍛磽閸屾瑧顦︽い锔诲灦椤㈡岸顢橀姀鐘靛姦濡炪倖宸婚崑鎾寸節閳ь剟鏌嗗鍛紱闂佺粯姊婚崢褔寮告笟鈧弻鏇㈠醇濠垫劖效闂佺ǹ楠哥粔褰掑蓟濞戙垹鍗抽柕濞垮劚椤晠姊烘导娆戠暠缂傚秴锕獮鍐ㄎ旈崘鈺佹瀭闂佸憡娲﹂崣搴ㄥ汲閿熺姵鈷戦柛婵嗗椤ユ垿鏌涚€n偅宕屾慨濠冩そ瀹曨偊宕熼崹顐嵮囨⒑閸涘﹥鈷愰柣妤冨█楠炲啴鏁撻悩铏珫闂佸憡娲﹂崜娆撴偟娴煎瓨鈷戦梻鍫熺〒缁犳岸鏌涢幘瀵哥疄闁诡喒鈧枼鏋庨柟閭﹀枤椤旀洘绻濋姀锝嗙【妞ゆ垵妫涚划鍫ュ焵椤掑嫭鍊垫繛鍫濈仢濞呭秹鏌¢埀顒勫础閻戝棗娈梺鍛婃处閸嬫帡宕ョ€n喗鐓曢柡鍥ュ妼楠炴ɑ淇婇崣澶婄伌婵﹥妞藉畷顐﹀礋椤愮喎浜惧ù鐘差儜缂嶆牕顭跨捄鐑樻拱闁稿繑绮撻弻娑㈩敃閿濆棛顦ㄩ梺鍝勬媼閸撶喖骞冨鈧幃娆撴濞戞顥氱紓鍌欒兌婵數鏁垾鎰佹綎濠电姵鑹鹃悙濠囨煥濠靛棙鍣稿瑙勬礋濮婃椽鎳¢妶鍛€惧┑鐐插级閸ㄥ潡骞婂Δ鍐╁磯閻炴稈鍓濋悘渚€姊虹涵鍛涧闂傚嫬瀚板畷鏇㈠箣閿旇棄鈧敻鏌ㄥ┑鍡涱€楁鐐瘁缚缁辨帡鎮╁畷鍥р拰闂佸搫澶囬崜婵嗩嚗閸曨偀妲堟繛鍡楁禋娴硷拷