首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 31 毫秒
1.
An all-sky CCD imager designed to measure wave structure of the OH, O2b(0,1) and OI557.7 nm airglow emission layers in the mesosphere and lower thermosphere (MLT) region has been operated near the equatorial region at São João do Cariri (Cariri), Brazil, (7.5°S, 36.5°W). A large number of gravity wave was observed from September 2000 to September 2002 and among them 64 wave events were identified as mesospheric bores. The bore front shows a horizontal extension greater than 1000 km, and observed in the airglow layers as a complementary brilliance between the three emissions. At the first time mesospheric bore events were observed and analyzed in the equatorial region. Their predominant characteristics, occurrence, local time dependency, morphology and propagation direction will be presented and discussed.  相似文献   

2.
Lidar measurements of upper atmospheric sodium were made at Alcântara (2°20′ S, 44°30′ W), Brazil, on 8 days during the month of November 1996. These are the first sodium measurements to be reported from an almost equatorial location. The average sodium distribution at Alcântara was found to be very similar to the annual mean for São José dos Campos (23° S), but about 1 km higher than the November mean value of 92 km for the latter site. Only one, weak, sporadic sodium layer was observed during more than 26 h of observations, an occurrence frequency much less than that seen at low latitude sites in the northern hemisphere. A survey of the observations of sporadic layers at various locations suggests that both magnetic dip and declination angles might influence their rate of occurrence.  相似文献   

3.
Least-squares collocation technique was used to process regional gravity data of the SE South American lithospheric plate in order to map intermediate (10–2000 km) wavelength geoid anomalies. The area between 35–10° S and 60–25° W includes the Paraná CFB Province, the Southern São Francisco Craton and its marginal fold/thrust belts, the Brazilian continental margin and oceanic basins. The main features in the geoid anomaly map are: (a) Paraná CFB Province is characterized by a 1000 km long and 500 km wide, NE-trending, 9 m-amplitude negative anomaly which correlates with the distribution of sediments and basalts within the Paraná basin. (b) A circular (600–800 km in diameter) positive, 8 m-amplitude geoid anomaly is located in the southern S. Francisco craton and extends into the northeastern border of the Paraná CFB Province. This anomaly partially correlates with Alto Paranaíba Igneous Province (APIP), where alkalic volcanism and tholeiitic dikes of ages younger than 80 Ma are found and where a low-velocity zone in the mantle has been mapped using seismic tomography. This positive geoid anomaly extends towards the continental margin at latitude 21° S and joins a linear sequence of short wavelength positive geoid anomalies associated with Vitoria–Trindade seamounts. (c) A NE-trending, 1000 km long and 800 km wide, 4 m-amplitude, positive geoid anomaly, which is located along the southeastern coast of Brazil, from latitude 24 to 35° S. The northern part of this anomaly correlates with the Ponta Grossa Arch and Florianopolis dyke swarm provinces. The age of this intrusive volcanism is 130–120 Ma. (d) A circular positive anomaly with 9 m of amplitude, located over the Rio Grande and Uruguay shields and offshore Pelotas basin. Few alkaline intrusives with ages between 65 and 80 Ma are found in the region and apatite fission track ages in basement rocks indicates cooling at around 30 Ma. A semi-quantitative analysis of the observed geoid anomalies using isostatic considerations suggests that the mechanism which generated Paraná CFB Province did not change, in a significant manner, the lithospheric thermal structure, since the same geoid pattern observed within this province continues northward over the Neoproterozoic fold/thrust belts systems separating the São Francisco and Amazon cratons. Therefore, this observation favours Anderson’s idea of rapid basaltic outpouring through a pull-apart mechanism along a major suture zone. A thermal component may still be present in the Southern São Francisco Craton and in the Rio Grande Shield and contiguous continental margins, sites of Tertiary thermal and magmatic reactivations.  相似文献   

4.
Continuous MF and meteor radar observations allow detailed studies of winds in the mesosphere and lower thermosphere (MLT) as well as temperatures around the mesopause. This height region is characterized by a strong variability in winter due to enhanced planetary wave activity and related stratospheric warming events, which are distinct coupling processes between lower, middle and upper atmosphere. Here the variability of mesospheric winds and temperatures is discussed in relation with major and minor stratospheric warmings as observed during winter 2005/06 in comparison with results during winter 1998/99.Our studies are based on MF radar wind measurements at Andenes (69°N, 16°E), Poker Flat (65°N, 147°W) and Juliusruh (55°N, 13°E) as well as on meteor radar observations of winds and temperatures at Resolute Bay (75°N, 95°W), Andenes (69°N, 16°E) and Kühlungsborn (54°N, 12°E). Additionally, energy dissipation rates have been estimated from spectral width measurements using a 3 MHz Doppler radar near Andenes. Particular attention is directed to the changes of winds, turbulence and the gravity wave activity in the mesosphere in relation to the planetary wave activity in the stratosphere.Observations indicate an enhancement of planetary wave 1 activity in the mesosphere at high latitudes during major stratospheric warmings. Daily mean temperatures derived from meteor decay times indicate that strong warming events are connected with a cooling of the 90 km region by about 10–20 K. The onset of these cooling processes and the reversals of the mesospheric circulation to easterly winds occur some days before the changes of the zonal circulation in the stratosphere start indicating a downward propagation of the circulation disturbances from the MLT region to the stratosphere and troposphere during the stratospheric warming events. The short-term reversal of the mesospheric winds is followed by a period of strong westerly winds connected with enhanced turbulence rates and an increase of gravity wave activity in the altitude range 70–85 km.  相似文献   

5.
Whatever its origin is, a floating particle at the sea surface is advected by ocean currents. Surface currents could be derived from in situ observations or combined with satellite data. For a better resolution in time and space, we use satellite-derived sea-surface height and wind stress fields with a 1/3° grid from 1993 to 2001 to determine the surface circulation of the South Pacific Ocean. Surface currents are then used to compute the Lagrangian trajectories of floating debris. Results show an accumulation of the debris in the eastern-centre region of the South Pacific subtropical gyre ([120°W; 80°W]–[20°S; 40°S]), resulting from a three-step process: in the first two years, mostly forced by Ekman drift, the debris drift towards the tropical convergence zone (30°S). Then they are advected eastward mostly forced by geostrophic currents. They finally reach the eastern-centre region of the South Pacific subtropical gyre from where they could not escape.  相似文献   

6.
A study of emplacement temperatures was carried out for the largest of the 22 November 1994 nuée ardente deposits at Merapi Volcano, based mainly on the response of plastic and woody materials subjected to the hot pyroclastic current and the deposits, and to some extent on eyewitness observations. The study emphasizes the Turgo–Kaliurang area in the distal part of the area affected by the nuée ardente, where nearly 100 casualties occurred. The term nuée ardente as used here includes channeled block-and-ash flows, and associated ash-clouds of surge and fallout origins. The emplacement temperature of the 8 m thick channeled block-and-ash deposit was relatively high, 550°C, based mainly on eyewitness reports of visual thermal radiance. Emplacement temperatures for ash-cloud deposits a few cm thick were deduced from polymer objects collected at Turgo and Kaliurang. Most polymers do not display a sharp melting range, but polyethylene terephthalate used in water bottles melts between 245 and 265°C, and parts of the bottles that had been deformed during fabrication molding turn a milky color at 200°C. The experimental evidence suggests that deposits in the Turgo area briefly achieved a maximum temperature near 300°C, whereas those near Kaliurang were <200°C. Maximum ash deposit temperatures occurred in fallout with a local source in the channeled block-and-ash flow of the Boyong river valley; the surge deposit was cooler (180°C) due to entrainment of cool air and soils, and tree singe-zone temperatures were around 100°C.  相似文献   

7.
An in-situ experimental program in a jointed rock mass was designed and implemented to investigate the rock joint effects on stress wave propagation. Accelerometers were placed on rock surface along three lines at 0, 45 and 90° to the predominant rock joint strike direction. Eight blasts were detonated in a charge hole drilled in the rock mass. The equivalent TNT blast weight ranged from 2.5 to 50.0 kg and the loading density (charge weight divided by the charge chamber volume) varied from 1 to 20.0 kg/m3. A fully coupled detonation with a charge weight of 25.0 kg was also exploded to study the coupling effect. The recorded data are presented and analyzed in this paper. The effects of rock joints on characteristics of stress wave propagation such as peak value attenuation, spectrum, and spatial variations are discussed.  相似文献   

8.
Santo Antão, the northernmost island of the Cape Verde Archipelago, consists entirely of silica-undersaturated volcanic products and minor intrusions. 40Ar–39Ar incremental heating experiments have been carried out on 24 samples that cover the entire exposed chronological sequence. The oldest lavas (7.57±0.56 Ma), representing an older volcanic basement, are exposed about 620 m above mean sea level. After an interval of quiescence of up to 4.3 Ma the volcanic activity resumed and continued at low eruption rates. The older basement is unconformably overlain by a ca. 810-m-thick lava sequence that spans an age range from 2.93±0.03 to 1.18±0.01 Ma. This sequence is cut by many dykes and sills. Simultaneous volcanic activity occurred in the northeastern, central and eastern part of the island. A phonolitic pumice deposit that forms a noteworthy feature over most of the island has an estimated age of 0.20 Ma. This predates volcanic activity that formed the highest point of the island (Tope de Coroa) which has an age of 0.17±0.02 Ma. The most recent eruption on the island formed nephelinitic lavas in the Porto Novo region at 0.09±0.03 Ma. The oldest volcanism exposed on Santo Antão, which took place about 7.6 Ma ago, was simultaneous with waning activity on Maio at the eastern end of the Cape Verde Archipelago.  相似文献   

9.
The sea-surface height anomalies derived from Simple Ocean Data Assimilation (SODA) during 1958–2001, Topex/Poseidon satellite during 1993–2001 and the SODA heat content anomalies (125 m depth) during 1958–2001 are filtered into annual and biennial Rossby wave components using a two-dimensional Finite Impulse Response filter. The filtered Rossby wave components (both annual and biennial) in the southern Pacific and Indian Oceans have considerable strength and variability. The propagation of annual and biennial Rossby waves in the Indonesian through-flow region [12.5°S–7.5°S] of the Indian Ocean is in phase with the southern Pacific Ocean waves. So it is speculated that the Pacific Ocean influences the Indian Ocean, especially through the region 17.5°S to 7.5°S and thus the southern Pacific Rossby waves could be an unexplored contributor to the Indian Ocean Rossby waves. We also carried out Fast Fourier Transform (FFT) and wavelet analysis on the tide gauge sea-level data along the Australian coast to support our claim. Filtered annual and biennial components of SODA heat content anomalies (125 m depth) also support these findings.  相似文献   

10.
Seventy sites of sills, flows and dikes from Northeastern Paraná Magmatic Province (PMP), were submitted to paleomagnetic, chemical and radiometric analyses. The rocks are high in TiO2 content, and similar in composition to the rocks from the northern region of PMP. The sills intrude mainly Paleozoic sediments, and can be subdivided into two domains; the northern being characterized by sills showing reversed polarities, and the southern essentially by sills of normal polarities. 40Ar/39Ar dating of three distinct sills gave plateau ages (129.9 ± 0.1, 130.3 ± 0.1 and 131.9 ± 0.4 Ma) that are similar to surface-outcropping flows of the Northern Paraná Basin, and the Ponta Grossa dikes. The new paleomagnetic data combined with existing data from the northern PMP allowed the calculation of a paleomagnetic pole at 71.4° E and 83.0° S (N = 92; α95=2.4°; k = 39). This pole is in good agreement with poles for central and southern PMP, which are slightly older than the northern PMP, as well as for the contemporaneous Central Alkaline Province (Paraguay) on the western side of PMP. In contrast, the coeval pole for the Ponta Grossa dikes (eastern border of PMP), however, is slightly displaced from that group of poles, suggesting that dikes in that area may have undergone some tectonic tilting.  相似文献   

11.
We assessed leaf breakdown of five native riparian species from Brazilian Cerrado (Myrcia guyanensis, Ocotea sp., Miconia chartacea, Protium brasiliense, and Protium heptaphyllum), incubated in single and mixed species packs in two headwater streams with different physico-chemical properties in the Espinhaço Mountain range (Southeastern Brazil). Leaves were placed in plastic litter bags (15 cm×20 cm, 10 mm mesh size) and the experiments were carried out during the dry seasons of 2003 and 2004. Leaf nitrogen and phosphorus contents were similar in all species, but polyphenolic contents were different (P<0.001). M. guyanensis showed higher polyphenolics content (8.48% g−1 dry mass) and leaf toughness. Individually, higher breakdown rates were found in M. guyanensis at Indaiá stream (k=0.0063±0.0005 d−1) and in Ocotea sp. at Garcia stream (k=0.0088±0.0006 d−1). However, P. brasiliense and P. heptaphyllum showed lower breakdown rates at Indaiá and Garcia streams (Indaiá: k=0.0020±0.0002 and 0.0019±0.0001 d−1; Garcia: k=0.0042±0.0001 and 0.0040±0.0002 d−1). Single and mixed breakdown processes of each species were not statistically different on both streams. However, all species showed higher breakdown rates at Garcia stream (P<0.01). These results suggest that leaf breakdown is not altered when litter benthic patches are composed by a mixture of species in the same proportions that they occur on riparian leaf falls.  相似文献   

12.
In this paper we try to explain the empirical Na---K---Ca geothermometer introduced by Fournier and Truesdell (1973).At high temperatures (100°C < t < 370°C), the geothermometer can be related to an equilibrium between plagioclase, K-feldspar and water. If this is true, plagioclase must have a nearly constant composition, rich in albite. At low temperatures, the relation between temperature and the Fournier and Truesdell formula can be explained by an increase of the total concentration of dissolved salts with temperature. This increase and the increase of temperature calculated by the Fournier and Truesdell relation can be related to other effects, among which one of the most important is that we deal with a quasi-closed system.

Résumé

Dans cet article on essaye d'expliquer le géothermométre empirique de Fournier et Truesdell (1973) utilisant les concentrations en sodium, potassium et calcium dissous dans les eaux.A haute température (100°C < t < 370°C), le thermométre peut être expliqué par un équilibre entre plagioclase, feldspath potassique et eau, á condition que le plagioclase ait une composition á peu prés constante et trés proche du pôle albite. A basse température, la liaison entre la température et la formule de Fournier et Truesdell peut s'expliquer par une augmentation de la concentration globale en sels dissous, á composition relative constante, avec la température. Mais l'augmentation de la concentration globale et par suite l'augmentation de la température calculée peut être également due á d'autres facteurs dont le principal est le confinement de la solution.  相似文献   

13.
An earthquake of magnitude Ms = 6.0 (CSEM, Strasbourg) occurred at Constantine (Algeria) on 27 October 1985. This seismic event is the strongest felt in the Tellian Atlas since the El Asnam seismic crisis of October 10, 1980. A team from the Centre de Recherche d'Astronomie, d'Astrophysique et de Géophysique (CRAAG, Algeria), utilising 8 portable stations, registered the activity a few days after the main shock. The aftershocks follow a N045° direction, and show the existence of three ruptured segments. Cross sections display a remarkable vertical fault plane and suggest asperities in the rupture process. Surface breaks were found affecting Quaternary deposits. The principal segment is about 3.8 km long showing “enéchelon” cracks with left-lateral displacement while the main direction of the rupture is N055°. Although the vertical motion is small, the northwestern block shows a normal component of the main surface faulting, while the left-lateral displacement is about 10 cm. The strike-slip focal mechanism solution determined from the global seismic network and field observations are in good agreement.  相似文献   

14.
A sodium resonance lidar at 589 nm has been operated in São José dos Campos, Brazil (23°S, 46°W) since 1972 mainly for studies related to the origin, chemistry and dynamics of the mesospheric sodium layer. Beginning in 1993, the improved laser capability has also enabled the processing of the Rayleigh signal from which the temperatures from 35 to 65 km are retrieved on a nightly mean basis. We used these nightly profiles to determine the monthly temperature profiles from 1993 to 2006. The mean temperature characteristics for each year and for the whole period are obtained. Seasonal thermal amplitude is small (6 K peak to peak at 40 and 60 km). Compared with the MSISE-90 model, a large difference is noted, with temperature lower than the model below the stratopause and higher above. Also the seasonal variation has a large difference with better agreement occurring around local winter, but with temperatures higher by 8–10 K at the equinoxes. The semiannual component is dominant over the annual at all altitudes. Linear trends with decreasing temperature of 1.09, 2.29 and 1.42 K/decade are observed at 40, 50 and 60 km, respectively.  相似文献   

15.
Basic properties of the mid-latitude traveling ionospheric disturbances (TIDs) during the maximum phase of a major magnetic storm of 6–8 April 2000 are shown. Total electron content (TEC) variations were studied by using data from GPS receivers located in Russia and Central Asia. The nightglow response to this storm at mesopause and termospheric altitudes was also measured by optical instruments FENIX located at the observatory of the Institute of Solar-Terrestrial Physics (51.9°N,103.0°E), and MORTI located at the observatory of the Institute of Ionosphere (43.2°N, 77.0°E). Observations of the O (557.7 and 630.0 nm) emissions originating from atmospheric layers centered at altitudes of 90 and 250 km were carried out at Irkutsk and of the O2(b1∑g+−X3∑g−) (0-1) emission originating from an atmospheric layer centered at altitude of 94 km was carried out at Almaty. Our radio and optical measurement network observed a storm-induced solitary large-scale wave with duration of 1 h and a wave front width of no less than 5000 km, while it traveled equatorward with a velocity of 200 m/s from 62°N to 38°N geographic latitude. The TEC disturbance, basically displaying an electron content depression in the maximum of the F2 region, reveals a good correlation with growing nightglow emission, the temporal shift between the TEC and emission variation maxima being different for different altitudes. A comparison of the auroral oval parameters with dynamic spectra of TEC variations and optical 630 nm emissions in the frequency range 0.4–4 mHz (250–2500 s periods) showed that as the auroral oval expands into mid-latitudes, also does the region with a developed medium-sale and small-scale TEC structure.  相似文献   

16.
We report here crustal shear-wave anisotropy, ranging from 1% to 10.76% with an average of 2.4% in the aftershock zone of the 2001 Bhuj earthquake, Gujarat, India, from a study of leading shear-wave polarization directions (LPSDs), which vary on average from NNW–SSE to E–W with a delay of 0.07–0.14 s. The delays in the NNW–SSE to NE–SW directions observed at seven stations, near the seismogenic fault, suggest cracks parallel to the direction of the maximum horizontal regional compressional stress prevailing in the region, suggesting a dilatancy-induced anisotropy resulting from approximately stress-aligned parallel vertical micro-cracks. In contrast, the LPSDs at Ramvav, Rapar and Vondh stations, away from the seismogenic fault, are fault parallel, approximately E–W and almost orthogonal to the stress-aligned polarizations inferred elsewhere. The maximum average time delay of 0.14 s is observed at Lodai, where the fast polarization direction is found to be N338°W. This has been observed from anisotropic poro-elastic (APE) modelling and observations that these are 90° flips in shear-wave polarization, resulting from propagation through micro-cracks containing fluids at critically high pore-fluid pressure surrounding the hypocenter of the 2001 mainshock. The presence of high pore-fluid pressure in the seismogenic fault zone could also explain the observed scatter in shear-wave time delays. Further, the coincidence of the N–S trending intrusive bodies (as inferred from tomographic studies in the area) with the N–S direction of regional maximum horizontal compressional stress supports the interpretation of stress-aligned vertical extensive-dilatant anisotropic (EDA) cracks. The depth distribution of the estimated anisotropy (1–10.76%), b-values and stress drop values suggests an increase at 18–30 km depths, which could be attributed to high pore-fluid pressures resulting from a fluid-filled fractured rock matrix or open micro-cracks (characterized by high crack density and high porosity) coinciding with a low velocity zone (at 18–30 km depths) as delineated from tomographic studies in the area.  相似文献   

17.
Paleomagnetic data from 46 sites (674 specimens) of the Westcoast Crystalline Gneiss Complex on the west coast of Vancouver Island using AF and thermal demagnetization methods yields a high blocking temperature WCB component (> 560°C) with a pole at 335°W, 66°N (δp = 4°, δm = 6°) and a lower coercivity WCA component ( 25 mT, < 500°C) with a pole at 52°W, 79°N (δp = 7°, δm = 8°). Further thermal demagnetization data from 24 sites in the Jurassic Island Intrusions also defines two high blocking temperature components. The IIA component pole is at 59°W, 79°N (δp = 7°, δm = 8°) and IIB pole at 130°W, 73°N (δp = 12°, δm = 13°). Combined with previous data from the Karmutsen Basalts and mid-Tertiary units on Vancouver Island and from the adjacent Coast Plutonic Complex, the geotectonic motions are examined for the Vancouver Island segment of the Wrangellian Subterrane of composite Terrane II of the Cordillera. The simplest hypothesis invokes relatively uniform translation for Terrane II from Upper Triassic to Eocene time producing 39° ± 6° of northward motion relative to the North American craton, combined with 40° of clockwise rotation during the Lower Tertiary.  相似文献   

18.
We have constructed a new bathymetric chart of part of the Chile transform system, based mainly on an R/V “Endeavor” survey from 100°W to its intersection with the East Ridge of the Juan Fernandez microplate at 34°30′S, 109°15′W. A generally continuous lineated trend can be followed through the entire region, with the transform valley being relatively narrow and well-defined from 109°W to approximately 104°30′W. The fracture zone then widens to the east, with at least two probable en echelon offsets to the south at 104° and 102°W. Six new strike-slip mechanisms along the Chile Transform and one normal fault mechanism near the northern end of the Chile Rise, inverted together with other plate motion data from the eastern portion of the boundary, produce a new best fit Euler pole for the Nazca-Antarctic plate pair, providing tighter constraints on the relative plate motions.  相似文献   

19.
On July 8, 1976, eruptive activity broke out at la Soufrière de Guadeloupe (F.W.I) after about one year of increasing seismic activity. Seismic activity continued to increase until August 1976, reaching more than 1500 events (a 200-fold increase over the preceding quiet period of a few years) and an energy output of about 1017 ergs in a day. A total of 26 major phreatic eruptions similar to the July 8 outburst took place during an eight-months period. The steam blasts that characterized the eruptions gave rise to particle- and sometimes block-charged plumes that deposited an estimated 106 m3 of solids. The H2O-rich gases emitted during the blasts presumably contained other gases (H2S, SO2, CO2...) that were partly adsorbed on solid particles. All material was erupted at temperatures of the order of 100° to 200°C.The observation of vertical migration of earthquake foci in less than a few hours and over about 6 km depth, and of abnormal variations of the geomagnetic field, indicate a deep energy source for the phreatic eruptions. A small proportion of the gases adsorbed on solid particles had a magmatic origin. However, most of the steam and the tephra seemed to originate from superficial levels of a hydrothermal system. Similar phreatic eruptions have occurred several times in recorded history. In the case of la Soufrière, the origin of the phreatic eruptions is best described by an abnormal energy input (versus steady-state) from a crustal magma chamber. The occurrence of truly magmatic eruptions is presumably inhibited by an extensive hydrothermal system. The abrupt release of more power from the magma chamber could have resulted in an explosive pyroclastic eruption.Substantial improvement of the Guadeloupe volcano observatory has followed the 1975–1977 crisis. Permanent telemetered geophysical networks and regular geochemical observations have provided a five year data base of the volcano behavior in its noneruptive state which can be compared to crisis situations.  相似文献   

20.
Paleomagnetic measurements have been carried out on welded tuffs ranging in age between 58 Ma and 112 Ma from Yamaguchi and Go¯river areas in the central part of Southwest Japan. The new data, together with those of younger igneous rocks published previously, define the change of paleomagnetic field direction during the late Mesozoic/ Cenozoic period for Southwest Japan. The paleomagnetic direction from this area has pointed 56 ± 3° clockwise from the expected field direction estimated from APWP (apparent polar wandering path) of the whole of Eurasia during the period between 100 and 20 Ma. In comparison with the expected one from the eastern margin of Eurasia (Korea, China, Siberia), the Cretaceous field direction of Southwest Japan shows the clockwise deflection by 44–49°. These results establish that while the eastern margin of Eurasia, including Southwest Japan, was rotated more or less with respect to the main part of Eurasia during last 100 Ma, Southwest Japan was rotated clockwise through more than 40° with respect to the eastern margin of Eurasia since 20 Ma. The large amount of rotation for Southwest Japan implies that it is rotated by an opening of the southwestern part of the Japan Sea, which widens northeastward (fan shape opening). The tectonic feature of Southwest Japan and the Japan Sea is analogous to that of Corso-Sardinia and the Ligurian Sea in the Mediterranean, indicating that the fan shape opening is a specific feature of the rifting of the continental sliver at the continental rim.  相似文献   

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

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