首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 609 毫秒
1.
Daily sea level variability in the Adriatic Sea is studied from different data sets using Empirical Orthogonal Functions, in connection with atmospheric pressure and wind stress. The first mode explains 56–69% of total variance and consists of uniform sea level variability all over the basin, correlated with atmospheric pressure through the inverse barometer effect. The second mode explains 13–16% of variance and accounts for an along-basin sea level gradient, which is correlated with the meridional wind stress component. The first two Principal Components are used as proxies to pressure- and wind-induced components of storm surges in the northern Adriatic. The analysis of the frequency of the most intense events in the 1957–2005 period shows that the wind contribution to storm surges has decreased, while no significant trends are found in the contribution of atmospheric pressure.  相似文献   

2.
Two mathematical models of surface gravity wave refraction were applied to a coastal area of the northern Ligurian Sea. A series of refraction patterns was computed by means of wave data recorded at a moored platform operating 1.5 miles off the harbour of Genoa. The numerical calculations provided an effect of wave energy concentration near Chiavari and Lavagna where coastal structures and beaches are often damaged during storm surges.  相似文献   

3.
《Continental Shelf Research》2007,27(3-4):400-416
During February and June of 2003, selected physical and biological sediment properties were measured at nine mud-bottom sites (11–22 m) along the western margin of the Adriatic Sea. Seabed properties were compared with shipboard measurements of sediment erodibility made at the same sites to gain insight into the physical and biological controls on sediment erodibility operating in the western Adriatic. In addition, spatial patterns of erodibility were compared with long-term records of sediment accumulation in the region to determine if patterns of erodibility were responsible for the discrepancy between sediment sources and sinks on the shelf. Results indicate that sediment erodibility along the western Adriatic Sea varied both in time and space. In the wintertime, sediment in the vicinity of the Po delta (in the northern region) was less easily eroded than to the south. In the summertime, the pattern reversed. The physical characteristics of the seabed including porosity and grain-size were important factors controlling erodibility during the winter, although in a manner opposite to expectations (i.e., higher porosity sediment was less erodible). No relationship between the physical characteristics of the seabed and erodibility was observed in the summer, when it was likely that erodibility was influenced by benthic organisms, especially microphytobenthos. Results of this study suggest long-term sediment accumulation patterns along the western Adriatic Sea are unlikely to be controlled by patterns of erodibility. However, due to several complicating factors, especially along-margin variability in wave energy, the impact of sediment erodibility on long-term patterns of accumulation remains unclear. Further studies that resolve both finer and larger-scale spatial and temporal variability in sediment erodibility and that incorporate year-to-year forcing of the coastal ocean are needed to more accurately resolve the relationship between erodibility and sediment accumulation in the western Adriatic, as well as other coastal environments.  相似文献   

4.
Spatially averaged temperature and salinity profiles from individual cruises between 1990 and 2009 were analysed to outline the temporal evolution of water mass properties in the deep convection site in the South Adriatic Pit (SAP). The long-term variability in thermohaline conditions has been explained and related to a close feedback mechanism between the Ionian and the Adriatic Sea. Prominent influences of the Eastern Mediterranean Transient are manifested in changes in the vertical temperature and salinity patterns in the South Adriatic, and the whole studied period was divided into three stages according to the main thermohaline characteristics: 1990-1995, 1995-2004 and the last period from 2005 onwards.Particular attention was given to data collected during 2006-2009, which permitted us to situate the actual thermohaline properties in the context of the decadal variability. This last period was characterised by a very low production of dense water in the northern basin during 2007, while from winter 2008 high production of North Adriatic Deep Water (NAdDW) and Adriatic Deep Water (AdDW) in the northern and southern basins, respectively, was observed. Finally, we used the Optimum Multiparameter Analysis (OMP) to identify the percentages of the different water masses contained in the SAP, and this highlighted some differences between two recent periods studied (2007 and 2008) and the production of dense waters.  相似文献   

5.
This paper documents the occurrence of strong fortnightly oscillations in the Adriatic Sea frequently observed in the current measurements. For that purpose, we analyzed half-decadal sea level series and long-term currents collected within different parts and layers of the Adriatic Sea. Harmonic analysis and band-pass filter with cutoff periods at 10 and 20 days have been applied to the series. The sea level Mf tide is found to be the only significant tidal constituent over periods between the diurnal and semi-annual ones, having no significant phase variations in the Adriatic Sea. Moreover, the currents on periods between 10 and 20 days could only partially be explained by atmospheric forcing. The fact that they are in-phase with the Mf tide gives a hint on a possible mechanism responsible for the generation of these oscillations. Additionally, these currents are usually strongly amplified during a weakly stratified season, indicating a baroclinic mechanism responsible for the multiplication of the fortnightly currents.  相似文献   

6.
Modeling studies of future changes in coastal hydrodynamics, in terms of storm surges and wave climate, need appropriate wind and atmospheric forcings, a necessary requirement for the realistic reproduction of the statistics and the resolution of small scale features. This work compares meteorological results from different climate models in the Mediterranean area, with a focus on the Adriatic Sea, in order to assess their capability to reproduce coastal meteorological features and their possibility to be used as forcings for hydrodynamic simulations. Five meteorological datasets are considered. They are obtained from two regional climate models, implemented with different spatial resolutions and setups and are downscaled from two different global climate models. Wind and atmospheric pressure fields are compared with measurements at four stations along the Italian Adriatic coast. The analysis is carried out both on simulations of the control period 1960–1990 and on the A1B Intergovernmental Panel for Climate Change scenario projections (2070–2100), highlighting the ability of each model in reproducing the statistical coastal meteorological behavior and possible changes. The importance of simulated global- and regional-scale meteorological processes, in terms of correct spatial resolution of the phenomena, is also discussed. Within the Adriatic Sea, the meteorological climate is influenced by the local orography that controls the strengthening of north-eastern katabatic winds like Bora. Results show indeed that the increase in spatial resolution provides a more realistic wind forcing for the hydrodynamic simulations. Moreover, the chosen setup and the global climate models that drive the regional downscalings appear to play an important role in reproducing correct atmospheric pressure fields. The comparison between scenario and control simulations shows a small increase in the mean atmospheric pressure values, while a decrease in mean wind speed and in extreme wind events is observed, particularly for the datasets with higher spatial resolution. Finally, results suggest that an ensemble of downscaled climate models is likely to provide the most suitable climatic forcings (wind and atmospheric pressure fields) for coastal hydrodynamic modeling.  相似文献   

7.
Tide gauge (TG) data along the northern Mediterranean and Black Sea coasts are compared to the sea-surface height (SSH) anomaly obtained from ocean altimetry (TOPEX/Poseidon and ERS-1/2) for a period of nine years (1993–2001). The TG measures the SSH relative to the ground whereas the altimetry does so with respect to the geocentric reference frame; therefore their difference would be in principle a vertical ground motion of the TG sites, though there are different error sources for this estimate as is discussed in the paper. In this study we estimate such vertical ground motion, for each TG site, from the slope of the SSH time series of the (non-seasonal) difference between the TG record and the altimetry measurement at a point closest to the TG. Where possible, these estimates are further compared with those derived from nearby continuous Global Positioning System (GPS) data series. These results on vertical ground motion along the Mediterranean and Black Sea coasts provide useful source data for studying, contrasting, and constraining tectonic models of the region. For example, in the eastern coast of the Adriatic Sea and in the western coast of Greece, a general subsidence is observed which may be related to the Adriatic lithosphere subducting beneath the Eurasian plate along the Dinarides fault.  相似文献   

8.
The flow of groundwater beneath barrier islands has been cited as a possible pathway for salt water and chemical exchange between a protected embayment and the open sea. Evidence is presented that identifies an exchange of groundwater through a highly permeable paleoinlet along the barrier beach of Cavallino, which separates the northern Venice Lagoon from the Adriatic Sea. We utilized both point measurements of submarine groundwater discharge (SGD) and a geophysical investigation of the subsurface resistivity to analyze the movement of saline groundwater. Discharge of groundwater and associated nutrients, was higher at the site of a former inlet than at a similar site along the barrier and modulated by the difference in tidal water level between the lagoon and Adriatic Sea. If the measured conditions are typical, storm surge barriers could potentially result in a saline groundwater flow of up to 1.5×106 m3 d−1 into the lagoon.  相似文献   

9.
《Continental Shelf Research》2007,27(3-4):296-308
This paper focuses on the delivery of water and sediment to the northern Adriatic to better understand the short-term evolution of continental margin sedimentation under natural and human impact. For that reason, the Po and six Apennine rivers (Metauro, Musone, Potenza, Tronto, Chienti and Pescara) are investigated. The climate-driven hydrological model HydroTrend is used to simulate discharge and sediment loads where observational data are limited. The northern Apennine hinterland has a significant impact on the sediment flux leaving the Po River, contributing 56% of the sediment it delivers to the Adriatic Sea. The Po River experienced a strong decrease in its sediment load (17.2–6.4 Mt/yr) across 1933–1987, in contrast to a small increase in its water discharge. The rivers draining the southern Apennine hinterland contribute more than 50% of the sediment load entering the Adriatic Sea, and this is in spite of human modification of their discharge through numerous small reservoirs that invariably reduce a river's sediment load. As a result, hyperpycnal flows, which historically carried 20–40% of the sediment flux from these Apennine rivers, become rare. Sediment load reduction is also reflected by retreat of the Apennine coastline. Based on the ART model (used in HydroTrend), the total sediment load to the northern Adriatic is 43 MT/yr where the northern Alpine rivers contribute 8 MT/yr, the Po River 13 MT/yr and the Apennine rivers contribute 22 MT/yr.  相似文献   

10.
The impact on a large-scale sea level pressure field to the regional mean sea level changes of the German Bight is analysed. A multiple linear regression together with an empirical orthogonal function analysis is used to describe the relationship between the sea level pressure and the regional mean sea level considering the time period 1924–2001. Both, the part of the variability and of the long-term trend that can be associated with changes in the sea level pressure, are investigated. Considering the whole time period, this regression explains 58?% of the variance and 33?% of the long-term trend of the regional mean sea level. The index of agreement between the regression result and the observed time series is 0.82. As a proxy for large-scale mean sea level changes, the mean sea level of the North East Atlantic is subsequently introduced as an additional predictor. This further improves the results. For that case, the regression explains 74?% of the variance and 87?% of the linear trend. The index of agreement rises to 0.92. These results suggest that the sea level pressure mainly accounts for the inter-annual variability and parts of the long-term trend of regional mean sea level in the German Bight while large-scale sea level changes in the North East Atlantic account for another considerable fraction of the observed long-term trend. Sea level pressure effects and the mean sea level of the North East Atlantic provide thus significant contributions to regional sea level rise and variability. When future developments are considered, scenarios for their future long-term trends thus need to be comprised in order to provide reliable estimates of potential future long-term changes of mean sea level in the German Bight.  相似文献   

11.
This work describes a study of GPS heights, gravity and hydrological time series collected by stations located in northeastern Italy. During the last 12 years, changes in the long-term behaviors of the GPS heights and gravity time series are observed. In particular, starting in 2004–2005, a height increase is observed over the whole area. The temporal and spatial variability of these parameters has been studied as well as those of key hydrological variables, namely precipitation, hydrological balance and water table by using the Empirical Orthogonal Functions (EOF) analysis. The coupled variability between the GPS heights and the hydrological balance and precipitation data has been investigated by means of the Singular Value Decomposition (SVD) approach. Significant common patterns in the spatial and temporal variability of these parameters have been recognized. In particular, hydrology-induced variations are clearly observable starting in 2002–2003 in the southern part of the Po Plain for the longest time series, and from 2004–2005 over the whole area. These findings, obtained by means of purely mathematical approaches, are supported by sound physical interpretation suggesting that the climate-related fluctuations in the regional/local hydrological regime are one of the main contributors to the observed variations. A regional scale signal has been identified in the GPS station heights; it is characterized by the opposite behavior of the southern and northern stations in response to the hydrological forcing. At Medicina, in the southern Po Plain, the EOF analysis has shown a marked common signal between the GPS heights and the Superconducting Gravimeter (SG) data both over the long and the short period.  相似文献   

12.
A shallow water hydrostatic 2D hydrodynamic numerical model, based on the boundary conforming coordinate system, was used to simulate aspects of both general and small scale oceanic features occurring in the composite system constituted by the Adriatic Sea and the Lagoon of Venice (Italy), under the influence of tide and realistic atmospheric forcing. Due to a specific technique for the treatment of movable lateral boundaries, the model is able to simulate efficiently dry up and flooding processes within the lagoon. Firstly, a model calibration was performed by comparing the results of the model, forced using tides and ECMWF atmospheric pressure and wind fields, with observations collected for a set of 33 mareographic stations uniformly distributed in the Adriatic Sea and in the Lagoon of Venice. A second numerical experiment was then carried out by considering only the tidal forcing. Through a comparison between the results obtained in the two experiments it was possible to assess the reliability of the estimated parameter through the composite forcing. Model results were then verified by comparing simulated amplitude and phase of each tidal constituent as well as tidal velocities simulated at the inlets of the lagoon and in the Northern Adriatic Sea with the corresponding observed values. The model accurately reproduces the observed harmonics: mean amplitude differences rarely exceed 1 cm, while phase errors are commonly confined below 15°. Semidiurnal and diurnal currents were correctly reproduced in the northern basin and a good agreement was obtained with measurements carried out at the lagoon inlets. On this basis, the outcomes of the hydrodynamic model were analyzed in order to investigate: (i) small-scale coastal circulation features observed at the interface between the adjoining basins, which consist often of vortical dipoles connected with the tidal flow of Adriatic water entering and leaving the Lagoon of Venice and with along-shore current fields connected with specific wind patterns; (ii) residual oscillations, which are often connected to meteorological forcing over the basin. In particular, it emerges that small-scale vortical features generated near the lagoon inlet can be efficiently transported toward the open sea, thus contributing to the water exchange between the two marine regions, and a realistic representation of observed residual oscillations in the area would require a very detailed knowledge of atmospheric as well as remote oceanic forcing.  相似文献   

13.
The coastal zones are facing the prospect of changing storm surge statistics due to anthropogenic climate change. In the present study, we examine these prospects for the North Sea based on numerical modelling. The main tool is the barotropic tide-surge model TRIMGEO (Tidal Residual and Intertidal Mudflat Model) to derive storm surge climate and extremes from atmospheric conditions. The analysis is carried out by using an ensemble of four 30-year atmospheric regional simulations under present-day and possible future-enhanced greenhouse gas conditions. The atmospheric regional simulations were prepared within the EU project PRUDENCE (Prediction of Regional scenarios and Uncertainties for Defining EuropeaN Climate change risks and Effects). The research strategy of PRUDENCE is to compare simulations of different regional models driven by the same global control and climate change simulations. These global conditions, representative for 1961–1990 and 2071–2100 were prepared by the Hadley Center based on the IPCC A2 SRES scenario. The results suggest that under future climatic conditions, storm surge extremes may increase along the North Sea coast towards the end of this century. Based on a comparison between the results of the different ensemble members as well as on the variability estimated from a high-resolution storm surge reconstruction of the recent decades it is found that this increase is significantly different from zero at the 95% confidence level for most of the North Sea coast. An exception represents the East coast of the UK which is not affected by this increase of storm surge extremes.  相似文献   

14.
《Continental Shelf Research》2007,27(3-4):542-559
Flood deposition and storm reworking of sediments on the inner shelf can change the mixture of grain sizes on the seabed and thus its porosity, bulk density, bulk compressional velocity and reflectivity. Whether these changes are significant enough to be detectable by repeat sub-bottom sonar surveys, however, is uncertain. Here the question is addressed through numerical modeling. Episodic flooding of a large versus small river over the course of a century are modeled with HYDROTREND using the drainage basin characteristics of the Po and Pescara Rivers (respectively). A similarly long stochastic record of storms offshore of both rivers is simulated from the statistics of a long-term mooring recording of waves in the western Adriatic Sea. These time series are then input to the stratigraphic model SEDFLUX2D, which simulates flood deposition and storm reworking on the inner shelf beyond the river mouths. Finally, annual changes in seabed reflectivity across these shelf regions are computed from bulk densities output by SEDFLUX2D and compressional sound speeds computed from mean seafloor grain size using the analytical model of Buckingham [1997. Theory of acoustic attenuation, dispersion, and pulse propagation in unconsolidated granular materials including marine sediments. Journal of the Acoustical Society of America 102, 2579–2596; 1998. Theory of compressional and shear waves in fluidlike marine sediments. Journal of the Acoustical Society of America 103, 288–299; 2000. Wave propagation, stress relaxation, and grain-tograin shearing in saturated, unconsolidated marine sediments. Journal of the Acoustical Society of America 108, 2796–2815]. The modeling predicts reflectivities that change from <12 dB for sands on the innermost shelf to >9 dB for muds farther offshore, values that agree with reflectivity measurements for these sediment types. On local scales of ∼100 m, however, maximum changes in reflectivity are <0.5 dB. So are most annual changes in reflectivity over all water depths modeled (i.e., 0–35 m). Given that signal differences need to be ⩾2–3 dB to be resolved, the results suggest that grain-size induced changes in reflectivity caused by floods and storms will rarely be detectable by most current sub-bottom sonars.  相似文献   

15.
This study analyzes how the stochastically generated rainfall time series accounting for the inter-annual variability of rainfall statistics can improve the prediction of watershed response variables such as peak flow and runoff depth. The modified Bartlett–Lewis rectangular pulse (MBLRP) rainfall generation model was improved such that it can account for the inter-annual variability of the observed rainfall statistics. Then, the synthetic rainfall time series was generated using the MBLRP model, which was used as input rainfall data for SCS hydrologic models to produce runoff depth and peak flow in a virtual watershed. These values were compared to the ones derived from the synthetic rainfall time series that is generated from the traditional MBLRP rainfall modeling. The result of the comparison indicates that the rainfall time series reflecting the inter-annual variability of rainfall statistics reduces the biasness residing in the predicted peak flow values derived from the synthetic rainfall time series generated using the traditional MBLRP approach by 26–47 %. In addition, it was observed that the overall variability of the peak flow and run off depth distribution was better represented when the inter-annual variability of rainfall statistics are considered.  相似文献   

16.
This paper presents the results of the statistical analysis of a set of physico-chemical and biological water quality parameters, monthly collected from 2000 to 2007 in the Genoa Harbour area (Ligurian Sea). We applied multivariate methods, such as principal component analysis (PCA) and dynamic factor analysis (DFA) for investigating the spatial and temporal variability and for providing important background information on pollution problems in the region. PCA evidenced the role of the sewage and river discharges and of the exchanges with the open sea in determining the harbour water quality. DFA was used to estimate underlying common trends in the time series. The DFA results partly show a general improvement of water quality over the 8-years period. However, in other areas, we found inter-annual variations but no significant multi-annual trend. Furthermore, we included meteorological variables in our statistical analyses because of their potential influence on the water quality parameters. These natural forcings explain part of the variability in water quality parameters that are superimposed on the dominating anthropogenic pollution factors.  相似文献   

17.
—In the last decade temporary teleseismic transects have become a powerful tool for investigating the crustal and upper mantle structure. In order to gain a clearer picture of the lithosphere-asthenosphere structure in peninsular Italy, between 1994 and 1996, we have deployed three teleseismic transects in northern, central, and southern Apennines, in the framework of the project GeoModAp (European Community contract EV5V-CT94–0464). Some hundreds of teleseisms were recorded at each deployment which lasted between 3 and 4 months. Although many analyses are still in progress, the availability of this high quality data allowed us to refine tomographic images of the lithosphere-asthenosphere structure with an improved resolution in the northern and central Apennines, and to study the deformation of the upper mantle looking at seismic anisotropy through shear-wave splitting analysis. Also, a study of the depth and geometry of the Moho through the receiver function technique is in progress. Tomographic results from the northernmost 1994 and the central 1995 teleseismic experiments confirm that a high-velocity anomaly (HVA) does exist in the upper 200–250 km and is confined to the northern Apenninic arc. This HVA, already interpreted as a fragment of subducted lithosphere is better defined by the new temporary data, compared to previous works, based only on data from permanent stations. No clear high-velocity anomalies are detected in the upper 250 km below the central Apennines, suggesting either a slab window due to a detachment below southern peninsular Italy, or a thinner, perhaps continental slab of Adriatic lithosphere not detectable by standard tomography. We found clear evidence of seismic anisotropy in the uppermost mantle, related to the main tectonic processes which affected the studied regions, either NE–SW compressional deformation of the lithosphere beneath the mountain belt, or arc-parallel asthenospheric flow (both giving NW–SE fast polarization direction), and successive extensional deformation (～E–W trending) in the back-arc basin of northern Tyrrhenian and Tuscany. Preliminary results of receiver function studies in the northern Apennines show that the Moho depth is well defined in the Tyrrhenian and Adriatic regions while its geometry underneath the mountain belt is not yet well constrained, due to the observed high complexity.  相似文献   

18.
The Adriatic Sea general circulation model coupled to a third generation wave model SWAN and a sediment transport model was implemented in the Adriatic Sea to study the dynamics of the sediment transport and resuspension in the northern Adriatic Sea (NAS) during the Bora event in January 2001. The bottom boundary layer (BBL) was resolved by the coupled model with high vertical resolution, and the mechanism of the wave–current interaction in the BBL was also represented in the model. The study found that, during the Bora event of 13–17 January 2001, large waves with significant wave height 2 m and period of 5 s were generated by strong winds in the northwestern shelf of the Adriatic where the direction of wave propagation was orthogonal to the current. The combined motion of the wave and current in the BBL increased the bottom stress over the western Adriatic shelf, resulting in stronger sediment resuspension there. Combining stronger bottom resuspension and strong upward vertical flux of resuspended sediments due to turbulent mixing, the model predicted that sediment concentration near the Po River was much higher than that predicted by the model run without wave forcing. The study also shows that wave–current interaction in the BBL reduced the western Adriatic Coastal Currents (WACCs) in the shallower north. It is concluded that wave forcing significantly changed the sediment distributions and increased the total horizontal fluxes over the western shelf. These results signified wave effect on sediment flux and distribution in the NAS, and suggested that waves cannot be neglected in the study of dynamics of sediment transport and resuspension in the shallow coastal seas. By including the tidal forcing in the coupled model, we also examined the effect of tides on the sediment transport dynamics in the NAS.  相似文献   

19.
We argue that seismotectonic activity in the Central Mediterranean area and the Aegean–Balkan zone is driven by the NNE-ward motion of Africa and westward motion of Anatolia with respect to Eurasia. These boundary conditions can plausibly and coherently account for E–W shortening and roughly S–N extension in the Aegean domain, thrusting and uplift at the boundary between the Aegean–Balkan system and the Adriatic/Ionian domain (Hellenic trench, Cephalonia fault, Epirus, Albanides and Southern Dinarides), the kinematics of the Adria plate (a large block encompassing the Adriatic continental domain, the northern Ionian zone and Hyblean-Adventure block) and consequently, the complex pattern of deformation recognized at its boundaries. Furthermore, the fact that in our scheme Adria moves almost in connection with Africa is consistent with the lack of an active decoupling zone between Adria and Africa, an evidence that can hardly be reconciled with the kinematics so far proposed for these two plates. The reasons why we adopt an Africa–Eurasia relative motion different from that implied by the popular NUVEL-1 global solution are discussed in detail. Finally, we make some considerations about the possible implications of the presently available geodetic data on the long-term plate kinematics.  相似文献   

20.
With global warming and sea level rise, many coastal systems will experience increased levels of inundation and storm flooding, especially along sandy lowland coastal areas, such as the Northern Adriatic coast (Italy). Understanding how extreme events may directly affect groundwater hydrology in shallow unconfined coastal aquifers is important to assess coastal vulnerability and quantify freshwater resources. This study investigates shallow coastal aquifer response to storm events. The transitory and permanent effects of storm waves are evaluated through the real time monitoring of groundwater and soil parameters, in order to characterize both the saturated and unsaturated portions of the coastal aquifer of Ravenna and Ferrara (southern Po Delta, Italy). Results highlight a general increase in hydraulic head and soil moisture, along with a decrease in groundwater salinity and pore water salinity due to rainfall infiltration during the 2 days storm event. The only exceptions are represented by the observation wells in proximity to the coastline (within 100 m), which recorded a temporary increase in soil and water salinity caused by the exceptional high waves, which persist on top of the dune crest during the storm event. This generates a saline plume that infiltrates through the vadose zone down to the saturated portion of the aquifer causing a temporary disappearance of the freshwater lens generally present, although limited in size, below the coastal dunes. Despite the high hydraulic conductivity, the aquifer system does not quickly recover the pre‐storm equilibrium and the storm effects are evident in groundwater and soil parameters after 10 days past the storm overwash recess.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号