Fact-checked by Grok 2 weeks ago

Subsurface scattering

Subsurface scattering is a fundamental light transport phenomenon in computer graphics where incident light penetrates the surface of a translucent or semi-translucent material, interacts with its internal microstructure through multiple scattering events, and exits at a different location, producing a soft, diffused glow that enhances visual realism.^[1]^[2]^[3] This process contrasts with surface reflection models, as it accounts for light diffusion within the material rather than immediate bounce-off, making it essential for accurately simulating the appearance of organic and inorganic substances like human skin, marble, wax, milk, and plant leaves.^[2]^[4] The effect is particularly prominent in materials with high translucency and low absorption, where light can travel significant distances subsurface before re-emerging, softening shadows and reducing the visibility of fine surface details to mimic natural subsurface diffusion.^[3]^[4] Without subsurface scattering, rendered objects appear unnaturally opaque or waxy, failing to capture the characteristic translucency observed in real-world scenes, such as the reddish tint of earlobes under backlighting or the veined patterns in marble.^[1]^[2] This phenomenon is computationally intensive, as full simulation via path tracing or Monte Carlo methods can increase rendering times by factors of two or more, necessitating approximations for real-time applications in video games and interactive visualizations.^[4] In rendering pipelines, subsurface scattering is typically modeled using the Bidirectional Scattering-Surface Reflectance Distribution Function (BSSRDF), which extends traditional surface reflectance models like the Bidirectional Reflectance Distribution Function (BRDF) to handle light transport across material boundaries.^[2] Pioneering work in the late 1990s and early 2000s, including dipole diffusion approximations and photon mapping techniques, enabled efficient computation of these effects, earning recognition such as a 2004 Academy Award for advancements in rendering translucent materials.^[2] Modern implementations in engines like Unreal Engine and Blender incorporate real-time approximations, such as texture-space diffusion or precomputed radiance transfer, to balance fidelity and performance in film production, animation, and virtual reality.^[1]^[3]

Fundamentals

Definition

Subsurface scattering (SSS) is a light transport phenomenon in which light penetrates the surface of a translucent or semi-translucent material, undergoes multiple internal scattering events, and re-emerges at a different location on the surface, distinct from direct surface reflection modeled by bidirectional reflectance distribution functions (BRDFs).^[5] This process contrasts with surface reflection, where light bounces off the exterior without entering the material, allowing SSS to capture the diffusion of light within the volume.^[3] Common materials exhibiting subsurface scattering include human skin, marble, wax, milk, and plant leaves, where the internal structure causes light to interact extensively before exiting.^[5] In these materials, light absorption and scattering depend on wavelength and material properties, leading to realistic appearances that pure surface models cannot replicate.^[5] Visually, subsurface scattering produces effects such as softened edges around illuminated areas, color bleeding where one region's hue diffuses into adjacent shadowed zones, and reduced sharpness in shadows due to light diffusion across the surface.^[5] These characteristics contribute to the lifelike quality of rendered scenes, enhancing photorealism in translucent objects.^[3] In computer graphics, subsurface scattering gained recognition in the early 2000s as essential for photorealistic rendering, particularly through practical models that enabled efficient simulation of these effects in translucent materials.^[5]

Physical Principles

Subsurface scattering occurs in participating media, which are translucent materials such as skin, marble, or plant tissues that allow light to penetrate, interact internally through absorption, scattering, and transmission, and exit at different points on the surface.^[6] In these media, incoming light is partially transmitted across the surface interface, then undergoes multiple internal interactions where photons are absorbed (converting energy to heat or other forms), scattered (redirected by particles), or transmitted further without redirection, leading to a diffused appearance distinct from surface reflection.^[6] This internal propagation is governed by the radiative transfer equation, which balances these interactions to describe light transport.^[5] The primary scattering mechanisms in participating media depend on particle size relative to the light wavelength. For small particles (much smaller than the wavelength, e.g., molecules in clear materials), Rayleigh scattering dominates, producing isotropic diffusion with intensity inversely proportional to the fourth power of the wavelength, resulting in blue-tinted scattering, as observed in the Earth's atmosphere (the blue sky).^[7] For larger particles (comparable to or greater than the wavelength, e.g., cells in biological tissues), Mie scattering prevails, often forward-directed and less wavelength-dependent, enhancing light diffusion over longer paths in denser media like skin or leaves.^[8] These mechanisms collectively cause the soft, blurred glow characteristic of subsurface scattering by randomizing light directions multiple times.^[8] Material properties are quantified by the absorption coefficient \sigma_a, which measures the probability of photon absorption per unit distance, and the scattering coefficient \sigma_s, which indicates the likelihood of scattering events.^[5] The effective penetration depth of light is inversely related to these coefficients: high \sigma_a limits depth by rapid energy loss, while high \sigma_s promotes diffusion but can also trap light if absorption is low; the reduced scattering coefficient \sigma_s' = \sigma_s (1 - g), where g is the average cosine of the scattering angle, further refines this by accounting for forward bias in anisotropic media.^[5] At material interfaces, boundary conditions govern light entry and exit, with Fresnel reflection determining the fraction of light reflected based on the incident angle and refractive index difference (e.g., air to material).^[5] The Fresnel transmittance allows partial entry, while diffuse reflectance approximations (e.g., F_{dr} \approx -1.44/\eta^2 + 0.71/\eta + 0.668 + 0.066/\eta) model the average reflection for scattered light exiting the medium.^[5] In real-world observations, subsurface scattering intensifies at grazing angles due to longer internal path lengths, increasing interactions and diffusion; for example, dusty surfaces appear brighter and more diffuse when viewed obliquely compared to normal incidence, as scattering from subsurface particles dominates over direct reflection.^[8]

Theoretical Models

Diffusion Approximation

The diffusion approximation provides a mathematical simplification of the radiative transfer equation (RTE) for modeling multiple scattering of light in participating media, particularly where photons undergo numerous scattering events that randomize their directions. This approximation is derived by expanding the radiance in a low-order spherical harmonics basis (P1 approximation), assuming the angular distribution is nearly isotropic due to frequent scattering. Specifically, the radiance L(\mathbf{r}, \boldsymbol{\omega}) is approximated as L(\mathbf{r}, \boldsymbol{\omega}) \approx \frac{\phi(\mathbf{r})}{4\pi} + \frac{3\boldsymbol{\omega} \cdot \mathbf{J}(\mathbf{r})}{4\pi}, where \phi(\mathbf{r}) is the fluence and \mathbf{J}(\mathbf{r}) is the net flux vector. Substituting this expansion into the RTE and performing moment integrals over directions (zeroth and first moments) eliminates the angular dependence, yielding the diffusion equation under the assumptions of isotropic scattering (or mildly anisotropic with asymmetry factor g \approx 0) and low absorption relative to scattering (high albedo \omega = \sigma_s / (\sigma_a + \sigma_s) \approx 1).^[9]^[10] The resulting steady-state diffusion equation is

\nabla \cdot (-[D](/page/D*) \nabla \phi) + \sigma_a \phi = S,

where \phi is the scalar fluence (integrated radiance over all directions), \sigma_a is the absorption coefficient, S is the isotropic source term (typically representing light after the first scattering event plus any internal emission), and D is the diffusion coefficient defined as D = \frac{1}{3(\sigma_a + \sigma_s'(1 - g))}, with \sigma_s' the reduced scattering coefficient and g the asymmetry factor of the phase function.^[9] This model is valid in highly scattering media where the optical thickness \tau = (\sigma_a + \sigma_s) L \gg 1 (with L the characteristic medium size), ensuring the mean free path \ell = 1/(\sigma_a + \sigma_s) is much smaller than L (\ell / L < 10^{-4}) and light directions become effectively randomized after multiple interactions.^[9] The diffusion approximation was introduced to computer graphics by Jos Stam in 1995, who adapted it from radiative transfer theory to simulate multiple scattering in heterogeneous participating media, solving the equation numerically via finite element methods (including a blob basis approach) and multi-grid finite differences for efficiency.^[9] Its limitations include inaccuracy in low-scattering regimes (small \tau) where ballistic or single-scattered light dominates, or in highly directional media (low g) where the isotropic radiance assumption breaks down, leading to errors near boundaries or light entry points.^[9]

Dipole Model

The dipole model provides an analytical approximation to the diffusion equation for modeling subsurface scattering in semi-infinite homogeneous media, employing virtual point sources to enforce boundary conditions at the surface. It consists of a positive point source placed inside the medium to represent the inward light injection and a negative point source positioned outside to account for the extrapolated boundary, ensuring no flux escapes the domain unrealistically. This dipole configuration, originally derived in the context of biomedical optics, was adapted for computer graphics to enable efficient computation of light transport without numerical solving of partial differential equations.^[11] The formulation, introduced by Jensen et al., expresses the diffuse BSSRDF S_d(x_i, \omega_i; x_o, \omega_o) as the product of incoming and outgoing Fresnel transmittance terms with a radial profile R_d(r):

S_d(x_i, \omega_i; x_o, \omega_o) = \frac{1}{\pi} F_t(\eta, \omega_i) \, R_d(\|x_i - x_o\|) \, F_t(\eta, \omega_o)

where F_t is the Fresnel transmittance for light entering (\omega_i) and exiting (\omega_o) the medium with relative index of refraction \eta, and the profile is given by

R_d(r) = \frac{\alpha'}{4\pi} \left[ \frac{(\sigma_{tr} d_r + 1) e^{-\sigma_{tr} d_r}}{\sigma'_t d_r^3} + \frac{z_v (\sigma_{tr} d_v + 1) e^{-\sigma_{tr} d_v}}{\sigma'_t d_v^3} \right].

Here, d_r = \| x_o - x_r \| and d_v = \| x_o - x_v \| are distances from the exit point x_o to the real source x_r = (x_i, 0, z_r) at depth z_r = 1 / \sigma'_t and to the virtual source x_v = (x_i, 0, z_v) at z_v = z_r + 4 A D above the surface, with A = (1 + F_{dr}) / (1 - F_{dr}) as the dipole factor using diffuse reflectance F_{dr} and D = 1/(3 \sigma'_t) the mean diffusion length. The reduced albedo is \alpha' = \sigma_a / \sigma'_t, the effective transport extinction is \sigma_{tr} = \sqrt{3 \sigma_a \sigma'_t}, and \sigma'_t = \sigma_a + \sigma'_s with reduced scattering \sigma'_s = \sigma_s (1 - g), where \sigma_a, \sigma_s, and g are absorption, scattering coefficients, and the Henyey-Greenstein phase function asymmetry. The outgoing radiance is then L_o(x_o, \omega_o) = \int_{x_i, \omega_i} S(x_i, \omega_i; x_o, \omega_o) L_i(x_i, \omega_i) \cos \theta_i \, d\omega_i \, dx_i. This model primarily captures multiple scattering effects, though extensions can isolate single scattering for improved accuracy in low-albedo materials.^[5] Key profile parameters include the subsurface albedo \alpha', which controls the overall brightness of scattered light, and an effective scattering radius R \approx 1 / \sigma_{tr}, representing the typical distance light travels subsurface before absorption or escape, thus shaping the blur of incident light patterns like shadows or textures. These parameters allow fitting the model to measured scattering profiles of materials such as marble or skin using a few wavelengths for spectral rendering.^[5] Extensions of the dipole model include separable approximations for real-time rendering, where the BSSRDF is decomposed into directional and radial components evaluated via 1D convolutions on curved surfaces, achieving sub-millisecond performance for translucent materials. The model's closed-form nature facilitates its integration into production renderers, offering a balance of physical fidelity and computational efficiency over more general diffusion solvers.^[12]

Rendering Techniques

Monte Carlo Methods

Monte Carlo methods provide an unbiased approach to simulating subsurface scattering by stochastically tracing light paths through participating media, enabling physically accurate rendering of translucent materials. These techniques model the complex interplay of absorption and scattering events within a volume, integrating over all possible photon trajectories to compute outgoing radiance. By relying on random sampling, they avoid approximations inherent in analytical models, though at the expense of computational intensity.^[8]^[13] Random walk algorithms form the core of these simulations, where photon paths are traced probabilistically inside the medium. Starting from an entry point x_i with incident direction \omega_i, a random walk generates steps of length d = -\ln(r)/\sigma_t, with r uniformly sampled from [0,1] and \sigma_t = \sigma_a + \sigma_s the total extinction coefficient. At each scattering event, the direction is sampled according to a phase function, such as the Henyey-Greenstein function p(\theta) = \frac{1 - g^2}{4\pi (1 + g^2 - 2g \cos\theta)^{3/2}}, where g is the asymmetry parameter controlling forward or backward scattering bias; the path weight is updated as w \leftarrow w \cdot \sigma_s / \sigma_t. These walks continue until termination, accumulating contributions to estimate radiance at exit points x_o in direction \omega_o. This process directly simulates the volume scattering integral, capturing multiple bounces and spatial light diffusion.^[8]^[13] The bidirectional scattering-surface reflectance distribution function (BSSRDF), defined as S(x_i, \omega_i, x_o, \omega_o), encapsulates subsurface effects by relating incoming light at x_i, \omega_i to outgoing light at x_o, \omega_o, integrating over all internal paths connecting the points. In Monte Carlo rendering, the BSSRDF is not precomputed but evaluated on-the-fly through ensembles of random walks or path tracing from both entry and exit points, ensuring the full subsurface transport is accounted for without factorization assumptions. This integration allows Monte Carlo methods to handle arbitrary geometries and material properties, producing reference solutions for validating faster techniques.^[8]^[13] To mitigate high variance from rare scattering events, variance reduction techniques are essential. Importance sampling directs walks toward high-contribution directions by sampling from the phase function or an approximation of the scattering distribution, adjusting weights to maintain unbiasedness and reduce the number of samples needed for convergence. Russian roulette probabilistically terminates low-throughput paths with probability $1 - \min(1, w / w_{\max}), where w is the current weight and w_{\max} a threshold, reweighting survivors to preserve expectations and prevent infinite paths in low-albedo media. These methods significantly improve efficiency in subsurface simulations, where unbiased estimation otherwise demands extensive sampling.^[13]^[14] Monte Carlo methods for subsurface scattering saw early adoption in offline rendering during the 1990s, primarily for validating theoretical models of light transport in layered materials like skin and marble. Pioneering work integrated these techniques into ray tracers to compute BSSRDFs without closed-form solutions, establishing them as a gold standard for accuracy despite limitations in computational resources at the time.^[8]^[13] The primary drawback is high computational cost due to noisy estimates, as subsurface scattering involves diffuse, low-probability paths requiring millions of samples per pixel for visual convergence, especially in highly scattering materials with low absorption. Rendering times can span hours or days on modern hardware for complex scenes, limiting their use to production validation rather than interactive applications.^[13]^[8]

Analytical Approximations

Analytical approximations for subsurface scattering offer deterministic solutions to model light diffusion in translucent materials, providing efficiency over stochastic methods while maintaining reasonable fidelity. These techniques typically derive from diffusion theory, approximating the bidirectional scattering-surface reflectance distribution function (BSSRDF) using analytical expressions based on virtual point or ray sources. By solving simplified transport equations, they enable fast evaluation in rendering pipelines without requiring extensive precomputation or sampling. Seminal work in this area builds on the dipole approximation, extending it to handle layered structures and directional effects. One prominent approach is depth map-based subsurface scattering, which leverages screen-space depth buffers to estimate material thickness and accumulate irradiance. A depth map is rendered from the incident light's viewpoint, capturing the distance light travels through the geometry at entry and exit points. The effective path length s is computed as the difference between exit and entry depths, and the transmitted irradiance is approximated using an exponential attenuation: I(s) = I_0 \exp(-\sigma_t s), where \sigma_t is the total extinction coefficient and I_0 is the incident irradiance. This method approximates diffusion by integrating over nearby samples in the depth buffer, weighted by the falloff, and is particularly effective for capturing geometry-dependent scattering in opaque-translucent hybrids. The technique originates from early implementations in production rendering systems, balancing geometric awareness with computational simplicity. Texture space diffusion provides another efficient strategy by precomputing scattering profiles directly in the material's UV parameterization. Dipole-derived diffusion kernels, often fitted as sums of separable Gaussians (typically 4–8 terms for accuracy), are convolved with albedo textures in texture space to simulate local light redistribution. The convolution is performed separably along U and V directions using Gaussian filters scaled by scattering parameters, yielding a diffused irradiance map that can be sampled per-fragment or per-vertex. This approach excels for deformable surfaces like skin, as it avoids costly 3D convolutions and leverages texture hardware acceleration, producing results close to full dipole solutions at reduced cost. It was advanced in work on efficient skin rendering, where global scattering limitations are mitigated by hybrid combinations with other approximations.^[15] Multi-pole models extend the basic dipole by incorporating multiple virtual sources—such as quadrupoles or higher-order configurations—to better capture directionality and layered diffusion, improving accuracy for thin or anisotropic materials. For instance, the directional dipole model places positive and negative ray sources along the plane of incidence and emergence, deriving an analytical BSSRDF that accounts for the incident light direction relative to the normal. This formulation, S(\mathbf{x}_i, \omega_i, \mathbf{x}_o, \omega_o) \approx \sum_k c_k G(d_k, \sigma), where G is a Gaussian-like profile and d_k are source displacements, enhances translucency rendering without precomputation, outperforming standard multipoles in directional fidelity. Such extensions, rooted in multi-layered diffusion theory, allow fitting complex profiles like those in skin or marble.^[16]^[17] Integration of subsurface scattering with bidirectional reflectance distribution functions (BRDFs) creates hybrid models that unify surface and volume effects in a single shading framework. In these approaches, the subsurface term is added to the BRDF evaluation, often by replacing the Lambertian diffuse component with a normalized diffusion profile, R_d(r) = \frac{\exp(-r/d) + \exp(-r/(3d))}{8\pi d r (1 + g)}, scaled by material thickness and scattering distance d. This ensures energy conservation across reflection and transmission, as seen in production BSDFs where specular lobes handle surface interactions while subsurface modulates the base color. The resulting model supports path tracing and real-time evaluation, with parameters like subsurface transmission blending the contributions seamlessly.^[18] These analytical methods are well-suited for moderately complex scenes, offering a balance between visual accuracy and performance; for example, texture space diffusion achieves 30 frames per second on early GPUs for high-resolution models, while multi-pole evaluations add minimal overhead to standard shading. They prioritize efficiency in offline and interactive rendering, though they may require calibration for highly heterogeneous materials.^[15]

Real-Time Methods

Real-time methods for subsurface scattering (SSS) prioritize interactive performance in applications such as video games and virtual reality, leveraging GPU acceleration to approximate light transport within translucent materials like skin and marble. These techniques often employ screen-space processing or precomputation to achieve frame rates above 60 FPS on consumer hardware, balancing visual fidelity with computational efficiency.^[3] Screen-space subsurface scattering (SSSS) is a widely adopted approach that performs blurring operations in image space after the initial scene render. It begins by identifying material pixels (e.g., skin) using a stencil buffer, then generates albedo and thickness maps from the depth buffer and material properties. Separable Gaussian blurs are applied in multiple passes—typically three horizontal and three vertical—to convolve these maps, simulating light diffusion with varying kernel sizes corresponding to scattering distances. Curvature-aware falloff, derived from depth discontinuities, reduces scattering across surface edges to preserve geometry. This method, introduced by Jimenez et al., enables realistic translucency at interactive speeds on GPUs like the GeForce 8800 series, with costs scaling to screen coverage rather than geometry complexity.^[19]^[19] Precomputed radiance transfer (PRT) offers another GPU-friendly strategy by baking SSS effects into low-dimensional representations during a preprocessing phase. Light transport, including subsurface diffusion, is precomputed and stored as vertex colors or spherical harmonics coefficients, allowing real-time evaluation under dynamic lighting via dot products. For translucent materials, extensions to PRT incorporate BSSRDF models to capture view-dependent scattering, enabling relighting of complex scenes without per-frame ray tracing. This approach, as implemented in Direct3D, suits static or low-deformation geometry like character models, providing global illumination hints at minimal runtime cost.^[20]^[21] Recent advances integrate path tracing with resampling for real-time SSS. ReSTIR SSS (2024) adapts the ReSTIR framework to handle subsurface paths using reservoir-based sampling, combining spatial and temporal resampling to reduce noise while supporting dynamic scenes. This GPU-accelerated method achieves path-traced quality at interactive rates on modern hardware, extending beyond screen-space limitations for materials with high scattering anisotropy.^[22] A 2025 development presented at SIGGRAPH Advances in Real-Time Rendering introduces a hybrid real-time SSS method combining volumetric path tracing for single scattering with a physically based diffusion profile for multiple scattering, integrated with ReSTIR for efficient sample reuse and noise reduction. Developed by NVIDIA, this approach achieves near path-traced quality at interactive frame rates, particularly enhancing fidelity in thin or curved translucent materials such as digital humans.^[23] Modern GPUs facilitate these methods through parallel multi-pass rendering pipelines. NVIDIA's implementations in GPU Gems demonstrate efficient blurring via fragment shaders, exploiting texture units for separable convolutions and achieving 10-20 ms per frame for SSS on early programmable hardware.^[3] Despite their efficiency, real-time SSS techniques introduce trade-offs, including artifacts like light leaking where illumination bleeds across occluded boundaries due to incomplete geometry awareness in screen space. Thickness maps mitigate this by weighting scattering contributions based on estimated material depth, attenuating effects on thin regions such as ears or edges; curvature falloff further curbs haloing around silhouettes. These approximations, while not physically exact, yield compelling results for interactive rendering when tuned for specific materials.^[19]^[24]^[25]

Applications

Materials Simulation

Subsurface scattering parameters are tuned to replicate the optical properties of specific materials by adjusting the scattering coefficient \sigma_s, absorption coefficient \sigma_a, and phase function anisotropy factor g. For human skin, \sigma_s' (the reduced scattering coefficient) typically ranges from 0.74 mm^{-1} (red) to 1.01 mm^{-1} (blue), with \sigma_a values such as 0.032 mm^{-1} for red wavelengths to allow deeper penetration and produce reddish tones, while g \approx 0.85 accounts for forward-biased scattering in dermal tissues.^[26] In contrast, marble exhibits higher \sigma_s' values of 2.19–3.00 mm^{-1} across RGB channels and very low \sigma_a (0.0021–0.0071 mm^{-1}), enabling deeper light penetration and a more uniform translucency compared to skin's selective absorption.^[26] Measured data for these parameters are obtained through inverse rendering techniques applied to photographs or controlled illumination setups, where radial reflectance profiles are fitted to diffusion models to estimate \sigma_s and \sigma_a.^[26] For instance, Jensen et al. used a laser beam and high-dynamic-range imaging on skin and marble samples to derive these coefficients by matching simulated BSSRDFs to captured data.^[26] Goniophotometers complement this by measuring the scattering phase function on thin tissue sections, yielding the anisotropy g through angular distribution analysis of backscattered light.^[27] Material-specific effects arise from layered or heterogeneous structures that influence scattering patterns. In skin, subdermal blood layers in the dermis, modeled with hemoglobin concentrations of 0.001–0.1 and 75% oxygenation, cause color shifts toward red due to stronger absorption of shorter wavelengths, resulting in a warm subsurface glow.^[28] Marble's veining stems from mineral inclusions that create localized variations in \sigma_a and \sigma_s, simulating irregular translucency and color mottling as light scatters through calcite veins.^[26] Similarly, subcutaneous fat layers in skin feature lower \sigma_s relative to dermal values and minimal absorption, enhancing overall translucency and softening edges in simulations.^[28] Hyperspectral imaging facilitates acquisition of wavelength-dependent parameters by capturing full spectral reflectance from subsurface interactions, enabling precise fitting of \sigma_s(\lambda) and \sigma_a(\lambda) for materials across 400–1000 nm.^[29] This technique reveals variations such as power-law scattering \sigma_s'(\lambda) = 14.74\lambda^{-0.22} + 2.2 \times 10^{11}\lambda^{-4} mm^{-1} in dermal tissues, supporting accurate multi-spectral simulations.^[28]

Use in Production

Subsurface scattering (SSS) has evolved significantly in computer graphics production pipelines, transitioning from computationally intensive offline rendering techniques in early 2000s films to efficient real-time approximations in 2010s video games, enabling broader adoption in interactive media.^[30]^[19] In film production, Pixar employed SSS extensively in Monsters University (2013) to render realistic skin and subsurface effects on characters, including point-based pre-pass solutions for translucent materials like Dean Hardscrabble’s wings, which captured light penetration and scattering for enhanced visual depth.^[31] Disney Animation Studios advanced SSS through path-traced methods integrated with their physically based BSDF model, replacing traditional diffusion approximations to better handle geometric details in translucent surfaces such as skin, as demonstrated in production rendering workflows.^[32] Video games have leveraged real-time SSS for organic materials, with Naughty Dog implementing it in The Last of Us Part II (2020) to achieve high-fidelity skin rendering through standard subsurface effects combined with detailed textures and models, contributing to photorealistic character appearances.^[33] Similarly, Uncharted 4: A Thief's End (2016) utilized SSS on characters to simulate light refraction and penetration in skin, allowing realistic bending and movement during animations while maintaining performance on console hardware.^[34] Engine-specific shaders, such as Unreal Engine 5's Subsurface Shading Model, facilitate these implementations for materials like skin and wax, enabling artists to control scattering profiles for translucent effects in real-time game environments.^[35] Production challenges with SSS include balancing computational budgets against visual fidelity, as simulating light interactions increases render times in offline pipelines and strains GPU resources in real-time scenarios.^[36] Artist workflows often involve creating specialized SSS maps and selecting appropriate shaders to tune scattering distances and colors, requiring iterative testing to avoid artifacts like unnatural translucency while preserving material realism.^[36] Recent tools incorporate AI-assisted SSS, such as neural network models trained via optimal transport to approximate scattering in granular materials, streamlining synthesis for complex scenes in production rendering.^[37] These advancements, including neural representations for subsurface effects, further reduce computational overhead while maintaining high-quality results in both film and game pipelines.^[38] As of 2025, further progress includes real-time path-traced SSS techniques like ReSTIR and hybrid volumetric methods, presented at SIGGRAPH 2025, enhancing interactive rendering in games and virtual production.^[39]

References

[1]
What Is Subsurface Scattering? - ITU Online IT Training
Subsurface scattering (SSS) is a light interaction phenomenon where light penetrates a translucent material, scatters internally, and exits at a different point ...
[2]
Images with Subsurface Scattering by Henrik Wann Jensen
Subsurface scattering is when light scatters inside a translucent material before being absorbed or leaving at a different location.Missing: computer definition
[3]
Chapter 16. Real-Time Approximations to Subsurface Scattering
This chapter describes several methods of approximating the look of translucent materials, such as skin and marble, using programmable graphics hardware.Chapter 16. Real-Time... · 16.3 Simulating Absorption... · 16.4 Texture-Space DiffusionMissing: definition | Show results with:definition
[4]
[PDF] Real-Time Approximate Subsurface Scattering - UMD CS
Abstract— Subsurface scattering is a complex physical process, which, in many cases, significantly affects the appearance of certain materials.
[5]
[PDF] A Practical Model for Subsurface Light Transport
In this paper we have presented a new practical BSSRDF model for computer graphics. The model combines a dipole diffusion ap- proximation with an accurate ...Missing: seminal | Show results with:seminal
[6]
[PDF] A Survey on Participating Media Rendering Techniques - Hal-Inria
Oct 13, 2010 · It reviews both appearance-based and physically-based media methods, including the single scattering and the more general multiple-scattering ...<|control11|><|separator|>
[7]
[PDF] Reflection from Layered Surfaces due to Subsurface Scattering
The scattering of light from particles small compared to the wavelength of light is given by the Rayleigh scatter- ing formula, and the scattering due to ...
[8]
[PDF] Multiple scattering as a diffusion process - Dynamic Graphics Project
We implemented both a multi-grid finite differences scheme and a finite-element ... In this situation, a diffusion-type equation is obtained for the scattered.
[9]
A diffusion theory model of spatially resolved, steady-state ... - PubMed
A diffusion theory model of spatially resolved, steady-state diffuse reflectance for the noninvasive determination of tissue optical properties in vivo. Med ...
[10]
Separable Subsurface Scattering - Jimenez - Wiley Online Library
Feb 9, 2015 · In this paper, we propose two real-time models for simulating subsurface scattering for a large variety of translucent materials, which need under 0.5 ms per ...
[11]
[PDF] monte carlo solution of scattering equations for computer graphics a ...
We show how Monte Carlo techniques can be applied to solve these adding equations and apply them to the problem of rendering subsurface scattering. v. Page 6 ...
[12]
https://onlinelibrary.wiley.com/doi/10.1111/cgf.12529
[13]
[PDF] Efficient Rendering of Human Skin - Eugene d'Eon
Texture-space diffusion captures highly diffuse local scatter- ing typical in skin, but can miss global scattering between regions that are close in ...
[14]
Light diffusion in multi-layered translucent materials
This paper introduces a shading model for light diffusion in multi-layered translucent materials. ... https://doi.org/10.26599/CVM.2025.9450360. Show More ...
[15]
Directional Dipole Model for Subsurface Scattering
We present an improved analytical model for subsurface scattering that captures translucency effects present in the reference solutions but remaining absent ...
[16]
[PDF] Extending the Disney BRDF to a BSDF with Integrated Subsurface ...
We introduced a new physically based shading model on Wreck-It Ralph [Bur12], and we used this single, general-purpose BRDF on all materials (except hair).
[17]
Screen-Space Subsurface Scattering - Jorge Jimenez
Apr 6, 2011 · Subsurface scattering calculations are done at the required resolution, which can be seen as a SSS level of detail. The cost depends on the ...
[18]
Precomputed Radiance Transfer (Direct3D 9) - Win32 apps
Jan 6, 2021 · PRT with Subsurface Scattering. Subsurface scattering is a technique that models how light passes through certain materials. As an example ...
[19]
[PDF] Precomputed Radiance Transfer - Jan Kautz
The main trick we are going to use for precomputed radiance transfer (PRT) is to combine ... – Subsurface-Scattering: • View-dependent effects require transfer.
[20]
[PDF] ReSTIR Subsurface Scattering for Real-Time Path Tracing
ReSTIR SSS is a specialized technique to render paths with SSS. ... Efficient screen-space subsurface scattering using Burley's normalized diffusion in real-time.
[21]
Interactive Subsurface Scattering for Materials With High Scattering ...
Jul 26, 2020 · ... light bleeding can be noticed. ... Existing screen space subsurface scattering algorithms cannot be applied if high scattering distances are ...
[22]
Digital Humans | Unreal Engine 4.27 Documentation
Subsurface scattering is a screen space effect that would often cause halo-type artifacts in areas with other subsurface materials, similar to the example ...
[23]
[PDF] A Practical Model for Subsurface Light Transport - UCSD CSE
Only methods that consider subsurface scattering can capture the true appearance of translucent materials, such as mar- ble, cloth, paper, skin, milk, cheese, ...
[24]
[PDF] An Introduction to Light Interaction with Human Skin - CORE
In fact, Jacques et al., [43] showed through goniophotometric measurements that backscattered light from the dermis is diffuse. While Mie scattering produces ...
[25]
[PDF] A Spectral BSSRDF for Shading Human Skin - UC San Diego
The model uses four parameters: oil, melanin, and hemoglobin, to simulate light interaction with skin, accounting for subsurface and surface scattering.
[26]
[PDF] Spectral Subsurface Scattering for Material Classification
function p(θ), the total scattering coefficient σs, and the total extinction coeffi- cient σt [13]. All of these parameters are expressed as functions of ...
[27]
[PDF] Physically Based Shading at DreamWorks Animation - Self Shadow
Fiona was the first animated CG character to use subsurface scattering on her skin. Figure 3: Shrek 2 was the first PDI/DreamWorks film to use single bounce ...
[28]
Monsters University: rendering physically based monsters - fxguide
Jun 24, 2013 · How new developments in RenderMan helped Pixar deliver Monsters University ... subsurface scattering techniques. But as he had also received a ...
[29]
[PDF] Practical and Controllable Subsurface Scattering for Production Path ...
Jul 24, 2016 · We have integrated path-traced subsurface scattering with our BSDF [Burley 2015] by replacing the non-directional (i.e. Lam- bertian) diffuse ...
[30]
The Last of Us Part 2 Graphics Analysis – Naughty Dog Continues ...
Oct 3, 2019 · Subsurface scattering has become a standard feature of late, this, combined with the high quality texture work and high-poly models leads to ...
[31]
Tech Analysis: Uncharted 4: A Thief's End | Digital Foundry
Sep 13, 2016 · ... game giving extra depth · This scene demonstrates local light shafts as well as the use of subsurface scattering on characters · The majority of ...
[32]
Using Subsurface Scattering in Unreal Engine Materials
Subsurface scattering is light scattering through translucent surfaces. Unreal Engine uses a Subsurface Shading Model for materials like skin or wax. Enable it ...
[33]
Subsurface scattering explained - The Chaos Blog
Aug 20, 2025 · Subsurface scattering (SSS) is crucial for achieving photorealism in 3D rendering by simulating how light interacts with translucent ...
[34]
Learning subsurface scattering solutions of tightly-packed granular ...
In the Sandbox scene, the error is first dominated by Monte Carlo noise for lower sample counts and later converges to the bias that is very similar for both ...<|control11|><|separator|>
[35]
Subsurface Scattering for 3D Gaussian Splatting - arXiv
Figure 1: SSS GS – We propose photorealistic real-time relighting and novel view synthesis of subsurface scattering objects. We learn to reconstruct the ...