首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 578 毫秒
1.
This work investigates the spectral atmospheric transmittance due to aerosols in the urban environment of Athens during a period of one year. The spectral transmittance due to aerosols is derived using measurements of spectral direct-beam solar irradiance in the 310–575 nm spectral band. This derivation is accomplished by using a radiative transfer model for estimating the partial spectral atmospheric transmittance functions due to Rayleigh scattering, and absorption by ozone, nitrogen dioxide and water vapor. The seasonal and diurnal variation of the aerosol transmittance is investigated and the results are discussed with a view to air pollution sources, meteorological factors and topographic characteristics of the Athens basin.  相似文献   

2.

Photosynthetically active radiation (PAR) is essential for plant photosynthesis and carbon cycle, and is also important for meteorological and environmental monitoring. To advance China’s disaster and environmental monitoring capabilities, the HJ-1A/B satellites have been placed in Earth orbit. One of their environmental monitoring objectives is the study of PAR. We simulated direct solar, scattered and environment radiation between 400 and 700 nm under different atmospheric parameters (solar zenith angle, atmospheric water vapor, atmospheric ozone, aerosol optical thickness, surface elevation and surface albedo), and then established a look-up table between these input parameters and PAR. Based on the look-up table, we used HJ-1A/B aerosol and surface albedo outputs to derive the corresponding PAR. Validation of inversed instantaneous and observed PAR values using HJ-1 Heihe experimental data had a root mean square error of 25.2 W m−2, with a relative error of 5.9%. The root mean square error for accumulated daily PAR and observed values was 0.49 MJ m−2, with a relative error of 3.5%. Our approach improved significantly the computational efficiency, compared with using directly radiation transfer equations. We also studied the sensitivity of various input parameters to photosynthetically active radiation, and found that solar zenith angle and atmospheric aerosols were sensitive PAR parameters. Surface albedo had some effect on PAR, but water vapor and ozone had minimal impact on PAR.

  相似文献   

3.
Photosynthetically active radiation (PAR) is essential for plant photosynthesis and carbon cycle, and is also important for meteorological and environmental monitoring. To advance China’s disaster and environmental monitoring capabilities, the HJ-1A/B satellites have been placed in Earth orbit. One of their environmental monitoring objectives is the study of PAR. We simulated direct solar, scattered and environment radiation between 400 and 700 nm under different atmospheric parameters (solar zenith angle, atmospheric water vapor, atmospheric ozone, aerosol optical thickness, surface elevation and surface albedo), and then established a look-up table between these input parameters and PAR. Based on the look-up table, we used HJ-1A/B aerosol and surface albedo outputs to derive the corresponding PAR. Validation of inversed instantaneous and observed PAR values using HJ-1 Heihe experimental data had a root mean square error of 25.2 W m?2, with a relative error of 5.9%. The root mean square error for accumulated daily PAR and observed values was 0.49 MJ m?2, with a relative error of 3.5%. Our approach improved significantly the computational efficiency, compared with using directly radiation transfer equations. We also studied the sensitivity of various input parameters to photosynthetically active radiation, and found that solar zenith angle and atmospheric aerosols were sensitive PAR parameters. Surface albedo had some effect on PAR, but water vapor and ozone had minimal impact on PAR.  相似文献   

4.
Observing and Modeling Earth’s Energy Flows   总被引：1，自引：0，他引：1  
This article reviews, from the authors’ perspective, progress in observing and modeling energy flows in Earth’s climate system. Emphasis is placed on the state of understanding of Earth’s energy flows and their susceptibility to perturbations, with particular emphasis on the roles of clouds and aerosols. More accurate measurements of the total solar irradiance and the rate of change of ocean enthalpy help constrain individual components of the energy budget at the top of the atmosphere to within ±2 W m?2. The measurements demonstrate that Earth reflects substantially less solar radiation and emits more terrestrial radiation than was believed even a decade ago. Active remote sensing is helping to constrain the surface energy budget, but new estimates of downwelling surface irradiance that benefit from such methods are proving difficult to reconcile with existing precipitation climatologies. Overall, the energy budget at the surface is much more uncertain than at the top of the atmosphere. A decade of high-precision measurements of the energy budget at the top of the atmosphere is providing new opportunities to track Earth’s energy flows on timescales ranging from days to years, and at very high spatial resolution. The measurements show that the principal limitation in the estimate of secular trends now lies in the natural variability of the Earth system itself. The forcing-feedback-response framework, which has developed to understand how changes in Earth’s energy flows affect surface temperature, is reviewed in light of recent work that shows fast responses (adjustments) of the system are central to the definition of the effective forcing that results from a change in atmospheric composition. In many cases, the adjustment, rather than the characterization of the compositional perturbation (associated, for instance, with changing greenhouse gas concentrations, or aerosol burdens), limits accurate determination of the radiative forcing. Changes in clouds contribute importantly to this adjustment and thus contribute both to uncertainty in estimates of radiative forcing and to uncertainty in the response. Models are indispensable to calculation of the adjustment of the system to a compositional change but are known to be flawed in their representation of clouds. Advances in tracking Earth’s energy flows and compositional changes on daily through decadal timescales are shown to provide both a critical and constructive framework for advancing model development and evaluation.  相似文献   

5.
We examine the impact of the Atmospheric Brown Clouds on the direct radiative forcing of the Himalayan foothills and the Indo-Gangetic Plains (IGP) regions, home for over 500 million S. Asians. The NASA-Terra MODIS satellite data reveal an extensive layer of aerosols covering the entire IGP and Himalayan foothills region with seasonal mean AODs of about 0.4 to 0.5 in the visible wavelengths (0.55 micron), which fall among the largest seasonal mean dry season AODs for the tropics. We show new surface data which reveal the presence of strongly absorbing aerosols that lead to a large reduction in solar radiation fluxes at the surface during the October to May period. The three-year mean (2001 to 2003) October to May seasonal and diurnal average reduction in surface solar radiation for the IGP region is about 32 (±5) W m−2 (about 10% of TOA insolation or 20% of surface insolation). The forcing efficiency (forcing per unit optical depth) is as large as −27% (note that the forcing is negative) of top-of-atmosphere (TOA) solar insolation, and exceeds the forcing efficiency that has been observed for other polluted regions in America, Africa, East Asia, and Europe. General circulation model sensitivity studies suggest that both the local and remote influence of the aerosol induced radiative forcing is to strengthen the lower atmosphere inversion, stabilize the boundary layer, amplify the climatological tendency for a drier troposphere, and decrease evaporation. These aerosol-induced changes could potentially increase the life times of aerosols, make them more persistent, and decrease their single scattering albedos, thus potentially leading to a detrimental positive feedback between aerosol concentrations, aerosol forcing, and aerosol persistence. In addition, both the model studies and observations of pan evaporation suggest that the reduction in surface solar radiation may have led to a reduction in surface evaporation of moisture. These results suggest the vulnerability of this vital region to air pollution related direct and indirect (through climate changes) impacts on agricultural productivity of the region.  相似文献   

6.
测量夜间135.6 nm大气气辉辐射强度是目前有效的电离层探测方式之一,我国即将在风云三号卫星上搭载仪器,利用该波段夜气辉辐射测量来反演电子总含量(本文所指电子总含量表示卫星高度以下大气柱的电子含量)及峰值电子密度,因此非常有必要开展相关的气辉发光模型及反演研究.在介绍氧原子135.6 nm波段夜气辉激发机制基础上,考虑辐射在传输过程中受到大气氧原子的散射及氧气分子的吸收,采用迭代法求解包含多次散射及大气吸收衰减的辐射传输方程,得到该波段的体发射率,最终通过考虑包含辐射传输的路径积分计算得到135.6 nm气辉辐射强度值.对结果的分析表明:该气辉模型能较好地描述体发射率随高度的分布特征,计算得到的135.6 nm夜气辉辐射强度在不同时空及太阳活动的分布与相应条件下峰值电子密度(NmF2)及电子总含量(TEC)的分布基本一致.相同的时空及太阳活动输入条件下,模式计算的135.6 nm夜气辉辐射强度与国外同类模式结果的值平均偏差约为3%.文中最后介绍了通过135.6 nm夜气辉的辐射强度探测来反演电离层峰值电子密度NmF2及电子总含量TEC的反演方法.  相似文献   

7.
Fog is an atmospheric phenomenon that has important environmental consequences related to visibility, air quality and climate change on local and regional scales. The formation of radiation fog results from a complex balance between surface radiative cooling, turbulent mixing in the surface layer, aerosol growth by deliquescence and activation of fog droplets. During the ParisFog field experiment, out of 16 events forecasted for radiation fog, activated fog materialized in seven events, while in five other events the visibility dropped to 1–2 km but haze particle size remained below the critical size of activation. To better understand the conditions that lead to or do not lead to sustained fog droplet activation, we performed a comparative study of dynamic, thermal, radiative and microphysical processes occurring between sunset and fog (or quasi-fog) onset. We selected two radiation fog events and two quasi-radiation fog events that occurred under similar large-scale conditions for this comparative study. We identified that aerosol growth by deliquescence and droplet activation actually occurred in both quasi-fog events, but only during <1 h. Based on ParisFog measurements, we found that the main factors limiting sustained activation of droplets at fog onset in the Paris metropolitan area are (1) lack of mixing in the surface layer (typically wind speed <0.5 ms?1), (2) relative humidity exceeding 90 % throughout the residual layer, (3) low cooling rate in the surface layer (typically less than ?1 °C per hour on average) due to weak radiative cooling (0 to ?30 Wm?2) and near zero sensible heat fluxes, and (4) a combination of the three factors listed above during the critical phase of droplet activation preventing the transfer of cooling from the surface to the liquid layer. In addition, we found some evidence of contrasted aerosol growth by deliquescence under high relative humidity conditions in the four events, possibly associated with the chemical nature of the aerosols, which could be another factor impacting droplet activation.  相似文献   

8.
—This paper deals with a new broadband atmospheric model designed for predicting the total and diffuse solar radiation incident on the earth’s surface in medium or large-sized coastal or near-coastal cities, under a clear or cloudy sky. The revised solar spectrum is used. The atmospheric transmittance of each atmospheric parameter contributing to solar radiation depletion, water vapor, ozone, uniformly mixed gases, molecules and aerosols, is calculated using parameterized expressions resulting from integrated spectral transmittance functions. The beam and diffuse radiation components are obtained as a function of the specific atmospheric transmittances. The model requires the following parameters as inputs: total water vapor and ozone amount in a vertical column, sunshine duration and the surface albedo. The model has been used for validation purposes at two stations with slightly different characteristics (NOA and Penteli) in the Athens basin, where total and diffuse radiation measurements are available, for a period of 34 months for NOA and 23 for Penteli. The NOA station is located on a small hill (107 m a.m.s.l.) near the center of Athens, while the Penteli station (500 m a.m.s.l.) is situated in a relatively less polluted area in northern Athens. The clear sky part of the model was tested for 70 individual “clear” days with 2-minute intervals, while the whole model was checked with monthly “mean” days and mean hourly values. A close agreement between the calculated and the measured values of total and diffuse solar radiation is observed.  相似文献   

9.
This paper discusses the formation mechanisms of infrared radiation in the mesosphere and lower thermosphere (MLT), the energetic effects of the radiative absorption/emission processes, and the retrieval of atmospheric parameters from infrared radiation measurements. In the MLT and above, the vibrational levels of the molecules involved in radiative transitions are not in local thermodynamic equilibrium (LTE) with the surrounding medium, and this then requires specific theoretical treatment. The non-LTE models for CO2, O3, and H2O molecules are presented, and the radiative cooling/heating rates estimated for five typical atmospheric scenarios, from polar winter to polar summer, are shown. An optimization strategy for calculating the cooling/heating rates in general circulation models is proposed, and its accuracy is estimated for CO2. The sensitivity of the atmospheric quantities retrieved from infrared observations made from satellites to the non-LTE model parameters is shown.  相似文献   

10.
This study used realistic representations of cloudy atmospheres to assess errors in solar flux estimates associated with 1D radiative transfer models. A scene construction algorithm, developed for the EarthCARE mission, was applied to CloudSat, CALIPSO and MODIS satellite data thus producing 3D cloudy atmospheres measuring 61 km wide by 14,000 km long at 1 km grid-spacing. Broadband solar fluxes and radiances were then computed by a Monte Carlo photon transfer model run in both full 3D and 1D independent column approximation modes. Results were averaged into 1,303 (50 km)2 domains. For domains with total cloud fractions A c  < 0.7 top-of-atmosphere (TOA) albedos tend to be largest for 3D transfer with differences increasing with solar zenith angle. Differences are largest for A c  > 0.7 and characterized by small bias yet large random errors. Regardless of A c , differences between 3D and 1D transfer rarely exceed ±30 W m?2 for net TOA and surface fluxes and ±10 W m?2 for atmospheric absorption. Horizontal fluxes through domain sides depend on A c with ～20% of cases exceeding ±30 W m?2; the largest values occur for A c  > 0.7. Conversely, heating rate differences rarely exceed ±20%. As a cursory test of TOA radiative closure, fluxes produced by the 3D model were averaged up to (20 km)2 and compared to values measured by CERES. While relatively little attention was paid to optical properties of ice crystals and surfaces, and aerosols were neglected entirely, ～30% of the differences between 3D model estimates and measurements fall within ±10 W m?2; this is the target agreement set for EarthCARE. This, coupled with the aforementioned comparison between 3D and 1D transfer, leads to the recommendation that EarthCARE employ a 3D transport model when attempting TOA radiative closure.  相似文献   

11.
The effect of atmospheric aerosols on the temperatures of a zonal average climate model is investigated. This is done by introducing into the climate model the results of extensive calculations of the effect of aerosols on the partitioning of solar radiation. Calculations are performed for a non-absorbing and a sloghtly absorbing aerosol, for average and for heavy aerosol amounts. The results indicate that the presence of atmospheric aerosols causes a cooling of the earth's surface and atmosphere, at all latitudes, with aerosol amounts, especially for the absorbing aerosols. The results are compared with some of the previous estimates of aerosol effects on climate.  相似文献   

12.
For calculating photolysis rates and solar heating in the atmosphere, the radiation field has to be calculated very accurately. Previous investigations have shown that for large solar zenith angles a solution of the radiation equation which accounts for the Earth\’s curvature is needed. A new simplified version of the 3D radiation equation in spherical geometry allowing for anisotropic scattering is presented. The horizontal variation of physical quantities, the variation of the solar zenith angle with different longitude and latitude for the scattering calculation for one vertical column of air and any effects of refraction are neglected. A numerical model is introduced which efficiently solves this new 3D radiation equation accurately. The effects of anisotropic scattering are shown to be very important for the directional dependence of the scattered intensity. Anisotropic scattering by aerosols and air molecules can change the intensity in certain directions by up to 180% and 25%, respectively. However, most of these changes cancel each other out when averaged over all angles, so that the effect of anisotropic scattering for large solar zenith angles on the mean intensity (actinic flux) is much smaller, i.e. less than 10%. For the heating rates, the effect of anisotropic scattering for large solar zenith angles is even smaller, being less than a few percent. Generally, the effects of anisotropic scattering and the effects of including aerosols are the larger on higher altitudes the larger the solar zenith angle is. Results of the model are shown to compare well with results of previous investigations, including the independent work of Dahlback and Stamnes. The agreement is especially good in the case of isotropic scattering by air molecules and neglecting the effects of aerosols.  相似文献   

13.
辐射参数化对海南岛海风雷暴结构模拟的影响   总被引：4，自引：1，他引：3       下载免费PDF全文
本文利用高分辨率WRF模式探讨了两组短波、长波辐射参数化方案(Dudhia+RRTM、RRTMG+RRTMG)对海南岛一次海风雷暴模拟的影响及其可能的物理机制.结果表明,辐射参数化能影响大气的加热程度和近地面能量,决定海陆温差和气压差,改变海南岛的海风特征,最终影响海风雷暴的发生发展.Dudhia+RRTM方案模拟的短波、长波综合加热率、感热通量以及潜热通量都大于RRTMG+RRTMG方案,造成了前者模拟的近地面能量偏高,大气层结也表现得更加不稳定,进而使得该方案下的海陆温差和气压差相对较大,Dudhia+RRTM方案模拟的海风明显强于RRTMG+RRTMG方案,能提供更好的水汽输送和抬升条件,有利于海风雷暴的发生发展,因此其模拟的雷暴活动范围和对流中心强度都要大于RRTMG+RRTMG方案.  相似文献   

14.
The influence of atmospheric solar radiation absorption on the photodetachment, dissociative photodetachment, and photodissociation rate coefficients (photodestruction rate coefficients) of O?, Cl?, O2 ?, O3 ?, OH?, NO2 ?, NO3 ?, O4 ?, OH?(H2O), CO3 ?, CO4 ?, ONOO?, HCO3 ?, CO3 ?(H2O), NO3 ?(H2O), O2 +(H2O), O4 +, N4 +, NO+(H2O), NO+(H2O)2, H+(H2O) n for n = 2–4, NO+(N2), and NO+(CO2) at D-region altitudes of the ionosphere is studied. A numerical one-dimensional time-dependent neutral atmospheric composition model has been developed to estimate this influence. The model simulations are carried out for the geomagnetically quiet time period of 15 October 1998 at moderate solar activity over the Boulder ozonesonde. If the solar zenith angle is not more than 90° then the strongest influence of atmospheric solar radiation absorption on photodestruction of ions is found for photodissociation of CO4 ? ions when CO3 ? ions are formed. It follows from the calculations that decreases in the photodestruction rate coefficients of ions under consideration caused by this influence are less than 2 % at 70 km altitude and above this altitude if the solar zenith angle does not exceed 90°.  相似文献   

15.
Agriculture crop residue burning in tropics is an important source of atmospheric aerosols and monitoring their long-range transport is an important element in climate change studies. Synchronous measurements using micro-pulsed lidar, MICROTOPS-II sun photometer, multi-filter rotating shadow band radiometer (MFRSR) on aerosol optical depth and ground reaching solar irradiance were carried at an urban location in central region of India. Aerosol backscatter profiles obtained from micro-pulse lidar showed elevated aerosol layers up to ～3 km on certain days during October 2007. Satellite data observations on aerosol properties suggested transport of particles from agriculture crop residue burning in Indo-Gangetic Plains over large regions. Radiative forcing of aerosols estimated from SBDART model with input information on aerosol chemical properties, aerosol optical depth and single scattering albedo and broadband solar irradiance measurements using MFRSR showed good correlation (R=0.98).  相似文献   

16.
Estimating transpiration of the trees in agroforestry system is important in water management of the site. Sap flow of intercropped fast‐growing young poplar trees and microclimate factors in semiarid northeastern China was measured in two growing seasons (2008 and 2009). Sapwood growth and water storage of wood and leaf increment during the growing season were involved in the calculation of sap flow. The results showed that diurnal variation of sap flow followed to that of short wave solar radiation. Sap flows both in 10 min mean and daily gross values mainly depended on solar radiation and vapor pressure deficit, and the relations well fit hyperbolic function. The regression coefficients of monthly window data indicated that the seasonal variation of sap flow capacities decreased gradually from June to September. Moderate soil water stress of upper soil layer (0–50 cm) did not constrain the sap flow because the trees could use the water at deeper soil layer. The daily sap flow per tree ranged 0.8 to 18.1 and 3.7 to 23.8 kg d?1 tree?1, with averages of 8.7 and 14.3 kg d?1 tree?1 in 2008 and 2009 respectively. An empirical model was established to estimate the sap flow of the poplar trees by solar radiation, vapor pressure deficit, leaf area index and Julian days. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

17.
With the intensification of pollution and urbanization, the aerosol radiation effect continues to play an important role in the urban boundary layer. In this paper, a winter pollution process in Beijing has been taken as an example, and a new aerosol vertical profile in the radiative parameterization scheme within the Weather Forecast Research and Forecasting (WRF) model has been updated to study the effect of aerosols on radiation and the boundary layer. Furthermore, the interactions among aerosols, urbanization, and planetary boundary layer (PBL) meteorology were discussed through a series of numerical experiments. The results show the following: (1) The optimization improves the performance of the model in simulating the distribution features of air temperature, humidity, and wind in Beijing. (2) The aerosols reduce the surface temperature by reducing solar radiation and increasing the temperature in the upper layer by absorbing or backscattering solar radiation. The changes in the PBL temperature lead to more stable atmospheric stratification, reducing the energy transfer from the surface and the height of the boundary layer. (3) With the increase in the aerosol optical depth, the atmospheric stratification most likely becomes stable over rural areas, most likely becomes stable over suburb areas, and has great difficultly becoming stable over urban areas. Aerosol radiative forcing, underlying urban surfaces, and the interaction between them are the main factors that affect the changes in the meteorological elements in the PBL.  相似文献   

18.
The aim of this study was to validate evaporation models that can be used for palaeo‐reconstructions of large lake water levels. Lake Titicaca, located in a high‐altitude semi‐arid tropical area in the northern Andean Altiplano, was the object of this case study. As annual evaporation is about 90% of lake output, the lake water balance depends heavily on the yearly and monthly evaporation flux. At the interannual scale, evaporation estimation presents great variability, ranging from 1350 to 1900 mm year?1. It has been found that evaporation is closely related to lake rainfall by a decreasing relationship integrating the implicit effect of nebulosity and humidity. At the seasonal scale, two monthly evaporation data sets were used: pan observations and estimations derived from the lake energy budget. Comparison between these data sets shows that (i) there is one maximum per year for pan evaporation and two maxima per year for lake evaporation, and (ii) pan evaporation is greater than lake evaporation by about 100 mm year?1. These differences, mainly due to a water depth scale factor, have been simulated with a simple thermal model θw(h, t) of a free‐surface water column. This shows that pan evaporation (h = 0·20 m) is strongly correlated with direct solar radiation, whereas the additional maximum of lake evaporation (h = 40 m) is related to the heat restitution towards the atmosphere from the water body at the end of summer. Finally, five monthly evaporation models were tested in order to obtain the optimal efficiency/complexity ratio. When the forcing variables are limited to those that are most readily available in the past, i.e. air temperature and solar radiation, the best results are obtained with the radiative Abtew model (r = 0·70) and with the Makkink radiative/air temperature model (r = 0·67). Copyright © 2007 John Wiley & Sons, Ltd.  相似文献   

19.
Snow and ice present interesting challenges to hydrologists. Simulating the radiative balance over snow, which is an important part of surface–atmosphere interactions, is particularly challenging because of the decay in albedo over time and the difficulty in estimating surface temperature and incoming long-wave radiation fluxes. Few models are available that include a comprehensive energy and water balance for cold season conditions. The simultaneous heat and water model (SHAW) is a detailed, physical process model of a vertical, one-dimensional canopy–snow–residue–soil system which integrates the detailed physics of heat and water transfer through a plant canopy, snow, residue and soil into one simultaneous solution. Detailed provisions for metamorphosis of the snowpack are included. The SHAW model was applied to data for one winter/spring season (November to May) on a ploughed field in Minnesota without prior calibration to test the performance of the radiation components. Maximum snow depth during this period was 30 cm. For the nearly 100 days of snowcover, the model accounted for 69% of the variation in net solar radiation, 66% of the variation in incoming long-wave radiation, 87% of the variation in emitted long-wave radiation, 26% of the variation in net long-wave radiation and 55% of the variation in net radiation balance. Mean absolute error in simulated values ranged from 10 W m−2 for emitted long-wave radiation to 27 W m−2 for the entire net radiation balance. Mean bias error ranged from 8 W m−2 for emitted long-wave radiation to −16 W m−2 for the entire net radiation balance. When the entire 170 days of simulation, which included periods without snowcover, were included in the analysis, the variation in observed values increased greatly. As a result, the variation in observed values accounted for by the model increased to 97, 71, 93, 56 and 94%, respectively, while the mean absolute and mean bias errors in simulated values remained nearly the same. Model modifications and parameter adjustments necessary to improve winter-time simulation were investigated. Simulation results suggest that the SHAW model may be a useful tool in simulating the interactive influences of radiative transfer at the surface–atmosphere interface.  相似文献   

20.
G. Kockarts 《Annales Geophysicae》1994,12(12):1207-1217
The Schumann-Runge bands of molecular oxygen between 175 nm and 205 nm cover a spectral region of great importance for the effect of solar radiation in the middle atmosphere. The highly structured absorption cross sections that can be computed with good accuracy lead to prohibitive computer resources, however, when they have to be used in complex atmospheric models. Although various approximations have been developed, another approach is proposed to avoid spurious effects, even when the new approximations are used outside of their validity range, corresponding to a solar attenuation of the order of 10?10. Comparisons between exact computations and present approximations lead to satisfactory agreement.  相似文献   

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

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