Photon flux density prediction system and photon flux density prediction program

Abstract

To provide a photon flux density prediction system for predicting photon flux density from sunlight.SOLUTION: A photon flux density prediction system includes: a storage unit that associates and stores date and time information, location information, and sunlight intensity information; a communication unit for communication; a first calculation unit that calculates sunlight intensity information associated with sunlight inquiry information included in prediction information inquiry information input via the communication unit; a second calculation unit that calculates direction-specific sunlight intensity information using the calculation result of the first calculation unit; and a third calculation unit that calculates photon flux density prediction information for predicting the photon flux density corresponding to the prediction information inquiry information input via the communication unit, using the calculation result of the second calculation unit. The communication unit transmits the photon flux density prediction information calculated by the third calculation unit to a transmission source of the prediction information inquiry information.SELECTED DRAWING: Figure 1

Description

[Technical Field]

[0001] The present invention relates to a photon flux density prediction system and a photon flux density prediction program. [Background technology]

[0002] Recently, the effects of sunlight (e.g., ultraviolet light, visible light, and infrared light) on the human body and living environment have become known, and attention has begun to be paid to the amount of sunlight exposure (e.g., ultraviolet light exposure and solar radiation). For example, in the exposure amount estimation system described in Japanese Patent No. 5524741, location information indicating a location and environmental information indicating the amount of an exposure target, such as ultraviolet light, present at that location are associated and stored in an environmental information storage unit, the exposure amount is estimated based on the environmental information acquired from the environmental information storage unit and an exposure rate determined according to the user's behavior, etc., and the estimated exposure amount is provided to the user. [Prior art documents] [Patent documents]

[0003] [Patent Document 1] Patent No. 5524741 Summary of the Invention [Problem to be solved by the invention]

[0004] Incidentally, the exposure estimation system described in Patent Publication No. 5524741 estimates exposure using environmental information (pollen dispersion amount, ultraviolet radiation amount, aerosol amount) identified by location information stored in an environmental information storage unit. However, there is a problem in that the user receiving the information cannot easily obtain more detailed information using this environmental information identified by location information.

[0005] For example, in the real world, the parts and locations where one wishes to know the energy intensity of ultraviolet radiation, solar radiation, etc. are not necessarily horizontal or perpendicular to the sun's rays, but are various. In order to know the energy intensity of such various parts and locations, it is conceivable to calculate the energy intensity on the actual irradiated surface from the energy intensity on the horizontal or perpendicular to the sun's rays using trigonometric functions, etc. However, because the energy intensity from the real sky that is scattered and reflected by atmospheric components differs depending on the direction, there is a problem in that calculations using trigonometric functions, etc., do not provide sufficient accuracy.

[0006] It has been known that sunlight affects photon flux density. Furthermore, photon flux density affects plant photosynthesis, so there is a demand for predicting photon flux density. However, no system has been provided that can predict photon flux density from sunlight.

[0007] The present invention has been made to solve the above-mentioned problems, and has an object to provide a photon flux density prediction system that predicts photon flux density from sunlight. [Means for solving the problem]

[0008] In order to solve the above-mentioned problems, the present invention provides a photon flux density prediction system, comprising: a storage unit that stores date and time information, which is information about a date and time; location information, which is information about a location; and sunlight intensity information, which is information about sunlight intensity at the date and time indicated by the date and time information and the location indicated by the location information, in association with each other; a communication unit that performs communication; a first calculation unit that calculates the sunlight intensity information associated with sunlight inquiry information included in prediction information inquiry information input via the communication unit; a second calculation unit that calculates direction-specific sunlight intensity information using a calculation result of the first calculation unit; and photon flux density prediction information that predicts a photon flux density for the prediction information inquiry information input via the communication unit using a calculation result of the second calculation unit. , by the number 5and a third calculation unit that calculates the sunlight ray inquiry information, wherein the sunlight ray inquiry information includes the date and time information, the location information, and direction information that indicates the direction in which an irradiated surface that is irradiated with sunlight faces, and the direction-specific sunlight ray intensity information is information that includes the intensity of sunlight received by the irradiated surface that is indicated by the direction information at the date and time indicated by the date and time information and the location that is indicated by the location information, and the communication unit transmits the calculation result of the third calculation unit to the sender of the forecast information inquiry information. [Effects of the Invention]

[0009] According to the present invention, it is possible to provide a photon flux density prediction system that predicts photon flux density from sunlight. [Brief explanation of the drawings]

[0010] [Figure 1] 1 is a block diagram showing a configuration of a photon flux density prediction system according to a first embodiment of the present invention. [Figure 2] FIG. 2 is a diagram illustrating the irradiation direction of sunlight. [Figure 3] 1 is a diagram illustrating sunlight irradiating an irradiation surface A that is a surface parallel to the ground. [Figure 4] 10 is a diagram illustrating sunlight irradiating an irradiation surface B, which is a surface that forms an angle of 30° with the ground. [Figure 5] 10 is a diagram illustrating sunlight irradiating an irradiation surface C that is a surface that forms an angle of 90° with the ground. [Figure 6] 10 is a diagram illustrating sunlight reflected by irradiation surface E and irradiated onto irradiation surface D. FIG. [Figure 7] 7A and 7B are diagrams illustrating an example of how to determine the spectral reflectance of another irradiated surface, in which FIG. 7A shows how sunlight directly received by irradiated surface E is measured, and FIG. 7B shows how reflected light from irradiated surface E is measured. [Figure 8] 10A and 10B are diagrams illustrating light reflected from other irradiation surfaces that enters the irradiation surface. [Figure 9]2 is a diagram showing an example of information stored in a storage unit 13 shown in FIG. 1. FIG. [Figure 10] 2 is a flowchart showing the operation of the photon flux density prediction system 10 shown in FIG. [Figure 11] FIG. 10 is a block diagram showing the configuration of a photon flux density prediction system according to a second embodiment of the present invention. [Figure 12] 12 is a flowchart showing the operation of the photon flux density prediction system 100 shown in FIG. [Figure 13] 1 is a graph showing spectral irradiance. [Figure 14] FIG. 10 is a diagram showing the conditions for calculating the temperature rise of an irradiated material. [Figure 15] FIG. 10 is a diagram showing an example of a trial calculation of the temperature rise after one hour depending on the material. [Figure 16] Figure 16(a) is a graph showing an example of spectral illumination intensity from each direction on a sidewalk, and Figure 16(b) is a graph showing an example of spectral reflectance on a sidewalk. [Figure 17] Fig. 17(a) is a graph showing an example of spectral irradiance from each direction on a grassland, and Fig. 17(b) is a graph showing an example of spectral reflectance on a grassland. [Figure 18] Figure 18(a) is a graph showing an example of spectral irradiation intensity from each direction on asphalt, and Figure 18(b) is a graph showing an example of spectral reflectance on asphalt. [Figure 19] FIG. 1 shows the relationship between CO2 respiration rate and light intensity. [Figure 20] FIG. 1 is a diagram showing the relationship between photosynthetic rate and light intensity for each plant. [Figure 21] FIG. 1 is a diagram showing the light energy required for photosynthesis. [Figure 22] FIG. 1 shows the dependence of photosynthetic rate on light intensity and temperature. [Figure 23] FIG. 1 shows the dependence of photosynthetic rate on light intensity and temperature. [Figure 24] FIG. 10 is a diagram showing conditions in an example of prediction by the photon flux density prediction system of the present invention. [Figure 25]25 is a graph showing solar radiation energy calculated by the photon flux density prediction system of the present invention under the conditions of FIG. 24. [Figure 26] 25 is a graph showing the photon flux density calculated (calculated using Equation 3) by the photon flux density prediction system of the present invention under the conditions of FIG. 24. [Figure 27] 25 is a graph showing the photosynthetic photon flux density calculated (calculated using Equation 5) by the photon flux density prediction system of the present invention under the conditions of FIG. 24. [Figure 28] 25 is a graph showing the measured values ​​of photosynthetic photon flux density under the conditions of FIG. 24. [Figure 29] 25 is a graph showing a combination of solar radiation energy (taking into account clear skies and cloud cover), photon flux density (taking into account clear skies and cloud cover), and photosynthetic photon flux density (taking into account clear skies and cloud cover) calculated by the photon flux density prediction system of the present invention under the conditions of FIG. 24, as well as actual measured values ​​of photosynthetic photon flux density. [Figure 30] FIG. 25 is a diagram showing the amount of carbon fixation predicted on a fine day by the photon flux density prediction system of the present invention under the conditions of FIG. 24. DETAILED DESCRIPTION OF THE INVENTION

[0011] The photon flux density prediction system according to the present invention will be described in detail below with reference to the drawings. The embodiments shown below are preferred examples of the system according to the present invention, and may include various limitations based on typical hardware and software configurations. However, the technical scope of the present invention is not limited to these aspects unless otherwise specified. Furthermore, the components in the embodiments shown below can be appropriately replaced with existing components, and various variations, including combinations with other existing components, are possible. Therefore, the description of the embodiments shown below does not limit the content of the invention described in the claims.

[0012] In the following examples, the present invention is applied to a photon flux density prediction system to provide a user with information about sunlight, such as solar radiation, but the present invention may also provide a user with individual information about ultraviolet rays, infrared rays, visible light, or other electromagnetic waves contained in sunlight. In the following examples, sunlight intensity is also referred to as solar radiation intensity. [Example]

[0013] 1 is a block diagram showing the configuration of a photon flux density prediction system according to a first embodiment of the present invention. A photon flux density prediction system 10 of this embodiment is, for example, a server machine made up of a computer.

[0014] As shown in FIG. 1 , the photon flux density prediction system 10 includes a storage unit 13 that stores sunlight intensity information 16 (described in detail below) and various other information; a communication unit 17 that performs communication; a first calculation unit 11 that calculates sunlight intensity information 16 associated with sunlight inquiry information included in prediction information inquiry information received via the communication unit 17; a second calculation unit 12 that calculates direction-specific sunlight intensity information using the calculation result of the first calculation unit 11; and a third calculation unit 21 that calculates photon flux density prediction information that predicts the photon flux density corresponding to the prediction information inquiry information input via the communication unit 17 using the calculation result of the second calculation unit 12. The communication unit 17 communicates with, for example, an external device. The prediction information inquiry information includes sunlight inquiry information. The sunlight inquiry information includes date and time information 14, location information 15, and direction information 16a that indicates the direction in which the irradiated surface is irradiated by sunlight, i.e., the sun's radiation. The direction-specific sunlight intensity information is the intensity of sunlight received by the irradiated surface indicated by direction information 16a at the date and time indicated by date and time information 14 and the location indicated by location information 15, and the communication unit 17 transmits the calculation result of the third calculation unit 21 to the sender of the prediction information inquiry information. Each component shown in FIG. 1 may be configured as hardware. Furthermore, each component shown in FIG. 1 may be realized by the photon flux density prediction system 10 executing a program, and the memory unit 13 may store the program to be executed by the photon flux density prediction system 10. The memory unit 13 has a volatile storage device or a non-volatile storage device depending on the use of the data.

[0015] The communication unit 17 may be configured to transmit direction-specific sunlight intensity information indicating the direction-specific sunlight intensity calculated by the second calculation unit 12 to the sender of the forecast information inquiry information.

[0016] In the photon flux density prediction system 10, the communication unit 17 may be configured to transmit the sunlight intensity information calculated by the first calculation unit 11 to the sender of the prediction information inquiry information. In this case, the client machine 2 or 3 may have a configuration corresponding to the function of the second calculation unit 12. That is, the present invention is a photon flux density prediction system in which a server machine and a client machine are connected via a network, and the server machine includes: a storage unit 13 that stores date and time information 14, which is information about a date and time; location information 15, which is information about a location; and sunlight intensity information 16, which is information about the date and time indicated by the date and time information and the sunlight intensity at the location indicated by the location information, in association with each other; a first communication unit that communicates with the client machine; and a first calculation unit 11 that calculates the sunlight intensity information 16 associated with the sunlight inquiry information included in the prediction information inquiry information input via the first communication unit. The first communication unit transmits sunlight intensity information 16 calculated by the first calculation unit 11 to the client machine. The client machine includes a second communication unit that communicates with the server machine 1 and a second calculation unit (corresponding to the function of the second calculation unit 12) that calculates direction-specific sunlight intensity information using the sunlight intensity information 16 input via the second communication unit. The second communication unit (corresponding to the function of the second calculation unit 12) transmits sunlight inquiry information to the server machine. The sunlight inquiry information may include date and time information 14, the location information, and direction information 16a indicating the direction in which the irradiated surface facing the sunlight is facing. Note that the client machine and the server machine may be located in the same terminal device without a network. In this case, the client machine and the server machine may be the same device, or may be connected via a bus, for example. When the client machine and the server machine are the same device, connecting the client machine and the server machine may refer to data exchange between a program that realizes the functions of the client machine and a program that realizes the functions of the server machine via a memory or the like.

[0017] The memory unit 13 stores date and time information 14, which is information about the date and time, location information 15, which is information about the location, and sunlight intensity information 16, which is information about the sunlight intensity at the date and time indicated by the date and time information 14 and the location indicated by the location information 15, in association with each other.

[0018] The communication unit 17 of the photon flux density prediction system 10 is connected to a network 4 such as the Internet. Client machines 2 and 3 used by users are connected to the network 4, and the client machines 2 and 3 communicate with the photon flux density prediction system 10 via the network 4.

[0019] The client machine 2 or 3 transmits prediction information inquiry information to the photon flux density prediction system 10 via the network 4. The prediction information inquiry information includes sunlight inquiry information. The sunlight inquiry information includes date and time information 14, location information 15, and direction information 16a indicating the direction of the irradiated surface irradiated with sunlight. The direction indicated by the direction information 16a is a direction perpendicular to the plane in which the irradiated surface irradiated with sunlight extends. The direction information 16a will be described later with reference to FIGS. 2 to 5. The calculation result by the first calculation unit 11 is transmitted from the photon flux density prediction system 10 to the client machine 2 or 3, which is the sender of the prediction information inquiry information, via the network 4.

[0020] The first calculation unit 11 calculates sunlight ray intensity information 16, which will be described in detail later, using the date and time information and location information included in the sunlight ray inquiry information. The second calculation unit 12 calculates direction-specific sunlight ray intensity information, which is information including the intensity of sunlight received by the irradiated surface indicated by the direction information 16a, using the calculation result by the first calculation unit 11. The communication unit 17 transmits the direction-specific sunlight ray intensity information, which indicates the direction-specific sunlight ray intensity calculated by the second calculation unit 12, via the network 4 to the client machine 2 or the client machine 3, which is the sender of the forecast information inquiry information.

[0021] The direction information 16a will be described below. Fig. 2 is a diagram illustrating the irradiation direction of sunlight. Fig. 3 is a diagram illustrating sunlight irradiating illumination surface A, which is a surface parallel to the ground. Fig. 4 is a diagram illustrating sunlight irradiating illumination surface B, which is a surface that forms an angle of 30° with the ground. Fig. 5 is a diagram illustrating sunlight irradiating illumination surface C, which is a surface that forms an angle of 90° with the ground. In Figs. 3, 4, and 5, the direction in which the illumination surface that is irradiated with sunlight faces, i.e., the direction indicated by the direction information 16a, is referred to as the "direction of the illumination surface."

[0022] As shown in Figure 2, the sunlight irradiating the Earth consists of direct solar radiation irradiated from the direction of the sun, as well as diffuse solar radiation irradiated from directions other than the direction of the sun. Diffuse solar radiation is irradiated onto the irradiated surface from the entire sky. Irradiated surface A, which is a surface parallel to the ground, is irradiated with both diffuse solar radiation and direct solar radiation from the entire sky, as shown in Figure 3.

[0023] Furthermore, for irradiated surface B, which forms an angle of 30° with the ground, diffuse solar radiation and direct solar radiation are irradiated from the direction of the sky toward which irradiated surface B faces, as shown in Figure 4. For irradiated surface C, which forms an angle of 90° with the ground, the sun is facing the backside of irradiated surface C, so direct solar radiation is not irradiated, but diffuse solar radiation is irradiated from the direction of the sky toward which irradiated surface C faces, as shown in Figure 5. In addition to the diffuse solar radiation and direct solar radiation that directly irradiates the irradiated surface, sunlight that is reflected by the ground and then irradiates the irradiated surface.

[0024] 3, 4, and 5, the intensity of the irradiated sunlight varies greatly depending on the direction the irradiated surface faces. Therefore, in this embodiment, different direction-specific sunlight intensity information is calculated for each direction the irradiated surface faces, and this information is provided to client machine 2 or client machine 3, which is the sender of the forecast information inquiry information.

[0025] In this embodiment, in addition to diffuse solar radiation and direct solar radiation, sunlight reflected by other irradiation surfaces such as the ground or a wall surface and then irradiated onto the irradiation surface can also be taken into consideration. This point will be explained with reference to Fig. 6. Fig. 6 is a diagram explaining sunlight reflected by irradiation surface E and then irradiated onto irradiation surface D.

[0026] In an actual environment, the amount of sunlight received by the irradiated surface D includes not only the sunlight (direct sunlight and diffuse sunlight) directly received by the irradiated surface D, but also sunlight (reflected sunlight) that is sunlight (direct sunlight and diffuse sunlight) reflected by other irradiated surfaces (irradiated surface E) such as the ground or a wall surface. In this embodiment, by including the intensity of the reflected sunlight reflected by these other irradiated surfaces in the intensity of the sunlight received by the irradiated surface, the accuracy of determining the intensity of the sunlight received by the irradiated surface can be further improved.

[0027] Normally, the reflectance of the other irradiated surface is used to calculate the reflected light from the other irradiated surface. The reflectance of a material surface is generally calculated using the reflectance at a specific representative wavelength, but in order to calculate the energy intensity accurately, it is desirable to use the spectral reflectance. Here, how to calculate the spectral reflectance of the other irradiated surface will be explained with reference to Figure 7.

[0028] FIG. 7 is a diagram illustrating an example of how to determine the spectral reflectance of another irradiated surface, where FIG. 7(a) illustrates measuring sunlight directly incident on irradiated surface E, and FIG. 7(b) illustrates measuring reflected light from irradiated surface E. As shown in FIGS. 7(a) and 7(b), a measuring instrument 50 is used here. The measuring instrument 50 functions as a spectrophotometer. First, as shown in FIG. 7(a), the measuring instrument 50 is used to measure the spectral irradiance from above irradiated surface E, which is the other irradiated surface (the spectral irradiance of sunlight irradiating irradiated surface E). The measuring instrument 50 is also used to measure the spectral irradiance reflected by irradiated surface E (the spectral irradiance of light reflected from irradiated surface E). The spectral reflectance of irradiated surface E is calculated using Equation 1 using the calculated spectral irradiance from above irradiated surface E and the spectral irradiance reflected by irradiated surface E.

number

[0029] Next, calculation of the irradiation energy received by the irradiation surface using the spectral reflectance obtained by Equation 1 will be described. Fig. 8 is a diagram explaining the light that enters the irradiation surface out of the light reflected from another irradiation surface. Fig. 8 shows a case where irradiation surface G, which is the other irradiation surface, is the ground, and takes into consideration the proportion of the light reflected from irradiation surface G that enters irradiation surface F. The energy of light that enters irradiation surface F from irradiation surface G (ground) can be calculated by Equation 2.

number

[0030] The irradiation energy received by the irradiated surface G in Equation 2 can be obtained, for example, by the method shown in FIG. 7(a). As the reflectance in Equation 2, for example, the spectral reflectance calculated by Equation 1 can be used. The area ratio H in Equation 2 can be calculated as described below with reference to FIG. 8. In FIG. 8, the field of view area of ​​the ground (irradiated surface G) at the irradiated surface F is determined by the angle θ corresponding to the orientation of the irradiated surface F. The angle θ that determines the field of view area of ​​the irradiated surface G, which is the ground, is geometrically the angle between the irradiated surface F and the horizontal direction if the irradiated surface G is an ideal ground surface (horizontal). However, in practice, it is preferable to take into account changes in topography and the fact that the more distant the position on the irradiated surface G from the irradiated surface F, the less reflected light enters the irradiated surface F and the attenuated contribution, and therefore it is preferable to determine the angle θ according to the required accuracy.

[0031] The ratio of the field of view area of ​​the ground on the irradiated surface F to the hemispherical area of ​​the irradiated surface F is the area ratio H. In this embodiment, when calculating the irradiation energy received by the irradiated surface F, it is possible to calculate with higher accuracy by adding the energy of light incident on the irradiated surface F calculated by Equation 2 to the irradiation energy due to diffuse solar radiation and direct solar radiation that directly irradiates the irradiated surface F. Furthermore, although Equation 2 takes the ground into account as another irradiated surface, it is possible to calculate with higher accuracy the irradiation energy received by the irradiated surface F by taking into account not only the light reflected by the ground but also all other irradiated surfaces onto which reflected light may enter the irradiated surface F.

[0032] FIG. 9 is a diagram showing an example of information stored in the storage unit 13 shown in FIG. 1. The storage unit 13 stores the information shown in FIG. 9 in, for example, a database format. The storage unit 13 stores date and time information 14 as a first primary key. The date and time information 14 may include the year, month, day, and time. The storage unit 13 stores location information 15 as a second primary key. The location information 15 is information that identifies a location on the earth using, for example, east longitude and north latitude.

[0033] The storage unit 13 stores the solar solid angle for the first primary key and the second primary key as the stored value of the sunlight intensity information 16. The storage unit 13 stores the direct sunlight intensity for the first primary key and the second primary key as the stored value of the sunlight intensity information 16. The storage unit 13 stores the scattered sunlight intensity for the first primary key and the second primary key as the stored value of the sunlight intensity information 16. The storage unit 13 stores the albedo value for the first primary key and the second primary key as the stored value of the sunlight intensity information 16. The albedo value is the ratio of the reflected sunlight intensity to the irradiated sunlight intensity. For example, in a certain area where the location information 15 is located, when the date and time information 14 indicates summer, the ground is soil and the albedo value is low, and when the date and time information 14 indicates winter, the ground is snowy and the albedo value is high. The values ​​shown in FIG. 9 are, for example, collected measured values ​​and stored in the storage unit 13. The direction information 16a included in the sunlight ray intensity information 16 is included in the stored values ​​shown in Fig. 9. The information included in the sunlight ray intensity information 16 includes, for example, a value calculated by solving the radiative transfer equation, or a value calculated in the process of solving the radiative transfer equation.

[0034] Fig. 10 is a flowchart showing the operation of the photon flux density prediction system 10 shown in Fig. 1. In step S71, it is determined whether or not prediction information inquiry information has been received from the client machine 2 or the client machine 3 via the network 4. If prediction information inquiry information has been received (step S71: Yes), the process proceeds to step S72, and if prediction information inquiry information has not been received (step S71: No), the process returns to step S71.

[0035] In step S72, based on the sunlight inquiry information included in the forecast information inquiry information received in step S71, sunlight intensity information 16 is calculated, which is information including the intensity of sunlight received at the location indicated by location information 15 included in the sunlight inquiry information on the date and time indicated by date and time information 14 included in the sunlight inquiry information, and stored in the memory unit 13. The sunlight intensity information 16 is calculated using, for example, a radiative transfer equation.

[0036] In step S73, a mathematical calculation is performed on the sunlight intensity information 16 calculated in step S72 using the direction information 16 as an input value to calculate direction-specific sunlight intensity information, which is sunlight intensity information received by the irradiated surface indicated by the direction information 16a at the date and time indicated by the date and time information 14, at the location indicated by the location information 15 included in the sunlight inquiry information.

[0037] In step S74, the direction-specific sunlight intensity information calculated in step S73 is used to calculate photon flux density prediction information that predicts the photon flux density for the prediction information inquiry information received in step S71. The photon flux density prediction information includes, for example, information on at least one of the photon flux density for the prediction information inquiry information and the photosynthetic photon flux density.

[0038] In step S75, the photon flux density prediction information calculated in step S74 is transmitted via network 4 to client machine 2 or client machine 3, which is the sender of the current prediction information inquiry information. Client machine 2 or client machine 3 can obtain photon flux density prediction information and easily obtain more detailed photon flux densities simply by transmitting prediction information inquiry information including sunlight inquiry information including date / time information 14, location information 15, and direction information 16a to photon flux density prediction system 10. Client machine 2 or client machine 3 can provide various applications to end users using the photon flux density prediction information obtained from photon flux density prediction system 10. [Example]

[0039] 11 is a block diagram showing the configuration of a photon flux density prediction system according to Example 2 of the present invention. A photon flux density prediction system 100 of this example is, for example, a server machine made up of a computer.

[0040] 11 , the photon flux density prediction system 100 includes a first calculation unit 110 that calculates sunlight intensity information 116 (described in detail below) in advance using date and time information and location information; a storage unit 113 that stores the sunlight intensity information 116 calculated by the first calculation unit 110 and various other information; a communication unit 117 that performs communication; an extraction unit 111 that extracts, from the storage unit 113, sunlight intensity information 116 associated with sunlight inquiry information included in prediction information inquiry information input via the communication unit 117; a second calculation unit 112 that calculates direction-specific sunlight intensity information using the sunlight intensity information 116 extracted by the extraction unit 111; and a third calculation unit 121 that calculates photon flux density prediction information that predicts the photon flux density corresponding to the prediction information inquiry information input via the communication unit 117 using the calculation result of the second calculation unit 112. The communication unit 117 communicates with, for example, an external device. The prediction information inquiry information includes sunlight inquiry information. The sunlight inquiry information includes date and time information 114, location information 115, and direction information 116a. The direction-specific sunlight intensity information is information including the intensity of sunlight received by the irradiated surface indicated by the direction information 116a at the date and time indicated by the date and time information 114 and the location indicated by the location information 115. The communication unit 117 transmits the calculation results of the third calculation unit 121 to the sender of the prediction information inquiry information. Each component shown in FIG. 11 may be configured as hardware. Furthermore, each component shown in FIG. 11 can also be realized by the photon flux density prediction system 100 executing a program, and the memory unit 113 may store the program to be executed by the photon flux density prediction system 100. The memory unit 113 has a volatile storage device or a non-volatile storage device depending on the use of the data.

[0041] The communication unit 117 may be configured to transmit direction-specific sunlight intensity information indicating the direction-specific sunlight intensity calculated by the second calculation unit 112 to the sender of the forecast information inquiry information.

[0042] In the photon flux density prediction system 100, the communication unit 117 may be configured to transmit the sunlight intensity information extracted by the extraction unit 111 to the sender of the prediction information inquiry information. In this case, the client machine 2 or 3 may have a configuration corresponding to the function of the second calculation unit 112. That is, the present invention is a photon flux density prediction system in which a server machine and a client machine are connected via a network, and the server machine includes a first calculation unit 110 that calculates sunlight intensity information 116, which is information on sunlight intensity at a date and time indicated by date and time information 114, which is information on the date and time, and a location indicated by location information 115, which is information on the location; a storage unit 113 that stores the date and time information 114, the location information 115, and the sunlight intensity information 116 calculated by the first calculation unit 110 in association with each other; a first communication unit that communicates with the client machine; and a sunlight intensity information associated with the sunlight inquiry information included in the prediction information inquiry information input via the first communication unit. The server machine may be configured to include an extraction unit 111 that extracts sunlight intensity information 116 from a memory unit 113, a first communication unit that transmits the sunlight intensity information 116 extracted by the extraction unit 111 to the client machine, and the client machine may be configured to include a second communication unit that communicates with the server machine and a second calculation unit (a configuration corresponding to the function of the second calculation unit 112) that calculates direction-specific sunlight intensity information using the sunlight intensity information 116 input via the second communication unit, and the second communication unit transmits sunlight inquiry information to the server machine, the sunlight inquiry information including date and time information 114, location information 115, and direction information 116a that indicates the direction in which the irradiated surface that is irradiated by sunlight faces.

[0043] The memory unit 113 stores date and time information 114, which is information about date and time, location information 115, which is information about location, and sunlight intensity information 116, which is information about sunlight intensity at the date and time indicated by the date and time information 114 and the location indicated by the location information 115 (sunlight intensity information 116 calculated by the first calculation unit 110), in association with each other.

[0044] The communication unit 117 of the photon flux density prediction system 100 is connected to a network 4 such as the Internet. Client machines 2 and 3 used by users are connected to the network 4, and the client machines 2 and 3 communicate with the photon flux density prediction system 100 via the network 4.

[0045] Prediction information inquiry information is transmitted from the client machine 2 or the client machine 3 to the photon flux density prediction system 100 via the network 4. The sunlight inquiry information includes date and time information 114, location information 115, and direction information 116a. Direction-specific sunlight intensity information, which is the extraction result by the extraction unit 112, is transmitted from the photon flux density prediction system 100 via the network 4 to the client machine 2 or the client machine 3 that has sent the prediction information inquiry information.

[0046] In this embodiment, an example of information stored in the storage unit 113 is the information shown in FIG.

[0047] 12 is a flowchart showing the operation of the photon flux density prediction system 100 shown in FIG. 11. In step S91, information contained in the sunlight intensity information 116 for all locations at all dates and times is calculated using date and time information and location information. The information contained in the sunlight intensity information 116 is calculated using, for example, a radiative transfer equation. That is, the information contained in the sunlight intensity information 116 includes, for example, values ​​calculated by solving the radiative transfer equation and values ​​calculated in the process of solving the radiative transfer equation. In step S92, the information contained in the sunlight intensity information 116 calculated in step S91 is stored in the storage unit 113.

[0048] Next, if forecast information inquiry information is received from client machine 2 or client machine 3 via network 4 (step S93: Yes), proceed to step S94; if forecast information inquiry information is not received (step S93: No), return to step S93. The forecast information inquiry information includes sunlight inquiry information. Note that the calculation of the information included in sunlight intensity information 116 and the storage of the calculated information in memory unit 13 may be completed before the process of receiving sunlight inquiry information from client machine 2 or 3 is performed, or the calculated information may be updated whenever new data becomes available, for example, for an area that was previously unavailable.

[0049] In step S94, sunlight intensity information 116 is extracted from storage unit 113 based on sunlight inquiry information included in the data received from client machine 2 or client machine 3. That is, sunlight intensity information 116 corresponding to date and time information 114 and location information 115 included in the sunlight inquiry information is extracted from storage unit 113. In step S94, a mathematical calculation is further performed on the extracted sunlight intensity information 116 using direction information 116a as an input value to calculate direction-specific sunlight intensity information, which is information including the intensity of sunlight received by the irradiated surface indicated by direction information 116a at the date and time indicated by date and time information 114 and at the location indicated by location information 15 included in the sunlight inquiry information.

[0050] In step S95, the direction-specific sunlight intensity information calculated in step S94 is used to calculate photon flux density prediction information that predicts the photon flux density corresponding to the prediction information inquiry information received in step S93. The photon flux density prediction information includes, for example, information on at least one of the photon flux density corresponding to the prediction information inquiry information and the photosynthetic photon flux density. The calculation of the photon flux density prediction information will be described in detail later.

[0051] In step S96, the photon flux density prediction information calculated in step S95 is transmitted via network 4 to client machine 2 or client machine 3, which is the sender of the current prediction information inquiry information. Client machine 2 or client machine 3 can obtain photon flux density prediction information and easily obtain more detailed photon flux density prediction information simply by transmitting prediction information inquiry information including sunlight inquiry information including date and time information 14, location information 15, and direction information 16a to photon flux density prediction system 100. Client machine 2 or client machine 3 can provide various applications to end users using the photon flux density prediction information obtained from photon flux density prediction system 100.

[0052] Furthermore, in this embodiment, the sunlight intensity information 116 is calculated in advance, so that the response is faster and direction-specific sunlight intensity information can be provided with greater immediacy than when the sunlight inquiry information is received from the client machine 2 or the client machine 3 and then the information is calculated. [Example]

[0053] <Calculation of irradiated heat amount> In this example, the photon flux density prediction system 10 can calculate the amount of irradiated heat in the configuration shown in FIG. 1. FIG. 13 is a graph showing spectral irradiance. In FIG. 13, the horizontal axis represents the wavelength of light, and the vertical axis represents the spectral irradiance. FIG. 13 shows the results of calculating the amount of heat irradiated to the irradiated surface after calculating the energy intensity of the sky relative to the irradiated surface through a simulation based on the above-described Example 1. FIG. 13 also shows an example of calculation of the amount of heat received by an irradiated surface located on Miyako Island from 12:00 to 13:00 on June 20, 2016. According to this example, the integrated value of the energy irradiated to the irradiated surface is determined to be 1,029.8 [W / m2], and the amount of heat received by the irradiated surface is determined to be 3,707,358 [J / m2]. This embodiment makes it possible to predict the amount of heat radiation received by high-rise structures and forests, which have been difficult to measure, and can be used in the protection design of structures and forest plantation plans. [Example]

[0054] <Calculation of the amount of radiation irradiated onto a structure> In this embodiment, in the configuration shown in FIG. 1, the photon flux density prediction system 10 can predict the temperature rise according to the irradiated material (the material of the irradiated material having the irradiated surface) from the amount of heat received by the irradiated surface, which is the result obtained in Example 3. In this embodiment, the prediction of the temperature rise of this irradiated surface will be explained. FIG. 14 is a diagram showing the calculation conditions for the temperature rise of the irradiated material. As shown in FIG. 14, in this embodiment, it is assumed that sunlight is incident on the irradiated surface J2 of the irradiated material J1, and the reflectance of the irradiated surface J2 is 30%. It is also assumed that the thermal radiation from the irradiated material J1 is zero. It is also assumed that the heat transfer and heat conduction from the irradiated material J1 are zero.

[0055] FIG. 15 shows an example of a calculation of the temperature rise after one hour depending on the material. When sunlight with a heat quantity of 3,707,358 [J / m2] is incident on irradiation material J1, whose irradiation surface J2 has a reflectance of 30%, the reflected heat quantity at irradiation surface J2 is 1,112,207 [J / m2], and the heat quantity absorbed by irradiation material J1 is 2,595,151 [J / m2]. Considering the thermal characteristics of each material shown in FIG. 15, the temperature rise after one hour for each material can be calculated as shown in FIG. 15. For example, if irradiation material J1 is made of steel, its heat capacity is 18,864 [J / K], so the heat quantity received in one hour is 2,595.151 [J / m2] × 1 [m2] / 18,864 [J / K] = 137.6 [K]. Therefore, the temperature of the steel, initially at 20°C, rises to 20°C + 137.6°C = 157.6°C after one hour. Similarly, the temperature of the glass wool insulation plate, initially at 20°C, rises to 20°C + 25745.5°C = 25765.5°C. Note that the shape of the material model of the irradiated material J1 is 1m x 1m x 5mm thick. This calculation assumes that there is no heat radiation (thermal radiation, thermal conduction, heat transfer, etc.) from the irradiated material J1. In reality, the temperature would not rise this much because of heat radiation from the irradiated material J1. Heat radiation can be taken into account in the calculation as needed. This calculation makes it possible to predict temperature rises on irradiated surfaces, such as structures and land, that are difficult to measure because of their high elevations or because people cannot approach them. This example can be used to predict structural deterioration and for forest plantation planning. [Example]

[0056] <Inclusion of reflections from other illuminated surfaces (example of a sidewalk)> In this embodiment, in the configuration shown in FIG. 1, the photon flux density prediction system 10 calculates the reflected energy based on the reflectance of the sidewalk and can calculate the amount of irradiation irradiated onto the irradiation surface. The reflectance of the sidewalk can be determined, for example, by the method described with reference to FIGS. 7(a) and 7(b). FIG. 16(a) is a graph showing an example of the spectral irradiation intensity from each direction on the sidewalk. In FIG. 16(a), the horizontal axis represents the wavelength of light, and the vertical axis represents the irradiation intensity. FIG. 16(b) is a graph showing an example of the spectral reflectance on the sidewalk. In FIG. 16(b), the horizontal axis represents the wavelength of light, and the vertical axis represents the spectral reflectance. Referring to FIG. 16(b), it can be seen that the spectral reflectance of the sidewalk is 10 to 20%. [Example]

[0057] <Inclusion of reflections from other illuminated surfaces (e.g. grass)> In this embodiment, in the configuration shown in FIG. 1, the photon flux density prediction system 10 calculates the reflected energy based on the reflectance of the grassland and can calculate the amount of irradiation irradiated onto the irradiation surface. The reflectance of the grassland can be calculated, for example, using the method described with reference to FIGS. 7(a) and 7(b). FIG. 17(a) is a graph showing an example of the spectral irradiation intensity from each direction on a grassland. In FIG. 17(a), the horizontal axis represents the wavelength of light, and the vertical axis represents the irradiation intensity. FIG. 17(b) is a graph showing an example of the spectral reflectance of a grassland. In FIG. 17(b), the horizontal axis represents the wavelength of light, and the vertical axis represents the spectral reflectance. Referring to FIG. 17(b), it can be seen that the spectral reflectance of the grassland is 5 to 10%. Also referring to FIG. 17(b), it can be seen that the reflectance of the grassland increases sharply in the near-infrared region. [Example]

[0058] <Inclusion of reflections from other irradiated surfaces (example of asphalt)> In this embodiment, in the configuration shown in FIG. 1, the photon flux density prediction system 10 calculates the reflected energy based on the reflectance of the asphalt and can calculate the amount of irradiation irradiated onto the irradiation surface. The reflectance of asphalt can be determined, for example, by the method described with reference to FIGS. 7(a) and 7(b). FIG. 18(a) is a graph showing an example of the spectral irradiation intensity from each direction on asphalt. In FIG. 18(a), the horizontal axis represents the wavelength of light, and the vertical axis represents the irradiation intensity. FIG. 18(b) is a graph showing an example of the spectral reflectance of asphalt. In FIG. 18(b), the horizontal axis represents the wavelength of light, and the vertical axis represents the spectral reflectance. Referring to FIG. 18(b), it can be seen that the spectral reflectance of asphalt is 5%. [Example]

[0059] <Explanation of Photon Flux Density Prediction Information> The photon flux density prediction information will be further explained below. The above-mentioned third calculation units 21 and 121 calculate photon flux density prediction information. The photon flux density prediction information includes, for example, information on at least one of the photon flux density and the photosynthetic photon flux density in response to the prediction information inquiry information.

[0060] (Prediction of photon flux density) The third calculation units 21 and 121 calculate the photon flux density (PFD) based on the direction-specific solar ray intensity information calculated by the second calculation units 12 and 112. In this embodiment, the third calculation units 21 and 121 obtain the photon flux density from the energy intensity of the spectral solar radiation (spectral radiance) using Equation 3.

number

[0061] The irradiance is correlated with the photon flux density and can be calculated using Equation 4. The third calculation units 21 and 121 can calculate the spectral irradiance using Equation 4.

number

[0062] (Prediction of photosynthetic photon flux density) The third calculation units 21 and 121 calculate the photosynthetic photon flux density (PPFD) based on the direction-specific solar ray intensity information calculated by the second calculation units 12 and 112. In this embodiment, the third calculation units 21 and 121 obtain the photosynthetic photon flux density from the energy intensity in the wavelength range of 400 to 700 nm of the spectral solar radiation (spectral radiance) using Equation 5.

number

[0063] (Predicting photosynthetic rate) Fig. 19 is a diagram showing the relationship between CO2 respiration rate and light intensity, and Fig. 20 is a diagram showing the relationship between photosynthetic rate and light intensity for each plant.

[0064] The third calculation units 21 and 121 calculate the photosynthetic rate from Equation 6 using the photosynthetic photon flux density obtained from Equation 5 based on the relationship between the photosynthetic rate and light intensity according to the plant species in FIG. 20. A1 Ask for.

number

[0065] FIG. 21 is a diagram showing the light energy required for photosynthesis. FIG. 21 shows two types of photoreceptors, chlorophyll a and chlorophyll b. The third calculation units 21 and 121 use the relationship of the light energy required for photosynthesis according to the photoreceptors (chlorophyll, etc.) of the plant species shown in FIG. 21 and the spectral irradiance to calculate a more precise photosynthetic rate using Equation 7. B1 Ask for.

number

[0066] Figures 22 and 23 are diagrams showing the dependence of the photosynthetic rate on light intensity and temperature. As shown in Figures 22 and 23, when the light intensity is A (weak light), the light intensity is the limiting factor for the photosynthetic rate, and when the light intensity is B (strong light), the temperature is the limiting factor for the photosynthetic rate.

[0067] 22 and 23, Equation 6 is changed to Equation 8, and Equation 7 is changed to Equation 9. The third calculation units 21 and 121 calculate the photosynthetic rate by Equation 8. A2 The third calculation units 21 and 121 can obtain the photosynthetic rate by using Equation 9. B2 For example, the leaf surface temperature can be obtained by the third calculation units 21 and 121 according to the invention described in Patent Application No. 2021-041495 previously filed by the applicant.

number

number

[0068] In terrestrial plants, the amount of CO2 absorbed per photon in the photosynthetic reaction shown in equation 10 is called the quantum yield.

number

[0069] The theoretical maximum value of this quantum yield is given by Equation 11.

number

[0070] For example, the PPFD (daily cumulative value) in Colombia on February 1, 2020 was 47.8 mol / m -2 ], so the leaf area is 1.0 [m 2 ] This means that approximately 12 units of CO2 are fixed.

number

[0071] In Colombia, it would take approximately 3,800 m to fix 1 ton of CO2 per day. -2 ] = 1,000,000 / 262.9. In this way, according to this embodiment, the third calculation units 21 and 121 estimate the amount of fixed CO2 from the predicted integrated value of PPFD at the point of interest, thereby making it possible to quantitatively calculate carbon credits.

[0072] <Example of prediction results from the photon flux density prediction system> An example of the prediction results obtained by the photon flux density prediction system of the present invention will now be described. FIG. 24 is a diagram showing conditions in an example of prediction by the photon flux density prediction system of the present invention. In this example, the location is Colombia. The dates and times are from 0:00 on February 1, 2020 to 24:00 on February 3, 2020. As described above, according to this example, the third calculation units 21 and 121 calculate solar radiation energy, photon flux density, and photosynthetic photon flux density.

[0073] Fig. 25 is a graph showing solar radiation energy calculated by the photon flux density prediction system of the present invention under the conditions of Fig. 24. In Fig. 25, the horizontal axis represents time, and the vertical axis represents energy.

[0074] Fig. 26 is a graph showing the photon flux density calculated (calculated using Equation 3) by the photon flux density prediction system of the present invention under the conditions of Fig. 24. In Fig. 26, the horizontal axis represents time, and the vertical axis represents photon flux density.

[0075] Fig. 27 is a graph showing the photosynthetic photon flux density calculated (calculated using Equation 5) by the photon flux density prediction system of the present invention under the conditions of Fig. 24. In Fig. 27, the horizontal axis represents time, and the vertical axis represents photon flux density.

[0076] Fig. 28 is a graph showing the measured values ​​of photosynthetic photon flux density under the conditions of Fig. 24. In Fig. 28, the horizontal axis represents time and the vertical axis represents photon flux density.

[0077] Figure 29 is a graph showing a combination of solar radiation energy (taking into account clear skies and cloud cover), photon flux density (taking into account clear skies and cloud cover), and photosynthetic photon flux density (taking into account clear skies and cloud cover) calculated by the photon flux density prediction system of the present invention under the conditions of Figure 24, as well as the measured value of photosynthetic photon flux density. As can be seen from the graph in Figure 29, the photosynthetic photon flux density (taking into account cloud cover) calculated by the photon flux density prediction system agrees well with the measured value of photosynthetic photon flux density, demonstrating that precise prediction is possible.

[0078] Fig. 30 is a diagram showing the prediction of the amount of carbon fixation on a clear day by the photon flux density prediction system of the present invention under the conditions of Fig. 24. The photon flux density prediction system of the present invention can predict the amount of carbon fixation by applying Equation 11 to the predicted results of photosynthetic photon flux density.

[0079] As described above, according to the present invention, it is possible to predict the photon flux density and photosynthetic photon flux density in a specified area at a specified time.

[0080] Although preferred embodiments of the present invention have been described above, the present invention is not limited to these embodiments. The object of the present invention can also be achieved by providing a storage medium storing program code (computer program) that realizes the functions of the above-described embodiments to a system or device, and having a computer in the system or device read and execute the program code stored in the storage medium. In this case, the program code read from the storage medium itself realizes the functions of the above-described embodiments, and the storage medium storing the program code constitutes the present invention. Furthermore, in the above-described embodiments, a computer executes a program to function as each processing unit, but some or all of the processing may be implemented using dedicated electronic circuits (hardware). The present invention is not limited to the specific embodiments described, and various modifications and variations are possible within the spirit and scope of the present invention as defined by the claims. [Explanation of symbols]

[0081] 2, 3 client machines 4 Network 10 Photon flux density prediction system 11 First calculation section 12 Second calculation section 21 Third Calculation Department 13 Storage section 17 Communications Department

Claims

1. a storage unit that stores date and time information, which is information about a date and time; location information, which is information about a location; and sunlight intensity information, which is information about the date and time indicated by the date and time information and the sunlight intensity at the location indicated by the location information, in association with each other; a communication unit for performing communication; a first calculation unit that calculates the sunlight ray intensity information associated with sunlight ray inquiry information included in the forecast information inquiry information input via the communication unit; a second calculation unit that calculates direction-specific solar ray intensity information using the calculation result of the first calculation unit; a third calculation unit that uses the calculation result of the second calculation unit to calculate photon flux density prediction information that predicts a photon flux density in response to the prediction information inquiry information input via the communication unit using Equation 1; Equipped with the sunlight inquiry information includes the date and time information, the location information, and direction information indicating a direction in which an irradiation surface irradiated with sunlight faces, the direction-specific sunlight intensity information is information including the intensity of sunlight received by an irradiation surface indicated by the direction information at the date and time indicated by the date and time information and the location indicated by the location information, the communication unit transmits the calculation result of the third calculation unit to a sender of the prediction information inquiry information. A photon flux density prediction system characterized by: [Equation 1]

2. a first calculation unit that calculates sunlight intensity information that is information about sunlight intensity at a date and time indicated by date and time information that is information about date and time and at a location indicated by location information that is information about a location; a storage unit that stores the date and time information, the location information, and the sunlight intensity information calculated by the first calculation unit in association with each other; a communication unit for performing communication; an extracting unit that extracts, from the storage unit, the sunlight ray intensity information associated with sunlight ray inquiry information included in the forecast information inquiry information input via the communication unit; a second calculation unit that calculates direction-specific sunlight intensity information using the sunlight intensity information extracted by the extraction unit; a third calculation unit that calculates photon flux density prediction information that predicts a photon flux density in response to the prediction information inquiry information input via the communication unit using a calculation result of the second calculation unit, using Equation 2; Equipped with the sunlight inquiry information includes the date and time information, the location information, and direction information indicating a direction in which an irradiation surface irradiated with sunlight faces, the direction-specific sunlight intensity information is information including the intensity of sunlight received by an irradiation surface indicated by the direction information at the date and time indicated by the date and time information and the location indicated by the location information, the communication unit transmits the calculation result of the third calculation unit to a sender of the prediction information inquiry information. A photon flux density prediction system characterized by: [Equation 2]

3. 3. The photon flux density prediction system according to claim 1 or 2, The third calculation unit a fourth calculation unit that calculates a photon flux density in response to the prediction information inquiry information input via the communication unit; a fifth calculation unit that calculates a photosynthetic photon flux density for the prediction information inquiry information input via the communication unit using the calculation result of the fourth calculation unit; a sixth calculation unit that calculates a fixed carbon amount for the prediction information inquiry information input via the communication unit using a calculation result of the fourth calculation unit; Including, the communication unit transmits a calculation result of at least one of the fourth calculation unit, the fifth calculation unit, and the sixth calculation unit to a sender of the prediction information inquiry information. A photon flux density prediction system characterized by:

4. 4. The photon flux density prediction system according to claim 1, The direction indicated by the direction information is a direction perpendicular to the surface on which the irradiation surface receiving the sunlight spreads. A photon flux density prediction system characterized by:

5. 5. The photon flux density prediction system according to claim 1, the sunlight intensity information includes the intensity of sunlight reflected by a surface other than the irradiation surface facing the direction indicated by the direction information; A photon flux density prediction system characterized by:

6. Computer a storage unit that stores date and time information, which is information about a date and time; location information, which is information about a location; and sunlight intensity information, which is information about the date and time indicated by the date and time information and the sunlight intensity at the location indicated by the location information, in association with each other; a communication unit for performing communication; a first calculation unit that calculates the sunlight ray intensity information associated with sunlight ray inquiry information included in the forecast information inquiry information input via the communication unit; a second calculation unit that calculates direction-specific solar ray intensity information using the calculation result of the first calculation unit; a third calculation unit that uses the calculation result of the second calculation unit to calculate photon flux density prediction information that predicts a photon flux density in response to the prediction information inquiry information input via the communication unit using Equation 3; It functions as the sunlight inquiry information includes the date and time information, the location information, and direction information indicating a direction in which an irradiation surface irradiated with sunlight faces, the direction-specific sunlight intensity information is information including the intensity of sunlight received by an irradiation surface indicated by the direction information at the date and time indicated by the date and time information and the location indicated by the location information, the communication unit transmits the calculation result of the third calculation unit to a sender of the prediction information inquiry information. A photon flux density prediction program characterized by: [Equation 3]

7. Computer a first calculation unit that calculates sunlight intensity information that is information about sunlight intensity at a date and time indicated by date and time information that is information about date and time and at a location indicated by location information that is information about a location; a storage unit that stores the date and time information, the location information, and the sunlight intensity information calculated by the first calculation unit in association with each other; a communication unit for performing communication; an extracting unit that extracts, from the storage unit, the sunlight ray intensity information associated with sunlight ray inquiry information included in the forecast information inquiry information input via the communication unit; a second calculation unit that calculates direction-specific sunlight intensity information using the sunlight intensity information extracted by the extraction unit; a third calculation unit that calculates photon flux density prediction information that predicts a photon flux density in response to the prediction information inquiry information input via the communication unit using a calculation result of the second calculation unit, using Equation 4; It functions as the sunlight inquiry information includes the date and time information, the location information, and direction information indicating a direction in which an irradiation surface irradiated with sunlight faces, the direction-specific sunlight intensity information is information including the intensity of sunlight received by an irradiation surface indicated by the direction information at the date and time indicated by the date and time information and the location indicated by the location information, the communication unit transmits the calculation result of the third calculation unit to a sender of the prediction information inquiry information. A photon flux density prediction program characterized by: [Equation 4] 8. The photon flux density prediction system according to claim 1, comprising: The third calculation unit uses the calculation result of the second calculation unit to calculate the photosynthetic rate for the prediction information inquiry information input via the communication unit using Equation 5. A photon flux density prediction system characterized by: [Equation 5]

9. A photon flux density prediction program according to claim 6 or 7, The third calculation unit uses the calculation result of the second calculation unit to calculate the photosynthetic rate for the prediction information inquiry information input via the communication unit using Equation 6. A photon flux density prediction program characterized by: [Equation 6]

Images