首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 171 毫秒
1.
On 21 September 2012, we carried out spectral observations of a solar facula in the Si?i 10827 Å, He?i 10830 Å, and H\(\upalpha\) spectral lines. Later, in the process of analyzing the data, we found a small-scale flare in the middle of the time series. Based on the anomalous increase in the absorption of the He?i 10830 Å line, we identified this flare as a negative flare.The aim of this article is to study the influence of the negative flare on the oscillation characteristics in the facular photosphere and chromosphere.We measured the line-of-sight (LOS) velocity and intensity of all the three lines as well as the half-width of the chromospheric lines. We also used the Helioseismic and Magnetic Imager (HMI) magnetic field data. The flare caused a modulation of all these parameters. In the location of the negative flare, the amplitude of the oscillations increased four times on average. In the adjacent magnetic field local maxima, the chromospheric LOS velocity oscillations appreciably decreased during the flare. The facular region oscillated as a whole with a 5-minute period before the flare, and this synchronicity was disrupted after the flare. The flare changed the spectral composition of the LOS magnetic field oscillations, causing an increase in the low-frequency oscillation power.  相似文献   

2.
We present an investigation of line-of-sight (LOS) velocity oscillations in solar faculae and sunspots. To study the phase relations between chromospheric and photospheric oscillations of the LOS velocity, we measured the time lag of the chromospheric signal relative to the photospheric one for several faculae and sunspots in a set of spectral line pairs. The measured time lags are different for different objects. The mean measured delay between the oscillations in the five-minute band in faculae is 50?s for the Si?i 10?827?Å?–?He?i 10?830?Å pair; for the pair Fe?i 6569?Å?–?Hα 6563?Å the mean delay is 20?s; for the pair Fe?i 4551?Å?–?Ba?ii 4554?Å the mean delay is 7?s; for the pair Si?i 8536?Å?–?Ca?ii 8542?Å the mean delay is 20?s. For the oscillations in the three-minute band in sunspot umbrae the mean delay is 55?s for the Si?i 10?827?Å?–?He?i 10?830?Å pair; for the Fe?i 6569?Å?–?Hα 6563?Å pair it was not possible to determine the delay; for the Fe?i 4551?Å?–?Ba?ii 4554?Å pair the mean delay is 6?s; for the Si?i 8536?Å?–?Ca?ii 8542?Å pair the mean delay is 21?s. Measured delays correspond to the wave propagation speed, which significantly exceeds the generally adopted speed of sound in the photosphere. This raises the question of the origin of these oscillations. The possibility that we deal with slow MHD waves is not ruled out.  相似文献   

3.
The differences between physical conditions in solar faculae and those in sunspots and quiet photosphere (increased temperature and different magnetic field topology) suggest that oscillation characteristics in facula areas may also have different properties. The analysis of 28 time series of simultaneous spectropolarimetric observations in facula photosphere (Fe?i 6569 Å, 8538 Å) and chromosphere (Hα, Ca?ii 8542 Å) yields the following results. The amplitude of five-minute oscillations of line-of-sight (LOS) velocity decreases by 20?–?40% in facula photosphere. There are only some cases revealing the inverse effect. The amplitude of four- to five-minute LOS velocity oscillations increases significantly in the chromosphere above faculae, and power spectra fairly often show pronounced peaks in a frequency range of 1.3?–?2.5 mHz. Evidence of propagating oscillations can be seen from space?–?time diagrams. We have found oscillations of the longitudinal magnetic field (1.5?–?2 mHz and 5.2 mHz) inside faculae.  相似文献   

4.
The Multi-Application Solar Telescope is a 50 cm off-axis Gregorian telescope recently installed at the Udaipur Solar Observatory, India. In order to obtain near-simultaneous observations at photospheric and chromospheric heights, an imager optimized for two or more wavelengths is being integrated with the telescope. Two voltage-tuneable lithium-niobate Fabry–Perot etalons along with a set of interference blocking filters have been used for developing the imager. Both of the etalons are used in tandem for photospheric observations in Fe i 6173 Å and chromospheric observation in H\(\alpha\) 6563 Å spectral lines, whereas only one of the etalons is used for the chromospheric Ca II line at 8542 Å. The imager is also being used for spectropolarimetric observations. We discuss the characterization of the etalons at the above wavelengths, detail the integration of the imager with the telescope, and present a few sets of observations taken with the imager set-up.  相似文献   

5.
We analyze the observations of a quiescent prominence acquired by the Téléscope Heliographique pour l’Étude du Magnetisme et des Instabilités Solaires (THEMIS) in the He?i 5876 Å (He?i D3) multiplet aiming to measure the spectral characteristics of the He?i D3 profiles and to find for them an adequate fitting model. The component characteristics of the He?i D3 Stokes I profiles are measured by the fitting system by approximating them with a double Gaussian. This model yields an He?i D3 component peak intensity ratio of \(5.5\pm0.4\), which differs from the value of 8 expected in the optically thin limit. Most of the measured Doppler velocities lie in the interval ±?5 km?s?1, with a standard deviation of ±?1.7 km?s?1 around the peak value of 0.4 km?s?1. The wide distribution of the full-width at half maximum has two maxima at 0.25 Å and 0.30 Å for the He?i D3 blue component and two maxima at 0.22 Å and 0.31 Å for the red component. The width ratio of the components is \(1.04\pm0.18\). We show that the double-Gaussian model systematically underestimates the blue wing intensities. To solve this problem, we invoke a two-temperature multi-Gaussian model, consisting of two double-Gaussians, which provides a better representation of He?i D3 that is free of the wing intensity deficit. This model suggests temperatures of 11.5 kK and 91 kK, respectively, for the cool and the hot component of the target prominence. The cool and hot components of a typical He?i D3 profile have component peak intensity ratios of 6.6 and 8, implying a prominence geometrical width of 17 Mm and an optical thickness of 0.3 for the cool component, while the optical thickness of the hot component is negligible. These prominence parameters seem to be realistic, suggesting the physical adequacy of the multi-Gaussian model with important implications for interpreting He?i D3 spectropolarimetry by current inversion codes.  相似文献   

6.
J. Koza 《Solar physics》2010,266(2):261-275
We examine the sensitivity of selected Ba?ii, Fe?i, Fe?ii, and Cr?i spectral lines to changes of the line-of-sight velocity by sharpness of their line profiles and response functions to line-of-sight velocity evaluated by the 1-D model of the quiet solar atmosphere in the LTE approximation. The set of selected lines includes the Ba?ii 4554 Å line, generally considered to be an excellent Doppler mapper. Our findings confirm earlier results showing that the sensitivity increases not only with wavelength, as anticipated from the Doppler relation, but mainly with the sharpness of line profiles given by the ratio of their depths and widths. The line Fe?i 5247 Å is the most sensitive in our set, whereas the Fe?i and Fe?ii infrared lines show very low sensitivity because of their large thermal widths. The line Ba?ii 4554 Å shows only moderate sensitivity due to its large width, given by a broad hyperfine structure and isotopic split. For the first time we identify a very promising and so far unknown Doppler mapper of the solar photosphere and low chromosphere, which is the line Ba?ii 6497 Å. Its sensitivity is comparable with the sensitivity of Fe?i 5247 Å and clearly surpasses the sensitivity of Ba?ii 4554 Å. The line Ba?ii 6497 Å offers many advantages, making it a highly recommendable choice for future studies of line-of-sight velocities in the photosphere and low chromosphere.  相似文献   

7.
This study based on longitudinal Zeeman effect magnetograms and spectral line scans investigates the dependence of solar surface magnetic fields on the spectral line used and the way the line is sampled to estimate the magnetic flux emerging above the solar atmosphere and penetrating to the corona from magnetograms of the Mt. Wilson 150-foot tower synoptic program (MWO). We have compared the synoptic program λ5250 Å line of Fe?i to the line of Fe?i at λ5233 Å since this latter line has a broad shape with a profile that is nearly linear over a large portion of its wings. The present study uses five pairs of sampling points on the λ5233 Å line. Line profile observations show that the determination of the field strength from the Stokes V parameter or from line bisectors in the circularly polarized line profiles lead to similar dependencies on the spectral sampling of the lines, with the bisector method being the less sensitive. We recommend adoption of the field determined with the line bisector method as the best estimate of the emergent photospheric flux and further recommend the use of a sampling point as close to the line core as is practical. The combination of the line profile measurements and the cross-correlation of fields measured simultaneously with λ5250 Å and λ5233 Å yields a formula for the scale factor δ ?1 that multiplies the MWO synoptic magnetic fields. By using ρ as the center-to-limb angle (CLA), a fit to this scale factor is δ ?1=4.15?2.82sin?2(ρ). Previously δ ?1=4.5?2.5sin?2(ρ) had been used. The new calibration shows that magnetic fields measured by the MDI system on the SOHO spacecraft are equal to 0.619±0.018 times the true value at a center-to-limb position 30°. Berger and Lites (2003, Solar Phys. 213, 213) found this factor to be 0.64±0.013 based on a comparison using the Advanced Stokes Polarimeter.  相似文献   

8.
We report observations of a white-light solar flare (SOL2010-06-12T00:57, M2.0) observed by the Helioseismic Magnetic Imager (HMI) on the Solar Dynamics Observatory (SDO) and the Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI). The HMI data give us the first space-based high-resolution imaging spectroscopy of a white-light flare, including continuum, Doppler, and magnetic signatures for the photospheric Fe i line at 6173.34 Å and its neighboring continuum. In the impulsive phase of the flare, a bright white-light kernel appears in each of the two magnetic footpoints. When the flare occurred, the spectral coverage of the HMI filtergrams (six equidistant samples spanning ±172 mÅ around nominal line center) encompassed the line core and the blue continuum sufficiently far from the core to eliminate significant Doppler crosstalk in the latter, which is otherwise a possibility for the extreme conditions in a white-light flare. RHESSI obtained complete hard X-ray and γ-ray spectra (this was the first γ-ray flare of Cycle 24). The Fe i line appears to be shifted to the blue during the flare but does not go into emission; the contrast is nearly constant across the line profile. We did not detect a seismic wave from this event. The HMI data suggest stepwise changes of the line-of-sight magnetic field in the white-light footpoints.  相似文献   

9.
We determine the optimal width and shape of the narrow-band filter centered on the He?i D3 line for prominence and coronal mass ejection (CME) observations with the ASPIICS (Association of Spacecraft for Polarimetric and Imaging Investigation of the Corona of the Sun) coronagraph onboard the PROBA-3 (Project for On-board Autonomy) satellite, to be launched in 2020. We analyze He?i D3 line intensities for three representative non-local thermal equilibrium prominence models at temperatures 8, 30, and 100 kK computed with a radiative transfer code and the prominence visible-light (VL) emission due to Thomson scattering on the prominence electrons. We compute various useful relations at prominence line-of-sight velocities of 0, 100, and 300 km?s?1 for 20 Å wide flat filter and three Gaussian filters with a full-width at half-maximum (FWHM) equal to 5, 10, and 20 Å to show the relative brightness contribution of the He?i D3 line and the prominence VL to the visibility in a given narrow-band filter. We also discuss possible signal contamination by Na?i D1 and D2 lines, which otherwise may be useful to detect comets. Our results mainly show that i) an optimal narrow-band filter should be flat or somewhere between flat and Gaussian with an FWHM of 20 Å in order to detect fast-moving prominence structures, ii) the maximum emission in the He?i D3 line is at 30 kK and the minimal at 100 kK, and iii) the ratio of emission in the He?i D3 line to the VL emission can provide a useful diagnostic for the temperature of prominence structures. This ratio is up to 10 for hot prominence structures, up to 100 for cool structures, and up to 1000 for warm structures.  相似文献   

10.
K. P. Raju 《Solar physics》2016,291(12):3519-3526
The transition region network seen in solar extreme ultraviolet (EUV) lines is the extension of the chromospheric network. The network appears as an irregular web-like pattern over the solar surface outside active regions. The average width of transition region network boundaries is obtained from the two-dimensional autocorrelation function of SOlar and Heliospheric Observatory (SOHO)/Coronal Diagnostic Spectrometer (CDS) synoptic images of the Sun in two emission lines, He i 586 Å and O v 630 Å during 1996?–?2012. The width of the network boundaries is found to be roughly correlated with the solar cycle variation with a lag of about ten months. A comparison of the widths in the two emission lines shows that they are larger for the He i line. The SOHO/CDS data also show large asymmetry in boundary widths in the horizontal (x) and vertical (y) image directions, which is shown to be caused by image distortions that are due to instrumental effects. Since the network boundary widths are related to the magnetic flux concentration along the boundaries, the results are expected to have implications on the flux transport on the solar surface, solar cycle, and the mass and energy budget of network loops and jets.  相似文献   

11.
Using the spectral data of representative solar flares observed with the infrared detector system of the solar spectrograph at Purple Mountain Observatory, we study the spectroscopic characteristics of solar flares in the Hα, the Ca i i 8?542 Å, and the He i 10?830 Å lines in different phases and various locations of flares and discuss their possible implications coupled with space observations. Our results show that in the initial phase of a flare the Hα line displays a red shift only with no wide wing. Large broadenings of the Hα line are observed a few minutes after the flare onset within small regions of 3?–?5′′ in both disk and limb flares with and without nonthermal processes. Far wings similar to those of damping broadening appear not only in the Hα line but in the He i 10?830 Å line as well in flares with nonthermal processes. Sometimes we even detect weak far-wing emission in the Ca i i 8?542 Å line in disk flares. Such large broadenings are observed in both the footpoints and the flare loop-top regions and possibly result from strong turbulence and/or macroscopic motions. Therefore, the so-called nonthermal wing of the Hα line profile is not a sufficient condition to distinguish whether nonthermal electrons are accelerated or not in a flare. The Ca i i 8?542 Å line shows lower intensity in the loop-top regions and higher intensity in the parts close to the solar surface. Emissions larger than nearby continuum in the He i 10?830 Å line are detected only in small regions with strong X-ray emissions and avoid sunspot umbrae.  相似文献   

12.
At a spatial resolution of 1″ and with spectral passbands of 0.25–0.47 Å, the chromospheric fine structure was studied with an array of 512 silicon photodiodes. The high quantum efficiency of the diodes in the near infrared allowed low noise spectroheliograms to be constructed from observations in the lines of H i, Ca ii, and He i. Magnetograms of the underlying photosphere were obtained simultaneously. Tachograms in the He 10830 Å line revealed 1″ points and elongated features that are interpreted as spicules seen against the disk. Active regions and filaments at 10830 are compared with Hα and Ca ii (8542 Å) features. Filament contrast increases with the proximity of bright plage. Twelve subflares were studied and in eight cases, 3–5″ kernels of He I emission appeared over small, growing pores or over 5″ patches of magnetic field emerging through the photosphere. All the subflares showed 10830 emission, contrary to established belief that the 10830 line goes into emission only in the largest flares. All the subflares included at least one emission kernel over regions where the photospheric magnetic field, as seen with a resolution of 1–2″, broke down into a mosaic of both polarities with 3–5″ diam. elements.  相似文献   

13.
An empirical model of solar UV spectral irradiance has been developed that is based on observed spectral radiance measurements and full disk Ca ii K images. The Mg ii index is then calculated from the estimated spectra in a narrow wavelength range (180 Å) near the Mg ii doublet at 2800 Å. Our long term goal is to expand this wavelength range from 10 to 4000 Å in continuing studies based on spectral data covering this wavelength range (e.g. Skylab, UARS/SUSIM, TIMED/SEE, etc.). Our previous modeling effort produced spectra in this 180 Å range and the resulting Mg ii index values for the period from 1991 through 1995 and we have used observations during this time period to validate the model results. The current paper presents results from this model based on a 21-year portion of the recently digitized Ca ii K images from the Mt Wilson Observatory (MWO) film archive. Here we present details of the model, the required model modifications, and the resulting Mg ii index from 1961 through 1981. Since the NOAA Mg ii index did not begin until 1978, the present model results are compared to a Mg ii index estimated from the F10.7 radio flux over this 21-year period. The NOAA Mg ii index, which is derived from measured UV spectra, is also included for comparison from late 1978 through 1981.  相似文献   

14.
Nanoflare-heating theory predicts steady hot-plasma emission in non-flaring active regions. It is hard to find this emission with conventional non-monochromatic imagers (such as the Atmospheric Imaging Assembly or the X-Ray Telescope), because their images contain a cool-temperature background. In this work, we search for hot plasma in non-flaring active regions using the Mg?xii spectroheliograph onboard the Complex Orbital Observations Near-Earth of Activity on the Sun (CORONAS)-F/SPectroheliographIc X-ray Imaging Telescope (SPIRIT). This instrument acquired monochromatic images of the solar corona in the Mg?xii 8.42 Å line, which emits only at temperatures higher than 4 MK. The Mg?xii images only contain the signal from hot plasma, without any low-temperature background. We studied the hot plasma in active regions using SPIRIT data from 18?–?28 February 2002. During this period, the Mg?xii spectroheliograph worked with a 105-second cadence almost without data gaps. Hot plasma was observed only in the flaring active regions. We did not observe any hot plasma in non-flaring active regions. The hot-plasma column emission measure in the non-flaring active region is not expected to exceed \(3 \times10^{24}\) cm?5. The hot differential emission measure is lower than 0.01% of the DEM of the main temperature component. The absence of Mg?xii emission in the non-flaring active regions can be explained by weak and frequent nanoflares (with a delay of less than 500 seconds) or by very short and intense nanoflares that lead to non-equilibrium ionization.  相似文献   

15.
We estimate the electron density, \(n_{\mathrm{e}}\), and its spatial variation in quiescent prominences from the observed emission ratio of the resonance lines Na?i?5890 Å (D2) and Sr?ii?4078 Å. For a bright prominence (\(\tau_{\alpha}\approx25\)) we obtain a mean \(n_{\mathrm{e}}\approx2\times10^{10}~\mbox{cm}^{-3}\); for a faint one (\(\tau _{\alpha }\approx4\)) \(n_{\mathrm{e}}\approx4\times10^{10}~\mbox{cm}^{-3}\) on two consecutive days with moderate internal fluctuation and no systematic variation with height above the solar limb. The thermal and non-thermal contributions to the line broadening, \(T_{\mathrm{kin}}\) and \(V_{\mathrm{nth}}\), required to deduce \(n_{\mathrm{e}}\) from the emission ratio Na?i/Sr?ii cannot be unambiguously determined from observed widths of lines from atoms of different mass. The reduced widths, \(\Delta\lambda_{\mathrm{D}}/\lambda_{0}\), of Sr?ii?4078 Å show an excess over those from Na?D2 and \(\mbox{H}\delta\,4101\) Å, assuming the same \(T_{\mathrm{kin}}\) and \(V_{\mathrm{nth}}\). We attribute this excess broadening to higher non-thermal broadening induced by interaction of ions with the prominence magnetic field. This is suggested by the finding of higher macro-shifts of Sr?ii?4078 Å as compared to those from Na?D2.  相似文献   

16.
K. P. Raju 《Solar physics》2009,255(1):119-129
Relative Doppler velocities and spectral linewidths in a coronal hole and in the quiet Sun region outside have been obtained from Solar and Heliospheric Observatory (SOHO)/Coronal Diagnostic Spectrometer (CDS) observations. Five strong emission lines in the CDS wavelength range (namely, O? iii 599 Å, O?v 630 Å, Ne?vi 562.8 Å, He?ii 304 Å, and Mg?ix 368 Å), whose formation temperatures represent different heights in the solar atmosphere from the lower transition region to the inner corona, have been used in the study. As reported earlier, relative velocities in the coronal hole are generally blueshifted with respect to the quiet Sun, and the magnitude of the blueshifts increases with height. It has been found that the polar coronal hole has larger relative velocities than the equatorial extension in the inner corona. Several localized velocity contours have been found mainly on network brightenings and in the vicinity of the coronal hole boundary. The presence of velocity contours on the network may represent network outflows whereas the latter could be due to localized jets probably arising from magnetic reconnection at the boundary. All spectral lines have larger widths in the coronal hole than in the quiet Sun. In O?v 630 Å an extended low-linewidth region is seen in the coronal hole?–?quiet Sun boundary, which may indicate fresh mass transfer across the boundary. Also polar coronal holes have larger linewidths in comparison with the equatorial extension. Together with larger relative velocities, this suggests that the solar wind emanating from polar hole regions is faster than that from equatorial hole regions.  相似文献   

17.
We present moderate-resolution (<5 Å) long-slit optical spectra of 51 nebular objects in the nearby Sculptor Group galaxy NGC 300 obtained with the 2.3 meter Advanced Technology Telescope at Siding Spring Observatory, Australia. Adopting the criterion of [S?ii]Total:Hα≥0.4 to confirm supernova remnants (SNRs) from optical spectra, we find that of 28 objects previously proposed as SNRs from optical observations, 22 meet this criterion with six showing [S?ii]Total:Hα of less than 0.4. Of 27 objects suggested as SNRs from radio data, four are associated with the 28 previously proposed SNRs. Of these four, three (included in the 22 above) meet the criterion. In all, 22 of the 51 nebular objects meet the [S ii]Total:Hα criterion as SNRs while the nature of the remaining 29 objects remains undetermined by these observations.  相似文献   

18.
The Helioseismic and Magnetic Imager (HMI) instrument onboard the Solar Dynamics Observatory produces line-of-sight (LOS) observables (Doppler velocity, magnetic-field strength, Fe i line width, line depth, and continuum intensity) as well as vector magnetic-field maps at the solar surface. The accuracy of LOS observables is dependent on the algorithm used to translate a sequence of HMI filtergrams into the corresponding observables. Using one hour of high-cadence imaging spectropolarimetric observations of a sunspot in the Fe i line at 6173 Å through the Interferometric Bidimensional Spectrometer installed at the Dunn Solar Telescope, and the Milne–Eddington inversion of the corresponding Stokes vectors, we test the accuracy of the observables algorithm currently implemented in the HMI data-analysis pipeline: the MDI-like algorithm. In an attempt to improve the accuracy of HMI observables, we also compare this algorithm to others that may be implemented in the future: a least-squares fit with a Gaussian profile, a least-squares fit with a Voigt profile, and the use of second Fourier coefficients in the MDI-like algorithm.  相似文献   

19.
As a coronal mass ejection (CME) passes, the flank and wake regions are typically strongly disturbed. Various instruments, including the Large Angle and Spectroscopic Coronagraph (LASCO), the Atmospheric Imaging Assembly (AIA), and the Coronal Multi-channel Polarimeter (CoMP), observed a CME close to the east limb on 26 October 2013. A hot (\({\approx}\,10~\mbox{MK}\)) rising blob was detected on the east limb, with an initial ejection flow speed of \({\approx}\, 330~\mbox{km}\,\mbox{s}^{-1}\). The magnetic structures on both sides and in the wake of the CME were strongly distorted, showing initiation of turbulent motions with Doppler-shift oscillations enhanced from \({\approx}\, \pm 3~\mbox{km}\,\mbox{s}^{-1}\) to \({\approx}\, \pm 15~\mbox{km}\,\mbox{s}^{-1}\) and effective thermal velocities from \({\approx}\,30~\mbox{km}\,\mbox{s}^{-1}\) to \({\approx}\,60~\mbox{km}\,\mbox{s}^{-1}\), according to the CoMP observations at the Fe?xiii line. The CoMP Doppler-shift maps suggest that the turbulence behaved differently at various heights; it showed clear wave-like torsional oscillations at lower altitudes, which are interpreted as the antiphase oscillation of an alternating red/blue Doppler shift across the strands at the flank. The turbulence seems to appear differently in the channels of different temperatures. Its turnover time was \({\approx}\,1000\) seconds for the Fe 171 Å channel, while it was \({\approx}\,500\) seconds for the Fe 193 Å channel. Mainly horizontal swaying rotations were observed in the Fe 171 Å channel, while more vertical vortices were seen in the Fe 193 Å channel. The differential-emission-measure profiles in the flank and wake regions have two components that evolve differently: the cool component decreased over time, evidently indicating a drop-out of cool materials due to ejection, while the hot component increased dramatically, probably because of the heating process, which is suspected to be a result of magnetic reconnection and turbulence dissipation. These results suggest a new turbulence-heating scenario of the solar corona and solar wind.  相似文献   

20.
Using mainly the 1600 Å continuum channel and also the 1216 Å Lyman-α channel (which includes some UV continuum and C iv emission) aboard the TRACE satellite, we observed the complete lifetime of a transient, bright chromospheric loop. Simultaneous observations with the SUMER instrument aboard the SOHO spacecraft revealed interesting material velocities through the Doppler effect existing above the chromospheric loop imaged with TRACE, possibly corresponding to extended nonvisible loops, or the base of an X-ray jet.  相似文献   

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

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