首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 95 毫秒
1.
The ionospheric total electron content (TEC) in the northern hemispheric equatorial ionization anomaly (EIA) crest region is investigated by using dual-frequency signals of the Global Positioning System (GPS) acquired from Rajkot (Geog. Lat. 22.29°N, Geog. Long. 70.74°E; Geom. Lat. 14.21°N, Geom. Long. 144.90°E), India. The day-to-day variability of EIA characteristics is examined during low solar activity period (F10.7∼83 sfu). It is found that the daily maximum TEC at EIA crest exhibits a day-to-day and strong semi-annual variability. The seasonal anomaly and equinoctial asymmetry in TEC at EIA is found non-existent and weaker, respectively. We found a moderate and positive correlation of daily magnitude of crest, Ic with daily F10.7 and EUV fluxes with a correlation coefficient of 0.43 and 0.33, respectively indicating an existence of a short-term relation between TEC at EIA and the solar radiation even during low solar activity period. The correlation of daily Ic with Dst index is also moderate (r = −0.35), whereas no correlation is found with the daily Kp index (r = 0.14) respectively. We found that the magnitude of EIA crest is moderately correlated with solar flux in all seasons except winter where it is weakly related (0.27). The magnitude of EIA crest is also found highly related with EEJ strength in spring (r = 0.69) and summer (r = 0.65) than autumn (0.5) and winter (r = 0.47), though EEJ is stronger in autumn than spring.  相似文献   

2.
We show that the higher range of the heliolongitudinal asymmetry of the solar wind speed in the positive polarity period (A > 0) than in the negative polarity period (A < 0) is one of the important reasons of the larger amplitudes of the 27-day variation of the galactic cosmic ray (GCR) intensity in the period of 1995–1997 (A > 0) than in 1985–1987 (A < 0). Subsequently, different ranges of the heliolongitudinal asymmetry of the solar wind speed jointly with equally important corresponding drift effect are general causes of the polarity dependence of the amplitudes of the 27-day variation of the GCR intensity. At the same time, we show that the polarity dependence is feeble for the last unusual minimum epoch of solar activity 2007–2009 (A < 0); the amplitude of the 27-day variation of the GCR intensity shows only a tendency of the polarity dependence. We present a three dimensional (3-D) model of the 27-day variation of GCR based on the Parker’s transport equation. In the 3-D model is implemented a longitudinal variation of the solar wind speed reproducing in situ measurements and corresponding divergence-free interplanetary magnetic field (IMF) derived from the Maxwell’s equations. We show that results of the proposed 3-D modeling of the 27-day variation of GCR intensity for different polarities of the solar magnetic cycle are in good agreement with the neutron monitors experimental data. To reach a compatibility of the theoretical modeling with observations for the last minimum epoch of solar activity 2007–2009 (A < 0) a parallel diffusion coefficient was increased by ∼40%.  相似文献   

3.
The ionospheric total electron content (TEC), derived by analyzing dual frequency signals from the Global Positioning System (GPS) recorded near the Indian equatorial anomaly region, Varanasi (geomagnetic latitude 14°, 55′N, geomagnetic longitude 154°E) is studied. Specifically, we studied monthly, seasonal and annual variations as well as solar and geomagnetic effects on the equatorial ionospheric anomaly (EIA) during the solar minimum period from May 2007 to April 2008. It is found that the daily maximum TEC near equatorial anomaly crest yield their maximum values during the equinox months and their minimum values during the summer. Using monthly averaged peak magnitude of TEC, a clear semiannual variation is seen with two maxima occurring in both spring and autumn. Statistical studies indicate that the variation of EIA crest in TEC is poorly correlated with Dst-index (r = −0.03) but correlated well with Kp-index (r = 0.82). The EIA crest in TEC is found to be more developed around 12:30 LT.  相似文献   

4.
The ionospheric variability at equatorial and low latitude region is known to be extreme as compared to mid latitude region. In this study the ionospheric total electron content (TEC), is derived by analyzing dual frequency Global Positioning System (GPS) data recorded at two stations separated by 325 km near the Indian equatorial anomaly region, Varanasi (Geog latitude 25°, 16/ N, longitude 82°, 59/ E, Geomagnetic latitude 16°, 08/ N) and Kanpur (Geog latitude 26°, 18/ N, longitude 80°, 12/ E, Geomagnetic latitude 17°, 18/ N). Specifically, we studied monthly, seasonal and annual variations as well as solar and geomagnetic effects on the equatorial ionospheric anomaly (EIA) during the descending phase of solar activity from 2005 to 2009. It is found that the maximum TEC (EIA) near equatorial anomaly crest yield their maximum values during the equinox months and their minimum values during the summer. Using monthly averaged peak magnitude of TEC, a clear semi-annual variation is seen with two maxima occurring in both spring and autumn. Results also showed the presence of winter anomaly or seasonal anomaly in the EIA crest throughout the period 2005–2009 only except during the deep solar minimum year 2007–2008. The correlation analysis indicate that the variation of EIA crest is more affected by solar activity compared to geomagnetic activity with maximum dependence on the solar EUV flux, which is attributed to direct link of EUV flux on the formation of ionosphere and main agent of the ionization. The statistical mean occurrence of EIA crest in TEC during the year from 2005 to 2009 is found to around 12:54 LT hour and at 21.12° N geographic latitude. The crest of EIA shifts towards lower latitudes and the rate of shift of the crest latitude during this period is found to be 0.87° N/per year. The comparison between IRI models with observation during this period has been made and comparison is poor with increasing solar activity with maximum difference during the year 2005.  相似文献   

5.
We have analyzed 101 Coronal Mass Ejection (CME) events and their associated interplanetary CMEs (ICMEs) and interplanetary (IP) shocks observed during the period 1997–2005 from the list given by Mujiber Rahman et al. (2012). The aim of the present work is to correlate the interplanetary parameters such as, the speeds of IP shocks and ICMEs, CME transit time and their relation with CME parameters near the Sun. Mainly, a group of 10 faster CME events (VINT > 2200 km/s) are compared with a list of 91 normal events of Manoharan et al. (2004). From the distribution diagrams of CME, ICME and IP shock speeds, we note that a large number of events tends to narrow towards the ambient (i.e., background) solar wind speed (∼500 km/s) in agreement with the literature. Also, we found that the IP shock speed and the average ICME speed measured at 1 AU are well correlated. In addition, the IP shock speed is found to be slightly higher than the ICME speed. While the normal events show CME travel time in the range of ∼40–80 h with a mean value of 65 h, the faster events have lower transit time with a mean value of 40 h. The effect of solar wind drag is studied using the correlation of CME acceleration with interplanetary (IP) acceleration and with other parameters of ICMEs. While the mean acceleration values of normal and faster CMEs in the LASCO FOV are 1 m/s2, 18 m/s2, they are −1.5 m/s2 and −14 m/s2 in the interplanetary medium, respectively. The relation between CME speed and IP acceleration for normal and faster events are found to agree with that of  and  except slight deviations for the faster events. It is also seen that the faster events with less travel time face higher negative acceleration (>−10 m/s2) in the interplanetary medium up to 1 AU.  相似文献   

6.
As an initial effort to study the evolution of the Venus atmosphere, the influence of the solar wind density and the interplanetary magnetic field (IMF) x component (the x-axis points from Venus towards the Sun) on the O+ ion escape rate from Venus is investigated using a three-dimensional quasi-neutral hybrid (HYB-Venus) model. The HYB-Venus model is first applied to a case of the high-density (100 cm−3) solar wind interaction with Venus selected from the Pioneer Venus Orbiter observations to demonstrate its capability for the study. Two sets of simulations with a wide range of solar wind densities and different IMF x components are then performed. It is found that the O+ ion escape rate increases with increasing solar wind density. The O+ ion escape rate saturates when the solar wind density becomes high (above 100 cm−3). The results also suggest that the IMF x component enhances the O+ ion escape rate, given a fixed IMF component perpendicular to the x-axis. Finally, the results imply a higher ion loss rate for early-Venus, when solar conditions were dramatically different.  相似文献   

7.
The 20th century temperature anomaly record is reproduced using an energy balance model, with a diffusive deep ocean. The model takes into account all the standard radiative forcings, and in addition the possibility of a non-thermal solar component. The model is parameterized and then optimized to produce the most likely values for the climate parameters and radiative forcings which reproduce the 20th century global warming. We find that the best fit is obtained with a negligible net feedback. We also show that a non-thermal solar component is necessarily present, indicating that the total solar contribution to the 20th century global warming, of ΔTsolar = 0.27 ± 0.07 °C, is much larger than can be expected from variation in the total solar irradiance alone. However, we also find that the largest contribution to the 20th century warming comes from anthropogenic sources, with ΔTman = 0.42 ± 0.11 °C.  相似文献   

8.
In this article, we study fast shocks at CIR boundaries during an extended interval of 15 consecutive major high speed solar wind streams in 1992–1993. Ulysses was 4–5 AU from the sun. The Abraham-Schrauner shock normal method and the Rankine-Hugoniot relations were used to determine fast shock directions and speeds. Out of 33 potential CIR shocks, 14 were determined to be fast forward shocks (FSs) and 14 were fast reverse shocks (RSs). Of the remaining 5 events, 2 were forward waves and 3 were reverse waves. CIR edges at latitudes below ∼30o were, for the most part, bounded by fast magnetosonic shocks. The forward shocks were generally quasi-perpendicular (average θnBo = 67o). The reverse shocks were more oblique (average θnBo = 52o), but they extended to all angles. Both FSs and RSs had magnetosonic Mach numbers ranging from 1 to 5 or 6. The average Mach numbers were 2.4 and 2.6 for FSs and RSs, respectively. The shock Mach numbers were noted to generally decrease with increasing latitude. The non-shock events or waves were noted to occur preferentially at high (∼−30° to −35°) heliolatitudes where stream-stream interactions were presumably weaker. These results are consistent with expectations, indicating the general accuracy of the Abraham-Schrauner technique.  相似文献   

9.
The aim of the study is to explore whether age at death from cardiovascular diseases depends on solar and geomagnetic activities. The data were collected for 1970–1978 in Novosibirsk, West Siberia, for industrial workers of Siberian origin. The Spearman correlations are computed between linearly detrended lifespan and daily or monthly physical variables to establish immediate (lag, L = 0), delayed (L = 1–3 days) and cumulative (L = ±30 days) influences. Significant correlations ranging from r = −0.26 to r = −0.30 for L from 0 to 3, respectively, are found for men between solar radio flux at wavelength 10.7 cm and age at death from acute myocardial infarction (AMI) but not from acute heart failure, ischemic heart disease and stroke. For AMI, women’s longevity displays an opposite (direct) association with the average solar character occurred at the calendar month of death. The index of geomagnetic activity, Ap, exhibits inverse association with longevity for the AMI stratum for both sexes. GLM univariate procedure revealed higher contribution of Ap to the variance of lifespan compared to season of death. The individual age at death susceptibility to cosmic influences is found to depend upon solar activity at year of birth. It is concluded that associations between the lifespan for cardiovascular decedents and the indices of solar and geomagnetic activities at time of death and of birth are cause-of-death- and sex-specific.  相似文献   

10.
Ionogram observations from the ionosonde at Fuke (9.5°N geomagnetic latitude), a Chinese low latitude station, in 2010–2012 are analyzed to present the features of F3 layer under low and moderate solar activity conditions. Structure of the ionogram, displaying the F3 layer, was more distinct and clear during MSA than LSA periods especially during spring to summer. Start time of occurrence of the F3 layer is about at 0830–0900 LT and is approximately the same for LSA and MSA conditions. The average duration time of the F3 layer occurrence was 181 min per day under F10.7 = 75 condition, 263 min in F10.7 = 99 and 358 min in F10.7 = 125, respectively. The differences of h′F2 and h′F3 exhibited obvious semiannual variation observed at Fuke from March 2010 to June 2012 and increased with increasing solar activity. The difference of foF2 and foF3 in the months February, March, September, October and November is less evident in the middle solar activity period 2011–2012 than the low solar activity 2010 and in the other period it shows a slight increase (0.5 MHz) or keeps constant. The results show that the solar activity dependence of the F3 layer occurrence at low latitude away from the magnetic equator is different from that at near the magnetic equator.  相似文献   

11.
Observations made with the two Voyager spacecraft confirmed that the solar wind decelerates to form the heliospheric termination shock and that it has begun its merger with the local interstellar medium. The compression ratio of this shock affects galactic cosmic rays when they enter the heliosphere. Hydrodynamic (HD) models show that the compression ratio can have a significant latitude dependence; with the largest value in the nose direction of the heliosphere, becoming significantly less towards the polar regions. The modulation effects of such large latitude dependence are studied, using a well-established numerical drift and shock modulation model. We focus on computing the modulated spectra for galactic protons with emphasis on the radial and polar gradients in the equatorial plane and at a polar angle of θ = 55°, corresponding to the heliolatitude of Voyager 1. Two sets of solutions are computed and compared each time; with and without a latitude dependence for the compression ratio. All computations are done for the two magnetic field polarity cycles assuming solar minimum conditions. Including the termination shock in the model allows the study of the re-acceleration of galactic protons in the outer heliosphere. We find that for the A < 0 polarity cycle the intensity between ∼200 MeV and ∼1 GeV in the vicinity of the shock in the heliospheric equatorial plane may exceed the local interstellar value specified at the heliopause. Unfortunately, at θ = 55°, the effect is reduced. This seems not possible during an A > 0 cycle because significant modulation is then predicted between the heliopause and the termination shock, depending on how strong global gradient and curvature drifts are in the heliosheath. The overall effect of the shock on galactic protons in the equatorial plane is to reduce the total modulation as a function of radial distance with respect to the interstellar spectrum. Making the compression ratio latitude dependent enhances these effects at energies E < 200 MeV in the equatorial plane. At larger heliolatitudes these effects are even more significant. The differences in the modulation between the two drift cycles are compelling when the compression ratio is made latitude dependent but at Earth this effect is insignificant. A general result is that the computed radial gradient changes for galactic protons at and close to the TS and that these changes are polarity dependent. In line with previous work, large polarity dependent effects are predicted for the inner heliosphere and also close to the shock’s position in the equatorial plane. In contrast, at θ = 55°, the largest polarity effect occurs in the middle heliosphere (50 AU), enhanced by the latitude dependence of the compression ratio. At this latitude, the amount of proton modulation between the heliopause and the termination shock is much reduced. If galactic cosmic rays were to experience some diffusive shock acceleration over the 100–1000 MeV range at the shock, the radial gradient should change its sign in the vicinity of the shock, how large, depends on the compression ratio and the amount of drifts taking place in the outer heliosphere. The effective polar gradient shows a strong polarity dependence at Earth but this dissipates at θ = 55°, especially with increasing radial distance. This tendency is enhanced by making the compression ratio latitude dependent.  相似文献   

12.
A new set of data obtained at low solar activity from Ilorin, Nigeria (geog. latitude 8.5°N, geog longitude, 4.6°E, dip 4.1°S) is used to validate the IRI 2001 model at low solar activity. The results show in general a good agreement between model and observed B0 at night but an over estimation during daytime. The overestimation is greatest during the morning period (0600LT–1000LT). The model prediction for B1 is fairly good at night and during the day. A dependence of B0 on solar zenith angle χ is observed during the daytime. A formulation of the form B0 = A[cos(χ)n] is therefore proposed. Values of the constants n and A were determined for the period of low solar activity for this station.  相似文献   

13.
We use hourly monthly median values of propagation factor M(3000)F2 data observed at Ouagadougou Ionospheric Observatory (geographic12.4°N, 1.5°W; 5.9o dip), Burkina Faso (West Africa) during the years Januar1987–December1988 (average F10.7 < 130 × 10−22 W/m2/Hz, representative of low solar flux conditions) and for January 1989–December1990 (average F10.7 ? 130 × 10−22 W/m2/Hz, representative of high solar epoch) for magnetically quiet conditions to describe local time, seasonal and solar cycle variations of equatorial ionospheric propagation factor M(3000)F2 in the African region. We show that that seasonal trend between solar maximum and solar minimum curves display simple patterns for all seasons and exhibits reasonable disparity with root mean square error (RMSE) of about 0.31, 0.29 and 0.26 for December solstice, June solstice and equinox, respectively. Variability Σ defined by the percentage ratio of the absolute standard deviation to the mean indicates significant dissimilarity for the two solar flux levels. Solar maximum day (10–14 LT) and night (22–02 LT) values show considerable variations than the solar minimum day and night values. We compare our observations with those of the IRI 2007 to validate the prediction capacity of the empirical model. We find that the IRI model tends to underestimate and overestimate the observed values of M(3000)F2, in particular, during June solstice season. There are large discrepancies, mainly during high solar flux equinox and December solstice between dawn and local midnight. On the other hand, IRI provides a slightly better predictions for M(3000)F2 between 0900 and 1500 LT during equinox low and high solar activity and equinox high sunspot number. Our data are of great importance in the area of short-wave telecommunication and ionospheric modeling.  相似文献   

14.
A multi-station analysis of geomagnetic field measurements conducted for a remarkable case event, at a European and an American array shows that, although several aspects of the geomagnetic field observations show a clear latitudinal and local time dependence, simultaneous oscillations at discrete frequencies (f ≈ 1.0, 1.3, 2.2 and 3.2 mHz) are almost ubiquitously detected, from low to high latitudes both in the light and dark sector. They are driven by fluctuations of the solar wind density and dynamic pressure at the same frequencies, via the modulation of the magnetopause current. We also report clear evidence for the occurrence of resonant coupling (at f ≈ 2.2 and 3.2 mHz) between such modes and high latitude field lines. Due to the variable length of the field line through the day, oscillation modes at the same frequencies resonate at different latitudes in the daytime and nighttime region, respectively.  相似文献   

15.
Response of the D-region of the ionosphere to the total solar eclipse of 22 July 2009 at low latitude, Varanasi (Geog. lat., 25.27° N; Geog. long., 82.98° E; Geomag. lat. = 14° 55’ N) was investigated using ELF/VLF radio signal. Tweeks, a naturally occurring VLF signal and radio signals from various VLF navigational transmitters are first time used simultaneously to study the effect of total solar eclipse (TSE). Tweeks occurrence is a nighttime phenomena but the obscuration of solar disc during TSE in early morning leads to tweek occurrence. The changes in D-region ionospheric VLF reflection heights (h) and electron density (ne: 22.6–24.6 cm−3) during eclipse have been estimated from tweek analysis. The reflection height increased from ∼89 km from the first occurrence of tweek to about ∼93 km at the totality and then decreased to ∼88 km at the end of the eclipse, suggesting significant increase in tweek reflection height of about 5.5 km during the eclipse. The reflection heights at the time of totality during TSE are found to be less by 2–3 km as compared to the usual nighttime tweek reflection heights. This is due to partial nighttime condition created by TSE. A significant increase of 3 dB in the strength of the amplitude of VLF signal of 22.2 kHz transmitted from JJI-Japan is observed around the time of the total solar eclipse (TSE) as compared to a normal day. The modeled electron density height profile of the lower ionosphere depicts linear variation in the electron density with respect to solar radiation as observed by tweek analysis also. These low latitude ionospheric perturbations on the eclipse day are discussed and compared with other normal days.  相似文献   

16.
A numerical model, based on Parker’s transport equation, describing the modulation of anomalous cosmic rays and containing diffusive shock acceleration is applied. The role of radial perpendicular diffusion at the solar wind termination shock, and as the dominant diffusion coefficient in the outer heliosphere, is studied, in particular the role it plays in the effectiveness of the acceleration of anomalous protons and helium when its latitude dependence is changed. It is found that the latitudinal enhancement of radial perpendicular diffusion towards the heliospheric poles and along the termination shock has a prominent effect on the acceleration of these particles. It results in a ‘break’ in the energy spectrum for anomalous protons at ∼6.0 MeV, causing the spectral index to change from E−1.38 to E−2.23, but for anomalous helium at ∼3.0 MeV, changing the spectral index from E−1.38 to E−2.30. When approaching the simulated TS, the changes in the modulated spectra as they unfold to a ‘steady’ power law shape at energies below 50 MeV are much less prominent as a function of radial distances when radial perpendicular diffusion is increased with heliolatitude.  相似文献   

17.
This work is a continuation of the previous article and it focuses on low solar activity and modeling effort. NeQuick model uses Epstein layer formalism to model each part of the profile. We study the diurnal and seasonal variations of B2bot, ΔB2 (B2best − B2NeQuick2) and R (B2best/B2NeQuick 2) at Hainan station during low solar activity. The results show it is possible to improve the B2bot parameter of the NeQuick model at that region during low solar activity. Then, we use a function ?(t) with LT in different seasons to correct the B2bot formula of NeQuick 2. The correction shows that (1) By the correction formula, the B2bot of NeQuick is improved. The maximum standard deviation is improved for 9 km. (2) The correction formula is more effective in summer than in equinox and winter and performs better during early morning hours than during the rest of the day.  相似文献   

18.
A current serious limitation on the studies of solar energetic particle (SEP) events is that their properties in the inner heliosphere are studied only through in situ spacecraft observations. Our understanding of spatial distributions and temporal variations of SEP events has come through statistical studies of many such events over several solar cycles. In contrast, flare SEPs in the solar corona can be imaged through their radiative and collisional interactions with solar fields and particles. We suggest that the heliospheric SEPs may also interact with heliospheric particles and fields to produce signatures which can be remotely observed and imaged. A challenge with any such candidate signature is to separate it from that of flare SEPs. The optimum case for imaging high-energy (E > 100 MeV) heliospheric protons may be the emission of π0-decay γ-rays following proton collisions with solar wind (SW) ions. In the case of E > 1 MeV electrons, gyrosynchrotron radio emission may be the most readily detectible remote signal. In both cases we may already have observed one or two such events. Another radiative signature from nonthermal particles may be resonant transition radiation, which has likely already been observed from solar flare electrons. We discuss energetic neutrons as another possible remote signature, but we rule out γ-ray line and 0.511 MeV positron annihilation emission as observable signatures of heliospheric energetic ions. We are already acquiring global signatures of large inner-heliospheric SW density features and of heliosheath interactions between the SW and interstellar neutral ions. By finding an appropriate observable signature of remote heliospheric SEPs, we could supplement the in situ observations with global maps of energetic SEP events to provide a comprehensive view of SEP events.  相似文献   

19.
The polarization pattern of ULF pulsations (f ≈ 1–100 mHz) at Terra Nova Bay (Antarctica, CGM λ ∼ 80°) has been determined for the entire 2003, soon after the solar maximum. A comparison with the results of previous investigations, conducted at the same station close to the solar minimum (1994–96), allows to focus common elements and major differences among different frequency bands which persist through the entire solar cycle. Basically, between f ∼ 1.5 and 5 mHz, the day can be divided into four sectors with alternate polarizations. The local time and latitudinal dependence of the observed pattern can be tentatively interpreted in terms of a latitude of resonant field lines reaching λ ∼ 80° in the noon sector; on the other hand, resonance effects of lower latitude field lines can be clearly identified also far from the noon meridian when the station moves into the deep polar cap. Moreover, in the morning sector the resonance region would extend to lower latitudes than in the evening sector. The proposed profile of the resonant region can interpret also the results obtained at other cusp/auroral stations and appears consistent with that one inferred in the northern hemisphere at smaller latitudes. The resonance region progressively shifts toward lower latitude with increasing frequency; correspondingly, the four-sector pattern progressively disappears at TNB. Above f ∼ 20 mHz, the experimental observations might suggest an additional contribution from Sunward propagating waves, possibly via the magnetotail lobes.  相似文献   

20.
We demonstrate that the general features of the radial and azimuthal components of the anisotropy of galactic cosmic rays can be studied by the harmonic analysis method using data from an individual neutron monitor with cut off rigidity <5 GV. In particular, we study the characteristics of the 27-day (solar rotation period) variations of the galactic cosmic ray intensity and anisotropy, solar wind velocity, interplanetary magnetic field strength and sunspot number. The amplitudes of the 27-day variations of the galactic cosmic ray anisotropy are greater, and the phases more clearly established, in A > 0 polarity periods than in A < 0 polarity periods at times of minimum solar activity. The phases of the 27-day variations of the galactic cosmic rays intensity and anisotropy are opposite with respect to the similar changes of the solar wind velocity in A > 0 polarity periods. No significant dependence of the amplitude of the 27-day variation of the galactic cosmic ray anisotropy on the tilt angle of the heliospheric neutral sheet is found. Daily epicyclegrams obtained by Chree’s method show that the 27-day variations of the galactic cosmic ray anisotropy during A > 0 polarity periods follow elliptical paths with the major axes oriented approximately along the interplanetary magnetic field. The paths are more irregular during A < 0 polarity periods.  相似文献   

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

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