首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 62 毫秒
1.
Some 8000 images obtained with the Solar Eclipse Coronal Imaging System (SECIS) fast-frame CCD camera instrument located at Lusaka, Zambia, during the total eclipse of 21 June 2001 have been analysed to search for short-period oscillations in intensity that could be a signature of solar coronal heating mechanisms by MHD wave dissipation. Images were taken in white-light and Fe xiv green-line (5303 ?) channels over 205 seconds (frame rate 39 s−1), approximately the length of eclipse totality at this location, with a pixel size of four arcseconds square. The data are of considerably better quality than those that we obtained during the 11 August 1999 total eclipse (Rudawy et al.: Astron. Astrophys. 416, 1179, 2004), in that the images are much better exposed and enhancements in the drive system of the heliostat used gave a much improved image stability. Classical Fourier and wavelet techniques have been used to analyse the emission at 29 518 locations, of which 10 714 had emission at reasonably high levels, searching for periodic fluctuations with periods in the range 0.1 – 17 seconds (frequencies 0.06 – 10 Hz). While a number of possible periodicities were apparent in the wavelet analysis, none of the spatially and time-limited periodicities in the local brightness curves was found to be physically important. This implies that the pervasive Alfvén wave-like phenomena (Tomczyk et al.: Science 317, 1192, 2007) using polarimetric observations with the Coronal Multi-Channel Polarimeter (CoMP) instrument do not give rise to significant oscillatory intensity fluctuations.  相似文献   

2.
Since 4 December 2006, the SECCHI instrument suites onboard the two STEREO A and B probes have been imaging the solar corona and the heliosphere on a wide range of angular scales. The EUVI telescopes have a plate scale of 1.7 arcseconds pixel−1, while that of the HI2 wide-angle cameras is 2.15 arcminutes pixel−1, i.e. 75 times larger, with the COR1 and COR2 coronagraphs having intermediate plate scales. These very different instruments, aimed at studying Coronal Mass Ejections and their propagation in the heliosphere, create a data visualization challenge. This paper presents FESTIVAL, a SolarSoftware package originally developed to be able to map the SECCHI data into dynamic composite images of the sky as seen by the STEREO and SOHO probes. Data from other imaging instruments can also be displayed. Using the mouse, the user can quickly and easily zoom in and out and pan through these composite images to explore all spatial scales from EUVI to HI2 while keeping the native resolution of the original data. A large variety of numerical filters can be applied, and additional data (i.e. coordinate grids, stars catalogs, etc.) can be overlaid on the images. The architecture of FESTIVAL is such that it is easy to add support for other instruments and these new data immediately benefit from the already existing capabilities. Also, because its mapping engine is fully 3D, FESTIVAL provides a convenient environment to display images from future out-of-the-Ecliptic solar missions, such as Solar Orbiter or Solar Probe.  相似文献   

3.
DuneXpress     
The DuneXpress observatory will characterize interstellar and interplanetary dust in-situ, in order to provide crucial information not achievable with remote sensing astronomical methods. Galactic interstellar dust constitutes the solid phase of matter from which stars and planetary systems form. Interplanetary dust, from comets and asteroids, represents remnant material from bodies at different stages of early solar system evolution. Thus, studies of interstellar and interplanetary dust with DuneXpress in Earth orbit will provide a comparison between the composition of the interstellar medium and primitive planetary objects. Hence DuneXpress will provide insights into the physical conditions during planetary system formation. This comparison of interstellar and interplanetary dust addresses directly themes of highest priority in astrophysics and solar system science, which are described in ESA’s Cosmic Vision. The discoveries of interstellar dust in the outer and inner solar system during the last decade suggest an innovative approach to the characterization of cosmic dust. DuneXpress establishes the next logical step beyond NASA’s Stardust mission, with four major advancements in cosmic dust research: (1) analysis of the elemental and isotopic composition of individual interstellar grains passing through the solar system, (2) determination of the size distribution of interstellar dust at 1 AU from 10 − 14 to 10 − 9 g, (3) characterization of the interstellar dust flow through the planetary system, (4) establish the interrelation of interplanetary dust with comets and asteroids. Additionally, in supporting the dust science objectives, DuneXpress will characterize dust charging in the solar wind and in the Earth’s magnetotail. The science payload consists of two dust telescopes of a total of 0.1 m2 sensitive area, three dust cameras totaling 0.4 m2 sensitive area, and a nano-dust detector. The dust telescopes measure high-resolution mass spectra of both positive and negative ions released upon impact of dust particles. The dust cameras employ different detection methods and are optimized for (1) large area impact detection and trajectory analysis of submicron sized and larger dust grains, (2) the determination of physical properties, such as flux, mass, speed, and electrical charge. A nano-dust detector searches for nanometer-sized dust particles in interplanetary space. A plasma monitor supports the dust charge measurements, thereby, providing additional information on the dust particles. About 1,000 grains are expected to be recorded by this payload every year, with 20% of these grains providing elemental composition. During the mission submicron to micron-sized interstellar grains are expected to be recorded in statistically significant numbers. DuneXpress will open a new window to dusty universe that will provide unprecedented information on cosmic dust and on the objects from which it is derived.  相似文献   

4.
We studied the characteristics of Coronal Mass Ejections (CMEs) associated with solar flares and Deca-Hectometric (DH) type II radio bursts, based on source position during 23rd solar cycle (1997–2007). We classified these CME events into three groups using solar flare locations as, (i) disk events (0–30); (ii) intermediate events (31–60) and (iii) limb events (61–90). Main results from this studies are, (i) the number of CMEs associated with solar flares and DH-type IIs decreases as the source position approaches from disk to limb, (ii) most of the DH CMEs are halo (72%) in disk events and the number of occurrence of halo CMEs decreases from disk to limb, (iii) the average width and speed of limb events (164 and 1447 km s−1) are higher than those of disk events (134 and 1035 km s−1) and intermediate events (146 and 1170 km s−1) and (iv) the average accelerations for disk, intermediate and limb events are −8.2 m s−2, −10.3 m s−2 and −4.5 m s−2 respectively. These analysis of CMEs properties show more dependency on longitude and it gives strong evidence for projection effect.  相似文献   

5.
The Hinode Solar Optical Telescope (SOT) is the first space-borne visible-light telescope that enables us to observe magnetic-field dynamics in the solar lower atmosphere with 0.2 – 0.3 arcsec spatial resolution under extremely stable (seeing-free) conditions. To achieve precise measurements of the polarization with diffraction-limited images, stable pointing of the telescope (<0.09 arcsec, 3σ) is required for solar images exposed on the focal plane CCD detectors. SOT has an image stabilization system that uses image displacements calculated from correlation tracking of solar granules to control a piezo-driven tip-tilt mirror. The system minimizes the motions of images for frequencies lower than 14 Hz while the satellite and telescope structural design damps microvibration in higher frequency ranges. It has been confirmed from the data taken on orbit that the remaining jitter is less than 0.03 arcsec (3σ) on the Sun. This excellent performance makes a major contribution to successful precise polarimetric measurements with 0.2 – 0.3 arcsec resolution. K. Kobayashi now at NASA/Marshall Space Flight Center, Huntsville, AL 35812, USA.  相似文献   

6.
We investigate the properties of acoustic events (AEs), defined as spatially concentrated and short duration energy flux, in the quiet Sun, using observations of a 2D field of view (FOV) with high spatial and temporal resolution provided by the Solar Optical Telescope (SOT) onboard Hinode. Line profiles of Fe i 557.6 nm were recorded by the Narrow-band Filter Imager (NFI) on a 82″×82″ FOV during 75 min with a time step of 28.75 s and 0.08″ pixel size. Vertical velocities were computed at three atmospheric levels (80, 130, and 180 km) using the bisector technique, allowing the determination of energy flux to be made in the range 3 – 10 mHz using two complementary methods (Hilbert transform and Fourier power spectrum). Horizontal velocities were computed using local correlation tracking (LCT) of continuum intensities providing divergences. We found that the net energy flux is upward. In the range 3 – 10 mHz, a full FOV space and time averaged flux of 2700 W m−2 (lower layer 80 – 130 km) and 2000 W m−2 (upper layer 130 – 180 km) is concentrated in less than 1 % of the solar surface in the form of narrow (0.3″) AE. Their total duration (including rise and decay) is of the order of 103 s. Inside each AE, the mean flux is 1.6×105 W m−2 (lower layer) and 1.2×105 W m−2 (upper). Each event carries an average energy (flux integrated over space and time) of 2.5×1019 J (lower layer) to 1.9×1019 J (upper). More than 106 events could exist permanently on the Sun, with a birth and decay rate of 3500 s−1. Most events occur in intergranular lanes, downward velocity regions, and areas of converging motions.  相似文献   

7.
J. Javaraiah 《Solar physics》2011,270(2):463-483
Using the combined Greenwich (1874 – 1976) and Solar Optical Observatories Network (1977 – 2009) data on sunspot groups, we study the long-term variations in the mean daily rates of growth and decay of sunspot groups. We find that the minimum and the maximum values of the annually averaged daily mean growth rates are ≈ 52% day−1 and ≈ 183% day−1, respectively, whereas the corresponding values of the annually averaged daily mean decay rates are ≈ 21% day−1 and ≈ 44% day−1, respectively. The average value (over the period 1874 – 2009) of the growth rate is about 70% more than that of the decay rate. The growth and the decay rates vary by about 35% and 13%, respectively, on a 60-year time scale. From the beginning of Cycle 23 the growth rate is substantially decreased and near the end (2007 – 2008) the growth rate is lowest in the past about 100 years. In the extended part of the declining phase of this cycle, the decay rate steeply increased and it is largest in the beginning of the current Cycle 24. These unusual properties of the growth and the decay rates during Cycle 23 may be related to cause of the very long declining phase of this cycle with the unusually deep and prolonged current minimum. A ≈ 11-year periodicity in the growth and the decay rates is found to be highly latitude and time dependent and seems to exist mainly in the 0° – 10° latitude interval of the southern hemisphere. The strength of the known approximate 33 – 44-year modulation in the solar activity seems to be related to the north-south asymmetry in the growth rate. Decreasing and increasing trends in the growth and the decay rates indicate that the next 2 – 3 solar cycles will be relatively weak.  相似文献   

8.
The properties of solar magnetic fields on scales less than the spatial resolution of solar telescopes are studied. A synthetic infrared spectropolarimetric diagnostic based on a 2D MHD simulation of magnetoconvection is used for this. Analyzed are two time sequences of snapshots that likely represent two regions of the network fields with their immediate surroundings on the solar surface with unsigned magnetic flux densities of 300 and 140 G. In the first region from the probability density functions of the magnetic field strength it is found that the most probable field strength at log τ 5=0 is equal to 250 G. Weak fields (B<500 G) occupy about 70% of the surface, whereas stronger fields (B>1000 G) occupy only 9.7% of the surface. The magnetic flux is −28 G and its imbalance is −0.04. In the second region, these parameters are correspondingly equal to 150 G, 93.3%, 0.3%, −40 G, and −0.10. The distribution of line-of-sight velocities on the surface of log τ 5=−1 is estimated. The mean velocity is equal to 0.4 km s−1 in the first simulated region. The average velocity in the granules is −1.2 km s−1 and in the intergranules it is 2.5 km s−1. In the second region, the corresponding values of the mean velocities are equal to 0, −1.8, and 1.5 km s−1. In addition the asymmetry of synthetic Stokes V profiles of the Fe i 1564.8 nm line is analyzed. The mean values of the amplitude and area asymmetry do not exceed 1%. The spatially smoothed amplitude asymmetry is increased to 10% whereas the area asymmetry is only slightly varied.  相似文献   

9.
10.
F. Cavallini 《Solar physics》2006,236(2):415-439
A new instrument for solar bi-dimensional spectroscopy, the Interferometric BIdimensional Spectrometer (IBIS), has been successfully installed at the Dunn Solar Telescope of the National Solar Observatory (USA-NM) in June 2003. This instrument is essentially composed of a series of two Fabry-Perot interferometers and a set of narrow-band interference filters, used in a classic mount and in axial-mode. It has been designed to take monochromatic images of the solar surface with high spectral (R ≥ 200 000), spatial ≃ 0.2″), and temporal resolution (several frames s−1). IBIS has a circular field of view, 80″ in diameter and, with suitable interference filters, it can be used in the wavelength range 580 – 860 nm. The wavelength stability of the instrumental profile is very high, the maximum drift in 10 hours amounting to ≃10 m s−1. In this paper the criteria used in the design and the expected instrumental characteristics are described.  相似文献   

11.
We performed two-dimensional spectroscopic observations of the preceding sunspot of NOAA 10905 located off disk center (S8 E36, μ≈0.81) by using the Interferometric BI-dimensional Spectrometer (IBIS) operated at the Dunn Solar Telescope (DST) of the National Solar Observatory, New Mexico. The magnetically insensitive Fe I line at 709.04 nm was scanned in wavelength repetitively at an interval of 37 s to calculate sequences of maps of the line-wing and line-core intensity, and the line-of-sight Doppler velocity at different line depths (3% to 80%). Visual inspection of movies based on speckle reconstructions computed from simultaneous broadband data and the local continuum intensity at 709.04 nm revealed an umbral dot (UD) intruding rapidly from the umbral boundary to the center of the umbra. The apparent motion of this object was particularly fast (1.3 km s−1) when compared to typical UDs. The lifetime and size of the UD was 8.7 min and 240 km, respectively. The rapid UD was visible even in the line-core intensity map of Fe I 709.04 nm and was accompanied by a persistent blueshift of about 0.06 km s−1.  相似文献   

12.
As part of a program to estimate the solar spectrum back to the early twentieth century, we have generated fits to UV spectral irradiance measurements from 1 – 410 nm. The longer wavelength spectra (150 – 410 nm) were fit as a function of two solar activity proxies, the Mg ii core-to-wing ratio, or Mg ii index, and the total Ca ii K disk activity derived from ground based observations. Irradiance spectra at shorter wavelengths (1 – 150 nm) where used to generate fits to the Mg ii core-to-wing ratio alone. Two sets of spectra were used in these fitting procedures. The fits at longer wavelengths (150 to 410 nm) were derived from the high-resolution spectra taken by the Solar Ultraviolet Spectral Irradiance Monitor (SUSIM) on the Upper Atmospheric Research Satellite (UARS). Spectra measured by the Solar EUV Experiment (SEE) instrument on the Thermosphere Ionosphere Mesosphere Energetics and Dynamics (TIMED) satellite were used for the fits at wavelengths from 1 to 150 nm. To generate fits between solar irradiance and solar proxies, this study uses the above irradiance data, the NOAA composite Mg ii index, and daily Ca ii K disk activity determined from images measured by Big Bear Solar Observatory (BBSO). In addition to the fitting coefficients between irradiance and solar proxies, other results from this study include an estimated relationship between the fraction of the disk with enhanced Ca ii K activity and the Mg ii index, an upper bound of the average solar UV spectral irradiance during periods where the solar disk contains only regions of the quiet Sun, as was believed to be present during the Maunder Minimum, as well as results indicating that slightly more than 60% of the total solar irradiance (TSI) variability occurs between 150 and 400 nm.  相似文献   

13.
We report very high temporal and spectral resolution interferometric observations of some unusual solar radio bursts near 1420 MHz. These bursts were observed on 13 September 2005, 22 minutes after the peak of a GOES class X flare from the NOAA region 10808. Our observations show 11 episodes of narrow-band intermittent emission within a span of ≈ 8 s. Each episode shows a heavily frequency-modulated band of emission with a spectral slope of about −245.5 MHz s−1, comprising up to 8 individual blobs of emission and lasts for 10 – 15 ms. The blobs themselves have a spectral slope of ≈ 0 MHz s−1, are ≈ 200 – 250 kHz wide, appear every ≈ 400 kHz and last for ≈ 4 – 5 ms. These bursts show brightness temperatures in the range 1012 K, which suggests a coherent emission mechanism. We believe these are the first high temporal and spectral resolution interferometric observations of such rapid and narrow-bandwidth solar bursts close to 1420 MHz and present an analysis of their temporal and spectral characteristics.  相似文献   

14.
The Extreme-ultraviolet Variability Experiment (EVE; see Woods et al., 2009) obtains continuous EUV spectra of the Sun viewed as a star. Its primary objective is the characterization of solar spectral irradiance, but its sensitivity and stability make it extremely interesting for observations of variability on time scales down to the limit imposed by its basic 10 s sample interval. In this paper we characterize the Doppler sensitivity of the EVE data. We find that the 30.4 nm line of He ii has a random Doppler error below 0.001 nm (1 pm, better than 10 km s−1 as a redshift), with ample stability to detect the orbital motion of its satellite, the Solar Dynamics Observatory (SDO). Solar flares also displace the spectrum, both because of Doppler shifts and because of EVE’s optical layout, which (as with a slitless spectrograph) confuses position and wavelength. As a flare develops, the centroid of the line displays variations that reflect Doppler shifts and therefore flare dynamics. For the impulsive phase of the flare SOL2010-06-12, we find the line centroid to have a redshift of 16.8 ± 5.9 km s−1 relative to that of the flare gradual phase (statistical errors only). We find also that high-temperature lines, such as Fe xxiv 19.2 nm, have well-determined Doppler components for major flares, with decreasing apparent blueshifts as expected from chromospheric evaporation flows.  相似文献   

15.
Khabarova  O.  Zastenker  G. 《Solar physics》2011,270(1):311-329
Analysis of the Interball-1 spacecraft data (1995 – 2000) has shown that the solar wind ion flux sometimes increases or decreases abruptly by more than 20% over a time period of several seconds or minutes. Typically, the amplitude of such sharp changes in the solar wind ion flux (SCIFs) is larger than 0.5×108 cm−2 s−1. These sudden changes of the ion flux were also observed by the Solar Wind Experiment (SWE), on board the Wind spacecraft, as the solar wind density increases and decreases with negligible changes in the solar wind velocity. SCIFs occur irregularly at 1 AU, when plasma flows with specific properties come to the Earth’s orbit. SCIFs are usually observed in slow, turbulent solar wind with increased density and interplanetary magnetic field strength. The number of times SCIFs occur during a day is simulated using the solar wind density, magnetic field, and their standard deviations as input parameters for a period of five years. A correlation coefficient of ∼0.7 is obtained between the modelled and the experimental data. It is found that SCIFs are not associated with coronal mass ejections (CMEs), corotating interaction regions (CIRs), or interplanetary shocks; however, 85% of the sector boundaries are surrounded by SCIFs. The properties of the solar wind plasma for days with five or more SCIF observations are the same as those of the solar wind plasma at the sector boundaries. One possible explanation for the occurrence of SCIFs (near sector boundaries) is magnetic reconnection at the heliospheric current sheet or local current sheets. Other probable causes of SCIFs (inside sectors) are turbulent processes in the slow solar wind and at the crossings of flux tubes.  相似文献   

16.
The gamma-ray imager (GRI) is a novel mission concept that will provide an unprecedented sensitivity leap in the soft gamma-ray domain by using for the first time a focusing lens built of Laue diffracting crystals. The lens will cover an energy band from 200–1,300 keV with an effective area reaching 600 cm2. It will be complemented by a single reflection multilayer coated mirror, extending the GRI energy band into the hard X-ray regime, down to ∼10 keV. The concentrated photons will be collected by a position sensitive pixelised CZT stack detector. We estimate continuum sensitivities of better than 10 − 7 ph cm − 2s − 1keV − 1 for a 100 ks exposure; the narrow line sensitivity will be better than 3 × 10 − 6 ph cm − 2s − 1 for the same integration time. As focusing instrument, GRI will have an angular resolution of better than 30 arcsec within a field of view of roughly 5 arcmin—an unprecedented achievement in the gamma-ray domain. Owing to the large focal length of 100 m of the lens and the mirror, the optics and detector will be placed on two separate spacecrafts flying in formation in a high elliptical orbit. R&D work to enable the lens focusing technology and to develop the required focal plane detector is currently underway, financed by ASI, CNES, ESA, and the Spanish Ministery of Education and Science. The GRI mission has been proposed as class M mission for ESAs Cosmic Vision 2015–2025 program. GRI will allow studies of particle acceleration processes and explosion physics in unprecedented detail, providing essential clues on the innermost nature of the most violent and most energetic processes in the universe. All authors are on behalf of a large international collaboration The GRI mission has been proposed as an international collaboration between (in alphabetical order) Belgium (CSR), China (IHEP, Tsinghua Univ.), Denmark (DNSC, Southern Univ.), France (CESR, APC, ILL, CSNSM, IAP, LAM), Germany (MPE), Ireland (UCD School of Physics), Italy (INAF/IASF Rome, Bologna, Milano, Palermo; INAF/OA Brera, Roma; UNIFE, CNR/IMEM), Poland (NCAC), Portugal (Combra Univ., Evora Univ.), Russia (SINP, MSU, Ioffe Inst.), Spain (IEEC-CSIC-IFAE, CNM-IMB), the Netherlands (SRON, Utrecht Univ.), Turkey (Sabanci Univ.), United Kingdom (Univ. of Southampton, MSSL, RAL, Edinburgh Univ.), and the United States of America (SSL UC Berkeley, Argonne National Lab., MSFC, GSFC, US NRL).  相似文献   

17.
Solar rotation rate has been measured using the sunspot positions recorded by W.C. Bond during the period 1847 – 1849 at the Harvard College Observatory. From the drawings carried out by Bond we have selected the sunspots and groups of sunspots with more reliable positions presented in three or more drawings on successive days. We have calculated from the positions of the selected sunspots (41 in total) a synodic rotation rate of ω=[(12.92±0.08)−(1.5±1.0)sin 2 φ] degrees/day, where φ is the heliographic latitude. This rate, although slightly lower, is similar to the actual solar rotation rate, confirming no important changes in the solar rotation during the last 160 years.  相似文献   

18.
An experiment was conducted in conjunction with the total solar eclipse on 29 March 2006 in Libya to measure both the electron temperature and its flow speed simultaneously at multiple locations in the low solar corona by measuring the visible K-coronal spectrum. Coronal model spectra incorporating the effects of electron temperature and its flow speed were matched with the measured K-coronal spectra to interpret the observations. Results show electron temperatures of (1.10±0.05) MK, (0.70±0.08) MK, and (0.98±0.12) MK, at 1.1 R from Sun center in the solar north, east and west, respectively, and (0.93±0.12) MK, at 1.2 R from Sun center in the solar west. The corresponding outflow speeds obtained from the spectral fit are (103±92) km s−1, (0+10) km s−1, (0+10) km s−1, and (0+10) km s−1. Since the observations were taken only at 1.1 R and 1.2 R from Sun center, these speeds, consistent with zero outflow, are in agreement with expectations and provide additional confirmation that the spectral fitting method is working. The electron temperature at 1.1 R from Sun center is larger at the north (polar region) than the east and west (equatorial region).  相似文献   

19.
The Solar Optical Telescope (SOT) onboard Hinode aims to obtain vector magnetic fields on the Sun through precise spectropolarimetry of solar spectral lines with a spatial resolution of 0.2 – 0.3 arcsec. A photometric accuracy of 10−3 is achieved and, after the polarization calibration, any artificial polarization from crosstalk among Stokes parameters is required to be suppressed below the level of the statistical noise over the SOT’s field of view. This goal was achieved by the highly optimized design of the SOT as a polarimeter, extensive analyses and testing of optical elements, and an end-to-end calibration test of the entire system. In this paper we review both the approach adopted to realize the high-precision polarimeter of the SOT and its final polarization characteristics.  相似文献   

20.
The solar wind quasi-invariant (QI) has been defined by Osherovich, Fainberg, and Stone (Geophys. Res. Lett. 26, 2597, 1999) as the ratio of magnetic energy density and the energy density of the solar wind flow. In the regular solar wind QI is a rather small number, since the energy of the flow is almost two orders of magnitude greater than the magnetic energy. However, in magnetic clouds, QI is the order of unity (less than 1) and thus magnetic clouds can be viewed as a great anomaly in comparison with its value in the background solar wind. We study the duration, extent, and amplitude of this anomaly for two groups of isolated magnetic clouds: slow clouds (360<v<450 km s−1) and fast clouds (450≤v<720 km s−1). By applying the technique of superposition of epochs to 12 slow and 12 fast clouds from the catalog of Richardson and Cane (Solar Phys. 264, 189, 2010), we create an average slow cloud and an average fast cloud observed at 1 AU. From our analysis of these average clouds, we obtain cloud boundaries in both time and space as well as differences in QI amplitude and other parameters characterizing the solar wind state. Interplanetary magnetic clouds are known to cause major magnetic storms at the Earth, especially those clouds which travel from the sun to the Earth at high speeds. Characterizing each magnetic cloud by its QI value and extent may help in understanding the role of those disturbances in producing geomagnetic activity.  相似文献   

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

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