首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
Topside sounding electron density profiles are analyzed to explore interrelations of the F2 layer critical frequency and the peak height for a representative set of conditions provided by ISIS1, ISIS2, IK19 and Cosmos-1809 satellites for the period of 1969–1987. The foF2 and hmF2 are delivered with exponential extrapolation of electron density profile to zero of its 1st derivative. It is shown that the linear regression exists between foF2 and hmF2 under different conditions. The linkage between the two parameters amended to the empirical model of the peak height [Gulyaeva, T.L., Bradley, P.A., Stanislawska, I., Juchnikowski, G. Towards a new reference model of hmF2 for IRI. Adv. Space Res. 42, 666–672, doi:10.1016/j.asr.2008.02.021, 2008] results in an empirical model of the both foF2 and hmF2 expressed by superposition of functions in terms of local-time, season, geodetic longitude, modified dip latitude and solar activity. For the solar activity we use a proxy Fsp index averaged from the mean solar radio flux F10.7s for the past 81 days (3 solar rotations) and F10.7 value for 1 day prior the day of observation. Impact of geomagnetic activity is not discernible with the topside sounding data due to mixed positive and negative storm-time effects. Appreciable differences have been revealed between IRI-CCIR predictions and outcome of the new model which might be attributed to the different techniques of the peak electron density and height derivation, different epochs and different global distribution of the source data as well as the different mathematical functions involved in the maps and the model presentation.  相似文献   

2.
Accuracy of IRI electron density profile depends on the F2 layer peak density and height converted by empirical formulae from the critical frequency and M3000F2 factor provided by the ITU-R (former CCIR). The CCIR/ITU-R maps generated from ground-based ionosonde measurements suffer from model assumptions, in particular, over the oceans where relatively few measurements are available due to a scarcity of ground-based ionosondes. In the present study a grid-point calibration of IRI/ITU-R maps for the foF2 and hmF2 over the oceans is proposed using modeling results based on the topside true-height profiles provided by ISIS1, ISIS2, IK-19 and Cosmos-1809 satellites for the period of 1969–1987. Topside soundings results are compared with IRI and the Russian standard model of ionosphere, SMI, and grouped to provide an empirical calibration coefficient to the peak density and height generated from ITU-R maps. The grid-point calibration coefficients maps are produced in terms of the solar activity, geodetic latitude and longitude, universal time and season allowing update of IRI–ITU-R predictions of the F2 layer peak parameters.  相似文献   

3.
To improve the accuracy of the real time topside electron density profiles given by the Digisonde software a new model-assisted technique is used. This technique uses the Topside Sounder Model (TSM), which provides the plasma scale height (Hs), O+–H+ transition height (HT), and their ratio Rt = Hs/HT, derived from topside sounder data of Alouette and ISIS satellites. The Topside Sounder Model Profiler (TSMP) incorporates TSM and uses the model quantities as anchor points in construction of topside density (Ne) profiles. For any particular location, TSMP calculates topside Ne profiles by specifying the values of foF2 and hmF2. In the present version, TSMP takes the F2 peak characteristics – foF2, hmF2, and the scale height at hmF2 – from the Digisonde measurements. The paper shows results for the Digisonde stations Athens and Juliusruh. It is found that the topside scale height used in Digisonde reconstruction is less than that extracted from topside sounder profiles. Rough comparison of their bulk distributions showed that they differ by an average factor of 1.25 for locations of Athens and Juliusruh. When the Digisonde scale heights are adjusted by this factor, the reconstructed topside profiles are close to those provided by TSM. Compared with CHAMP reconstruction profiles in two cases, TSMP/Digisonde profiles show lower density between 400 and 2000 km.  相似文献   

4.
This paper presents the observed ionospheric F-region critical frequency, foF2, and peak height, hmF2, at northern crest of equatorial ionization anomaly (EIA) area station, namely Chung-Li (24.9°N, 121.1°E, dip 35°), and to be compared with International Reference Ionosphere model (IRI-2001) predictions for the period from 1994 to 1999, corresponding to half of the 23rd solar cycle. The diurnal and seasonal variation of foF2 and hmF2 are analyzed for different solar phases, respectively. The result shows the largest discrepancies were observed during nighttime for foF2 and hmF2, respectively. The value of foF2 both CCIR and URSI selected in the IRI model produced a good agreement during the daytime and underestimated during the noon time for high solar activities. The underestimation at noon time is mainly caused by the fountain effect from equator. Further, the peak height hmF2 shows a larger variability around the midnight than daytime in the equinox and winter seasons and reserved in summer, respectively. The study shows that the monthly median values of observed hmF2 is somewhat lower than those predicated by the IRI model, at night time in all the seasons except the period of 04:00–06:00 LT and reverse at daytime in summer. In general the IRI model predictions with respect to the observed in hmF2 is much better than foF2. The percentage deviation of the observed foF2 (hmF2) values with respect to the IRI model varies from 5% to 80% (0–25%) during nighttime and 2–17% (0–20%) at daytime, respectively. In general, the model generates good results, although some improvements are still necessary to implement in order to obtain better simulations for ionospheric low-latitudes region.  相似文献   

5.
We describe a Parameterized Regional Ionospheric Model (PARIM) to calculate the spatial and temporal variations of the ionospheric electron density/plasma frequency over the Brazilian sector. The ionospheric plasma frequency values as calculated from an enhanced Sheffield University Plasmasphere–Ionosphere Model (SUPIM) were used to construct the model. PARIM is a time-independent 3D regional model (altitude, longitude/local time, latitude) used to reproduce SUPIM plasma frequencies for geomagnetic quiet condition, for any day of the year and for low to moderately high solar activity. The procedure to obtain the modeled representation uses finite Fourier series so that all plasma frequency dependencies can be represented by Fourier coefficients. PARIM presents very good results, except for the F region peak height (hmF2) near the geomagnetic equator during times of occurrence of the F3 layer. The plasma frequency calculated by IRI from E region to bottomside of the F region present latitudinal discontinuities during morning and evening times for both solar minimum and solar maximum conditions. Both the results of PARIM and the IRI for the E region peak density show excellent agreement with the observational values obtained during the conjugate point equatorial experiment (COPEX) campaign. The IRI representations significantly underestimate the foF2 and hmF2 compared to the observational results over the COPEX sites, mainly during the evening–nighttime period.  相似文献   

6.
In this paper, the F2-layer critical frequency (foF2) and peak height (hmF2) measured by the FM/CW ionosonde at Thailand equatorial latitude station, namely Chumphon (10.72°N, 99.37°E, dip 3.22) are presented. The measurement data during low solar activity from January 2004 to December 2006 are analyzed based on the diurnal, seasonal variation. The results are then compared with IRI-2001 model predictions. Our study shows that: (1) In general, both the URSI and CCIR options of the IRI model give foF2 close to the measured ones, but the CCIR option produces a smaller range of deviation than the URSI option. The agreement during daytime is generally better than during nighttime. Overestimation mostly occurs in 2004 and 2006, while underestimation is during pre-sunrise hours in June solstice in 2005. The peak foF2 around sunset is higher during March equinox and September equinox than the other seasons, with longer duration of maximum levels in March equinox than September equinox. Large coefficients of variability foF2 occur during pre-sunrise hours. Meanwhile, the best agreement between the observed foF2 and the IRI model is obtained in June solstice. (2) In general, The IRI (CCIR) model predicts the observed hmF2 well during daytime in June solstice from 2004–2006, but it overestimates during March equinox, September equinox and December solstice. For nighttime, the model overestimates hmF2 values for all seasons especially during March equinox and September equinox. However, the model underestimates hmF2 values during September equinox and for some cases during June solstice and December solstice at pre-sunrise. The agreement between the IRI model and the hmF2(M3000OBS) is worst around noontime, post-sunset and pre-sunrise hours. All comparative studies give feedback for new improvements of CCIR and URSI IRI models.  相似文献   

7.
The paper is focused on ionospheric response to occasional magnetic disturbances above selected ionospheric stations located at middle latitudes of the Northern and Southern Hemisphere under extremely low solar activity conditions of 2007–2009. We analyzed changes in the F2 layer critical frequency foF2 and the F2 layer peak height hmF2 against 27-days running mean obtained for different longitudinal sectors of both hemispheres for the initial, main and recovery phases of selected magnetic disturbances. Our analysis showed that the effects on the middle latitude ionosphere of weak-to-moderate CIR-related magnetic storms, which mostly occur around solar minimum period, could be comparable with the effects of strong magnetic storms. In general, both positive and negative deviations of foF2 and hmF2 have been observed independent on season and location. However positive effects on foF2 prevailed and were more significant. Observations of stormy ionosphere also showed large departures from the climatology within storm recovery phase, which are comparable with those usually observed during the storm main phase. The IRI STORM model gave no reliable corrections of foF2 for analyzed events.  相似文献   

8.
Monthly median values of foF2, hmF2 and M(3000)F2 parameters, with quarter-hourly time interval resolution for the diurnal variation, obtained with DPS4 digisonde at Hainan (19.5°N, 109.1°E; Geomagnetic coordinates: 178.95°E, 8.1°N) are used to investigate the low-latitude ionospheric variations and comparisons with the International Reference Ionosphere (IRI) model predictions. The data used for the present study covers the period from February 2002 to April 2007, which is characterized by a wide range of solar activity, ranging from high solar activity (2002) to low solar activity (2007). The results show that (1) Generally, IRI predictions follow well the diurnal and seasonal variation patterns of the experimental values of foF2, especially in the summer of 2002. However, there are systematic deviation between experimental values and IRI predictions with either CCIR or URSI coefficients. Generally IRI model greatly underestimate the values of foF2 from about noon to sunrise of next day, especially in the afternoon, and slightly overestimate them from sunrise to about noon. It seems that there are bigger deviations between IRI Model predictions and the experimental observations for the moderate solar activity. (2) Generally the IRI-predicted hmF2 values using CCIR M(3000)F2 option shows a poor agreement with the experimental results, but there is a relatively good agreement in summer at low solar activity. The deviation between the IRI-predicted hmF2 using CCIR M(3000)F2 and observed hmF2 is bigger from noon to sunset and around sunrise especially at high solar activity. The occurrence time of hmF2 peak (about 1200 LT) of the IRI model predictions is earlier than that of observations (around 1500 LT). The agreement between the IRI hmF2 obtained with the measured M(3000)F2 and the observed hmF2 is very good except that IRI overestimates slightly hmF2 in the daytime in summer at high solar activity and underestimates it in the nighttime with lower values near sunrise at low solar activity.  相似文献   

9.
Diurnal and seasonal variations of critical frequency of ionospheric F2-region ‘foF2’ and the height of peak density ‘hmF2’ are studied using modern digital ionosonde observations of equatorial ionization anomaly (EIA) crest region, Bhopal (23.2°N, 77.6°E, dip 18.5°N), during solar minimum period 2007. Median values of these parameters are obtained at each hour using manually scaled data during different seasons and compared with the International Reference Ionosphere-2001 model predictions. The observations suggest that on seasonal basis, the highest values of foF2 are observed during equinox months, whereas highest values of hmF2 are obtained in summer and lowest values of both foF2 and hmF2 are observed during winter. The observed median and IRI predicted values of foF2 and hmF2 are analyzed with upper and lower bound of inter-quartile range (IQR) and it is find out that the observed median values are well inside the inter-quartile range during the period of 2007. Comparison of the recorded foF2 and hmF2 values with the IRI-2001 output reveals that IRI predicted values exhibit better agreement with hmF2 as compared to foF2. In general, the IRI model predictions show some agreement with the observations during the year 2007. Therefore it is still necessary to implement improvements in order to obtain better predictions for EIA regions.  相似文献   

10.
The comparison of the IRI model with the foF2 distribution in the equatorial anomaly region obtained by topside sounding onboard the Interkosmos-19 satellite has been carried out. The global distribution of foF2 in terms of LT-maps was constructed by averaging Intercosmos-19 data for summer, winter, and equinox. These maps, in fact, represent an empirical model of the equatorial anomaly for high solar activity F10.7 ~ 200. The comparison is carried out for the latitudinal foF2 profiles in the characteristic longitudinal sectors of 30, 90, 210, 270, and 330°, as well as for the longitudinal variations in foF2 over the equator. The largest difference between the models (up to 60%) for any season was found in the Pacific longitudinal sector of 210°, where there are a few ground-based sounding stations. Considerable discrepancies, however, are sometimes observed in the longitudinal sectors, where there are many ground-based stations, for example, in the European or Indian sector. The discrepancies reach their maximum at 00 LT, since a decay of the equatorial anomaly begins before midnight in the IRI model and after midnight according to the Interkosmos-19 data. The discrepancies are also large in the morning at 06 LT, since in the IRI model, the foF2 growth begins long before sunrise. In the longitudinal variations in foF2 over the equator at noon, according to the satellite data, four harmonics are distinguished in the June solstice and at the equinox, and three harmonics in the December solstice, while in the IRI model only two and one harmonics respectively are revealed. In diurnal variations in foF2 and, accordingly, in the equatorial anomaly intensity, the IRI model does not adequately reproduce even the main, evening extremum.  相似文献   

11.
A method is proposed for reconstructing the electron density profiles N(h) of the IRI model from ionograms of topside satellite sounding of the ionosphere. An ionograms feature is the presence of traces of signal reflection from the Earth's surface. The profile reconstruction is carried out in two stages. At the first stage, the N(h) –profile is calculated from the lower boundary of the ionosphere to the satellite height (total profile) by the method presented in this paper using the ionogram. In this case, the monotonic profile of the topside ionosphere is calculated by the classical method. The profile of the inner ionosphere is represented by analytical functions, the parameters of which are calculated by optimization methods using traces of signal reflection, both from the topside ionosphere and from the Earth. At the second stage, the profile calculated from the ionogram is used to obtain the key parameters: the height of the maximum hmF2 of the F2 layer, the critical frequency foF2, the values of B0 and B1, which determine the profile shape in the F region in the IRI model. The input of key parameters, time of observation, and coordinates of sounding into the IRI model allows obtaining the IRI-profile corrected to real experimental conditions. The results of using the data of the ISIS-2 satellite show that the profiles calculated from the ionograms and the IRI profiles corrected from them are close to each other in the inner ionosphere and can differ significantly in the topside ionosphere. This indicates the possibility of obtaining a profile in the inner ionosphere close to the real distribution, which can significantly expand the information database useful for the IRTAM (IRI Realmax Assimilative Modeling) model. The calculated profiles can be used independently for local ionospheric research.  相似文献   

12.
In this work, the foF2 and hmF2 parameters at the conjugate points near the magnetic equator of Southeast Asia are studied and compared with the International Reference Ionosphere (IRI) model. Three ionosondes are installed nearly along the magnetic meridian of 100°E; one at the magnetic equator, namely Chumphon (10.72°N, 99.37°E, dip angle 3.0°N), and the other two at the magnetic conjugate points, namely Chiang Mai (18.76°N, 98.93°E, dip angle 12.7°N) and Kototabang (0.2°S, 100.30°E, dip angle 10.1°S). The monthly hourly medians of the foF2 and hmF2 parameters are calculated and compared with the predictions obtained from the IRI-2007 model from January 2004 to February 2007. Our results show that: the variations of foF2 and hmF2 predicted by the IRI-2007 model generally show the similar feature to the observed data. Both parameters generally show better agreement with the IRI predictions during daytime than during nighttime. For foF2, most of the results show that the IRI model overestimates the observed foF2 at the magnetic equator (Chumphon), underestimates at the northern crest (Chiang Mai) and is close to the measured ones at the southern crest of the EIA (Kototabang). For hmF2, the predicted hmF2 values are close to the hmF2(M3000F2OBS) during daytime. During nighttime, the IRI model gives the underestimation at the magnetic equator and the overestimation at both EIA crests. The results are important for the future improvements of the IRI model for foF2 and hmF2 over Southeast Asia region.  相似文献   

13.
We describe a new version of the Parameterized Regional Ionospheric Model (PARIM) which has been modified to include the longitudinal dependences. This model has been reconstructed using multidimensional Fourier series. To validate PARIM results, the South America maps of critical frequencies for the E (foE) and F (foF2) regions were compared with the values calculated by Sheffield Plasmasphere-Ionosphere Model (SUPIM) and IRI representations. PARIM presents very good results, the general characteristics of both regions, mainly the presence of the equatorial ionization anomaly, were well reproduced for equinoctial conditions of solar minimum and maximum. The values of foF2 and hmF2 recorded over Jicamarca (12°S; 77°W; dip lat. 1°N; mag. declination 0.3°) and sites of the conjugate point equatorial experiment (COPEX) campaign Boa Vista (2.8°N; 60.7°W; dip lat. 11.4°; mag. declination −13.1°), Cachimbo (9.5°S; 54.8°W; dip lat. −1.8°; mag. declination −15.5°), and Campo Grande (20.4°S; 54.6°W; dip lat. −11.1°; mag. declination −14.0°) have been used in this work. foF2 calculated by PARIM show good agreement with the observations, except during morning over Boa Vista and midnight-morning over Campo Grande. Some discrepancies were also found for the F-region peak height (hmF2) near the geomagnetic equator during times of F3 layer occurrences. IRI has underestimated both foF2 and hmF2 over equatorial and low latitude sectors during evening-nighttimes, except for Jicamarca where foF2 values were overestimated.  相似文献   

14.
The ionospheric characteristics of the F2 layer peak have been measured with ionosondes from the ground or with satellites from space. The most common characteristics are the F2-peak density NmF2 and peak height hmF2. In addition to these two parameters this paper studies the F2-peak scale height. Comparing the median values of hmF2 and NmF2 obtained from topside and bottomside sounding shows good agreement in general. The Chapman scale height values for the F2 layer peak derived from topside profiles, Hm,top, are generally several times larger than Hm,bot derived from bottomside profiles.  相似文献   

15.
The paper describes the technique that has been implemented to model the electron density distribution above and below the F2 peak making use of only the profiles obtained from the INTERCOSMOS-19 topside ionograms. Each single profile from the satellite height to the ionosphere peak has been fitted by a semi-Epstein layer function of the type used in the DGR model with shape factor variable with altitude. The topside above the satellite height has been extrapolated to match given values of plasmaspheric electron densities to obtain the full topside profile. The bottomside electron density has been calculated by using the maximum electron density and its altitude estimated from the topside ionogram as input for a modified version of the DGR derived profiler that uses model values for the foF1 and foE layers of the ionosphere. Total electron content has also been calculated. Longitudinal cross sections of vertical profiles from latitudes 50° N to 50° S latitude are shown for low and high geomagnetic activity. These cross sections indicate the equatorial anomaly effect and the changes of the shape of low latitude topside ionosphere during geomagnetic active periods. These results and the potentiality of the technique are discussed.  相似文献   

16.
Monthly median values of foF2, hmF2 and M(3000)F2 parameters, with hourly time interval resolution for the diurnal variation, obtained with DPS-4 digisonde observations at Hainan (19.4°N, 109.0°E) are used to study the low latitude ionospheric variation behavior. The observational results are compared with the International Reference Ionospheric Model (IRI) predictions. The time period coverage of the data used for the present study is from March 2002 to February 2005. Our present study showed that: (1) In general, IRI predictions using CCIR and URSI coefficients follow well the diurnal and seasonal variation patterns of the experimental values of foF2. However, CCIR foF2 and URSI foF2 IRI predictions systematically underestimate the observed results during most time period of the day, with the percentage difference ΔfoF2 (%) values changing between about −5% and −25%, whereas for a few hours around pre-sunrise, the IRI predictions generally overestimate the observational ones with ΔfoF2 (%) sometimes reaching as large as ∼30%. The agreement between the IRI results and the observational ones is better for the year 2002 than for the other years. The best agreement between the IRI results and the observational ones is obtained in summer when using URSI coefficients, with the seasonal average values of ΔfoF2 (%) being within the limits of ±10%. (2) In general, the IRI predicted hmF2 values using CCIR M(3000)F2 option shows a poor agreement with the observational results. However, when using the measured M(3000)F2 as input, the diurnal variation pattern of hmF2 given by IRI2001 has a much better agreement with the observational one with the detailed fine structures including the pre-sunrise and post-sunset peaks reproduced reasonably well. The agreement between the IRI predicted hmF2 values using CCIR M(30,000)F2 option and the observational ones is worst for the afternoon to post-midnight hours for the high solar activity year 2002. During daytime hours the agreement between the hmF2 values obtained with CCIR M(30,000)F2 option and the observational ones is best for summer season. The discrepancy between the observational hmF2 and that obtained with CCIR M(30,000)F2 option stem from the CCIR M(3000)F2 model, which does not produce the small scale structures observed in the measured M(3000)F2.  相似文献   

17.
The behaviour of the equatorial/low-latitude ionosphere and the transport processes during magnetic disturbed and quiet periods of a high solar activity year, 2014, in the American sector are investigated. Parameters used include vertical drift (Vz), transport term (W), NmF2, hmF2 and scale-height (H). The F2 plasma variations followed the diurnal local solar pattern, being higher at daytime. The sunset maximum and sunrise minimum peaks of hmF2 were directly opposite to the scale height (H) pattern. The plasma distribution was basically controlled by combined actions of the electrodynamic convection/thermospheric composition, which is geomagnetic activity dependent. The annual, semi-annual and winter-anomalies of the F2 parameters were higher at the dip equator in comparison with the low-latitude. The Vz pre-reversal peak magnitude coincided with hmF2 peak and the effects are more pronounced during geomagnetic disturbed conditions. The transport term pattern was similar to that of the scale height and it is suggested as a proxy parameter for quantifying low-latitude plasma irregularities and distribution of thermospheric composition.  相似文献   

18.
The method of updating the IRI-1990 model for the effects of large scale travelling ionospheric disturbances (LSTIDs) of the F-region have been established as a result of the statistical analysis of VI sounding data of European- Asian sector for 380 substorms. The effect of a LSTID occurring during a substorm is modelled as an additive correction to the quiet IRI foF2 and hmF2 values which depends on the time since the onset of the substorm and on the maximum AE-index value reached in it. In case of a sequence of substorms the additive corrections ΔfoF2 and ΔhmF2 are taking into account. The method for night-time ionosphere is described in detail, the listing of the FORTRAN-77 computer program is available on the floppy disk. The method is workable as correction to any ionospheric model.  相似文献   

19.
Monthly average electron density profiles have been calculated from hourly electron density N(h) recorded in 26 digisonde stations distributed worldwide encompassing the time interval 1998–2006. The ionospheric electron density peak height of the F2 region, hmF2, and the effective scale height at the hmF2, Hm, deduced from average profiles have been analyzed to obtain the quiet-time behavior and have been analytically modeled by the spherical harmonic analysis (SH) technique using the modip latitude as the coordinate of the reference system. The coefficients of the SH models of hmF2 and Hm are bounded to the solar activity, and the temporal and seasonal variations are considered by Fourier expansion of the coefficients. The SH models provide a tool to predict hmF2 and Hm located anywhere in the range of latitudes between of 70°N and 70°S and at any time. The SH analytical model for hmF2 improves the fit to the observations by 10% in average compared to the IRI prediction, and it might improve the IRI prediction of hmF2 by more than 30% at high and low latitudes. The analytical model for Hm predicts the quiet behavior of the effective scale height with accuracy better than 15% in average which enables to obtain a good estimation of vertical profiles. These results could be useful to estimate information for the topside profile formulation.  相似文献   

20.
An analysis of properties and peculiarities of the nighttime winter foF2 increases (NWI) in the East Siberia is made on data of ionospheric station Irkutsk in the periods 1958–1992 and 2002–2009 and the empirical model of the F2 layer critical frequency under the geomagnetic quiet conditions deduced from these data (model Q-F2). It is revealed, that the NWI is the stable regularity of the quiet ionosphere over Irkutsk. The amplitude of the NWI (the difference between maximum and minimum foF2 values at night hours) is the greatest in December–January and nearly the same at low and middle solar activity. It is a peculiarity of the quiet ionosphere in the East Siberia. Maximum in night foF2 under quiet geomagnetic conditions is observed mainly after midnight (02-04 LT) and is shifted to predawn hours as solar activity increases. At low solar activity the quiet ionosphere at ∼02–04 LT shows the following properties: (a) the fluctuations of foF2 and hmF2 are in the reverse correlation but this dependence is weak; (b) very strong fluctuations of foF2 (|δfoF2| > 30%) occur seldom (∼4% of events) and almost all of them are positive; an example of very strong fluctuations of foF2 up to 60% can be an extreme increase in the foF2 on 19.12.2008; (c) the very strong enhancements of foF2 in the NWI maximum can be observed at the low geomagnetic activity, they occur more often during substorms but very seldom during geomagnetic storms. Possible reasons of these properties of NWI are discussed.  相似文献   

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

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