Evaluation of Visual and Nonvisual Levels of Daylight from Spectral Power Distributions Considering Orientation and Seasonality

: The evaluation of both visual and nonvisual effects from the spectral power distribution (SPD) of outdoor light is critical in lighting design. The dome-light SPD characteristically changes continuously depending on the seasonality, orientation, altitude of the sun or hour of the day. Traditional photopic parameters, such as the illuminance, luminance or correlated colour temperature (CCT), have been widely studied, but presently, there is no melanopic measurement or evaluation method. This article discusses the processes involved in establishing a simple method to determine the SPD of daylight and solar radiation over the skydome in a location to accurately account for the effects of both photopic and circadian levels around a location. Once per month for one year, natural daylight was spectrally measured in the city of Zaragoza (Spain); radiometric and photometric characteristics were analysed by season; and circadian effects were calculated in terms of standard parameters described by the Commission International de l’Eclairage (CIE), factors recommended by normative and scientiﬁc backgrounds. Finally, we suggest that the best parameter is the melanopic versus photopic irradiance ratio, which achieves reliable results at simplifying and correlating calculations.


Introduction
The sun is the main source of natural light on Earth; and as natural light passes through the atmosphere, its energy is manifested in the different types of radiation along the electromagnetic spectrum, measured from the ground.The spectral power distribution (SPD) of the light over the skydome is rarely uniform.The amount of visible energy that reaches us depends on the atmosphere.Visible energy can be absorbed, reflected and scattered in all directions by gas and dust particles and clouds of water vapour.The characteristics of the entire sky can change within minutes, so overcast, partly cloudy and clear skies are highly dynamic due to winds.The Commission International de l'Eclairage (CIE) defines the Standard General Skies in the CIE S 011:2003/ISO 15469:2004 [1] as a mathematical model based on the variation in three characteristics: sky clearness and brightness, both derived from a comparison of direct normal and diffuse horizontal irradiance (W/m 2 ) and/or illuminance (lux); and solar altitude that can be either calculated for a specific date and time or given as a measured angle above the horizontal plane.This model combines five coefficients, which can then be calculated directly from these three values and the position of the sun within the sky, to describe the spatial distribution of 15 diverse sky types; and the luminance for any point on the skydome can be calculated relative to the zenith illuminance [1][2][3].
The use of irradiance data makes luminance calculations more accurate, as it has been described [4,5].These general sky types represent a range of different sky conditions from overcast to cloudless with varying levels of direct sunlight.They are approximately illustrative of common sky settings with varying levels of cloudiness and turbidity, and they are commonly used in indoor lighting design software.Knowledge of the irradiance distribution of the sky is essential to achieve a good natural lighting design in indoor areas.It is important to know how the characteristics of the sky change in different periods of the year, depending on the climate of the specific place in which we want to conduct the measurements, the types of clouds, the light distribution, the illuminance, the geographical location, the time, etc. Daylight greatly influences the energy efficiency of a building.Despite the advances in artificial lighting, it is still difficult to reduce the energy consumption of large buildings.According to the Office Energy Saving and Efficiency Guide, 40% of the world's energy consumption is produced by buildings.A good estimate of daylight will help reduce electric light consumption.Using the universal method provided by the CIE, the luminance of the sky and thus daylight can be partially predicted.Currently, both irradiance and illuminance levels, mainly at vertical and inclined planes, from the sky are needed for proper building design [3,[5][6][7], but studies seeking to resolve the problem of finding artificial lighting with similar behaviour to the solar SPD have been partially developed.
The correlated colour temperature (CCT) and the colour rendering index (CRI) are parameters used to visually characterize light and lighting strategies, but light reaching the eye level has two important objectives: first, this light is responsible for forming images; and second, this light is part of the exposome [8], meaning it is essential in the synchronization of the circadian system, the system in charge of setting our biological clock on time each day in accordance with the cycles of day and night (light-dark) [9].Regulating the amount of light that reaches the corneal level and improving environmental conditions could have positive impacts on quality of life, wellbeing, and ageing-related concerns.These lighting projects, which consider the potential effects of light on people and optimize them to achieve the greatest wellbeing, are called human-centric lighting (HCL) projects.However, the increasing number of such projects does not correspond to the importance that circadian light should have according to relevant indoor studies.A lack of knowledge has revealed how the characteristics of outdoor daylight are considered to perform indoor lighting designs.Beyond visual stimuli, the CIE S 026/E: 2018 standard [10] includes information concerning the effects of light on intrinsically photoreceptive retinal ganglion cells (ipRGCs) that are related to the SPD, temporal distribution or photoperiod: dynamic light, which is the right light at the right time [11,12].Even the International WELL Building Standard (Visual lighting design and Circadian lighting design) recommends lighting guidelines to provide appropriate photopic and scotopic light levels, minimize disruption to the body's circadian system, and support good sleep quality [13].Currently, lighting projects consider that the circadian rhythm in humans is described as an approximately 24 h cycle controlled by the lighting environment and one of the keys to regulating sleep, mood, heart rate and body temperature [14].Nowadays, people stay indoors with artificial lighting and with relatively limited sunlight most of their time awake, or under the influence of secondary sources such as screens; these environmental contributions should be globally evaluated in circadian terms [15].The CIE divides the visual from nonvisual effects of light and is developing integrative standards from this perspective [10].To fulfil this objective, special attention has to be paid to the SPD.Both the excitation of traditional photoreceptors and ipRGCs are highly wavelength-dependent and form two interdependent pathways for light in the brain [16].To evaluate and quantify the visual and nonvisual influences of light in humans belonging anatomically to these two ocular pathways, the CIE adopted the proposal of Lucas et al. defining that the influences are physically characterized by five equivalent ↵-opic illuminance approaches.This standard employs light measurement methods that quantify the effective irradiance for each of the photoreceptive inputs to the visual system independently.The CIE S 026/E:2018 [10] adopted melanopic equivalent daylight (D 65 ) illuminance, which the Equivalent Daylight Illuminance (EDI) defined as a light source-type D 65 that employs photopic illuminance E photopic,D65 to provide the same melanopic irradiance as a light source with the SPD and photopic illuminance E photopic,SPD .An important issue found for accurate visual and nonvisual daylight implications is the limited information available for this analysis because at least six factors triggering nonvisual light effects have been summarized, including the SPD, the irradiance level, the directionality, the timing, the duration, and the history, which are grouped into luminous and temporal categories [11].Nevertheless, there have been efforts to develop a clear method to evaluate the effects of lighting on circadian rhythm regulation from different perspectives [17,18].
This research hypothesizes that daylighting SPDs modify the photopic and melanopic levels inside buildings depending on the orientation, seasonality and characteristics of the skydome [17,18].The aim of this study is to analyse these real seasonal changes in natural sunlight to evaluate how changes in the SPD affect photometric quantities and calculate the variations in circadian magnitudes based on the CIE standard.Knowledge of the outdoor daylight is essential to achieve a good natural lighting design in indoor areas from a wellbeing perspective.Differences in daylight with geographical location of a building, façade orientations or seasonality should be analyzed and considered to evaluate the nonvisual potential of daylight for indoor lighting purposes.This analysis leads us to evaluate the main parameters that should be considered by indoor designers to simulate daylight conditions with artificial lighting promoting people's wellbeing.

Spectral Measurements
The methodology of this research was based on the analysis of the radiometric behaviour of the daylight at a location, which serves to model their circadian contribution at this same point.We conducted several spectral measurements of the skydome over one year from December 2018 to November 2019.We sought a sunny day around the 15th of each month, and measurements were performed from 10:00 to 11:00 a.m.The experiment was located at the roof of the Sciences Faculty at the University of Zaragoza in Spain (41 38 0 31 00 N, 0 53 0 60 00 W; 243 m above sea level).The SPD measurements were performed with a calibrated spectroradiometer (Avaspec-1024, Avantes, Apeldoorn, The Netherlands, NPL E01110063/DDK calibration and NIST traceability).A tripod was used to ensure accurate measurements of the spectra in the four cardinal directions (north, south, east and west) 15 , 30 , and 45 from the horizontal plane.To register information about the weather, photos of the sky were taken.The sky conditions were always clear except intermediate skies with visible sun and a few clouds.

Theoretical Considerations
The SPD (µW/cm 2 ) from 300 nm to 1000 nm, the CCT, and the colour coordinates x and y of the selected day were recorded for further processing.The photopic irradiance W for photopic vision, SPD(l) is the spectral power distribution of the lighting source and V(l) is the CIE photopic spectral luminous efficiency function) [10,19].The melanopic irradiance E e,melanopic W/m 2 = R 780 l=380 SPD(l)S mel (l)dl was analogously calculated (S mel (l) is the CIE melanopic spec- tral weighting function) [10].The ratio E e,melanopic E e,photopic is named the Melanopic Action Factor (MAF) [20], following the Circadian Action Factor (CAF) proposed by Gall [21], based on spectral irradiance behaviour and was used to propose a method to evaluate the nonvisual potential of daylight.In a previous publication by this research group, the conversion factors from photopic values of illuminance to the CIE melanopic standard and to other recommendations such as WELL were calculated in detail [20].To understand these conversions, a summary of the equations is added.
The normalization imposed by the CIE standard with melanopic illuminance was provided by the standard illuminant D 65 (daylight CCT = 6500 K) in accordance with the EDI.A light-source type D 65 furnishing photopic illuminance E e,photopic,D65 to provide the same melanopic illuminance E melanopic,D65 can be calculated.
If the photometric quantity is defined by E photopic,D65 = Km E e,photopic,D65 [19], the same melanopic E melanopic,D65 = Km melanopic,D65 E e,melanopic,D65 for light source D 65 can be described in terms of an equality as follows: E photopic,D65 = E melanopic,D65 .For any source: In addition, the relation of the EDI with the photometric illuminance can be deduced as follows [20]: The Equivalent Melanopic Lux (EML) is defined by the WELL recommendations [13] with reference to the equi-energy illuminance and can be calculated for any source, in terms of melanopic parameters; the EML is defined as follows [20]: In addition, the relation of the EML to photometric illuminance is described as follows: Using these conversion factors and photopic values of the illuminance, transformations from one melanopic metric to another are easily calculated by EML = 1.104EDI.

Results and Discussion
Outdoor SPD measurements in Zaragoza, Spain (latitude: 41.6370082 , longitude: 0.9078846 ) are provided facing the four main cardinal points (north, east, south and west) 15 , 30 , and 45 from the horizontal plane.The seasonal irradiance (µW/cm 2 ) from 380 to 780 nm, luminance (cd/m 2 ), illuminance (lux), and CCT (K) were calculated from the SPD measurements to estimate variations depending on altitude and orientation.Calculations of nonvisual effects using the previously described metrics were performed to analyse values, compare current melanopic evaluation methods and estimate the variability ranges of daylight.

Sky Characteristics and Photopic Descriptions
The outdoor solar conditions of the study during the experiment were registered (Table 1), and the SPDs were measured (Figure 1).
The highest irradiances can be found in the east direction at 30 -45 in practically every month except in December, January and October.In those three months, at 15 , the irradiance is higher than that at 45 .This is followed by the south direction, except in January and November.In those two months, at 15 , the irradiance is higher than at 45 .The lowest irradiances are in the west orientation except in five months (February, March, September, October and November), which slightly exceed the values in the north direction (Table 1 and Figure 2).From April to September, the sun is located in the east and at the highest latitudes (Table 1 and Figure 2).In the normalized SPD curves (Figure 1), spectral variations were found throughout the months of the year and depended on the cardinal points at which the measurements were taken.In most of the data, higher CCTs were found in the north and west orientations than in the east and south orientations and in the winter months (from December to March) and some of the autumn months, especially in October (Figure 3).In general, the CCTs are higher in measurements taken at 45 and lower in those taken at 15 .The average CCT was 6122 K; the maximum value was taken in March, facing north at 45 (12,380 K); and the minimum was taken in December, facing south at 15 (4892 K).In general, the CCTs where the sun is located (south and east) always comprised warmer restricted ranges irrespective of the season and the sky altitude.Conversely, in the opposite locations (north and west), the CCTs are altitude-dependent; and the higher the CCT is, the cooler the results that are found.The measured CCTs show variations as a function of the season, orientation and altitude in correspondence with CIE D illuminances, even with clear sky conditions [22,23].When the value of the sky is considered, daylight approaches 6500 K.The average value of the data taken in this study is 6122 K, similar to the dominant CCT of the study conducted in the city of Granada, Spain (5700 K).However, the recorded data show there were variations throughout the year, which highlights the importance of studying the parameters of the daylight in each type of sky and in different climatic conditions [24].Although the CCT measurement gives us information mainly on the appearance of colour and on aesthetic perception, its influence on circadian effects can also be considered.
Our experiment was performed in an urban environment, and it did not consider either the level or the components of the atmosphere.The presence of aerosol particles at certain times of the year could influence the presented SPD measures and, consequently, the subsequent calculations [25].Peyvandi et al. [26] calculated the 25% and 75% interquartile values of the CCT to find the range containing 50% of the CCTs for each condition, obtaining a range for clear sky without aerosol particles of 5889-6294 K, a range for clear sky with aerosol particles of 5712-6325 K, and a range for overcast sky of 6218-7757 K.The overcast sky and atmosphere with aerosol particles represented extreme conditions where the limits of the CCT were observed.Daylight measured in Granada (Spain), with the most observed CCTs of 5555-5715 K [24], and in Boulder, Colorado (United States), with obtained CCTs of 5500-6400 K, showed that lower CCTs represented atmospheres with higher aerosol particle levels [23,24].Similar results were found in Rochester, New York (United States), where daylight with CCT as low as 3530 K in a hazy sky was reported [27].In Zaragoza, the mean CCTs obtained at 45 were 7574 K to the north (range 5660-12,380 K), 5377 K to the south (range 4954-5671 K), 5294 K to the east (range 5085-5543 K), and 8025 K to the west (range 5883-9587 K), irrespective of the presence of clouds, mainly due to the season and showing a broad range of daylight CCTs.Comparing our results with those previous results show very low rates of aerosol elements during all the seasons of the evaluated year.
Luminance patterns of the 15 CIE standard skies have been described and used together with several skydome parameters including different illuminance levels or vertical sky components to classify daylight [1].Sky luminance values have a strong dependence on its particular conditions such as clouds, which can vary in minutes; or the locations of measurement points that can have different conditions, even when they are nearby [28].Therefore, the same categories of sky conditions, including clear sky, partly cloudy and overcast, cannot be considered to have similar luminance levels.We chose similar clear skies to perform our experiment, but the location of the sun, partially cloudy skies, and even the buildings of the surroundings highly influenced the results.At 15 , values elevated enough to be considered to have been influenced by the buildings located around the measurement station can be observed (Figure 4).Despite this result, the highest luminance values were found at 45 at the nearby positions where the sun was placed during the measurements.Clouds influenced the measurements, as reflected in the September data, due to the presence of some clouds during the measurement (Table 1).
Higher levels of photopic illuminance (Figure 5) are reached in the east orientation from March to August and in the south direction from September to March compared to the north and west orientations.The months with the highest photopic illuminance are generally in the summer months and in part of the spring months (from May to August), decreasing in the autumn and winter months in the east at the three elevations (15 , 30 and 45 ).However, in the south orientation, many differences in photopic illuminance are seen between the three elevations in the spring and summer months (from April to August), which become equal in the autumn and winter months (from September to March).These differences in the three elevations in spring and summer also occur in the north (from May to August).This does not happen in the west orientation, where the photopic illuminance at the three elevations is very equal during all months of the year.The minimum photopic illuminance was found in the month of January in the west orientation with an angle of 45 (5084 lux).The maximum value was found in the month of August in the east orientation with an angle of 45 (99,822 lux).In general, the lowest photopic illuminances were obtained for measurements taken at 15 , and the highest values were obtained for measurements taken at 45 .Illuminance measurements reveal the changing position of the sun with seasons.Regarding the measurement locations, Figures 1 and 5 show that the illuminance level reached higher values while it was positioned 45 in the southeast during autumn-winter and northeast in spring-summer.

Circadian Characteristics
To assess the circadian characteristics, we began by analysing the parameter MAF (Table 2).
The highest MAFs occur in the north and in the west with a 45 elevation from October to March.Specifically, the highest MAF is 1.156, which is located in the north in March at 45 .The most uniform MAFs occur in the east during all months of the year and at the three elevations.The lowest MAF is 0.759, located in the south in December at 15 .
The mean MAF of all calculated values for the north, south, east and west and for the three elevations (mean ± standard deviation) is 0.871 ± 0.089 with an annual variation of less than 10%.That is, the MAFs are very close to 1, which means that the evaluated daylights have almost equal photopic and melanopic intensities, which will generate the same image-forming and non-image-forming effects.
According to the CIE calculations, higher levels of melanopic EDI are reached (Figure 6) in the east and south orientations compared to the north and west orientations.The months with the highest EDIs were generally the summer months, and the minimum EDI was found in December facing west at an angle of 15 (5656 melanopic lux).The maximum value was taken in August in an east orientation with an angle of 45 (88,714 melanopic lux).In general, the lowest melanopic illuminances have been obtained in measurements taken at 15 , and the highest values have been obtained in measurements taken at 45 .Both the photopic illuminances and EMLs reach higher illuminances in the east and south orientations and in the summer months, as is the case with the EDIs that can be calculated according to the relation EML = 1.104EDI.The minimum illuminances, both photopic and melanopic, were taken in December at a west orientation with an angle of 15 .The maximum illuminances, both photopic and melanopic, were taken in August at an east orientation with an angle of 45 , as has been previously described for EDI due to their proportionality relationship.
Bella et al. [17] modelled the nonvisual potential of daylight in terms of the melanopic/ photopic ratio calculated from the SPD, finding values very close to 1 and similar to those calculated by using the D 40 to D 120 illuminances.Although this definition differs from our MAF, both indicate practically the same melanopic contribution as the photopic for daylight, and the spectra could be considered to be covering similar ranges to our CCTs.
Other authors have developed circadian theories.One of the most widespread is Rea et al. who defined the circadian stimulus (CS) and the CL A , parameters directly proportional to the levels of night-time light-induced melatonin suppression [29].This method is not based on any standard but has been adopted by most researchers to define the circadian effect of light [18].The highest CS, 0.699, is found in the south and east orientations than in the north and east orientations, especially in the winter months (Figure 7).The lowest CS found is in the month of December facing the west at 15 with a value of 0.685.Given the outdoor SPDs with higher short wavelength contributions than longer ones, the corresponding CSs are always very close to the 0.7 saturation value in the studied period during the evaluated hours.Our calculation shows a CS (mean ± standard deviation) of 0.695 ± 0.004 throughout the evaluated year for all orientations and elevations with a variation of <1%.There are higher CL A values in the east followed by the south orientation, and the lowest values are in the west followed by the north orientation (Figure 8).The measurements taken at 45 are generally the highest values compared to those measured at 30 and those measured at 15 , which have the lowest values.The highest values measured in the east occur in spring and summer (from March to August) at the three elevations, and the lowest occur in the west in winter (from December to March).Our calculations show a CL A (mean ± standard deviation) of 77,935 ± 62,347 circadian lux (range 11,665-215,042 circadian lux).It should be noted that our obtained values of CL A for daylight are proportional to the EDIs, as shown by the analysed data, with CL A = (2.3 ± 0.1) EDI.Therefore, based on this relationship and for cool daylight, CL A , EDI and EML are parameters proportional to each other.
According to the mathematical model based on spectral sensitivity data proposed by Rea et al. [29], CL A provides information on the SPD at the corneal level converted into circadian light [30,31]; and the CS, based on the CL A , is equal to the percentage of melatonin suppression [32], ranging from threshold (CS = 0.1) to saturation (CS = 0.7) [33].The CS seems to be necessary to promote circadian entrainment.As our mean CS obtained is 0.695 ± 0.004 with a variation of 1%, it could be considered sufficient for the described saturation.Therefore, CS could be used directly to quantify the circadian impact of lighting.The values calculated in our study show that there would be a high percentage of melatonin suppression independent of the seasonality, orientation and elevation.Studies show that a CS equal to or greater than 0.3 (equivalent to 30% melatonin suppression) in the first part of the day can improve the circadian cycle and the quality of sleep, reduce drowsiness and increase vitality and alertness [31,[33][34][35], which coincides with our data collected from 10:00 to 11:00 AM.
A value of CS = 0.35 was proposed as a sufficient and adequate value to promote the daily cycle [36].According to Rea, the daylight spectrum has been shown to be relatively efficient in providing circadian stimulation; and based on the CIE D 65 spectrum, a corneal illuminance of 233 lx corresponds to this 0.35 CS.Extrapolating our outdoor analysis to indoor situations, circadian efficiency should be treated with caution because artificial lighting seems to be less effective since for the same value of CS = 0.35, a corneal illuminance much higher than that necessary with natural light would be needed [37].All the outdoor SPDs collected for this paper are, regarding the CCT, above 4800 K, meaning that for indoor lighting purposes, it would be equivalent to cold light; and, assuming the same light intensity, higher CCT values indicate a greater perception of luminosity.This fact could improve visual performance and the relative effectiveness of light sources with higher melanopic intensity.This situation will be favourable for improving indoor performance in the morning hours; however, in the afternoon and at night, it will be more convenient to use warm lighting with low melanopic intensity [9,15].In an investigation of a method to evaluate the nonvisual potential of daylight in two European cities, CCTs between 4631 K and 11,871 K were recorded in Naples and CCTs between 3585 K and 9812 K were recorded in Bialystock [17].These results are similar to those of this study since our CCT measurements vary between 4892 and 12,380 K.However, due to the low correlation between the CCT and circadian effects, it is considered that the CCT is not a good indicator to establish the levels of circadian response [20,37].Although Rea et al. reported that LEDs with high CCTs above 5299 k could slightly stimulate the circadian system after one hour of exposure to outdoor lighting at 95 lux, when taking advantage of natural light in interior spaces, the reflectance of the walls must also be considered because they vary the CCT and the spectral composition of light that regulate the circadian system [38][39][40].Additionally, our measurements show that for indoor purposes where daylight reaching outdoors is important, new parameters, such as the seasonality, façade orientation and altitude of the sun, are substantial and have sufficient contributions that should always be considered.

Conclusions
The CIE Standard General Skies models the SPD of daylight according to the omnidirectional diffuse light reflected and refracted from the entire skydome and not just the highly directional rays from the sun.This standard makes it possible to derive coefficient values directly from measured instantaneous sky illuminance or irradiance or to adjust them independently to simulate and to match any sky conditions.Our benchmarking study reveals the differences in daylight SPDs, with the orientation, altitude and seasonality, on photopic and melanopic estimations, which could potentially modify occupants' responses in outdoor and indoor scenarios.Knowledge of the spectral and spatial patterns of daylight is needed to develop simple outdoor circadian models and procedures to apply them, especially when they should be robust to indoor lighting projects.We have found that based on the CIE recommendations, intercomparison among metrics with high correlation has potential.In addition, the results presented in our study demonstrate that despite the significant variations in daylight SPD in a specific location, it is possible to evaluate the nonvisual potential of daylight based on the analysis of these spectra using the MAF parameter and traditional photopic illuminance.
In this study, the skies seasonality of Zaragoza during a year in terms of the SPD was studied.Selected clear skies have the highest contribution to both photopic and melanopic illuminances, but data relating to seasonality, sun altitude and glazing orientation in buildings have been measured.The results show that these factors should be considered when outdoor or indoor lighting projects are implemented.
As a limitation, further long-term SPD measurements are required for the annual study of a specific local daylight contribution.This could contribute to the establishment of general daylight indoor design rules based on real scenarios that consider the total information on the light that exists in the physical environment of each city.Information is lacking on how the daylight SPD may vary depending on geographic location and weather conditions; and another important aspect to consider is that due to the position of the sun vault in the skydome, the SPD is different on different vertical planes depending on both season and orientation.Having access to this information will allow the development of better guidelines not only for interior lighting design but also for the management and control of some biomedical questions such as dermatological [41,42] and psychiatric pathologies [43] or the improvement of some ocular dysfunctions in the fields of optometry and ophthalmology [44], among others.Finally, it was demonstrated that despite the dynamic variations in daylight and the corresponding variations in circadian values, for each orientation, altitude and season, it is possible to easily calculate a value describing the nonvisual potential of the daylight incidence at the building façade from photopic illuminance and knowing the SPD of the incident light.

Figure 8 .
Figure 8. CL A values calculated depending on the season (from December 2018 to November 2019); orientation of north (top), east (right), south (bottom), and west (left); and elevation from the horizontal level of 15 (black), 30 (striped), and 45 (grey).

Table 2 .
Melanopic Action Factor (MAF) values calculated depending on the season (from December 2018 to November 2019); orientation of north, east, south, and west; and elevation from the horizontal level of 15 , 30 , and 45 .