首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 15 毫秒
1.
The structure of electric current and magnetic helicity in the solar corona is closely linked to solar activity over the 11-year cycle, yet is poorly understood. As an alternative to traditional current-free “potential-field” extrapolations, we investigate a model for the global coronal magnetic field which is non-potential and time-dependent, following the build-up and transport of magnetic helicity due to flux emergence and large-scale photospheric motions. This helicity concentrates into twisted magnetic flux ropes, which may lose equilibrium and be ejected. Here, we consider how the magnetic structure predicted by this model – in particular the flux ropes – varies over the solar activity cycle, based on photospheric input data from six periods of cycle 23. The number of flux ropes doubles from minimum to maximum, following the total length of photospheric polarity inversion lines. However, the number of flux rope ejections increases by a factor of eight, following the emergence rate of active regions. This is broadly consistent with the observed cycle modulation of coronal mass ejections, although the actual rate of ejections in the simulation is about a fifth of the rate of observed events. The model predicts that, even at minimum, differential rotation will produce sheared, non-potential, magnetic structure at all latitudes.  相似文献   

2.
Studying the evolution of magnetic clouds entrained in coronal mass ejections using in-situ data is a difficult task, since only a limited number of observational points is available at large heliocentric distances. Remote sensing observations can, however, provide important information for events close to the Sun. In this work we estimate the flux rope orientation first in the close vicinity of the Sun (2?–?20 R ) using forward modeling of STEREO/SECCHI and SOHO/LASCO coronagraph images of coronal mass ejections and then in situ using Grad–Shafranov reconstruction of the magnetic cloud. Thus, we are able to measure changes in the orientation of the erupted flux ropes as they propagate from the Sun to 1 AU. We present both techniques and use them to study 15 magnetic clouds observed during the minimum following Solar Cycle 23 and the rise of Solar Cycle 24. This is the first multievent study to compare the three-dimensional parameters of CMEs from imaging and in-situ reconstructions. The results of our analysis confirm earlier studies showing that the flux ropes tend to deflect towards the solar equatorial plane. We also find evidence of rotation on their travel from the Sun to 1 AU. In contrast to past studies, our method allows one to deduce the evolution of the three-dimensional orientation of individual flux ropes rather than on a statistical basis.  相似文献   

3.
The hemispheric pattern of solar filaments is considered in the context of the global magnetic field of the solar corona. In recent work Mackay and van Ballegooijen have shown how, for a pair of interacting magnetic bipoles, the observed chirality pattern could be explained by the dominant range of bipole tilt angles and helicity in each hemisphere. This study aims to test this earlier result through a direct comparison between theory and observations, using newly developed simulations of the actual surface and 3D coronal magnetic fields over a 6-month period, on a global scale. We consider two key components: (1) observations of filament chirality for the sample of 255 filaments and (2) our new simulations of the large-scale surface magnetic field. Based on a flux-transport model, these will be used as the lower boundary condition for the future 3D coronal simulations. Our technique differs significantly from those of other authors, where the coronal field is either assumed to be purely potential or has to be reset back to potential every 27 days for the photospheric field to remain accurate. In our case we ensure accuracy by the insertion of newly emerging bipolar active regions, based on observed photospheric synoptic magnetograms. The large-scale surface field is shown to remain accurate over the 6-month period, without any resetting. This new technique will enable future simulations to consider the long-term buildup and transport of helicity and shear in the coronal magnetic field over many months or years.  相似文献   

4.
Flux ropes are twisted magnetic structures that can be detected by in-situ measurements in the solar wind. However, different properties of detected flux ropes suggest different types of flux-rope populations. As such, are there different populations of flux ropes? The answer is positive and is the result of the analysis of four lists of flux ropes, including magnetic clouds (MCs), observed at 1 AU. The in-situ data for the four lists were fitted with the same cylindrical force-free field model, which provides an estimate of the local flux-rope parameters such as its radius and orientation. Since the flux-rope distributions have a broad dynamic range, we went beyond a simple histogram analysis by developing a partition technique that uniformly distributes the statistical fluctuations across the radius range. By doing so, we found that small flux ropes with radius R<0.1 AU have a steep power-law distribution in contrast to the larger flux ropes (identified as MCs), which have a Gaussian-like distribution. Next, from four CME catalogs, we estimated the expected flux-rope frequency per year at 1 AU. We found that the predicted numbers are similar to the frequencies of MCs observed in-situ. However, we also found that small flux ropes are at least ten times too abundant to correspond to CMEs, even to narrow ones. Investigating the different possible scenarios for the origin of these small flux ropes, we conclude that these twisted structures can be formed by blowout jets in the low corona or in coronal streamers.  相似文献   

5.
It is generally accepted that the presence of a hot magnetic corona provides the source of X-ray emission in cool stars. With this connection one could expect to see the variation of magnetic flux in the activity cycle of a star mirrored by a similar variation in the stars X-ray emission. Using magnetic maps produced from flux emergence and transport simulations and assuming a potential field for the corona, we can extrapolate the coronal magnetic field and hence calculate the variation of the X-ray emission. We consider three types of activity cycle that successfully reproduce the pattern of intermingled magnetic flux at high latitudes, a feature observed with Zeeman–Doppler imaging. The three different cycles take the form of (1) an enhanced butterfly pattern where flux emergence is extended to a latitude of 70°, (2) an extended emergence profile as before but with an overlap of 4 yr in the butterfly diagram and (3) where no butterfly diagram is used. The cyclic variation in the X-ray emission is around two orders of magnitude for cases (1) and (3), but less than one order of magnitude for case (2). For all three cases, the rotational modulation of the X-ray emission is greatest at cycle minimum, but the emission measure weighted density varies little over the cycle. For cases (1) and (2) the fraction of the total flux that is open (along which a wind can escape) varies little over the cycle, but for case (3) this is three times larger at cycle minimum than at maximum. Our results clearly show that although magnetic cycles may exist for stars they are not necessarily observable in the X-ray emission.  相似文献   

6.
G. de Toma 《Solar physics》2011,274(1-2):195-217
We analyze coronal holes present on the Sun during the extended minimum between Cycles 23 and 24, study their evolution, examine the consequences for the solar wind speed near the Earth, and compare it with the previous minimum in 1996. We identify coronal holes and determine their size and location using a combination of EUV observations from SOHO/EIT and STEREO/EUVI and magnetograms. We find that the long period of low solar activity from 2006 to 2009 was characterized by weak polar magnetic fields and polar coronal holes smaller than observed during the previous minimum. We also find that large, low-latitude coronal holes were present on the Sun until 2008 and remained important sources of recurrent high-speed solar wind streams. By the end of 2008, these low-latitude coronal holes started to close down, and finally disappeared in 2009, while smaller, mid-latitude coronal holes formed in the remnants of Cycle 24 active regions shifting the sources of the solar wind at the Earth to higher latitudes.  相似文献   

7.
The Mechanism involved in the Reversals of the Sun's Polar Magnetic Fields   总被引:2,自引:0,他引:2  
Durrant  C.J.  Turner  J.P.R.  Wilson  P.R. 《Solar physics》2004,222(2):345-362
Models of the polarity reversals of the Sun's polar magnetic fields based on the surface transport of flux are discussed and are tested using observations of the polar fields during Cycle 23 obtained by the National Solar Observatory at Kitt Peak. We have extended earlier measurements of the net radial flux polewards of ±60° and confirm that, despite fluctuations of 20%, there is a steady decline in the old polarity polar flux which begins shortly after sunspot minimum (although not at the same time in each hemisphere), crosses the zero level near sunspot maximum, and increases, with reversed polarity during the remainder of the cycle. We have also measured the net transport of the radial field by both meridional flow and diffusion across several latitude zones at various phases of the Cycle. We can confirm that there was a net transport of leader flux across the solar equator during Cycle 23 and have used statistical tests to show that it began during the rising phase of this cycle rather than after sunspot maximum. This may explain the early decrease of the mean polar flux after sunspot minimum. We also found an outward flow of net flux across latitudes ±60° which is consistent with the onset of the decline of the old polarity flux. Thus the polar polarity reversals during Cycle 23 are not inconsistent with the surface flux-transport models but the large empirical values required for the magnetic diffusivity require further investigation.  相似文献   

8.
We study the solar-cycle variation of subsurface flows from the surface to a depth of 16 Mm. We have analyzed Global Oscillation Network Group (GONG) Dopplergrams with a ring-diagram analysis covering about 15 years and Helioseismic and Magnetic Imager (HMI) Dopplergrams covering more than 6 years. After subtracting the average rotation rate and meridional flow, we have calculated the divergence of the horizontal residual flows from the maximum of Solar Cycle 23 through the declining phase of Cycle 24. The subsurface flows are mainly divergent at quiet regions and convergent at locations of high magnetic activity. The relationship is essentially linear between divergence and magnetic activity at all activity levels at depths shallower than about 10 Mm. At greater depths, the relationship changes sign at locations of high activity; the flows are increasingly divergent at locations with a magnetic activity index (MAI) greater than about 24 G. The flows are more convergent by about a factor of two during the rising phase of Cycle 24 than during the declining phase of Cycle 23 at locations of medium and high activity (about 10 to 40 G MAI) from the surface to at least 10 Mm. The subsurface divergence pattern of Solar Cycle 24 first appears during the declining phase of Cycle 23 and is present during the extended minimum. It appears several years before the magnetic pattern of the new cycle is noticeable in synoptic maps. Using linear regression, we estimate the amount of magnetic activity that would be required to generate the precursor pattern and find that it should be almost twice the amount of activity that is observed.  相似文献   

9.
Durrant  C.J.  Mccloughan  J. 《Solar physics》2004,219(1):55-78
We describe the application of the synoptic transport equation to simulate the temporal evolution of the magnetic flux over the solar surface. This provides a means of predicting each day both the synoptic maps for the Carrington rotation starting the next day and the instantaneous map of the solar flux over the whole solar surface for the next day. The reliability of the predicted synoptic maps is tested by comparing the locations of the zero-flux contour with those of the observed maps produced by the National Solar Observatory, Kitt Peak and with the locations of Hα filaments measured on filtergrams obtained by the Big Bear Solar Observatory. We conclude that the best match at high latitudes is obtained by long-term simulations (over 20 rotations) with flux updates each rotation between latitudes ± 60°. We illustrate the use of the simulations to describe the evolution of the polar fields at the time of the polarity reversals in Cycle 23. The reconstruction of the instantaneous maps is tested by comparison with full-disk magnetograms. The method provides a simple means of estimating the large-scale flux distribution over the whole surface. It does not take account of flux emerging after the central meridian passage each rotation so it is only approximate in the activity belts but provides a reliable map beyond those latitudes.  相似文献   

10.
Flux ropes ejected from the Sun may change their geometrical orientation during their evolution, which directly affects their geoeffectiveness. Therefore, it is crucial to understand how solar flux ropes evolve in the heliosphere to improve our space-weather forecasting tools. We present a follow-up study of the concepts described by Isavnin, Vourlidas, and Kilpua (Solar Phys. 284, 203, 2013). We analyze 14 coronal mass ejections (CMEs), with clear flux-rope signatures, observed during the decay of Solar Cycle 23 and rise of Solar Cycle 24. First, we estimate initial orientations of the flux ropes at the origin using extreme-ultraviolet observations of post-eruption arcades and/or eruptive prominences. Then we reconstruct multi-viewpoint coronagraph observations of the CMEs from ≈?2 to 30 R with a three-dimensional geometric representation of a flux rope to determine their geometrical parameters. Finally, we propagate the flux ropes from ≈?30 R to 1 AU through MHD-simulated background solar wind while using in-situ measurements at 1 AU of the associated magnetic cloud as a constraint for the propagation technique. This methodology allows us to estimate the flux-rope orientation all the way from the Sun to 1 AU. We find that while the flux-ropes’ deflection occurs predominantly below 30 R, a significant amount of deflection and rotation happens between 30 R and 1 AU. We compare the flux-rope orientation to the local orientation of the heliospheric current sheet (HCS). We find that slow flux ropes tend to align with the streams of slow solar wind in the inner heliosphere. During the solar-cycle minimum the slow solar-wind channel as well as the HCS usually occupy the area in the vicinity of the solar equatorial plane, which in the past led researchers to the hypothesis that flux ropes align with the HCS. Our results show that exceptions from this rule are explained by interaction with the Parker-spiraled background magnetic field, which dominates over the magnetic interaction with the HCS in the inner heliosphere at least during solar-minimum conditions.  相似文献   

11.
The solar photospheric magnetic flux distribution is key to structuring the global solar corona and heliosphere. Regular full-disk photospheric magnetogram data are therefore essential to our ability to model and forecast heliospheric phenomena such as space weather. However, our spatio-temporal coverage of the photospheric field is currently limited by our single vantage point at/near Earth. In particular, the polar fields play a leading role in structuring the large-scale corona and heliosphere, but each pole is unobservable for \({>}\,6\) months per year. Here we model the possible effect of full-disk magnetogram data from the Lagrange points \(L_{4}\) and \(L_{5}\), each extending longitude coverage by \(60^{\circ}\). Adding data also from the more distant point \(L_{3}\) extends the longitudinal coverage much further. The additional vantage points also improve the visibility of the globally influential polar fields. Using a flux-transport model for the solar photospheric field, we model full-disk observations from Earth/\(L_{1}\), \(L_{3}\), \(L_{4}\), and \(L_{5}\) over a solar cycle, construct synoptic maps using a novel weighting scheme adapted for merging magnetogram data from multiple viewpoints, and compute potential-field models for the global coronal field. Each additional viewpoint brings the maps and models into closer agreement with the reference field from the flux-transport simulation, with particular improvement at polar latitudes, the main source of the fast solar wind.  相似文献   

12.
Based on 11 years of SOHO/MDI observations from the cycle minimum in 1997 to the next minimum around 2008, we compare observed and modeled axial dipole moments to better understand the large-scale transport properties of magnetic flux in the solar photosphere. The absolute value of the axial dipole moment in 2008 is less than half that in the corresponding cycle-minimum phase in early 1997, both as measured from synoptic maps and as computed from an assimilation model based only on magnetogram data equatorward of 60° in latitude. This is incompatible with the statistical fluctuations expected from flux-dispersal modeling developed in earlier work at the level of 7 – 10 σ. We show how this decreased axial dipole moment can result from an increased strength of the diverging meridional flow near the Equator, which more effectively separates the two hemispheres for dispersing magnetic flux. Based on the combination of this work with earlier long-term simulations of the solar surface field, we conclude that the flux-transport properties across the solar surface have changed from preceding cycles to the most recent one. A plausible candidate for such a change is an increase of the gradient of the meridional-flow pattern near the Equator so that the two hemispheres are more effectively separated. The required profile as a function of latitude is consistent with helioseismic and cross-correlation measurements made over the past decade.  相似文献   

13.
Using in situ observations from the Advanced Composition Explorer (ACE), we have identified 70 Earth-affecting interplanetary coronal mass ejections (ICMEs) in Solar Cycle 24. Because of the unprecedented extent of heliospheric observations in Cycle 24 that has been achieved thanks to the Sun Earth Connection Coronal and Heliospheric Investigation (SECCHI) instruments onboard the Solar Terrestrial Relations Observatory (STEREO), we observe these events throughout the heliosphere from the Sun to the Earth, and we can relate these in situ signatures to remote sensing data. This allows us to completely track the event back to the source of the eruption in the low corona. We present a summary of the Earth-affecting CMEs in Solar Cycle 24 and a statistical study of the properties of these events including the source region. We examine the characteristics of CMEs that are more likely to be strongly geoeffective and examine the effect of the flare strength on in situ properties. We find that Earth-affecting CMEs in the first half of Cycle 24 are more likely to come from the northern hemisphere, but after April 2012, this reverses, and these events are more likely to originate in the southern hemisphere, following the observed magnetic asymmetry in the two hemispheres. We also find that as in past solar cycles, CMEs from the western hemisphere are more likely to reach Earth. We find that Cycle 24 lacks in events driving extreme geomagnetic storms compared to past solar cycles.  相似文献   

14.
It is a basic feature of the Babcock-Leighton model of the solar cycle that the polar field reversal is due to the diffusive decay and poleward drift of the active region fields. The flux from follower regions moves preferentially polewards in each hemisphere, where it cancels with, and then replaces, the previously existing polar fields. A number of workers have attempted to model this process by numerical solutions of the flux transport equation, which include the surface effects of supergranule diffusion, differential rotation and meridional flow, with conflicting results.Here we describe recent changes in the polar fields using synoptic magnetic data provided by the Mount Wilson Observatory, and compare them with simulations using the flux transport equation and based on the observed fields for Carrington rotation 1815. These changes include a part-reversal of the north polar field. It is shown that the evolution of the polar fields cannot be reproduced accurately by simulations of the diffusion and poleward drift of the emerging active regions at sunspot latitudes.Histograms of the distribution of the field intensities derived from the daily magnetograms obtained at the Kitt Peak Station of the National Solar Observatory provide independent evidence that flux is emerging at high latitudes and that this flux makes a contribution to the evolution of these patterns. This implies the presence of some form of sub-surface dynamo action at high latitudes.On leave from the School of Mathematics, University of Sydney.  相似文献   

15.
We studied the soft X-ray solar events that could be associated with the interplanetary magnetic flux ropes observed by the WIND satellite during 1995 through 1998. The timings of the launches of the magnetic flux ropes from the Sun were estimated using flux rope speeds derived by the fitting of a cylindrical model. In the reasonable time window, soft X-ray solar signatures were found in approximately 70% of the flux ropes. Parameters (e.g., axis direction, strength of magnetic field, radius, and helicity) of the magnetic flux ropes obtained by the model fitting were compared with the characteristics of the corresponding soft X-ray events observed by Yohkoh. According to the result of the comparison, the magnetic flux ropes with strong magnetic fields or high speeds were observed in association with higher soft X-ray solar activities.  相似文献   

16.
The solar-cycle oscillations of the toroidal and poloidal components of the solar magnetic field in the northern solar hemisphere have a persistent phase difference of about \(\pi \). We propose a symmetrical Kuramoto model with three coupled oscillators as a simple way to understand this anti-synchronization. We solve an inverse problem and reconstruct natural frequencies of the top and bottom oscillators under the conditions of a constant coupling strength and a non-delayed coupling. These natural frequencies are associated with angular velocities of the meridional flow circulation near the solar surface and in the deep layer of the solar convection zone. A relationship between our reconstructions of the shallow and the deep meridional flow speed during recent Solar Cycles 21?–?23 is in agreement with estimates obtained in helioseismology and flux-transport dynamo modeling. The reconstructed top oscillator speed presents significant solar-cycle like variations that agree with recent helioseismical reconstructions. The evolution of reconstructed natural frequencies strongly depends on the coupling strength. We find two stable regimes in the case of strong coupling with a change of regime during anomalous solar cycles. We see the onset of a new transition in Solar Cycle 24. We estimate the admitted range of coupling values and find evidence of cross-equatorial coupling between solar hemispheres not accounted for by the model.  相似文献   

17.
M. A. Weber  Y. Fan  M. S. Miesch 《Solar physics》2013,287(1-2):239-263
We study how active-region-scale flux tubes rise buoyantly from the base of the convection zone to near the solar surface by embedding a thin flux tube model in a rotating spherical shell of solar-like turbulent convection. These toroidal flux tubes that we simulate range in magnetic field strength from 15 kG to 100 kG at initial latitudes of 1° to 40° in both hemispheres. This article expands upon Weber, Fan, and Miesch (Astrophys. J. 741, 11, 2011) (Article 1) with the inclusion of tubes with magnetic flux of 1020 Mx and 1021 Mx, and more simulations of the previously investigated case of 1022 Mx, sampling more convective flows than the previous article, greatly improving statistics. Observed properties of active regions are compared to properties of the simulated emerging flux tubes, including: the tilt of active regions in accordance with Joy’s Law as in Article 1, and in addition the scatter of tilt angles about the Joy’s Law trend, the most commonly occurring tilt angle, the rotation rate of the emerging loops with respect to the surrounding plasma, and the nature of the magnetic field at the flux tube apex. We discuss how these diagnostic properties constrain the initial field strength of the active-region flux tubes at the bottom of the solar convection zone, and suggest that flux tubes of initial magnetic field strengths of ≥?40 kG are good candidates for the progenitors of large (1021 Mx to 1022 Mx) solar active regions, which agrees with the results from Article 1 for flux tubes of 1022 Mx. With the addition of more magnetic flux values and more simulations, we find that for all magnetic field strengths, the emerging tubes show a positive Joy’s Law trend, and that this trend does not show a statistically significant dependence on the magnetic flux.  相似文献   

18.
We used an automatic image-processing method to detect solar-activity features observed in white light at the Kislovodsk Solar Station. This technique was applied to automatically or semi-automatically detect sunspots and active regions. The results of this automated recognition were verified with statistical data available from other observatories and revealed a high detection accuracy. We also provide parameters of sunspot areas, of the umbra, and of faculae as observed in Solar Cycle 23 as well as the magnetic flux of these active elements, calculated at the Kislovodsk Solar Station, together with white-light images and magnetograms from the Michaelson Doppler Imager onboard the Solar and Heliospheric Observatory (SOHO/MDI). The ratio of umbral and total sunspot areas during Solar Cycle 23 is ≈?0.19. The area of sunspots of the leading polarity was approximately 2.5 times the area of sunspots of the trailing polarity.  相似文献   

19.
D. V. Erofeev 《Solar physics》1996,167(1-2):25-45
Discrete rigidly rotating components (modes) of the large-scale solar magnetic field have been investigated. We have used a specially calculated basic set of functions to resolve the observed magnetic field into discrete components. This adaptive set of functions, as well as the expansion coefficients, have been found by processing a series of digitized synoptic maps of the background magnetic field over a 20-year period. As a result, dependences have been obtained which describe the spatial structure and the temporal evolution of the 27-day and 28-day rigidly rotating modes of the Sun's magnetic field.The spatial structure of the modes has been compared with simulations based on the known flux-transport equation. In the simulations, the rigidly rotating modes were regarded as stationary states of the magnetic field whose rigid rotation and stability were maintained by a balance between the emergence of magnetic flux from stationary sources located at low latitudes and the horizontal transport of flux by turbulent diffusion and poleward directed meridional flow. Under these assumptions, the structure of the modes is determined solely by the horizontal velocity field of the plasma, except for the low-latitude zone where sources of magnetic flux concentrate. We have found a detailed agreement between the simulations and the results of the data analysis, provided that the amplitude of the meridional flow velocity and the diffusion constant are equal to 9.5 m s–1 and 600 km2 s–1, respectively.The analysis of the expansion coefficients has shown that the rigidly rotating modes undergo rapid step-like variations which occur quasi-periodically with a period of about two years. These variations are caused by separate surges of magnetic flux in the photosphere, so that each new surge gives rise to a rapid replacement of old large-scale magnetic structures by newly arisen ones.  相似文献   

20.
We propose a new model for the magnetic field at different distances from the Sun during different phases of the solar cycle. The model depends on the observed large-scale non-polar (\({\pm}\, 55^{\circ }\)) photospheric magnetic field and on the magnetic field measured at polar regions from \(55^{\circ }\) N to \(90^{\circ }\) N and from \(55^{\circ }\) S to \(90^{\circ }\) S, which are the visible manifestations of cyclic changes in the toroidal and poloidal components of the global magnetic field of the Sun. The modeled magnetic field is determined as the superposition of the non-polar and polar photospheric magnetic field and considers cycle variations. The agreement between the model predictions and magnetic fields derived from direct in situ measurements at different distances from the Sun, obtained with different methods and at different solar activity phases, is quite satisfactory. From a comparison of the magnetic fields as observed and calculated from the model at 1 AU, we conclude that the model magnetic field variations adequately explain the main features of the interplanetary magnetic field (IMF) radial, \(B_{\mathrm{x}}\), component cycle evolution at Earth’s orbit. The modeled magnetic field averaged over a Carrington rotation (CR) correlates with the IMF \(B_{\mathrm{x}}\) component also averaged over a CR at Earth’s orbit with a coefficient of 0.691, while for seven CR-averaged data, the correlation reaches 0.81. The radial profiles of the modeled magnetic field are compared with those of already existing models. In contrast to existing models, ours provides realistic magnetic-field radial distributions over a wide range of heliospheric distances at different cycle phases, taking into account the cycle variations of the solar toroidal and poloidal magnetic fields. The model is a good approximation of the cycle behavior of the magnetic field in the heliosphere. In addition, the decrease in the non-polar and polar photospheric magnetic fields is shown. Furthermore, the magnetic field during solar cycle maxima and minima decreased from Cycle 21 to Cycle 24. This implies that both the toroidal and poloidal components, and therefore the solar global magnetic field, decreased from Cycle 21 to Cycle 24.  相似文献   

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

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