Geog595 Ecological ModelingSpring 2010Due: 5pm Jan 26, 2010Lab 1 Calculating Incoming Solar Radiation with C1. Objectives(1) Learning programming in ecological modeling in C.(2) Development a module that calculates solar radiation arriving at the Earth surface at any place any day and any time.2. Input DataLatitude and longitudeMeasured total photosynthetically active radiation (PAR)Month, Date and Hour of measurementsAir Pressure3. TheorySolar radiation is the ultimate source of energy that drives all ecosystem process. It is essential to modelers to know how much solar radiation arriving at the Earth surface in order to model other ecological processes. At the very top of the atmosphere, the solar radiation flux density is usually considered a constant, “Solar Constant” (S0=1367.0 W/m2). However, for scientific computation, it can no longer be considered as constant, it actually varies as a result of Sun-Earth distance changes through the year. ))365/)10(2cos(033.00.1(*13670 JS,where J is Julian date. The amount of solar radiation arriving at the Earth surface changesbecause of the change in atmospheric conditions and position of the Sun in the sky. The total radiation at the top of the atmosphere is composed of ultraviolet (SU), visible (SV) and near-infrared (SN) as in the following:0034.0 SSU0439.0 SSV0527.0 SSNDue to the lack of account for SU in the literature, let assume ultraviolet has the same optical properties as the Visible and combine them into a single component of SUV. Then we will separate Visible light out from SVU when modeling photosynthesis. For convenience we will call SVU as SV.0473.0 SSVAccording to (Weiss and Norman, 1985), the potential U and V direct radiation is)cos())/(1556.0exp(0mPPSRVDVWhere P is air pressure, and P0 is standard air pressure (P0=101,325Pa). m is the atmospheric optical depth as)cos(/1mThe potential diffuse Visible reaching the horizontal surface is considered 40% of the scattered radiation:)cos()(4.0DVVUdVRSR The potential direct NIR reaching the horizontal surface is )cos())/(0086.0exp(0wmPPSRNDNAnd the potential diffuse NIR radiation on the horizontal surface is)cos())((6.0wRSRDNNdVω is the water absorption in the near infrared for 10 mm of precipitation water))(log0345.0log4459.0195.1(log2101010mmantiSwNThen, the fraction of direct beam in PAR and NIR are321BRATIOARRfVDVV321CRATIOCRRfNDNNWhere A=0.9, B=0.7, C=0.88, and D=0.68. RATIO is the measured to potential solar radiation, i.e. RATIO=RT/(RV+RN). Often we do not have RT measure, but RV only, we can convert RV to RT as RT=RV/0.439. Calculating solar zenith angle by time of day:)cos()cos()cos()sin()sin()sin()cos( hWhere τ is local time hour angle. τ=15×(t-12) degrees, t is in hours; φ is latitude; δ is sun declination angle. 365)284(2sin5.23JThe Goudriaan and Van Laar (1994) model, as given in Leuning et al. (1995) for fractions of diffuse radiation for total radiation based on atmospheric transmittance (τa):7.07.03.03.0 2.0)3.0(211aaaadfThe Spitters et al (1986) model for fractions of diffuse radiation for total radiation: 35.066.147.135.022.0)22.06.4(-122.01aaa2aadKRKfaWhere R=0.847-1.61cos(θ)+1.04cos2(θ), and K=(1.47-R)/1.66.4. Lab Report(1) Run your program to get the diffuse fraction of total radiation for the three approaches and make a time-series plot for them. Compare the difference.(2) Modify the model so that the program only runs the Weiss and Norman (1985)model with direct and diffuse radiation (not fractions) for PAR and NIR, delete all unnecessary statements and add programs as needed, save the new program as Weiss_Norman.c. (3) Plot the four components of radiation you produced in (2) above as a function of time of day with excel. Explain the seasonal and diurnal patterns you