• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    Use of Microwave Radiances from Metop-C and Fengyun-3 C/D Satellites for a Northern European Limited-area Data Assimilation System

    2021-07-08 09:29:24MagnusLINDSKOGAdamDYBBROEandRogerRANDRIAMAMPIANINA
    Advances in Atmospheric Sciences 2021年8期

    Magnus LINDSKOG, Adam DYBBROE, and Roger RANDRIAMAMPIANINA

    1Swedish Meteorological and Hydrological Institute, Folkborgsv?gen 17, 60361 Norrk?ping, Sweden

    2Norwegian Meteorological Institute, Henrik Mohns Plass 1, 0371 Oslo, Norway

    ABSTRACT MetCoOp is a Nordic collaboration on operational Numerical Weather Prediction based on a common limited-area km-scale ensemble system.The initial states are produced using a 3-dimensional variational data assimilation scheme utilizing a large amount of observations from conventional in-situ measurements, weather radars, global navigation satellite system, advanced scatterometer data and satellite radiances from various satellite platforms.A version of the forecasting system which is aimed for future operations has been prepared for an enhanced assimilation of microwave radiances.This enhanced data assimilation system will use radiances from the Microwave Humidity Sounder, the Advanced Microwave Sounding Unit-A and the Micro-Wave Humidity Sounder-2 instruments on-board the Metop-C and Fengyun-3 C/D polar orbiting satellites.The implementation process includes channel selection, set-up of an adaptive bias correction procedure,and careful monitoring of data usage and quality control of observations.The benefit of the additional microwave observations in terms of data coverage and impact on analyses, as derived using the degree of freedom of signal approach,is demonstrated.A positive impact on forecast quality is shown, and the effect on the precipitation for a case study is examined.Finally, the role of enhanced data assimilation techniques and adaptions towards nowcasting are discussed.

    Key words: data assimilation, regional numerical weather prediction, microwave radiances

    1.Introduction

    The aim of Numerical Weather Prediction (NWP) is to predict the future state of the atmosphere.This is achieved by using a numerical weather forecast model which starts from an initial state and integrates the atmospheric state forward in time.The initial state is prepared through a data assimilation (DA) process (Lorenc, 1986), which optimally blends observations with a forecast from the atmospheric model in order to produce the best possible initial state.The atmospheric state is represented by a number of variables defined on discrete grid points and the model is run forward in time to simulate the atmosphere and its evolution.Some processes take place on scales that cannot be explicitly resolved by the model and are accounted for through various parameterizations.It was realized early on by Lorenz (1965) that the forecast quality is strongly dependent on an accurate description of the initial state and hence on the capability of the DA system.An NWP model is said to be regional if the geographical domain over which the model is being resolved covers only a region of the globe.A regional model uses forecasts from global NWP as lateral boundary conditions.Usually the global systems have coarser resolution than the limited-area systems and focus on forecast ranges from a day or so up to typically twoweeks.A km-scale limited-area NWP system, on the other hand, is typically concerned with producing frequently updated detailed forecasts of the atmospheric state out to one or two days from initialization.

    An important category of instruments measuring the upwelling microwave radiance from the earth and atmosphere in various frequency bands is the passive microwave(PMW) sounder.Examples of PMW sounding instruments placed on-board presently flying satellites are the Microwave Humidity Sounder (MHS) and the Advanced Microwave Sounding Unit-A (AMSU-A), both on-board the Metop spacecrafts (Saunders, 1993; Bonsignori, 2007;Klaes et al., 2007), the Advanced Technology Microwave Sounder (ATMS) on-board the Joint Polar Satellite System(JPSS) satellites (Goldberg and Weng, 2006), and the Micro-Wave Humidity Sounder-2 (MWHS-2) on-board the Fengyun-3 (FY-3) C and D satellites (Zhang et al., 2019).These PMW instruments are all so-called cross-track scanners and placed on various sun-synchronous polar orbiting satellites, each with different Equator crossing times (ECT),so over a particular area of the globe the data coverage and distribution in time is dependent of the satellite ECTs and the geographical location of the area.The horizontal resolution of individual observations in terms of their instantaneous field of view (IFOV) vary with viewing angle, but is typically around 20—50 km at nadir.

    In addition to existing PMW sounding instruments onboard current orbiting satellites, there are instruments planned to be part of future satellite programs.One such example is a small cross-track PMW sounder planned onboard the Arctic Weather Satellite (AWS), a Swedish led initiativefor a prototype satellite under the Earth Watch Programme of the European Space Agency (ESA).The launch for the prototype AWS is planned for the first half of 2024,and is conceived as a demonstrator for a possible follow-on constellation of up to 20 small ( ≈ 120 kg) AWS satellites in various different orbital planes (ESA, 2021).Also a new polar orbiting FY-3 satellite, FY-3E, placed in an early morning orbit and carrying a MWHS-2 instrument, is scheduled for launch in 2021 (Zhang et al., 2019).

    Satellite data is an important source of information in NWP DA for accurately describing structures of atmospheric temperature, moisture, wind fields, surface pressure,and clouds.In particular, PMW radiances sensitive to atmospheric temperature and moisture have been demonstrated to be important observation types for global (Li and Liu,2016; Geer et al., 2017; Lawrence et al., 2018; Carminati et al., 2020; Jiang et al., 2020), as well as for regional NWP(Storto and Randriamampianina, 2010; Schwartz et al.,2012; Xu et al., 2016; Zou et al., 2017).The PMW radiances demonstrated to be useful in global NWP reside from a wide range of instruments, including AMSU-A and MHS(Geer et al., 2017; Lawrence et al., 2018) as well as MWHS-2 (Li and Liu, 2016; Lawrence et al., 2018; Carminati et al., 2020; Jiang et al., 2020).Geer et al.(2017)have reported on a growing impact from humidity sensitive PMW radiances.For regional models the importance of PMW radiances has been demonstrated using various approaches, including an energy-norm based method (Storto and Randriamampianina, 2010), case studies (Li and Liu,2016; Jiang et al., 2020), and a data denial procedure(Schwartz et al., 2012).

    PMW radiances are influenced by clouds and precipitation and until recently have been used for NWP DA only in clear-sky conditions.However, today all-sky PMW radiances are also assimilated by several operational centers(Geer et al., 2018).These improvements have been facilitated by enhanced radiative transfer models and representation of moist and cloud processes in the observation operators and have been obtained in global as well as regional model frameworks.In addition, research is being conducted aiming at an improved use of microwave sounding channels peaking low in the atmosphere and thus more sensitive to surface conditions (Karbou et al., 2005; English, 2008; Frolov et al., 2020).Key challenges are to handle surface emissivities and temperatures in combination with sometimes highly heterogeneous surface properties within the IFOV.

    Here, we use a km-scale limited-area NWP system over a northern European domain (see Fig.1).At high latitudes,radiances from geostationary satellites are in general much less useful than over areas closer to the Equator.This is due to the oblique viewing and larger IFOVs at high latitudes.Radiance observations from sensors on board polar orbiting satellites, on the other hand, are crucial over northern Europe.These observations become more and more frequent with increasing latitude, and can achieve higher horizontal resolution with their lower orbit and more nadir view.

    Fig.1.MetCoOp modelling domain used in this study.

    In this paper, the focus is on the current use of PMW radiances and on the effect of an increased use of PMW radiance sounding data by extending the usage to also include radiances from PMW instruments on board the European Metop-C and the Chinese FY-3C and FY-3D satellites.The components of the limited-area forecasting system are the subject of section 2, with special focus on the handling of PMW radiances.In section 3, the experimental design is described, followed by a presentation of results in section 4.Concluding remarks are presented in section 5 together with a future outlook.

    2.The limited-area forecasting system

    2.1.General description

    The national meteorological weather services of Sweden, Norway, Finland, and Estonia have joined forces in a collaboration, named MetCoOp, around a common operational km-scale forecasting system (Müller et al., 2017).It is a configuration of the shared Aire Limitée Adaptation dynamique Developpement InterNational (ALADIN)-High Resolution Limited Area Model (HIRLAM) NWP system.This system can be run with different configurations and in MetCoOp the HIRLAM-ALADIN Regional Meso-scale Operational NWP In the Europe Application of Research to Operations at Mesoscale (HARMONIE-AROME) is used (Bengtsson et al., 2017) and is run as an ensemble forecasting system.In our study we use the cy43 version of the MetCoOp forecasting system.The northern European model domain is illustrated in Fig.1.It has 900 × 960 horizontal grid points with a grid distance of 2.5 km, and 65 vertical model levels,extending from roughly 12 m above ground up to approximately 33 km (10 hPa).

    The three main components of this forecasting system are surface DA, upper-air DA, and the forecast model.The system is run with a 3-hourly assimilation cycle and launching forecasts at 0000, 0300, 0600, 0900, 1200, 1500, 1800,and 2100 UTC.Only the forecasts launched at the synoptic hours 0000, 0600, 1200, and 1800 UTC are being used by the duty forecasters.Therefore, these have a relatively strict observation cut-off time of 1 h and 15 min compared to the 3 h and 20 min for the asynoptic cycles 0300, 0900, 1500,and 2100 UTC.For the asynoptic cycles, the only time-constraint is to produce a 3 h forecasts serving as background state for the DA at the following synoptic cycle.Thus, this 3 h forecast can be produced just prior to the synoptic cycle DA.Due to the operational cut-off constraints mentioned above, in practice only observations within a time-range from -1 h and 30 min to +1 h and 15 min are used for the cycles at 0000, 0600, 1200, and 18 UTC.For the asynoptic cycles observations within the entire time-range -1 h and 30 min to + 1 h and 29 min are used.

    A detailed description of the forecast model setup is given in Seity et al.(2011) and Bengtsson et al.(2017).It is a non-hydrostatic model formulation with a spectral representation of the model state (Bubnová et al., 1995; Bénard et al.,2010).Sub-grid scale parameterization of clouds, including shallow convection is handled by the EDMF (Eddy Diffusitivity Mass Flux) originating from de Rooy and Siebesma(2008) and Neggers et al.(2009).Deep convection is resolved by the model.Turbulence and vertical diffusion is represented using the so-called HARATU scheme, which is based on a Turbulent Kinetic Energy scheme by Lenderink and Holtslag (2004).The radiative transfer is modelled as described by Fouquart and Bonnel (1980) and Mlawer et al.(1997) for short- and long-wave radiation processes, respectively.Surface processes are modeled using the SURFEX (Surface Externalisée) scheme (Masson et al., 2013).Global forecasts provided by the European Centre for Medium-Range Weather Forecasts (ECMWF) are used as lateral boundary conditions.These forecasts are launched every 6 h with a 1 h output frequency.Global model information is also used to replace larger-scale information in the background state with lateral boundary information (Müller et al., 2017).The idea is to make use of high-quality large-scale information from the ECMWF global fields, in the MetCoOp analysis.

    In the surface DA synoptic observations of two-meter temperature, two-meter relative humidity and snow cover are used to estimate the initial state of the soil temperature,soil moisture, and snow field.The DA is comprised of a horizontal optimal interpolation (Taillefer, 2002), which for soil moisture and temperature is followed by a vertical optimal interpolation procedure (Giard and Bazile, 2000).The upper-air DA is based on a 3-dimensional variational (3DVar) approach (Fischer et al., 2005).Many types of observations are assimilated including conventional in-situ measurements (pilot—balloon wind, radiosonde, aircraft, buoy, ship,and synop), Global Navigation Satellite System (GNSS)Zenith Total Delay (ZTD) data, weather radar reflectivity information, as well as infrared (IR) and PMW radiances from satellite-based instruments.The IR radiances are sensed by the Infrared Atmospheric Sounding Interferometer (IASI) placed on board the Metop satellites.PMW radiances are traditionally provided by the Advanced TIROS Operational Vertical Sounder (ATOVS) instrument family,including AMSU-A and MHS.Recently the DA system has been prepared to also utilize data from the MWHS-2 instrument on board the FY-3C and FY-3D polar orbiting satellites.Background error covariances are based on a climatological assumption and their representation is based on a multivariate formulation under the assumptions of horizontal homogeneity and isotrophy.They are calculated from an ensemble of forecast differences (Berre, 2000; Brousseau et al., 2012).These are produced by Ensemble Data Assimilation (EDA) experiments carried out with the HARMONIEAROME system.The HARMONIE-AROME EDA uses perturbed observations and ECMWF global EDA (Bonavita et al., 2012) forecasts as lateral boundary conditions.Scaling is applied to the derived statistics in order to be in agreement with the amplitude of HARMONIE-AROME + 3 h forecast errors (Brousseau et al., 2012).

    2.2.Handling of microwave radiances

    The NWP system uses atmospheric temperature and mois-ture information extracted from the PMW radiances sensed by satellite instruments summarized in Table 1.The AMSU-A instrument is primarily used to retrieve information on the vertical distribution of atmospheric temperature.The MHS instrument, on the other hand, provides information on the vertical structure of water vapour.Similarly to MHS, the MWHS-2 instrument has capability to retrieve information on atmospheric moisture.In addition, it has some capacity to retrieve information on temperature.To produce the model counterparts of the observed PMW radiances an observation operator,

    H

    , is applied.This operator is based on the radiative transfer model RTTOV (Radiative Transfer for TOVS) version 11.2.0 (Saunders et al., 2018)as developed under the EUMETSAT Satellite Application Facility to support NWP.At present, only clear-sky radiances with small contributions from surface can be efficiently handled by the observation operator.Furthermore,the observation operator requires that the measurement is sensitive mainly to atmospheric conditions below the model top.It implies that a clear-sky radiance with a major part(roughly 90%) of the integral of the Jacobians below the model top can be assimilated, if influence by the surface is small enough.The resulting channel usage from the different instruments is presented in Table 2.The AMSU-A channels used are the temperature sensitive ones located around 55 GHz while the MHS channels utilized are centered around 1 83 GHz.The used channels from MWHS-2 are mainly the moisture sensitive ones around the 1 83 GHz water vapor absorption line (channels 11 to 15).MWHS-2 channels 5 and 6, located around 1 1 8 GHz and sensitive also to temperature, are used with very low weight within the assimilation system.Some low-peaking channels (MHS chan-nel 5 and MWHS-2 channel 15) affected by the surface are used over sea only.The reason for this is that the surface emissivity and skin temperature contributions are better represented in the observation operator for radiances over sea than for radiances over land.Due to the problem of contamination by the surface, many of the PMW instrument channels (AMSU-A channel 6, MHS channels 3, 4, and MWHS-2 channels 5, 6, 11, 12, 13, and 14) are used over low level terrain only.The discrimination of low and high level terrain is done using the model orography over the field of view of the observation.

    Table 1.PMW radiance observation usage.

    To indicate whether radiance data from particular instruments and channels are affected by clouds, radiances in associated window channels within the instrument, capable of identifying clouds, are compared with the corresponding model state equivalents.If the window channel departures are larger than a particular threshold value, the corresponding nonwindow channel radiances are considered to be affected by clouds, and are therefore rejected from use in the DA.The window channel for AMSU-A channels 6 and 7 is AMSUA channel 4 (5 2.8 GHz).The window channel used for MHS channels 3, 4, and 5 is MHS channel 2 at 1 57 GHz.The window channel used for MWHS-2 channels 5 and 6 is MWHS-2 channel 1 (89 GHz) and for MWHS-2 channels 11, 12, 13, 14, and 15 the window channel is MWHS-2 channel 10 at 1 50 GHz.Note that for the rather high-peaking AMSU-A channels 8 and 9 no cloud detection is applied.For AMSU-A channel 6, a complementary check, aiming at identifying cloud and rain is applied, using liquid water path and scattering index.In case of too high liquid water path or enhanced scattering from precipitation-sized particles, radiance observations from AMSU-A channel 6 are considered contaminated by large hydrometeors and are therefore rejected.

    Systematic errors that might be present in the clear-sky radiances that have passed the cloud detection are handled by applying an adaptive Variational Bias Correction(VARBC) as proposed by Dee (2005) and further adapted by Auligné et al.(2007).A linear model of the following form is applied to describe the bias

    b

    in the PMW observations:

    The journey over the heath and the deep sand was only tooquickly at an end. They stopped before the house of mourning, wherethey found plenty of guests within and without. Waggon after waggonstood side by side, while the horses and oxen had been turned out tograze on the scanty pasture. Great sand-hills like those at home bythe North Sea rose behind the house and extended far and wide. How had they come here, so many miles inland? They were as large and high as those on the coast, and the wind had carried them there; there was also a legend attached to them.

    Table 2.PMW radiance channels used in the DA (window channels used in the data filtering are not included).

    Here

    p

    represents the predictors,

    x

    represents the model state, βrepresents the bias parameters, and

    N

    describes the number of predictors.Predictors are adopted to the channels and typically include atmospheric layer thicknesses and satellite instrument viewing angle.The bias parameters and the model initial state are simultaneously derived within the variational DA framework.A quality control procedure is applied to remove radiance observations that are considered to be of poor quality.Based on the operational monitoring experiences, we do not use AMSU-A and MHS radiances from field of views close to scan-line edges.Channels from instruments on specific satellite platforms might be temporarily or permanently blacklisted in the DA system.This blacklisting may occur due to known problems reported by the satellite agencies (NOAA-18 MHS channels 3-5, NOAA-19 MHS channel 3), by other collaborating NWP partners assimilating the same satellite radiances (NOAA-19 AMSU-A channels 7—8 and Metop-B channel 7), or due to quality limitations found during observation monitoring (i.e., where observed radiances are compared with corresponding model state equivalents over a longer period).It should be mentioned that the Metop-B AMSU-A channel 7 blacklisting is based on experiences with noisy radiances starting already in 2017.However, this noise is not presently observed, and the plan is to re-activate the assimilation of this channel.Further, as part of the quality control, a filtering is applied to get rid of radiance observations affected by gross errors.In this gross error check, a radiance observation,

    y

    , is rejected if it satisfies the following inequality:

    The bias corrected radiances that have passed the data selection and quality control are assumed to have a Gaussian error distribution.The associated observation error standard deviations are derived from long-term observation monitoring.The observation errors are comprised of instrument errors, representativity errors, persistence errors, and errors in the observation operator, and they are slightly inflated to account for the lack of representation of observation error correlations (Bormann and Bauer, 2010).For the AMSU-A channels, the estimated error standard deviations are approximately 0.2 K.For the MHS channels, the values are 1.8 K, for MWHS-2 channels 11—15 the values are 1.8-2.0 K, and for MWHS-2 channels 5 and 6 the estimated error standard deviations are 0.9 K.To further alleviate the effects on the initial state of spatially correlated observation errors, a thinning of radiances is applied.For radiances from the AMSUA and MHS instruments, a thinning distance of 8 0 km is used.This choice is based on previous experiences by Randriamampianina (2006).Due to lack of previous experience with the MWHS-2 instrument, as a first step we chose to apply the larger thinning distance of 1 6 0 km.This could probably be reduced at a later stage, based on a posteriori diagnosis of observation error correlations following Bormann and Bauer (2010).

    The rejections in the data selection and quality control are dominated by the cloud detection and the thinning procedure.Roughly 3% of the observations that have passed the data selection (except thinning, since it is applied after the quality control) are identified as gross errors and rejected by the gross error check described by Eq.2.

    3.Experimental design

    To evaluate the impact of an enhanced use of PMW radiances in the MetCoOp limited-area NWP system a parallel data assimilation and forecast experiment was designed.A version of the MetCoOp HARMONIE-AROME forecasting system was run over a domain illustrated in Fig.1.To avoid spin-up, the experiment was initiated 21 August 2019 from MetCoOp pre-operational initial state and bias-correction coefficients, cycled for many months.Thereafter, a passive assimilation of Metop-C (AMSU-A and MHS) radiances as well as FY3-C and FY3-D (MWHS-2) radiances was carried out in the experiment for the period of one month, until 21 September 2019.With passive assimilation of radiances from the newly introduced satellites and instruments, this means that the new data enters the DA in a way so that they are used to estimate and evolve the VARBC predictor coefficients only, and do not have an impact in the subsequent analyses.On 21 September 2019, the initial state files and the VARBC predictor coefficients from the one month spin-up experiment were used in two parallel experiments:

    REF:DA settings with reference system and observation usage as described in section 2.

    SAT: Like in REF, except that the new PMW type observations were also actively (both initial state and VARBC predictor coefficient were influenced by the new types of observations) assimilated in addition to the observations assimilated in REF.The additional new PMW data assimilated were from the AMSU-A and MHS instruments on board Metop-C and from the MWHS-2 instrument on board FY-3C and FY-3D.

    The parallel experiments REF and SAT were performed for a period extending from 21 September 2019 to 27 October 2019.During this period, forecasts up to a range of 36 h were launched four times a day, at 0000, 0600, 1200 and 1800 UTC.The first six days were excluded from the verification, giving some time for systems to do potential minor adjustments to the VARBC predictor coefficients,which was particularly important regarding the SAT experiment due to active assimilation of the newly introduced additional PMW radiances.

    4.Results

    4.1.Data coverage

    The REF satellite-based PMW radiance observation usage results in a rather uneven distribution of microwave observations between the assimilation cycles.Some cycles are, not at all, or only by a small fraction, covered by PMW satellite radiances.This is illustrated by the left part of Fig.2,which shows, for REF, the fraction of the model domain covered by satellite-based PMW radiances, assuming a 20 min latency (or timeliness, time from observation is made until being available for use in model), which is considered realistic.The right part of Fig.2 shows the corresponding fractions covered for SAT experiment, when satellite-based PMW radiances from Metop-C, FY-3C, and FY-3D are also being used.With these additional PMW observations included there is a more even distribution, between different assimilation cycles, of the fraction of the area covered by PMW radiances.In SAT, almost 80% of the domain is covered by PMW radiance observations for all assimilation cycles.In particular, for the 0000 UTC assimilation cycle, a large part of the domain is covered by PMW data in SAT,whereas there are no PMW observations at all in REF.The additional PMW satellite radiances thus have the potential to improve forecast quality by filling existing data gaps.

    4.2.Observation monitoring

    One way to evaluate the quality of the observations in the DA is to compare the observed values with the corresponding model state equivalents.This was done for all the PMW instruments and channels used.Results are shown for AMSU-A in Fig.3, MHS in Fig.4, and MWHS-2 in Fig.5.Channels for which statistics are lacking are permanently rejected by our DA system since they are considered to be of poor quality.Note that the innovation (observation minus background equivalent departure) statistics do basically have a Gaussian distribution so that the standard deviation of the innovations can be estimated from where the distribution has dropped to roughly one third of its maximum value.For unbiased data with observation errors uncorrelated with background errors, the square of the innovation standard deviations for a particular channel is further built up by the sum of the square of the observation error standard deviation and the square of the standard deviation of the error for the background observation equivalent.Taking into account the slight inflation of the observation error standard deviations,to compensate for the lack of representation of observation error correlations, the width of the innovations are consistent with the estimated observation error standard deviations(Figs.3, 4 and 5).Innovations are smaller for temperature sensitive channels (Fig.3) than for moisture sensitive channels (Figs.4 and 5).Statistics from AMSU-A on Metop satellites look very similar to that of AMSU-A statistics from NOAA-18 and NOAA-19 satellites (not shown).For MWHS-2, only the moisture sensitive channels are shown.From Fig.3 it is clear that the newly introduced Metop-C AMSU-A radiances are of at least similar quality as the corresponding AMSU-A measurements from instruments onboard the other Metop satellites already used in the system.

    The MHS Metop-C channel 3 and 4 radiances are of slightly worse quality than the corresponding radiances from the other Metop satellites and NOAA-19.It is clear from the width of the histograms that the observation error standard deviations for Metop-C MHS channels 3 and 4 should be increased by a few tenths of a Kelvin in the future due to larger instrument errors.As a starting point,however, the same observation error standard deviations are applied for Metop-C as for the other satellites.The MHS channel 5 distributions are skewed for all satellites.The reason is that only observations over sea are assimilated, but close to the coast, they are likely influenced by land surface, which causes an undesirable deviation between observation and model counterpart in monitoring and DA.This will be improved by refined data selection in future versions of the system, taking the satellite footprint into account.TheMWHS-2 instrument channels 11-15 are demonstrated to be of a similar quality to those of MHS.Again, skewed distributions appear for channel 15, for the same reason as for MHS channel 5.

    Fig.2.PMW radiance observation coverage over MetCoOp domain for different assimilation cycles with current(left) and enhanced (right) PMW radiance observation usage and with operational cut-off settings and a 20 min latency.

    Fig.3.Normalized histograms of AMSU-A brightness temperature innovation statistics (units: K) based on all data within the domain that have passed the quality control during the period 0000 UTC 29 September 2019 to 2100 UTC 9 October 2019.Different rows represents different channels and different columns represents different Metop satellites.

    4.3.Observation impact on analyses

    The impact of observations on the analysis system can be evaluated using the degrees of freedom for signal (DFS)(Randriamampianina et al., 2011).DFS is the derivative of the analysis increments in observation space with respect to the observations used in the analysis system.As proposed by Chapnik et al.(2006), DFS can be computed through a randomization technique, as follows:

    Here,

    y

    and

    y

    ~ represent the vector of unperturbed and perturbed observations, respectively.Similarly,

    x

    and

    x

    ~are the analyses produced when the unperturbed and perturbed observations are used in the DA, respectively.

    H

    represents a tangent-linear version of the observation operator

    H

    , linearised around the background state, while

    R

    is a matrix describing the observation errors.Randriamampianina et al.(2011)discussed the potential of both absolute and relative DFS values in evaluating the impact of observations on the DA.

    The DFS values vary depending on the assimilation cycles due to differences in observation coverage and also from day to day due to variations in the meteorological situation.We calculate the DFS subdivided into various observation types based on data from the three selected dates 17,20, and 24 October 2019.The days are separated well in time to provide independent weather situations.The DFS was calculated for each type of observation assimilated in REF and SAT.The calculation was done based on all eight assimilation cycles within each of the three days, as presented in Fig.6.In terms of absolute DFS, the contribution from satellite PMW radiances has clearly increased in SAT as compared to REF.The satellite PMW radiances can be seen to be the fourth most influential observation category on the analysis, after radar reflectivities, satellite IR radiances, and aircraft reports.One should keep in mind that at 0000 UTC there are very few aircraft observations available and no Metop satellite overpasses, and hence, no IR satel-lite data in neither of the parallel experiments.Satellite PMW data at 0000 UTC are available only in SAT, and together with radiosondes and radar reflectivities are the largest contribution to absolute DFS (not shown).Thus,there is an enhanced effect on the initial state of the additional observations, in particular at 0000 UTC.By including the additional PMW observations, the DFS of almost all the other observation types increases, indicating a consistency between the additional PMW radiances and the other observation types.It seems that the additional PMW observations also contribute by increasing the impact on the analysis of the other observation types.The exception is satellite IR radiances, for which the DFS slightly decreases when including additional PMW radiances.One potential explanation for the reduced impact of IR radiances when introducing additional PMW radiances is interactions through VARBC in combination with a relatively small number of anchoring observations.The reduction in DFS for IR radiances might also indicate some kind of deficiency in the handling of IR radiances, such as sub-optimal cloud-detection procedures.This is something that will be studied in more detail in the future.

    Fig.4.Normalized histograms of MHS brightness temperature innovation statistics (units: K) based on all data within the domain that have passed the quality control during the period 0000 UTC 1 October 2020 to 2100 UTC 9 October 2020.Different rows represents different satellites and different columns represents different channels.

    4.4.Observation impact on the forecasts

    To evaluate the quality of the forecasts from the parallel experiment, we verified them against radiosonde and synoptic weather observations within the model domain.The verification was carried out for surface and upper-air model variables.Special emphasis was put on verification of humidity,clouds, precipitation, and temperature.Results revealed that the largest differences in verification statistics between the two parallel experiments were found for forecast ranges up to +18 h.In Fig.7 the verification statistics for +12 h temperature and specific humidity forecasts are shown as a function of vertical level, for verification against radiosondes.

    In terms of both bias and standard deviation, the temperature forecasts of SAT and REF are of a similar quality.On the other hand, with respect to standard deviations below 700 hPa, SAT forecasts are better than REF forecasts.Above 700 hPa humidity forecast quality is rather similar for SAT and REF, in terms of bias as well as standard deviation.The largest impact on low-level humidity forecasts can be explained by most of the humidity sensitive Jacobians peaking between 600 and 800 hPa in combination with background error humidity standard deviation profiles (not shown) having the largest values between roughly 700 and 900 hPa.At 700 hPa the bias is negative for both SAT and REF but the magnitude (absolute value) is larger for SAT than for REF.Thus, compared to radiosonde observations,both SAT and REF are slightly drier around 700 hPa, and SAT is the driest.One reason for the slightly more dry SAT forecast as compared to REF could be the too high model background equivalents close to coastlines for low-peakingmoisture sensitive channels (Figs.4 and 5), causing negative moisture increments.However, the magnitude of the bias actually seems to be smaller for SAT below the 800 hPa level.

    Fig.5.Normalized histograms of MWHS-2 brightness temperature innovation statistics (units: K)based on all data within the domain that have passed the quality control during the period 0000 UTC 1 October 2020 to 2100 UTC 9 October 2020.Different rows represents different satellites and different columns represents different channels.

    Fig.6.DFS subdivided into various observation types for the two experiments REF and SAT.Results were based on data from three different dates/cases, including all day (i.e., 24 different data assimilation cycles).Where SYNOP—the surface weather; GNSS ZTD—ground-based zenith total delay; TEMP—radiosonde; AIREP—aircraft; DRIBU—drifting Buoy; SAT MW—PMW radiances; SAT IR—IR radiances; and RAD RFL—radar reflectivity observations.

    Fig.7.Verification statistics, in terms of bias (BIAS) and standard deviation (STDV), of +12 h forecasts against radiosonde observations of temperature (left, units: K) and specific humidity (right, units: g kg-1) averaged over all observations within the domain and over the one month period.Scores are shown as function of vertical level and dark grey solid lines are for REF while light grey solid lines represent SAT.The grey dashed line illustrates the number of observations used in the verification.

    Verification scores for forecasts of +12 h total cloud cover and 12 h accumulated precipitation between +6 h and+18 h are better for SAT than for REF, as shown in Fig.8,in terms of Kuiper skill score, when verifying against synop land weather stations in the domain over the one month period.Both the cloud cover and accumulated precipitation forecasts are slightly better for SAT than for REF for almost all thresholds.These results are consistent for cloud and precipitation and with better low-level specific humidity forecasts in SAT as compared to REF.

    Although low-level SAT moisture forecasts are better on average than REF forecasts during the month period, it should be noted that there is a considerable variation in verification scores between different days.In addition, one can identify some differences between forecasts launched from 0000 and 1200 UTC.Figure 9 illustrates a clear daily variation for +12 h 925 hPa relative humidity forecasts launched from 00 UTC (left) and 12 UTC (right).In terms of magnitude of standard deviations for forecasts launched at 00 UTC, SAT performs better than REF in 63% of thecases.For forecasts launched at 12 UTC, SAT forecasts are better than REF forecasts in 56% of the cases.The improved SAT forecast as compared with REF is particularly evident for forecasts launched at 00 UTC and with valid times between 17 and 24 October.One contributing factor for this difference between 00 and 12 UTC based forecasts is likely that at 00 UTC there are no humidity sensitive PMW observations assimilated in REF while there is a substantial amount of such observations available for assimilation in SAT (recall Fig.2).

    Fig.8.Kuiper skill score for +12 h total cloud cover forecasts (left) and 12 h accumulated precipitation forecasts for accumulation between +6 h and +18 h forecast ranges (right).Dark grey line is for REF and light grey for SAT.

    4.5.Case study

    As an example of impact of assimilation of the additional PMW observations on forecast quality, one particular case has been selected.We illustrate in Fig.10 the precipitation forecasts launched from 0000 UTC 18 October 2019 and valid at 19 October between 0000 and 0300 UTC.The forecasts are for 3 h accumulated precipitation (unit: mm).For comparison, also the corresponding weather radar derived 3 h accumulated precipitation is shown.

    The weather situation is characterised by a synoptic scale cyclone, with associated frontal systems, over the southern part of the model domain and moving towards northeast.The event resulted in substantial precipitation amounts in the southern parts of Norway as well as in southern and central parts of Sweden.The main precipitation pattern was rather well predicted, both in terms of amount and pattern.In front of this main precipitation pattern there are also twoother frontal structures, accompanied by precipitation.One of these is situated close to the eastern coast and over the Baltic sea in the middle part of Sweden, and one is situated over the northern parts of Sweden and Finland.Interestingly, the substantial precipitation amounts close to the coast in eastern Sweden to the north-east of the main frontal structure (and marked with a red dotted circle in Fig.10) are better predicted with the SAT-based forecast than with the REF-based forecast.It can be argued that the exact position of such small-scale precipitation structures is in general not predictable at a forecast range of 24 h.Nevertheless, the assimilation of additional PMW radiances in SAT as compared to REF might have provided improved lower-level atmospheric moisture fields, enabling improved precipitation forecasts.

    Fig.9.Time variability of the +12 h relative humidity forecast bias and standard deviation scores (unit: %) at the vertical level of 925 hPa during the one month period (27 September—27 October, 2019) for verification against radiosonde observations.The scores are for the experiments REF (dark grey solid line) and SAT (light grey solid line).The grey dashed curve illustrates the number of observations used within the verification.Left panel is for forecast launched from 0000 UTC, and right panel is for forecast launched from 1200 UTC.

    Fig.10.Prediction of +24 h to +27 h 3 h accumulated precipitation (units: mm) for SAT (middle) and REF (right).Predictions are launched from 0000 UTC 18 October 2019 and valid between 0000 and 0300 UTC 19 October 2019.The corresponding radar derived accumulated precipitation (units: mm) as shown in left panel is derived from the radar reflectivity observations from the weather radar network for the Baltic Sea Region (BALTRAD).The grey shaded area represents the coverage of the radar network.The red dotted circle highlights area of improvement with SAT as compared to that of REF.

    5.Concluding remarks and outlook

    The use of PMW radiances from the MWHS-2 instrument in a northern European regional km-scale forecasting system is presented.The handling of radiance observations is described, including several aspects, such as cloud detection, systematic observation errors, error correlations, and exclusion of gross errors.

    Extending the regional forecasting system with assimilation of MWHS-2 microwave radiance observations as well as Metop-C MHS and AMSU-A radiances was demonstrated to be beneficial in different ways.It was shown that the newly introduced observations fill the lack of PMW satellite data over the northern European domain during midnight and early morning, which leads to a more even availability of PMW radiances in DA over the course of the day.Observation monitoring results for humidity sensitive PMW channels point out that the quality of the MWHS-2 is comparable to the already assimilated MHS instruments, on board the NOAA-19, Metop-A, and Metop-B satellites.On the other hand, MHS channels 3 and 4 on Metop-C seem to be of slightly worse quality than the corresponding MHS observations from other satellites already used in operations.Through the DFS diagnostic we have shown that the additional PMW radiances also had a clear impact on the analyses.Finally, it was demonstrated through verification scores that the additional observations do have a positive impact on the short-range forecasts of humidity, clouds, and precipitation, in a statistical sense.The potential impact of the additional microwave radiances was demonstrated through a case study where an example of prediction of precipitation was discussed.Based on the work and findings presented here, radiances from MHS and AMSU-A from Metop-C as well as MWHS-2 on board FY-3C and FY-3D are now assimilated in a pre-operational version of the regional forecasting system.

    Despite the encouraging results obtained, there is room for further improvements regarding the handling of satellitebased humidity-sensitive PMW observations in our regional NWP model.As a first step we plan to investigate the horizontal correlations of observation errors for the MWHS-2 radiances using a posteriori diagnosis of observation error correlations following Bormann and Bauer (2010).From such a study, we hope to justify the application of a shorter horizontal thinning distance than the currently used 160 km for MWHS-2.The thinning procedures, when rejecting radiances affected by land surface from low-peaking humidity sensitive channels, needs to be improved to better take the actual satellite footprint into account.In the long term, important future improvements are the assimilation of all-sky radiances and introduction of flow-dependent data assimilation techniques, such as 4D-Var (Gustafsson et al., 2012) or a hybrid ensemble-variational data assimilation (Gustafsson et al., 2014).

    Within MetCoOp, the researchers in regional data assimilation are working towards the introduction of an NWP based nowcasting system, with data assimilation carried out each hour and with a 15 min observation cut-off time.Figure 11 shows the PMW radiance observation coverage and availability for such a nowcasting system.Clearly theextended observation usage will also benefit such a system,but there is still room for enhancements regarding the availability of PMW data.A system of many small polar satellites equipped with a PMW instrument would be sufficient for filling these many remaining gaps in a cost-effective way.The AWS concept has been born out of exactly such a need.The prototype AWS is expected to be launched in early 2024, and a follow-on constellation, of up to 20 small satellites providing observational coverage over the Arctic every 30 min, is considered for 2025 and onwards.

    Fig.11.PMW radiance observation coverage and availability over MetCoOp area for different assimilation cycles with current (left) and enhanced (right) PMW radiance observation usage and with nowcasting cut-off settings and a 20 min latency.

    Acknowledgements.

    The work has been carried out within MetCoOp.We acknowledge Philippe Chambon for support and discussions regarding handling of PMW radiances in HARMONIEAROME.We are grateful for technical assistance from Eoin WHELAN, Frank GUILLAUME, Ulf ANDR?, and Ole VIGNES.We are grateful to Susanna HAGELIN for support on improving the readability of the manuscript.We also thank the anonymous reviewers for useful comments.

    Open Access

    This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source,provide a link to the Creative Commons license, and indicate if changes were made.

    97碰自拍视频| 精品国产乱子伦一区二区三区| 757午夜福利合集在线观看| 亚洲精品美女久久久久99蜜臀| 99国产综合亚洲精品| 一级作爱视频免费观看| 国产精品野战在线观看| 夜夜看夜夜爽夜夜摸| 国产不卡一卡二| 亚洲无线在线观看| 非洲黑人性xxxx精品又粗又长| 国产蜜桃级精品一区二区三区| 色综合站精品国产| 制服诱惑二区| 久久午夜综合久久蜜桃| 999久久久精品免费观看国产| 婷婷六月久久综合丁香| 免费看十八禁软件| 久久热在线av| 免费在线观看影片大全网站| 老汉色av国产亚洲站长工具| 日本一区二区免费在线视频| 一级毛片高清免费大全| 2021天堂中文幕一二区在线观 | 色播在线永久视频| 18禁观看日本| 免费高清在线观看日韩| 国产视频内射| 一本综合久久免费| 国产一区在线观看成人免费| 母亲3免费完整高清在线观看| 精品久久蜜臀av无| 亚洲精品美女久久av网站| 黄片大片在线免费观看| 脱女人内裤的视频| 自线自在国产av| 国产黄a三级三级三级人| 欧美色视频一区免费| 观看免费一级毛片| 99久久久亚洲精品蜜臀av| 深夜精品福利| 亚洲色图 男人天堂 中文字幕| 久久午夜综合久久蜜桃| 成人av一区二区三区在线看| 日韩视频一区二区在线观看| 一级a爱视频在线免费观看| 黄色 视频免费看| 国产精品二区激情视频| 午夜激情福利司机影院| 女人爽到高潮嗷嗷叫在线视频| 自线自在国产av| 琪琪午夜伦伦电影理论片6080| 中文字幕人妻熟女乱码| 欧美日本视频| 美女 人体艺术 gogo| 夜夜躁狠狠躁天天躁| 久久久久久免费高清国产稀缺| 在线免费观看的www视频| 丝袜人妻中文字幕| 精品国内亚洲2022精品成人| 精品一区二区三区视频在线观看免费| 国产精品九九99| 一级作爱视频免费观看| av电影中文网址| 亚洲自拍偷在线| 91av网站免费观看| 国产精品 欧美亚洲| 国语自产精品视频在线第100页| 久久精品91无色码中文字幕| 黄色丝袜av网址大全| 在线观看66精品国产| 丰满的人妻完整版| 麻豆国产av国片精品| 中文字幕久久专区| 国产男靠女视频免费网站| 视频区欧美日本亚洲| 99久久国产精品久久久| 亚洲国产精品合色在线| 首页视频小说图片口味搜索| 日韩国内少妇激情av| 久久国产精品影院| 香蕉av资源在线| 久久中文字幕一级| 给我免费播放毛片高清在线观看| 一级毛片高清免费大全| 欧美日本亚洲视频在线播放| 最近最新免费中文字幕在线| 国产激情欧美一区二区| 午夜福利18| 在线永久观看黄色视频| 好男人在线观看高清免费视频 | 琪琪午夜伦伦电影理论片6080| 日本一区二区免费在线视频| 国产精品乱码一区二三区的特点| 夜夜看夜夜爽夜夜摸| 日韩欧美三级三区| 国产成+人综合+亚洲专区| 一进一出抽搐gif免费好疼| 久久香蕉国产精品| 免费看日本二区| www日本在线高清视频| 男人舔女人下体高潮全视频| 制服诱惑二区| 国产免费男女视频| 香蕉av资源在线| 精品欧美国产一区二区三| 中文字幕久久专区| 国产成人系列免费观看| 禁无遮挡网站| av电影中文网址| 久久人人精品亚洲av| 性色av乱码一区二区三区2| 老汉色∧v一级毛片| 99re在线观看精品视频| 国产精华一区二区三区| 9191精品国产免费久久| 久久精品aⅴ一区二区三区四区| 免费看十八禁软件| 色播在线永久视频| 少妇被粗大的猛进出69影院| 亚洲成人久久爱视频| 午夜日韩欧美国产| x7x7x7水蜜桃| 欧美国产精品va在线观看不卡| 男人舔女人下体高潮全视频| 久久久国产成人免费| 国产成人精品久久二区二区免费| 亚洲片人在线观看| 真人一进一出gif抽搐免费| 亚洲中文日韩欧美视频| 久久亚洲精品不卡| 女警被强在线播放| 成年版毛片免费区| 日韩三级视频一区二区三区| 欧美黄色淫秽网站| 午夜激情av网站| 国产91精品成人一区二区三区| 中文字幕精品免费在线观看视频| av中文乱码字幕在线| 搡老妇女老女人老熟妇| 国内精品久久久久精免费| 韩国精品一区二区三区| 欧美不卡视频在线免费观看 | 中文资源天堂在线| www.自偷自拍.com| 最近最新免费中文字幕在线| 热re99久久国产66热| 无遮挡黄片免费观看| 欧美黑人巨大hd| 国产成人av激情在线播放| 大型黄色视频在线免费观看| 最近最新免费中文字幕在线| 91在线观看av| 欧美激情 高清一区二区三区| 操出白浆在线播放| 国产黄片美女视频| 女生性感内裤真人,穿戴方法视频| 亚洲人成伊人成综合网2020| 在线观看免费视频日本深夜| 丝袜美腿诱惑在线| 亚洲人成网站高清观看| 久久天躁狠狠躁夜夜2o2o| 欧美中文综合在线视频| 国产爱豆传媒在线观看 | 亚洲成a人片在线一区二区| 手机成人av网站| 久久久久久久久久黄片| 免费在线观看亚洲国产| 欧美日韩精品网址| 久久伊人香网站| 亚洲色图av天堂| 亚洲国产欧洲综合997久久, | 两个人免费观看高清视频| 国产色视频综合| 一区二区三区国产精品乱码| 午夜福利一区二区在线看| 亚洲一卡2卡3卡4卡5卡精品中文| 婷婷亚洲欧美| 国产精品日韩av在线免费观看| 成年女人毛片免费观看观看9| 在线观看免费视频日本深夜| 午夜亚洲福利在线播放| 精品国内亚洲2022精品成人| 一区福利在线观看| 欧美不卡视频在线免费观看 | 日韩欧美免费精品| 正在播放国产对白刺激| 午夜福利18| 这个男人来自地球电影免费观看| 国产亚洲精品一区二区www| 一进一出好大好爽视频| 国产精品久久视频播放| 免费在线观看成人毛片| 成人特级黄色片久久久久久久| 亚洲av日韩精品久久久久久密| 亚洲成人久久性| 精品国产乱码久久久久久男人| 正在播放国产对白刺激| av中文乱码字幕在线| 中国美女看黄片| 波多野结衣巨乳人妻| 亚洲国产看品久久| 非洲黑人性xxxx精品又粗又长| 国语自产精品视频在线第100页| 国产精品九九99| 99热这里只有精品一区 | 午夜免费鲁丝| 国产主播在线观看一区二区| 99久久精品国产亚洲精品| 女性被躁到高潮视频| 日韩三级视频一区二区三区| 亚洲电影在线观看av| 美女高潮到喷水免费观看| 成人亚洲精品av一区二区| 老司机靠b影院| а√天堂www在线а√下载| xxx96com| 精华霜和精华液先用哪个| 国产色视频综合| 亚洲自偷自拍图片 自拍| 国产成年人精品一区二区| 日韩中文字幕欧美一区二区| 国产成人精品无人区| 欧美最黄视频在线播放免费| 18禁国产床啪视频网站| 一级毛片高清免费大全| av欧美777| 18禁黄网站禁片午夜丰满| 国产精品一区二区精品视频观看| 俄罗斯特黄特色一大片| 真人做人爱边吃奶动态| 亚洲,欧美精品.| 日本撒尿小便嘘嘘汇集6| 一边摸一边做爽爽视频免费| 日韩欧美免费精品| 人人妻人人看人人澡| 色综合亚洲欧美另类图片| e午夜精品久久久久久久| 亚洲av成人av| 精品无人区乱码1区二区| 变态另类成人亚洲欧美熟女| 亚洲精品中文字幕一二三四区| 黄色a级毛片大全视频| 日韩成人在线观看一区二区三区| 99久久国产精品久久久| 老司机靠b影院| 亚洲欧美精品综合一区二区三区| 国产伦人伦偷精品视频| 久久青草综合色| 日韩欧美一区视频在线观看| 亚洲三区欧美一区| 极品教师在线免费播放| 精品一区二区三区av网在线观看| 午夜a级毛片| 热99re8久久精品国产| 在线观看66精品国产| 中文字幕另类日韩欧美亚洲嫩草| 亚洲 欧美一区二区三区| 欧美黑人精品巨大| 最新在线观看一区二区三区| 欧美+亚洲+日韩+国产| 欧美成人一区二区免费高清观看 | 久久午夜综合久久蜜桃| 亚洲国产精品999在线| 亚洲人成电影免费在线| 日本成人三级电影网站| 99国产精品一区二区蜜桃av| 国产免费男女视频| 国产精品日韩av在线免费观看| 亚洲va日本ⅴa欧美va伊人久久| 99在线人妻在线中文字幕| 亚洲精品国产一区二区精华液| 欧美又色又爽又黄视频| 亚洲成人久久爱视频| 亚洲成人免费电影在线观看| av超薄肉色丝袜交足视频| 国产伦在线观看视频一区| 变态另类丝袜制服| 黑人欧美特级aaaaaa片| 色综合亚洲欧美另类图片| 成人av一区二区三区在线看| 99国产精品一区二区蜜桃av| 91九色精品人成在线观看| 亚洲国产欧洲综合997久久, | 深夜精品福利| 久久人妻av系列| 丁香欧美五月| 性欧美人与动物交配| 国产高清视频在线播放一区| 亚洲成av人片免费观看| 亚洲中文字幕日韩| 这个男人来自地球电影免费观看| av免费在线观看网站| 欧美激情高清一区二区三区| 好看av亚洲va欧美ⅴa在| 久久精品影院6| 中文字幕最新亚洲高清| 精品久久久久久久末码| 91在线观看av| 嫩草影院精品99| 脱女人内裤的视频| 久久精品影院6| 国产又黄又爽又无遮挡在线| 国产精品久久久人人做人人爽| 国产97色在线日韩免费| 欧美又色又爽又黄视频| 精品无人区乱码1区二区| 丝袜在线中文字幕| 黄片大片在线免费观看| 日韩中文字幕欧美一区二区| 嫁个100分男人电影在线观看| 免费无遮挡裸体视频| 免费观看人在逋| 人人妻人人看人人澡| 长腿黑丝高跟| xxx96com| 中文字幕人妻丝袜一区二区| 人人澡人人妻人| 精品国产美女av久久久久小说| 人妻久久中文字幕网| 国产野战对白在线观看| 久久久久久亚洲精品国产蜜桃av| 免费看a级黄色片| 国产精品美女特级片免费视频播放器 | 岛国视频午夜一区免费看| 久99久视频精品免费| 成人18禁在线播放| 国产又色又爽无遮挡免费看| 色综合欧美亚洲国产小说| 色在线成人网| 女人被狂操c到高潮| 嫩草影视91久久| 国产av又大| 久久香蕉精品热| 欧美成人一区二区免费高清观看 | 男男h啪啪无遮挡| 亚洲精品美女久久久久99蜜臀| 成人18禁高潮啪啪吃奶动态图| 久久99热这里只有精品18| 亚洲av电影不卡..在线观看| 精品国产一区二区三区四区第35| 欧美成人免费av一区二区三区| 亚洲全国av大片| 脱女人内裤的视频| 亚洲精品美女久久久久99蜜臀| 欧美大码av| 他把我摸到了高潮在线观看| 中文字幕最新亚洲高清| 18禁国产床啪视频网站| 精品久久久久久久久久免费视频| 久久久久久久久久黄片| 高潮久久久久久久久久久不卡| 国产精品免费一区二区三区在线| 在线播放国产精品三级| 久久中文看片网| 18禁黄网站禁片免费观看直播| 国产真实乱freesex| 免费在线观看影片大全网站| 精品熟女少妇八av免费久了| 亚洲专区中文字幕在线| 国产日本99.免费观看| 久久久久亚洲av毛片大全| 18禁黄网站禁片午夜丰满| 日本免费一区二区三区高清不卡| 一进一出抽搐动态| 国产男靠女视频免费网站| 99久久精品国产亚洲精品| 久久久久亚洲av毛片大全| 99热6这里只有精品| 亚洲专区中文字幕在线| 99久久精品国产亚洲精品| av中文乱码字幕在线| 97超级碰碰碰精品色视频在线观看| 午夜福利在线在线| 最新在线观看一区二区三区| 男人操女人黄网站| 日本熟妇午夜| 日韩av在线大香蕉| 婷婷精品国产亚洲av| 人成视频在线观看免费观看| 国产精品99久久99久久久不卡| 国产精品 国内视频| 一级a爱视频在线免费观看| 久久精品成人免费网站| 精品国产超薄肉色丝袜足j| 亚洲成av人片免费观看| 午夜福利高清视频| 亚洲人成伊人成综合网2020| 成人三级黄色视频| 成人永久免费在线观看视频| 欧美zozozo另类| 曰老女人黄片| 久久精品成人免费网站| 国产一区二区三区视频了| 国产熟女午夜一区二区三区| 国内少妇人妻偷人精品xxx网站 | 嫁个100分男人电影在线观看| 日本黄色视频三级网站网址| 1024手机看黄色片| 嫩草影视91久久| 欧美色欧美亚洲另类二区| 91成人精品电影| 十分钟在线观看高清视频www| 色在线成人网| 日韩欧美在线二视频| 精品国产美女av久久久久小说| 91大片在线观看| 欧美黑人欧美精品刺激| 国产激情偷乱视频一区二区| 久久人妻福利社区极品人妻图片| 国产三级黄色录像| 国产久久久一区二区三区| 一区福利在线观看| 九色国产91popny在线| 精品国产美女av久久久久小说| 天堂影院成人在线观看| 国产欧美日韩精品亚洲av| 免费一级毛片在线播放高清视频| 免费看美女性在线毛片视频| 老司机在亚洲福利影院| 两性夫妻黄色片| 日本免费一区二区三区高清不卡| 欧美丝袜亚洲另类 | www国产在线视频色| 天天躁狠狠躁夜夜躁狠狠躁| 欧美+亚洲+日韩+国产| 男人操女人黄网站| av天堂在线播放| 性欧美人与动物交配| 久久狼人影院| 中文在线观看免费www的网站 | 18美女黄网站色大片免费观看| 精品不卡国产一区二区三区| 最近在线观看免费完整版| 国内少妇人妻偷人精品xxx网站 | 可以在线观看的亚洲视频| 国产av在哪里看| 99久久综合精品五月天人人| 久久天堂一区二区三区四区| 窝窝影院91人妻| 欧美性猛交╳xxx乱大交人| www日本在线高清视频| 国产精品久久视频播放| 亚洲精品国产区一区二| 热99re8久久精品国产| 亚洲欧美日韩无卡精品| 老熟妇仑乱视频hdxx| 欧美成人一区二区免费高清观看 | or卡值多少钱| 日日爽夜夜爽网站| 国产精品av久久久久免费| 国产精品乱码一区二三区的特点| 日韩大码丰满熟妇| 12—13女人毛片做爰片一| 很黄的视频免费| 啦啦啦免费观看视频1| www.熟女人妻精品国产| 国产午夜精品久久久久久| 亚洲av成人一区二区三| 国产精品久久久久久人妻精品电影| 亚洲专区国产一区二区| 久久香蕉激情| 国产成人系列免费观看| 国语自产精品视频在线第100页| 精品久久久久久久末码| 一进一出抽搐gif免费好疼| 亚洲免费av在线视频| 男女下面进入的视频免费午夜 | 97人妻精品一区二区三区麻豆 | 99re在线观看精品视频| 黄片播放在线免费| 亚洲人成77777在线视频| 久久久久国产精品人妻aⅴ院| 久久久久国内视频| 男人舔女人的私密视频| 成熟少妇高潮喷水视频| 亚洲成av片中文字幕在线观看| 热re99久久国产66热| 成人午夜高清在线视频 | 99久久精品国产亚洲精品| 男女下面进入的视频免费午夜 | 婷婷六月久久综合丁香| 国产精品电影一区二区三区| 黄频高清免费视频| 日本黄色视频三级网站网址| 国产午夜福利久久久久久| 精品欧美国产一区二区三| 2021天堂中文幕一二区在线观 | 天堂√8在线中文| 一本综合久久免费| 高清毛片免费观看视频网站| 欧美又色又爽又黄视频| 色综合亚洲欧美另类图片| 99精品在免费线老司机午夜| 国产黄色小视频在线观看| 侵犯人妻中文字幕一二三四区| tocl精华| 午夜影院日韩av| 精品久久久久久久久久久久久 | 国产片内射在线| 黄色视频不卡| 色婷婷久久久亚洲欧美| 制服诱惑二区| 1024视频免费在线观看| 成人国产综合亚洲| 好男人在线观看高清免费视频 | 中文字幕最新亚洲高清| 精品久久久久久久毛片微露脸| 亚洲七黄色美女视频| 成人特级黄色片久久久久久久| 免费在线观看亚洲国产| 日韩欧美一区二区三区在线观看| 久久 成人 亚洲| 久久久久久久精品吃奶| 国产精品综合久久久久久久免费| 亚洲精品中文字幕在线视频| 欧美日韩一级在线毛片| avwww免费| 国产精品电影一区二区三区| 99国产精品一区二区蜜桃av| 午夜视频精品福利| 特大巨黑吊av在线直播 | 国产亚洲精品av在线| 啪啪无遮挡十八禁网站| 大型黄色视频在线免费观看| √禁漫天堂资源中文www| 女生性感内裤真人,穿戴方法视频| 欧美日韩瑟瑟在线播放| av福利片在线| 一区福利在线观看| 一区二区三区国产精品乱码| 精品久久久久久久末码| 人人妻,人人澡人人爽秒播| 国内久久婷婷六月综合欲色啪| 国产亚洲欧美在线一区二区| 国产精品98久久久久久宅男小说| 午夜老司机福利片| 黄色女人牲交| 91在线观看av| 高清毛片免费观看视频网站| 亚洲自拍偷在线| 国产午夜福利久久久久久| 麻豆国产av国片精品| 久久国产精品男人的天堂亚洲| 亚洲欧美日韩高清在线视频| 两个人视频免费观看高清| 久久婷婷成人综合色麻豆| 欧美一级毛片孕妇| 99在线人妻在线中文字幕| 成人国产综合亚洲| 国产麻豆成人av免费视频| 亚洲精品久久成人aⅴ小说| 麻豆一二三区av精品| 最近最新免费中文字幕在线| 午夜免费成人在线视频| 九色国产91popny在线| 少妇裸体淫交视频免费看高清 | 18禁观看日本| 欧美国产日韩亚洲一区| 久久久精品欧美日韩精品| 黄网站色视频无遮挡免费观看| 国产精品亚洲av一区麻豆| 亚洲国产欧洲综合997久久, | 香蕉国产在线看| 美女免费视频网站| 成人18禁高潮啪啪吃奶动态图| 少妇裸体淫交视频免费看高清 | 国产成人影院久久av| 人人妻人人澡欧美一区二区| 午夜久久久久精精品| 制服丝袜大香蕉在线| 悠悠久久av| 国产久久久一区二区三区| 日韩视频一区二区在线观看| 动漫黄色视频在线观看| 色综合欧美亚洲国产小说| 久久国产乱子伦精品免费另类| 最近最新免费中文字幕在线| 亚洲真实伦在线观看| 男女那种视频在线观看| 狠狠狠狠99中文字幕| 午夜福利高清视频| 国产高清videossex| 91av网站免费观看| 12—13女人毛片做爰片一| 日韩欧美在线二视频| 黄色视频,在线免费观看| 妹子高潮喷水视频| 老司机福利观看| a级毛片a级免费在线| 亚洲电影在线观看av| 又黄又粗又硬又大视频| 桃红色精品国产亚洲av| 巨乳人妻的诱惑在线观看| 久久 成人 亚洲| 高清毛片免费观看视频网站| 巨乳人妻的诱惑在线观看| 精品高清国产在线一区| 在线观看免费日韩欧美大片| 99re在线观看精品视频| 亚洲av片天天在线观看| 国产色视频综合| 久久精品夜夜夜夜夜久久蜜豆 | www日本在线高清视频| 热99re8久久精品国产| 久久久久久久久久黄片| 老熟妇乱子伦视频在线观看| 97超级碰碰碰精品色视频在线观看| 亚洲熟妇中文字幕五十中出| 国产精品98久久久久久宅男小说| 99久久久亚洲精品蜜臀av| 91麻豆av在线| 亚洲av电影在线进入| 久久久久久免费高清国产稀缺| 国产成人啪精品午夜网站|