Page updated: November 21, 2020
Author: Curtis Mobley
View PDF

Raman Scattering

First of all, the pronunciation is RAman, not raMAN. Raman scattering is a type of inelastic scattering of light by molecules. Compton scattering (see The Nature of Light) is inelastic scattering of X-rays by electrons. Upon learning of Compton’s discovery (for which Compton received the 1928 Nobel Prize in Physics), C. V. Raman thought there should be a similar effect at visible wavelengths and set out to find it. He succeeded and subsequently received the 1930 Nobel Prize for the discovery and explanation of what is now called Raman scattering or the Raman effect. (Incidentally, Raman received his bachelor’s degree in physics at 16, graduating at the head of his class, and he published his first paper at 18. The Raman family produced a number of scientists, included Raman’s nephew Subrahmanyan Chandrasekhar, who received the Nobel in 1983. The family seems to have had good genes for science.)

The physics of Raman scattering is quite complicated, but it can be conceptualized as following. Incident light can excite a molecule in its ground state to a higher “virtual” energy level, which then immediately decays back to a lower level accompanied by the emission of light. If the decay returns the molecule to its initial state, the scattering is elastic and is called Rayleigh scattering. If the decay is to a molecular vibrational level above the ground state, then then emitted light has a longer wavelength (lower energy) than the incident light: this is Raman scattering. (See The Physics of Absorption for a discussion of electronic, vibrational, and rotational energy levels.) The use of Raman scattering, often called Raman spectroscopy, is widely used in chemistry as a way to study the vibrational energy levels of molecules. In those applications, the incident light usually comes from a laser. In the ocean, the molecule of interest is water, and the exciting light can be either sunlight or a laser.

The wavelength shift for Raman scatter by water is exceptionally large, corresponding to a wavenumber (1/wavelength) shift of about 3400cm1, which at visible wavelengths is many tens to more than a hundred nanometers. The time scale for Raman scattering (i.e, the time for the molecular excitation and subsequent decay) is roughly 1013 to 1012 seconds. This is much faster than the time scale of fluorescence by chlorophyll or CDOM molecules, which is on the order of 109s or longer. Thus Raman scatter can be thought of as an almost instantaneous scattering process, rather than as the absorption of light followed much later by the emission of new light.

Quantities Defining Raman Scattering

The quantities needed to compute Raman scatter contributions to the radiance are

  • the Raman scattering coefficient bR(λ), with units of m1
  • the Raman wavelength distribution function fR(λ,λ), with units of nm1
  • the Raman scattering phase function beta β̃R(ψ), with units of sr1

The next sections discuss each of these quantities in turn.

The Raman scattering coefficient

The Raman scattering coefficient bR(λ) tells how much of the irradiance at the excitation wavelength λ scatters into all emission wavelengths λ > λ, per unit of distance traveled by the excitation irradiance. The most recently published values of bR(488nm) for water are (2.7 ± 0.2) × 104m1 (Bartlett et al. (1998)) and 2.4 × 104m1 (Desiderio (2000)). (The current version (6.0) of the HydroLight radiative transfer model uses bR(488nm) = 2.6 × 104m1 as the default value.)

Various values for the wavelength dependence of bR can be found in the literature. Bartlett et al. (1998) reviewed the wavelength dependence of the Raman scattering coefficient in detail and found, based on their measurements, a wavelength dependence of λ4.8±0.3 for calculations performed in energy units (as in HydroLight). In terms of the excitation wavelength, Bartlett et al. found bR(λ) = bR(488)(488λ)5.5±0.4 for energy computations. (The current version of HydroLight uses bR(λ) = bR(488)(488λ)5.5 as the default.) For calculations in terms of photon numbers (as in a Monte Carlo simulation), Bartlett et al. found wavelength dependencies of (λ)5.3±0.3 or (λ)4.6±0.3.

The Raman wavelength distribution function

The Raman wavelength distribution function fR(λ,λ) relates the excitation and emission wavelengths, i.e., what wavelengths λ receive the scattered spectral irradiance for a given excitation wavelength λ or, conversely, what wavelengths λ excite a given emission wavelength λ. (The term “wavelength distribution function” is non-standard, but descriptive, so that is what was used in Light and Water and again here.) The function fR(λ,λ) is most conveniently described in terms of the corresponding wavenumber distribution function fR(κ), where κ is the wavenumber shift, expressed in units of cm1. This follows because the Raman-scattered light undergoes a frequency shift that is determined by the type of molecule and is independent of the incident frequency. The wavenumber κ in cm1 is related to the wavelength λ in nm by κ = 107λ, and to the frequency ν by κ = νc, where c is the speed of light. (The 107 factor converts nanometers to centimeters.)

According to Walrafen (1967), the shape of fR(κ) for water is given by a sum of four Gaussian functions:

fR(κ) = π 4 ln 2 1 2 i=14A i 1 j=14A j 1 Δκj exp 4 ln 2(κ κj)2 Δκj2 (cm), (1)


  • κ is the wavenumber shift of the Raman-scattered light, relative to the wavenumber κ of the incident light, in cm1
  • κj is the center of the jth Gaussian function, in cm1
  • Δκj is the full width at half maximum of the jth Gaussian function, in cm1
  • Aj is the nondimensional weight of the jth Gaussian function.

The values of Aj, κj, and Δκj for pure water at a temperature of 25 deg C are given in Table 1. Figure 1 shows fR(κ) evaluated for the water parameter values of Table 1. The function shows a peak and a shoulder, which result from the sums of the four gaussians seen in Eq. (1). For water, the wavenumber shift is roughly 3400cm1.

Δκ j

1 0.41 3250 210
2 0.39 3425 175
3 0.10 3530 140
4 0.10 3625 140

Table 1: Parameter values for the Raman wavenumber redistribution function fR(κ) of Eq. (1), for pure water at a temperature of 25 C. Data from Walrafen (1967).


Figure 1: The fR(κ) function of Eq. (1) for pure water at 25 C.

Consider, for example, incident light with λ = 500nm, which corresponds to κ = 20000cm1. Figure 1 shows that this light, if Raman scattered, will be shifted by roughly 3400cm1 to κ = 16600cm1, which corresponds to λ 602nm.

The function fR(κ) can be interpreted as a probability density function giving the probability that a light of any incident wavenumber κ = 107λ, if Raman scattered, will be scattered to a wavenumber

κ = κ κ. (2)

The function fR(κ) satisfies the normalization condition

0κ fR(κ)dκ = 1, (3)

as is required of any probability distribution function. The integration limits above come from observing that as λ then the wavenumber κ κ, and as λ λ then κ 0. If A change of variables from κ to λ in Eq. (3) leads to the corresponding wavelength redistribution function fR(λ λ). Thus

0κ fR(κ)dκ = λf R 107 λ dκ dλdλ =λf R 107 1 λ 1 λ107 λ2 dλ λf R(λ λ)dλ = 1,

where the wavelengths are in nanometers. In the last equation, we have identified the function

fR(λ,λ) 107 λ2 fR 107 λ = 107 λ2 fR 107 1 λ 1 λ ifλ < λ 0ifλ λ (4)

as being the desired Raman wavelength redistribution function, with units of nm1.

In fR(λ,λ) we can fix the incident wavelength λ and plot the corresponding emission wavelengths. We can also fix the emission wavelength λ and use the function to see where the light emitted at λ comes from. Both options are seen in Fig. 2.


Figure 2: Raman emission function for excitation at 488 nm (red), and excitation function for emission at 583 nm (blue). Compare with similar figures in Ge et al. (1993) (their Fig 3) and Kattawar and Xu (1992) (their Fig. 2).

Figure 3 shows fR(λ,λ) for four values of the incident wavelength λ. These plots show that as the incident wavelength λ increases, the emission band becomes broader and the shift from λ to λ becomes larger. Excitation at 400 nm gives emission centered at roughly 463 nm, a shift of 63 nm, but excitation at 550 nm gives emission centered at round 677 nm, a shift of 127 nm.


Figure 3: The Raman emission functions for four excitation wavelengths, as computed from Eq. (4).

Equation (2) can be rewritten as

λ = 107 107 λ 3400

and used to compute the approximate center of the emission wavelength band for a given excitation wavelength. Figure 4 shows the result.


Figure 4: Approximate center λ of the Raman emission band for a given excitation wavelength λ.

The Raman phase function

The Raman phase function β̃R(ψ) gives the angular distribution of the Raman scattered radiance. This function (averaging over all polarization states) is given by

β̃R(ψ) = 3 16π1 + 3ρ 1 + 2ρ 1 + 1 ρ 1 + 3ρ cos 2ψ,

where ψ is the scattering angle between the direction of the incident and scattered radiance, and ρ is the depolarization factor. The value of ρ depends on the wavenumber shift κ (Ge et al. (1993), their Fig. 2). For a value of κ = 3400cm1, ρ 0.18, in which case the phase function is

β̃R(ψ) = 0.068(1 + 0.53 cos 2ψ).

This phase function is similar in shape to the phase function for elastic scattering by pure water.

Incorporation of Raman Scatter into Radiative Transfer Calculations

We now have have the pieces needed to define the volume scattering function for Raman scattering, βR(ξ ξ; λ λ), where ξ and ξ represent the incident and final directions of the light. This VSF specifies the strength of the Raman scattering via the Raman scattering coefficient bR(λ), the wavelength distribution of the scattered light via the wavelength redistribution function fR(λ,λ), and its angular distribution relative to the direction of the incident light via the Raman phase function β̃R(ψ). Thus we have

βR(ξ ξ; λ λ) = b R(λ)f R(λ,λ)β̃ R(ψ)(m1nm1sr1). (5)

Raman scatter is incorporated into unpolarized radiative transfer calculations as a source term (See Eq. 3 of The Scalar Radiative Transfer Equation):

cos 𝜃dL(z,𝜃,ϕ,λ) dz = c(z,λ)L(z,𝜃,ϕ,λ) (6) + 02π0πβ(z; 𝜃,ϕ 𝜃,ϕ; λ)L(z,𝜃,ϕ,λ) sin 𝜃d𝜃dϕ + 0λ02π0πβ R(𝜃,ϕ 𝜃,ϕ; λ λ)L(z,𝜃,ϕ,λ) sin 𝜃d𝜃dϕdλ.

The scattering angle ψ is given by ψ = cos 1(ξ ξ), which can be computed from 𝜃,ϕ,𝜃,ϕ via Eq. (5) of the Geometry page.

Note that this formidiable equation cannot be solved at just the wavelength λ of interest; it must be solved at all wavelengths λ < λ that contribute Raman-scattered radiance to the radiance at λ.

Examples of Raman Effects

Effect on remote-sensing reflectance

Figure 5 shows the effect of Raman scattered radiance on the remote-sensing reflectance for chlorophyll values of Chl = 0.02 and 2mgm3 as simulated by HydroLight. The HydroLight runs used a bio-optical model for Case 1 water for which the absorption and scattering properties of the water are determined only by the chlorophyll concentration. The water was homogeneous and infinitely deep. The Sun was at a 30 deg zenith angle in a clear sky; the wind speed was 5ms1. Identical runs were made with and without Raman scatter included in the runs. It is seen that for the very clear water with Chl = 0.02, Rrs is as much as 22% higher, but for the water with Chl = 2, Raman increases Rrs by at most 5%. The Raman effect decreases but still can be significant in higher chlorophyll waters or in turbid Case 2 waters. Although difficult to see in Fig. 5, the Raman effect does not “turn on” until about 340 nm even though the HydroLight run started at 300 nm. This is because the start of the emission band for excitation near 300 nm is around 340 nm (recall Fig. 4).


Figure 5: Effect of Raman scattering on remote-sensing reflectance Rrs. The red curves included Raman scattering and the blue curves did not.

Effect on upwelling plane irradiance

Although Raman scattering does not have a large effect on Rrs except in the clearest water, it can be the dominant source of light at red wavelengths at depths where absorption by water has removed most of the incoming sunlight. This is illustrated in Fig. 6. HydroLight was run for Case 1 water with a chlorophyll concentration of Chl = 0.5mgm3, which is typical of open ocean water. The Sun was at a 30 deg zenith angle in a clear sky. The run started at 300 nm, so that Raman effects would be present at wavelengths greater than 340 nm. At 400 and 500 nm, the contribution by Raman scattered light to the upwelling irradiance is almost unnoticeable. This is because at these wavelengths elastically scattered solar radiance is the main contributor to the upwelling irradiance. At 580 nm, water absorption (aw(580) = 0.09m1) is beginning to filter out enough of the solar radiance that the Raman contribution, which comes from wavelengths around 485 nm (see Fig.4), where the light penetrates well to depth, is becoming the main contribution to Eu(z, 580). At 600 nm, water absorption (aw(600) = 0.22m1) has removed almost all of the solar light below about 20 m. Below 20 m, almost all of Eu(z, 600) comes from Raman scattered light that originates from wavelengths around 500 nm, where sunlight penetrates well to depth. It should be noted that below about 40 m, the depth rate of decay of Eu(z, 600) (i.e., Ku(z, 600)) is almost the same as the rate of decay the irradiance at 500 nm. This is because the light at 500 nm is the source of the light at 600 nm.

It was the measurement of unexpected upwelling irradiance Eu(z,λ) at depths below 50 m and wavelengths greater than 520 nm that led Sugihara et al. (1984) to suggest that the unexpected upwelling irradiance came from downwelling irradiance a blue-green wavelengths being Raman scattered. A number of subsequent studies (e.g., Stavin and Weidemann (1988), Marshall and Smith (1990)) confirmed this hypothesis. Many studies since have studied Raman effects on ocean light fields. For example, Kattawar and Xu (1992) and Ge et al. (1993) studied the filling in of Fraunhofer lines in underwater light by Raman scattering.


Figure 6: HydroLight simulations of Eu(z,λ) with and without Raman scattering.

Interpretation of Raman emission profiles

Another interesting example of the contribution of Raman scattering to in-water radiances is seen in Fig. 7. This plot shows the downwelling (zenith-viewing) radiance generated from a HydroLight run using its option for simulating LIDAR illumination at one wavelength in an otherwise black sky. The inputs were as follows:

  • The incident (LIDAR) irradiance was Ed(direct beam) = 1Wm2 at 488 nm; the sky was otherwise black.
  • The chlorophyll concentration was Chl = 0.05mgm3 for Case 1 water
  • The water was infinitely deep water with Eq. (6) solved down to 50 m
  • The output was at 560 to 610 nm at 1 nm resolution


Figure 7: Downwelling Raman-scattered radiance Ld(z,λ) as a function of depth and wavelength over the Raman emission band for 488 nm excitation. The wavelength shape of the Ld(z,λ) emission band depends on depth because of wavelength-dependent absorption.

The curves in Fig. 7 are explained as follows. At depth 0, just below the sea surface, the only contribution to Ld(0,λ) is upwelling radiance that is reflected back downward by the sea surface. By 5 m depth, there is now enough water above the simulated measurement instrument that the water column is generating significant downwelling radiance. Note that the shape of the emission has the same shape as the emission function seen in Fig. 2, namely a peak near 585 nm with a shoulder at about 580 nm. As the depth increases to 10 and then 20 m, the magnitude of Ld(z,λ) increases, but the shape of the emission begins to flatten between 580 and 585. By 30 m, the magnitude of Ld(z,λ) has decreased because of the decrease in the radiance penetrating to this depth from the water between the surface and 30 m, and the magnitude continues to decrease as the depth becomes greater. However, the shape of the emission at 30 m shows almost the same magnitudes at 580 and 585 nm, and at 40 and 50 m, the peak emission is actually greater at 580 than at 585 nm. This may seem strange because the shape of the Raman emission function seen in Fig. 2 is the same for all depths.

This “reversal” of the “peak-shoulder” shape of the emission is a consequence of the difference in absorption across the emission wavelengths. The total absorption (water plus phytoplankton) increases by a factor of three (from 0.071 to 0.221m1) between 570 and 600 nm, and by 22% (from 0.091 to 0.111m1) between 580 and 585 nm. These rapidly increasing absorption values change the shape of the local (at each depth) emission function when integrated over depth to obtain the total Ld(z,λ), which has contributions from all depths. Simply stated, the higher absorption at the 585 nm peak lets relatively less of the radiance emitted above a given depth reach the measurement depth than for the shoulder at 580 nm, so the 585 nm peak appears smaller relative to the shoulder than what is seen for the emission function of Fig. 2.

The claim that the change is shape with depth of the Raman emission is due to wavelength-dependent absorption can be verified as follows. An “artificial water” IOP data file was created with the IOP values between 570 and 600 nm having the values at 570 nm. The water IOPs are then the same over the entire range of emission wavelengths seen in Fig. 7. The resulting Raman Ld(z,λ) spectra are seen in Fig. 8. Now, without the wavelength-dependent absorption, the shape of Ld(z,λ) does not change with depth and, indeed, looks exactly like the shape of the emission function seen in Fig. 2. The magnitude of the Ld(z,λ) curves is greater than before because the absorption is less. Of course, if the HydroLight run had been made at 5 or 10 nm resolution, then the shape of the emission band would not have been resolved. However, the total Raman-scattered power would have been the same, but spread over the wider bands.


Figure 8: Downwelling Raman-scattered radiance as a function of depth for the case of artificial water IOPs having the same values from 570 to 600 nm. The shape of Ld(z,λ) is now independent of depth and is the same as the Raman-emission function of Fig. 2.

These simulations show that the interpretation of Raman-scattered spectra can be complicated because of IOP effects at both the excitation and emission wavelengths, even in the simplest possible case of excitation at one wavelength in an otherwise black sky. The situation becomes even more complicated for solar-stimulated Raman scatter because many excitation wavelengths can contribute to a range of emission wavelengths, and everything blurs together in a non-obvious fashion.

Temperature and Salinity Dependence

The Walrafen data presented in Table 1 were determined on pure water at a temperature of 25 deg C. There is, however, a small but significant dependence on temperature and salinity of the shape of the emission curve seen in Fig. 1. This dependence is shown in Fig. 9. The upper panel shows the shape of the emission curve as a function of temperature for a salinity of 15 PSU, and the lower curve shows the dependence on salinity for a temperature of 25 deg C. Artlett and Pask (2017) have shown in laboratory measurements that these differences can be used to simultaneously determine temperature and salinity with an RMS errors of ± 0.7degC and ± 1.4PSU, and they present the design for a three-channel Raman spectrometer that excites at 532 nm and measures the emission at three bands. The emission bands would be used to form band ratios, from which temperature and salinity can be extracted.


Figure 9: Dependence of the Raman wavelength redistribution function on temperature and salinity. (Figure 2 from Artlett and Pask (2017) reproduced by permission of the Optical Society of America)

Comments for Raman Scattering:

Loading Conversation