Investigation of the 3-D actinic flux field in mountainous terrain.

During three field campaigns spectral actinic flux was measured from 290-500 nm under clear sky conditions in Alpine terrain and the associated O3- and NO2-photolysis frequencies were calculated and the measurement products were then compared with 1-D- and 3-D-model calculations. To do this 3-D-radiative transfer model was adapted for actinic flux calculations in mountainous terrain and the maps of the actinic flux field at the surface, calculated with the 3-D-radiative transfer model, are given. The differences between the 3-D- and 1-D-model results for selected days during the campaigns are shown, together with the ratios of the modeled actinic flux values to the measurements. In many cases the 1-D-model overestimates actinic flux by more than the measurement uncertainty of 10%. The results of using a 3-D-model generally show significantly lower values, and can underestimate the actinic flux by up to 30%. This case study attempts to quantify the impact of snow cover in combination with topography on spectral actinic flux. The impact of snow cover on the actinic flux was ~ 25% in narrow snow covered valleys, but for snow free areas there were no significant changes due snow cover in the surrounding area and it is found that the effect snow-cover at distances over 5 km from the point of interest was below 5%. Overall the 3-D-model can calculate actinic flux to the same accuracy as the 1-D-model for single points, but gives a much more realistic view of the surface actinic flux field in mountains as topography and obstruction of the horizon are taken into account.

Introduction

Ultraviolet spectroradiometry has grown in importance during the last few decades due to the increasing need for high quality measurements of UV radiation in areas such as medicine, e.g. Melamed et al. (2008), atmospheric chemistry, e.g. Raff and Finlayson-Pitts (2010), and material science, e.g. Kaczmarek et al. (2009), for high-quality measurements of UV radiation (UVR). There are numerous case studies and long-term investigations into global solar irradiance, defined as the sum of direct and diffuse irradiance coming from the upper hemisphere, received on a flat horizontal surface (e.g. McKenzie et al. (1993), Pribullova and Chmelik (2008)). While irradiance is the standard quantity for solar radiation, measurements of actinic flux have become increasingly popular enabled by the development of suitable input optics (Bais et al. (2001), Hofzumhaus et al. (2002)). Actinic flux represents the radiation which is received at a point with each photon being assigned with the same weight irrespective of incident direction. Recent investigations in the actinic flux have used airborne measurements and model calculations of actinic flux (Palancar et al. (2011), Bierwirth et al. (2010)) and for less complex situations (low surface albedo, clear sky conditions, simple topography) the agreement between models and measurements is within 5–10%, but for more complex situations (variable surface albedo, partial cloud cover, complex topography) deviations are much larger (Bais et al. (2003), Thiel et al. (2008)). Spectral actinic flux can be integrated to derive photolysis frequencies J, which describe specific photochemical dissociation processes. Therefore the actinic flux is an elementary quantity for atmospheric photochemistry. However, measuring actinic flux is more difficult than irrradiance, since a detector with uniform response characteristics over one hemisphere (2 π sr) is needed (Volz-Thomas et al. (1996) and Junkermann et al. (1989)).

Spectroradiometers with suitable input optics are the most accurate instruments for measuring the spectrally resolved UV actinic flux (Bais et al. (2001). As ground-based high quality spectral UV actinic flux measurements are only available for a small number of research sites worldwide, different calculation methods have been developed over the last three decades (e.g. Stamnes et al., 1988). 1-D radiative transfer models like DISORT (Stamnes et al., 1988) or POLRADTRAN (Evans and Stephens, 1991) have been popular for about two decades while 3-D radiative transfer models (e.g., Cahalan et al., 2005) are a more recent development. A comprehensive intercomparison of 3-D radiative transfer models is available at http://i3rc.gsfc.nasa.gov/I3RC-intro.html.

The UV actinic flux under clear sky conditions, is determined largely by solar zenith angle and total column ozone (Feister and Grewe, 1995a, Burrows, 1997, Kerr, 2003), but the aerosol optical depth and the ground albedo also have a significant impact. Investigations into the impact of surface albedo on the radiation field were conducted by Blumthaler and Ambach (1988), Kylling et al. (2000), Schwander et al. (1999) and Weihs et al. (2001). All these studies suggest that surface albedo, especially in the presence of snow, has a significant influence on the UVR field. Actinic flux can be enhanced significantly more strongly than irradiance by surface albedo because of its angle dependence (Madronich, 1987). Other studies have shown that the UVR is significantly affected by variations in atmospheric aerosols especially in urban areas (Jaroslawski , 2005, Kylling et al., 1998). Furthermore the impact of both aerosols and surface albedo on UVR for a particular mountain site (Sonnblick 3106 m) has been investigated (Weihs et al., 1999, Simic et al., 2008,Rieder et al., 2010). The interplay of mountainous topography and high surface albedo with respect to spectral irradiance at 340 nm was studied by Kylling and Mayer (2001) using measurements and 3-D-model simulations. Lenoble et al. (2004) also studied the enhancement of the UV irradiance due to snow cover on the ground in mountainous terrain. Both studies report an enhancement due to snow cover of about 25% for clear sky situations. General aspects of 3-D-radiative transfer are discussed by Bass et al. (2010) and Deutschmann et al. (2011). Degünther et al. (1998) used a 3-D-model to calculate UV irradiance under cloud-free conditions in mountainous, partly snow-covered terrain and reported a large horizontal impact of snow-covered surface with significant enhancements of up to 3% due to ground albedo affects from more than 40 km away from the measurement site. It is not straight forward to convert UV irradiance to UV actinic flux but Webb et al. (2002) and Kylling et al. (2003) have developed conversion algorithms. The conversion works better at certain low solar zenith angles and information on the sky situation (clear sky or cloudy) is also crucial. However modeling of UV spectral actinic flux in mountainous areas, especially in the presence of snow, is complex and still poses a scientific challenge. So far the influence of snow and topography on surface actinic flux has not been quantified, a problem adressed by this work.

We compare measurements with 1-D and 3-D-model calculations of spectral actinic flux at 305, 310, 324 and 380 nm. In addition photolysis frequencies for ozone (O3) and nitrogen dioxide (NO2) can be calculated from both modeled and measured spectra. These two reactions are regarded as the most important photodissociation reactions in the troposphere. Photodissociation of ozone (O3 + hv → O2 + O(1D)) generates an exited oxygen atom which can react with water to give the hydroxyl radical (O(1D) + H2O) → 2OH). The hydroxyl radical is central to atmospheric chemistry as it initiates the oxidation of hydrocarbons in the atmosphere and so acts as a detergent. Photodissociation of nitrogen dioxide (NO2 → NO + O) is a key reaction of tropospheric ozone (O + O2 → O3). Quantifying the impact of surface albedo and topography on actinic fluxes under cloud-free conditions, in Alpine terrain is thus an important contribution to a better understanding of the chemistry in such regions.

Material and methods

Ground-based measurements

Ground-based measurements were performed during two campaigns near Innsbruck, in late summer 2007 (2nd September to 22nd September) and in late winter 2008 (19th February to 6th March) and one measurement campaign in spring 2008 (29th April to 9th May) near Sonnblick. It was the intention to perform measurements first for snow-free conditions in September and then for snow-covered conditions in February and May. Unfortunately due to the presence in Innsbruck region of more than the usual snow in September and less snow in February, the resulting differences in snowline were only small (1900 m in September and between 1300 and 1700 m in February). Prior to each campaign an instrument intercomparison was performed. Differences of less than 5% were found between the reference instrument from the Medical University Innsbruck and the instruments from the two other participating parties, the University of Manchester and the University of Natural Resources and Life Sciences (BOKU), Vienna.

During the first campaign, spectral UV actinic flux was measured at two stations (see Table 1 and Fig. 1 for details) using two Bentham DTM 300 spectroradiometers (Bernhard and Seckmeyer, 1999). The entrance optic is designed such that photons incident from all angles are equally represented — there is no significant angular dependence (Hofzumahaus et al., 1999). The Innsbruck station (Inn) is situated on a roof in the city center of Innsbruck and represents a typical valley site. The second station, Hafelekar (Haf) is only 7 km away from the Innsbruck site, but on top of the mountain range North of Innsbruck (Nordkette). During the second campaign in Innsbruck, measurements using Bentham DTM 300s were carried out at Innsbruck and Hafelekar sites and using a Bentham DM 150 at an additional site in Lans. This site is about 6 km away from Innsbruck, situated about 250 m above the valley bottom. During the measurement campaign near Sonnblick, UV actinic flux was measured at two sites Bodenhaus (Bod) and Kolm-Saigurn (Kol) (see Table 1 and Fig. 2 for details) with two Bentham spectroradiometers. Bod and Kol are typical valley stations situated in a steep valley stretching from the Alpine main ridge northward with Kol beeing furthest into the valley and thus experiencing the greatest obstruction of the horizon.

Table 1

Information on ground-based stations used in this study.

Measurement site (abbreviation)	Instrument	Latitude	Longitude	Altitude
Innsbruck (Inn)	BENTHAM DTM300	47.26428∘N	11.38511∘E	616 m a.s.l.
Lans (Lan)	BENTHAM DM150	47.24544∘N	11.43175∘E	833 m a.s.l.
Hafelekar (Haf)	BENTHAM DTM300	47.31267∘N	11.38386∘E	2275 m a.s.l.
Bodenhaus (Bod)	BENTHAM DM150	47.09947∘N	12.99575∘E	1296 m a.s.l.
Kolm-Saigurn (Kol)	BENTHAM DTM300	47.06822∘N	12.98406∘E	1600 m a.s.l.

Fig. 1

The digital elevation map used within the 3-D-radiative transfer model is shown for Innsbruck section. The three measurement sites are also shown (see text for details).

Fig. 2

Same as 1 but for the Sonnblick section. The two measurement sites Kolm-Saigurn and Bodenhaus and additionally the location of the Sonnblick observatory is marked.

At all stations, Bentham spectroradiometers were used to measure actinic flux from the upper hemisphere (solid angle 2 π) spectrally from 290 nm to 500 nm in 30 min intervalls. The stepwidth is 0.5 nm and the bandwidth (FWHM) is 0.77 nm. The instruments are calibrated regularly to a NIST (National Institute of Standards and Technology) calibrated 1000 W lamp. For further details on the instruments see Erb (1996). We used the instant values at the chosen wavelength from spectra obtained by the Bentham instruments. As the slit function varies between the instruments this has to be taken into account when comparing measurements from different sites, or with the model. The measured data were transformed to represent those from the same instrument but a standard triangular slit function with a 1 nm full width at half maximum (FWHM) using the ShicRIVM algorithm (Slaper et al., 1995). This removes characteristic instrumental differences and standardizes the comparison with model data.

Aerosol vertical distribution was determined by continuous lidar ceilometer observations with a VAISALA LD-40 (Angelini et al., 2009), while sample profiles where acquired by a polarization lidar (Gobbi et al., 2003). Both measurements were performed at Inn station. In addition direct sun irradiance measurements with a sun filter photometer were used to determine the aerosol optical depth on clear sky days. The wavelengths of the filters are 368, 412, 501 and 862 nm. The Ångstrøm coefficient β and the exponent α (Ångstrøm (1929)) are calculated by regression for the wavelength range 368–501 nm. During the days analyzed, aerosol optical depth was rather low, with the highest values being seen on 16th September (Table 2). However, even on 16th September the aerosol effect is rather small and is considered to be almost negligible (~ 5% at 324 nm and less for other wavelengths).

Table 2

Input parameters for the the 3-D-models (altitude, solar zenith angle (SZA), time, total column ozone (ozone), albedo, Ångstrøm exponent and coefficient (α and β) and height of snowline (snowline). sza, azimuth and time are not dependent on the station (Inn - Innsbruck station, Haf - Hafelekar station, Lan - Lans station, Bod - Bodenhaus station, Kol - Kolm-Saigurn station),but on the modeled wavelength since the measurement of a spectrum lasts about 10 min (first value was used for 305 nm run, second value for 380 nm run). There are two values for surface albedo and snowline for the 3-D-model, first for snow-covered areas and second for snow-free areas and first value for the south-facing slopes and second for the north-facing slopes, respectively. Bold values indicate input parameter for both 1-D and 3-D-models, italics (in brackets) indicate input parameter for the 1-D-model.

Date	16.09.2007		24.02.2008			07.05.2008
Station	Inn	Haf	Inn	Lan	Haf	Bod	Kol
Altitude [m]	576(616)	2278(2275)	576(616)	829(833)	2278(2275)	1288(1296)	1595(1600)
SZA [∘]	44.60 - 44.63		56.92 - 56.91			29.66 - 29.72
Azimut [∘]	182.69 - 183.76		173.54 - 173.83			186.04 - 187.47
Time [UTC]	11:17–11:20		11:16–11:17			11:17–11:20
Ozone [DU]	276	271	265	265	260	349	349
Albedo	0.6/0.03(0.14)		0.7/0.03(0.32)			0.7/0.03(0.58)
α	1.4	1.4	1.1	1.1	1.1	1.08	1.08
β	0.060	0.024	0.020	0.016	0.004	0.022	0.028
Snowline [m]	1900/1900		1700/1300			1700/1400

Total column ozone was measured with a Brewer Spectrophotometer at Sonnblick observatory and compared with data from Arosa (Switzerland) and Hohenpeissenberg (Germany). There were no significant differences found in total column ozone between these three stations on the days analyzed. Therefore measurements taken at Sonnblick are considered accurate enough for use in the Sonnblick and Innsbruck area. This is supported by Schmalwieser et al. (2008) who estimated the average difference in total column ozone between Innsbruck and Sonnblick to be on average less then 7%.

To estimate the surface albedo, observations of the snowline were used. The snowline was observed by eye and also classed as either a north- or south facing slope. We also tried to use satellite data to determine snow-covered areas, but this approach did not work well, because of the coarse resolution of the satellite data (500 m) in comparison with the model resolution (50 m). The difference in snowline between south-facing and north-facing slopes is largest on 24th February due to the warm and sunny conditions on the previous days, whereas there is no difference on 16th September (Table 2). For the snow-covered areas a surface albedo of 0.7 (surface partly covered with fresh snow) is assigned. For 16th September we assigned a value of 0.6 to allow for the presence of older snow following suggestions by Simic et al. (2005). For the 7th May and 24th February in the presence of relatively fresh snow the albedo for snow-covered areas was estimated to be 0.7 and for snow-free areas, an albedo of 0.03 was used. The albedo values used come from the findings of Feister and Grewe (1995b). Also it should be noted that Feister and Grewe (1995b) reported only a small wavelength dependence of the surface albedo for polluted snow and therefore the same value is used for all wavelengths.

Modeling

1-D-radiative transfer model SDISORT

The libRadtran software package (http://libradtran.org/doku.php) was used to calculate spectral UV actinic flux and photolysis frequencies. This package consists of a set of programs to perform radiative transfer calculations using the UVSPEC interface (Mayer and Kylling, 2005) to run the SDISORT model (Dahlback and Stamnes, 1991). The SDISORT model is a well tested model used in the UV community being relatively easy to use and fast in terms of computation time. Although 1-D-models cannot take into account the effect of the local terrain, they have produced good results in several intercomparisons, of actinic flux in homogeneous areas (Kylling et al., 2005). These models use numerical methods to solve the radiative transfer equation and are mainly limited by knowledge of the properties of the Earth's atmosphere used in the model run. The SDISORT radiative transfer solver was run here using six streams. The spectral actinic flux was calculated using 1 nm steps and convolved with a triangular slit function with a FWHM = 1 nm. The ATLAS 3 extraterrestrial solar spectrum with 1 nm resolution was used and Sun-Earth distance corrections were applied. The AFGL (Air Force Geophysics Laboratory) mid-latitude profiles (Anderson et al. (1986)) were used for ozone, temperature and air pressure. The ozone profile was scaled according to total column ozone derived from Brewer measurements at Sonnblick. The 532 nm aerosol extinction profile was measured by the polarization lidar and scaled according to the 532 nm aerosol optical depth obtained by interpolation of consecutive sunphotometer aerosol optical depths (Langley method). We ran the 1-D-model using the input information indicated in Table 2. In the model the mountain top station (Hafelekar) was treated as a point above the surface. The surface height used in the 1-D-model is the average surface height of the 3-D-model domain (Fig. 1), e.g. 1386 m, and the height of the point is the actual height (Table 1, e.g. Hafelekar - 2275 m). This is the most realistic approach within a 1-D-model. To estimate 3-D-effects, another calculation with the 1-D-model was performed. We ran the 1-D-model individually for all 2.6 ∗ 105 pixels of the 3-D-model domain using the same surface height and surface albedo as in the 3-D-model. Finally the results were compared with the outcome of the 3-D-model (Fig. 9).

Fig. 9

Ratios of actinic flux calculated with 3-D-model and 1-D-model.

3-D-radiative transfer model GRIMALDI

In comparison to 1-D models, 3-D models reproduce the radiation field more realistically in mountainous terrain as they allow to consider shadowing effects and reflections from tilted surfaces. A disadvantage of 3-D radiative transfer codes is that the computation time is much longer than for 1-D models. In this study, radiative transfer calculations were performed using the 3-D Monte Carlo model GRIMALDI (Scheirer and Macke, 2001). A detailed model description is also given in Scheirer and Macke (2003). GRIMALDI uses a Monte Carlo method for its 3-D-radiative transfer model and in the rest of the paper will be referred as the 3-D-model. Control runs of the 3-D-model without topography show an agreement of ± 5% with the 1-D-model. Furthermore GRIMALDI (note, here referred to as UNIK) took part in an international Intercomparison of 3D Radiation Codes (Cahalan et al., 2005). We applied periodic conditions, with photons leaving the model domain at one side and entering the model domain at the opposite side. Two different model domains were used characterizing the area of investigation of the different field campaigns. The Innsbruck model domain has the size of x = 26 km, y = 25 km, z = 40 km and the Sonnblick model domain x = 21.75 km, y = 24.1 km, z = 40 km. The horizontal resolution was 50 m and the vertical resolution 200 m in both domains. Input data for total column ozone, the aerosol optical depth and vertical profiles of both are identical to the values used in the 1-D-model calculations. For this study two major changes to GRIMALDI model code were made.

First the model terrain was incorporated using a digital elevation map with a horizontal resolution of 50 m (Figs. 1 and 2) which also takes into account the inclination of the surface. Therefore it is possible to model the 3-D-radiation field much more realistically than in the 1-D-approach.

The second major change in the code of the GRIMALDI model is related to the actinic flux calculations. For calculating irradiance within the model it is adequate to count “photons” on a surface, but actinic flux is more complex. Therefore to calculate actinic flux, photons are counted on the surface of a cuboid. To convert photon numbers to actinic flux one has to take the geometry and the surface area of each pixel into account to calculate a “projected surface”. The conversion from photon number to actinic flux is described in detail in Appendix A.

For each wavelength (305, 310, 324 and 380 nm) ten individual 3-D model runs were performed with 108 “photons” using input parameters indicated in Table 2. The final result was calculated using the mean of the ten model runs. With this approach, it is assured, that the numerical uncertainty is well below 5%. As a rule of thumb the CPU-time for one model run with 108 “photons” is about 36 h, therefore we ran the simulations in parallel on a computer cluster with 15 processors and obtained a calculation using 109 “photons” within about 24 h.

Determination of the photolysis frequencies

To obtain photolysis frequencies from spectral actinic flux measurements and 1-D-model calculations the following approach was used. The O3-photolysis frequencies were calculated from spectral actinic flux between 290 nm and 340 nm, NO2-photolysis frequencies are calculated using spectral actinic flux between 290 nm and 420 nm. Photolysis frequencies J are derived by numerical integration of Eq. (1) using 1 nm steps.

Here σ(λ) denotes quantum yield, ϕ(λ) absorption cross section and F(λ) the photon flux. Absorption cross sections and quantum yield data for a temperature of 0 °C were taken from the NASA Panel for Data Evaluation (http://jpldataeval.jpl.nasa.gov, Paur and Bass, 1985). To reduce computation time, we applied a slightly different method for 3-D-model calculations. We did not perform single calculations for each wavelength but instead just one calculation each for photons from two different wavebands (between 290 and 340 nm (O3 - photolysis) and between 290 and 420 nm (NO2 - photolysis)) was performed. For this calculation the contribution of each wavelength to the photolysis frequencies has to be determined in advance. This is achieved by multiplication of the spectral extraterrestrial irradiance with spectral absorption cross section and quantum yield. According to this value (between 0 and 1), 109 “photons” are distributed in 1 nm steps. The maximum impact and thus “photon” number occurs at 290 nm for O3 – photolysis and at 378 nm for NO2 – photolysis. Now these “photons” equally contribute to photolysis frequencies and one can add “photons” from different wavelengths to obtain photolysis frequencies.

Results and discussion

Comparison of model results and measurements

Figs. 3 and 4 show the ratio of model calculations to the measurements. It is shown that the 3-D-model calculations (black filled squares) have a similar site to site variation as the 1-D-model calculations (black open triangles). Further there is a significant variation in the ratios for different stations and on different dates but spectrally the results for each site are self-consistent. This indicates that the correct wavelength dependent model input parameters (e.g. total ozone column) are used in the models. Therefore this implies that changes in the ratios on different days are mainly due to real changes in the snow cover which cause changes in the agreement between the actual and the model snow-cover. The local topography of the station (valley, mountain top) also has a significant impact. The simulations using the 3-D-model are significantly closer to the measurements on 5th May 2008 especially for Bod station. The results suggest for Bod and Kol (both located in a steep valley), that obstructions on the horizon results in a decrease in actinic flux at both stations. The 1-D-model cannot take this phenomenon into account whereas the 3-Dmodel includes the effect of topography and thus more accurately reproduce the measurements. However in the case of Kol, this effect is overestimated in the 3-D-model.

For 16th September 2007 the 1-D-model reproduces the measurements better than the 3-D-model. For this day the 3-D-model significantly underestimates actinic flux for station Inn and slightly overestimates actinic flux for the station Haf. The agreement between 1-D-model results and measurements is surprisingly good and interestingly consideration of topography within the 3-D-model does not improve the results.

On 24th February the 3-D-model generally reproduces the measurements better than the 1-D-model for the valley stations Inn and Lan. The ratios of both models for the mountain top station Haf are similar.

The results suggest that in cases with less snow cover (16th September 2007 and 24th February 2008) the 1-D-model is able to reproduce the measurements to within 20% (Figs. 3 and 4). This result was expected, since multiple reflections on inhomogeneous surfaces are reduced by low surface albedo and the 1-D approximation works quite well. In general the 1-D-model agrees better with measurements at the mountain station (Haf) where the situation is closest to a site with a clear horizon. In situations with significant snow cover and a steep valley (5th May 2008) the 1-D-model does not agree well with the measurements. Under such conditions one expects a large impact of surface scattered radiation, highly influenced by topography, which cannot be properly taken into account in the 1-D-model. For these cases the 3-Dmodel reproduces the interaction between high surface albedo and topography better. However, at Kol the 3-D-model tends to overestimate the obstruction of the horizon, possibly due to an error in estimation of the snowline which was between the two sites Bod and Kol.

Influence of topography and surface albedo

The topography used in the 3-D-model is shown in Figs. 1 and 2 and the spatial distribution of spectral actinic flux for 324 nm is shown for two different days (24th February 2008 and 5th May 2008) and areas (Figs. 7 and 5). Apart from minor boundary problems (e.g. mountain shadow at the Southern boundary in Fig. 7), the typical features of radiation in mountainous areas are well reproduced with the actinic flux increasing with altitude. We observed an increase of about 15% 1000 m at 305 nm and about 5% at 400 nm. However the increase can vary a lot depending both on location and solar zenith angle (see Figs. 1, 2, 5 and 7), with the highest values occuring on south-facing slopes. Due to the high solar zenith angle on 24th February, the steepest parts of the north-facing slopes are shaded, resulting in significantly lower actinic flux values (Fig. 7). This feature is not visible in Fig. 5 where the solar zenith angle is 29.67°. The greatest values in both cases (Figs. 5 and 7) occur on south-facing snow-covered slopes in narrow valleys.

Fig. 7

Same as Fig. 5 but for Innsbruck area on 24/02/2008 11:17 UTC, position of the sun is at 56.92  zenith angle and 198.63 degree azimuth angle.

Fig. 5

Calculated actinic flux F from the 3-D model on 05/07/2008 for 324 nm 11:18 UTC with realistic surface albedo (with snow), position of the sun is at 29.67  zenith angle and 186.52  azimuth angle.

To study the influence of snow in combination with topography on surface actinic flux we performed additional calculations without snow using the 3-D-model. The results show that snow cover in combination with topography can significantly enhance actinic flux in some areas (up to 25%, see Figs. 6 and 8). Surprisingly the effect of snow is almost negligible in the snow-free valley bottom. The snowline within the model was at 1700 m for the south-facing and at 1300 m for the north-facing slopes. Although the distance between valley and snow-covered mountain slopes is small, e.g. about 5 km between Innsbruck and the snow-covered Nordkette (mountains north of Innsbruck), the snow-covered mountains don't significantly increase the actinic flux at the valley site in Innsbruck.

To determine 3-D-effects we run the 1-D-model was run for each pixel of the Innsbruck area 3-D model domain with surface height and albedo being taken from the 3-D-model. The ratio of spectral actinic flux at 324 nm calculated with the 3-D-model to the 1-D-model results is shown in Fig. 9. Red and yellow indicate where the 3-D-model calculates higher values than the 1-Dmodel, green and blue mark regions where 3-D-model results are lower than those of the 1-D-model. The snowline is clearly visible (border between red and green), here surface albedo jumps from 0.03 to 0.7 with altitude resulting in a jump in actinic flux within the 1-D-model of about 40%. In the 3-D-model the change over from snow-free to snow-covered areas results in a much smoother increase of actinic flux. Therefore we observe ratios up to 130% in snow-free areas close to snow-covered areas and ratios of around 80% in snow-covered areas close to snow-free areas. Blue areas with ratios below 60% indicate shaded regions (only diffuse radiation in the 3-Dmodel, but direct and diffuse radiation in the 1-D-model). From these results one could estimate an effective albedo for each pixel. Running the 1-D-model with such an albedo value would result in the same value as the 3-D-model. In this way one could incorporate 3-D-effects in a 1-D-model (and benefit from 1-D-model advantages such as being user friendly and computationally faster). However, at present the differences between the 3-D-model results and the measurements (Figs. 3 and 4) are too large for such an approach.

Overall the results suggest that topography has a large impact on the actinic flux field. Generally a rugged topography tends to reduce actinic flux in comparison with a flat topography. On the other hand the impact of snow on the ground is less pronounced than expected. It seems that the steep mountain sides which have an albedo of 0.03 below 1300 m/1700 m (north-facing and south-facing slopes respectively) act as a barrier to light reflected from the surface.

Conclusions and outlook

In general 1-D-model calculations are not able to reproduce spectral actinic flux and photolysis frequencies in complex Alpine terrain under cloud-free conditions within the desired accuracy of about 10%. Furthermore it is not possible to derive realistic 2-D-surface actinic flux fields with 1-D-models. In this case study the 1-D-model simulations showed a significant overestimation of actinic fluxes. Since crucial atmospheric parameters like total column ozone, ozone profile, aerosol optical depth and aerosol profiles are known with reasonable accuracy, the observed deviation between model results and observations are most likely resulting from 3-D effects. However, the modified 3-D-model was also not able to reproduce actinic flux measurements to within the 10% accuracy. The 3-D-model gave significantly different results from the 1-D-model and in general underestimated actinic flux. Overall the accuracy of the 3-D-model was of the same order as the 1-D-model, we conclude that topography generally tends to significantly reduce actinic flux in situations with moderate snow cover (average surface albedo in our cases was between 0.14 and 0.59), which implies it is very important to take topography into account if one is interested in actinic flux in mountainous terrain. However, a big source of uncertainty remains in the determination of the surface albedo. In our study only two different values for surface albedo were used, one for snow-covered areas and one for snow free areas (Table 2). More detailed surface albedo information (e.g. from land use maps) is needed to improve the 3-D-model calculations. Furthermore one should address the problem of different geometries. The entrance optic used for the measurements is a hemispheric sampling head (2 π response), but within the 3-D-model the “photons” are counted on five surfaces of a cuboid. It would be useful to determine if any significant errors result from this difference. Despite these limitations, our 3-D-model calculations allow us to estimate the potential impact of surface albedo on actinic flux. The results indicate that the influence of surface albedo can vary greatly depending on the area considered from being a very important effect (e.g. in steep valleys), to only a minor influence (e.g. bottom of wide valleys).