Skip to main content
Physics LibreTexts

2.10: A, p and q for General Reflectance Rules

  • Page ID
  • \( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

    \( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)

    \( \newcommand{\id}{\mathrm{id}}\) \( \newcommand{\Span}{\mathrm{span}}\)

    ( \newcommand{\kernel}{\mathrm{null}\,}\) \( \newcommand{\range}{\mathrm{range}\,}\)

    \( \newcommand{\RealPart}{\mathrm{Re}}\) \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)

    \( \newcommand{\Argument}{\mathrm{Arg}}\) \( \newcommand{\norm}[1]{\| #1 \|}\)

    \( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)

    \( \newcommand{\Span}{\mathrm{span}}\)

    \( \newcommand{\id}{\mathrm{id}}\)

    \( \newcommand{\Span}{\mathrm{span}}\)

    \( \newcommand{\kernel}{\mathrm{null}\,}\)

    \( \newcommand{\range}{\mathrm{range}\,}\)

    \( \newcommand{\RealPart}{\mathrm{Re}}\)

    \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)

    \( \newcommand{\Argument}{\mathrm{Arg}}\)

    \( \newcommand{\norm}[1]{\| #1 \|}\)

    \( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)

    \( \newcommand{\Span}{\mathrm{span}}\) \( \newcommand{\AA}{\unicode[.8,0]{x212B}}\)

    \( \newcommand{\vectorA}[1]{\vec{#1}}      % arrow\)

    \( \newcommand{\vectorAt}[1]{\vec{\text{#1}}}      % arrow\)

    \( \newcommand{\vectorB}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

    \( \newcommand{\vectorC}[1]{\textbf{#1}} \)

    \( \newcommand{\vectorD}[1]{\overrightarrow{#1}} \)

    \( \newcommand{\vectorDt}[1]{\overrightarrow{\text{#1}}} \)

    \( \newcommand{\vectE}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash{\mathbf {#1}}}} \)

    \( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

    \( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)

    Again, consider a sphere of radius a centred in an Oxyz frame with corresponding directional spherical coordinates (Θ, Φ), and let the sphere be irradiated with flux density F from the z-direction. For the geometrical albedo the phase angle α is zero and the incident and reflected radiation are given by μ0 = μ = cos Θ, so that

    \[p=\int_{0}^{2 \pi} \int_{0}^{\pi / 2} f_{r}(\cos \theta, \cos \Theta, 0 ; \ldots) \cos ^{2} \Theta \sin \theta d \Theta d \Phi,\]

    resulting in

    \[p=2 \pi \int_{0}^{1} f_{r}(\mu, \mu, 0 ; \ldots) \mu^{2} d \mu.\]

    Using the same geometry for the Bond albedo, for each point on the irradiated hemisphere we have μ0 = cos Θ, so that the directional hemispherical reflectance is

    \[\rho\left(\mu_{0}\right)=\int_{0}^{2 \pi} \int_{0}^{1} f_{r}\left(\mu_{0}, \mu, \alpha ; \ldots\right) \mu d \mu d \phi_{r},\]

    where the phase angle is that between the incident and reflected radiation at each stage of the integral,

    \[\cos \alpha=\mu_{0} \mu+\sqrt{\left(1-\mu_{0}^{2}\right)\left(1-\mu^{2}\right)} \cos \phi_{r},\]

    where φr is the azimuth of the reflected radiation. The Bond albedo is then given by

    \[A=\frac{1}{\pi} \int_{0}^{2 \pi} \int_{0}^{\pi / 2} \rho(\cos \Theta) \cos \Theta \sin \Theta d \Theta d \Phi,\]

    which reduces to

    \[A=2 \int_{0}^{t} \rho(\mu) \mu d \mu.\]

    The phase integral, equation (20), may be computed from equations (14), (15) and (16), the factor a2 F disappearing in the process, so that we may write, for the purposes of computation, equation (16) as

    \[I(\alpha)=\int_{\alpha-\pi / 2}^{\pi / 2} \int_{0}^{\pi} f_{r}\left(\mu_{0}, \mu, \alpha ; \ldots\right) \mu_{0} \mu \sin \Theta d \Theta d \Phi,\]

    where it can be seen that for Φ the range of integration is from α - π/2, the limb, to π/2, the terminator.

    In these equations it can be seen that the geometrical albedo is just a single integral and thus may be quickly and accurately integrated numerically with just about any method. The Bond albedo and the phase integral are, however, triple integrals, so that a method which combines the advantages of speed and accuracy is required; for this reason Gaussian Quadrature is the chosen method. In the following section we present this method in algorithmic form and discuss its application to the integrals at hand.

    For the theory and examples of Gaussian Quadrature, its performance compared to other methods of integration as well as tabulations of the roots and coefficients needed, the reader is referred to Celestial Mechanics, Chap. 1.

    This page titled 2.10: A, p and q for General Reflectance Rules is shared under a CC BY-NC 4.0 license and was authored, remixed, and/or curated by Max Fairbairn & Jeremy Tatum via source content that was edited to the style and standards of the LibreTexts platform; a detailed edit history is available upon request.