--- title: Shadow Accrual Maps in Urban and Gravitational Contexts url: https://www.emergentmind.com/topics/shadow-accrual-maps type: topic --- # Shadow Accrual Maps in Urban and Gravitational Contexts Shadow Accrual Maps designate two distinct constructions in recent literature. In computer graphics and urban analysis, they provide an efficient way to accumulate, analyze, and visualize city-scale shadows over user-defined periods without brute-force sampling at dense time steps. In gravitational theory, the same term denotes a linear kernel transform that accrues compact, conserved thermodynamic or effective-medium perturbations of \(\Delta T^{\mu}{}_{\nu}\) into shifts of black-hole shadow observables such as the shadow radius and photon-sphere frequency. The shared expression therefore refers not to a single methodology, but to two mathematically structured accrual frameworks operating in different domains: temporal accumulation of urban direct-sun shadow occupancy, and linear-response propagation from matter susceptibilities to geometric optics observables [1907.04435] [2601.17613]. ## 1. Terminological scope and conceptual commonality The city-scale formulation was introduced in "Shadow Accrual Maps: Efficient Accumulation of City-Scale Shadows Over Time" and published in 2019. It addresses the problem of determining the effects of building shadows across different periods in a year, with applications including thermal comfort, vegetation growth, planning and zoning, right-to-light, and solar access. Its core idea is to exploit short-interval coherence in solar motion so that shadow movement can be tracked between interval endpoints rather than sampled minute by minute [1907.04435]. The black-hole formulation appeared in 2026 in "From Thermodynamic Criticality to Geometric Criticality: A Linear Kernel Map from Matter Susceptibilities to Black-Hole Shadows." There, a Shadow Accrual Map is an explicit linear map from compact, conserved thermodynamic or effective-medium perturbations of the stress-energy tensor to the metric response, and from there to the shadow radius and photon-sphere frequency. The response is written through \(L^{1}\)-bounded kernels in a piecewise "local + tail" form [2601.17613]. The common conceptual element is accrual. In the urban setting, accrual refers to integrating or summing direct-shadow occupancy over time on receivers such as ground planes, parks, and façades. In the gravitational setting, accrual refers to integrating radial source perturbations against kernels that encode sensitivity of the shadow observable. This suggests that the phrase names a structural pattern—aggregation by a map from elementary shadow-generating contributions to a cumulative observable—rather than a domain-specific algorithm. ## 2. Urban Shadow Accrual Maps: mathematical formulation and solar-geometry model In the urban literature, the central quantity is the direct-shadow indicator \(I(x,t)\in\{0,1\}\) at location \(x\) and time \(t\). Accrual over a time set \(T\) with temporal weights \(w(t)\) is defined by \[ S(x)=\int_{t\in T} w(t)\,I(x,t)\,dt, \] with discrete approximation \[ S(x)=\sum_{k=1}^{K} w_k\,I(x,t_k), \qquad \widehat{S}(x)=\frac{S(x)}{\sum_{k=1}^{K} w_k}. \] The paper distinguishes gross shadow, meaning total time or fraction in shadow over \(T\), from continuous shadow, meaning the longest contiguous duration in shadow over \(T\) [1907.04435]. The computational difficulty arises because a naïve baseline computes shadows for many time steps. The paper gives the example that 3 hours per day over a week at 1-minute resolution implies 1260 sun directions. At city scale, with thousands of buildings and millions of triangles, such dense temporal sampling is prohibitively expensive and impedes interactivity. The key observation is that over short intervals of approximately 1 hour, the sun’s direction changes nearly linearly, so shadow occupancy can be obtained by tracking boundary motion and effectively “sweeping” shadows across the receiver. The solar model is stated in conventional astronomical terms. For day-of-year \(n\), approximate declination is \[ \delta(n)\approx 23.44^\circ \sin\!\left(\frac{2\pi}{365}(284+n)\right), \] and with local solar time \(\mathrm{LST}\), the hour angle is \[ H=15^\circ(\mathrm{LST}-12). \] Given latitude \(\varphi\), solar zenith is \[ \theta_z=\arccos\!\big(\sin\varphi\,\sin\delta+\cos\varphi\,\cos\delta\,\cos H\big), \] elevation is \(\alpha=90^\circ-\theta_z\), and the sun direction vector is \[ \mathbf{s}(t)=[\cos\alpha\cos A,\ \cos\alpha\sin A,\ \sin\alpha]^\top. \] A notable empirical result is that, in New York City, the cosine similarity between actual and interpolated sun directions is approximately \(1.0\) for intervals up to 60 minutes, with mean approximately \(0.99999\) and \(\sigma \approx 4.9\times 10^{-6}\) at 60 minutes; the error grows beyond an hour. The method therefore segments analyses into short intervals and uses linear interpolation within each one. Geometrically, hard shadows are cast by directional light onto a receiver that may be a ground plane \(z=0\) or an arbitrary plane. For a scene point \(\mathbf{s}\) on an occluder, its shadow projection on the receiver along \(-\mathbf{s}(t)\) is obtained from the parametric line \[ \mathbf{r}(\lambda)=\mathbf{s}-\lambda\,\mathbf{s}(t), \qquad \lambda\ge 0. \] The swept region \(\bigcup_{t\in T} p(t)\) traces a curve that is approximately straight under linear sun motion. Rasterization of this sweep yields occupancy across the interval, and occlusions are resolved by depth tests that keep the farthest depth per texel per slice. ## 3. Extensions to shadow maps and ray tracing The 2019 paper develops two extensions: a shadow-map-based SAM and a ray-tracing-based Inverse Accrual Maps (IAM) method. In the shadow-map formulation, conventional shadow mapping would compute a separate shadow map per time step, whereas SAM computes a single 3D texture whose slices correspond to all steps in a short interval in one pass. For each scene point, the method projects to \(p_1\) and \(p_n\) on a shadow plane orthogonal to the start direction \(d_1\), and for intermediate slice \(i\) uses \[ p_i=p_1+(p_n-p_1)\,\frac{\tan(i'\theta/n)}{\tan(\theta)}, \] where \(\theta=\angle(\vec d_1,\vec d_n)\) and \(i'=i-1\). The paper states a theorem that, under linear sun movement, the 3D accrual texture is equivalent to computing \(n\) shadow maps at each time step; SAM consolidates \(n\) passes into one pass with no quality loss [1907.04435]. The shadow-map implementation uses image atomics in OpenGL 4.3 to store the maximum depth per texel per slice in one rendering pass. Trapezoidal shadow map warping improves texel density near the view frustum, and a depth bias mitigates shadow acne. Gross and continuous shadow are then computed from per-pixel bitvectors over slices. IAM addresses flat receivers such as parks and sidewalks. Its motivation is that ray tracing yields higher-quality results and avoids shadow map aliasing. The method traces, for each receiver pixel at \(t_1\), rays along \(-d_1\) to record up to \(l\) intersection points, each treated as a potential shadow source. The mapping from \(p_1\) to \(p_n\) is stored in a 3D texture whose slice \(j\) corresponds to the \(j\)-th closest source. To capture stretching and contracting shadow behavior, the method computes IAM both forward \((t_1\rightarrow t_n)\) and reverse \((t_n\rightarrow t_1)\). Accrual then proceeds by line rasterization: if a subsegment corresponding to time step \(i\) crosses a pixel, the bit \(B[i]\) is set to 1. The source-level tradeoff is explicit. Larger \(l\) increases accuracy at dawn and dusk, when many occluders may matter, but increases cost. The paper reports that \(l \le 3\) suffices with mean area error less than \(1\%\) and maximum error approximately \(2.6\%\), with the worst cases at dawn and dusk. This places IAM as a higher-fidelity method for flat receivers, while SAM remains the interactive method for broader urban exploration. ## 4. Implementation, validation, and urban-planning use The methods were implemented in the interactive visual analysis system Shadow Profiler, targeted at city planners and architects. The reported architecture combines C++, OpenGL 4.3 for GPU shadow accrual maps and image atomics, OpenCL 1.2 for bit operations and progressive accumulation, GLSL shaders for IAM ray traversal, and a 3D grid acceleration structure for ray tracing. The interface includes a 3D map widget, date and time selection, gross or continuous accumulation type, averaging over days, progressive computation and visualization, polygonal ROI selection, and a divergent colormap for positive or negative differences in minutes. Two modes are defined: exploration, based on SAM, and analysis, based on IAM [1907.04435]. The Manhattan case study uses OpenStreetMap geometry extruded to heights, with approximately 1.5 million triangles across 43k buildings. Experiments were conducted at resolutions from \(800\times 800\) to \(2048\times 2048\), with shadow-map resolutions up to \(2048\times 2048\), on a workstation with Intel Xeon E5-2620, 128 GB RAM, and an NVIDIA GTX 1080 8 GB. At \(1024\times 1024\), SAM uses approximately 240 MB and IAM approximately 56 MB; memory scales proportionally with resolution, reaching 960 MB for SAM and 224 MB for IAM at \(2048\times 2048\). The reported validation covers linearity of sun motion, clustering of directions, source truncation in IAM, and performance. With a direction graph using bins bounded by \(5^\circ\), mean direction similarity is \(0.9996\) with \(\sigma \approx 3\times 10^{-4}\). Shadow area computed with cluster representatives has mean absolute error approximately \(0.47\%\) of total area, median approximately \(0.35\%\), and approximately \(0.8\%\) of pixels misclassified as shadow or non-shadow on average. For performance over 6 hours per day, 10 random positions, and 20 random days, SAM is approximately \(10\times\) faster than minute-sampled shadow maps, and IAM is approximately \(5.3\times\) faster than minute-sampled ray tracing. For annual accumulation over 6 hours per day, the direction graph reduces hourly-interval computations from 2190 to 299 edges, more than \(7\times\) fewer, and end-to-end speedups over \(50\times\) are observed with increasing time periods. The planning-oriented case studies center on Manhattan. For Central Park skyscraper scenarios, the baseline impact of the Time Warner Center is described as localized, while new supertall towers cast widespread long shadows but keep the net increase above 30 minutes concentrated. An alternate scenario with shorter, wider towers of equal total area yields stronger near-base impact, especially under Boston-style thresholds of at least 60 minutes. The paper also reports that most Manhattan neighborhoods average more than 4 hours per day in shadow due to density, while wide streets, plazas, and parks fare better. Weighted “shadow score” maps are mostly negative overall in dense areas, with pockets of positive scores in lower-density regions. This suggests that the method supports comparisons between summer shade benefits and winter shadow penalties in planning practice. ## 5. Assumptions, limitations, and relation to adjacent urban methods The urban formulation makes several assumptions explicit. It models hard shadows only; penumbrae and soft-shadow integration are not included. It models direct sunlight only; diffuse skylight, atmospheric scattering, and refraction are not considered. IAM assumes flat receivers, so arbitrary topography and vertical receivers require extensions, including hybrid approaches or Monte Carlo. Accuracy declines at dawn and dusk because of rapid shadow motion and the presence of multiple occluders when \(l\) is small. The paper also notes sensitivity to user-defined temporal weights \(w(t)\), so planning conclusions must be contextualized by the weighting scheme [1907.04435]. The method is also conditioned by data quality and implementation choices. Building geometry errors, missing heights, and overhangs affect accuracy. Shadow-map artifacts require warping and bias tuning. The empirical validation of 60-minute segmentation is specific to New York City, and other latitudes or climates may require parameter tuning. Practical guidance in the paper therefore includes hourly intervals per short segment, direction clustering with at most \(5^\circ\) bins for multi-day or seasonal analyses, and spatial resolutions around \(800\times 800\) to \(1024\times 1024\) for city-scale exploration. For parks and streets, IAM at \(1024\times 1024\) with \(l=3\) is recommended in the reported experiments. Relative to earlier approaches, the contrast is methodological. Traditional shadow mapping requires separate shadow maps per time step, while SAM uses a single pass to populate all slices by exploiting sun-motion coherence. Minute-by-minute ray tracing is expensive, while IAM maps shadow-source events and rasterizes lines instead of sampling every minute. Voxel occupancy and precomputed occlusion intervals involve heavier precomputation and are less flexible under scene changes. Fernando’s occlusion interval maps are identified as prior art, but SAM and IAM are event-based and reuse direction bins across days. Deep shadow maps, adaptive or warped shadow maps, and shadow volumes remain per-light computations; the paper locates SAM’s novelty in leveraging time-coherent sun movement to accelerate temporal accumulation. A common misconception is that SAM is a general daylighting model. The cited formulation is narrower: it is a direct-sun, hard-shadow accumulation framework. Another misconception is that IAM is a universal replacement for shadow maps. In the paper’s own system design, SAM is the interactive approach for broad regions, whereas IAM is the higher-fidelity analysis mode for flat receivers. ## 6. Shadow Accrual Map in black-hole optics and critical phenomena In the 2026 gravitational literature, the term denotes a linear kernel functional rather than a rendering algorithm. The background is a static, spherically symmetric spacetime in areal-radius gauge with \[ ds^{2}=-f(r;\lambda)\,dt^{2}+h(r;\lambda)\,dr^{2}+r^{2}d\Omega^{2}, \qquad f=e^{2\Phi}A,\qquad h=A^{-1}, \] where \[ A(r;\lambda)=1-\frac{2\,m(r;\lambda)}{r}-\frac{\Lambda r^{2}}{3}. \] Linearization is performed around a vacuum background with \(m_0(r)=M\), \(\Phi_0(r)=0\), and \(f_0(r)=A_0(r)\). Matter perturbations are diagonal, \[ \Delta T^{\mu}{}_{\nu}=\mathrm{diag}\!\big(-\Delta\rho,\ \Delta p_r,\ \Delta p_t,\ \Delta p_t\big), \] with compact support in a shell \([a,b]\subset(2M,\infty)\), smooth profiles, and the conservation constraint \(\nabla_\mu \Delta T^\mu{}_\nu=0\) [2601.17613]. The optical observables are defined through the photon sphere. An unstable circular photon orbit satisfies \[ F(r):=r f'(r)-2f(r)=0,\qquad F'(r_{\rm ph})>0, \] and the corresponding shadow radius for a static observer at infinity is \[ R_{\rm sh}=\frac{r_{\rm ph}}{\sqrt{f(r_{\rm ph})}}. \] The photon-sphere frequency is \[ \Omega_{\rm ph}=\frac{\sqrt{f(r_{\rm ph})}}{r_{\rm ph}}=\frac{1}{R_{\rm sh}}. \] Linearizing Einstein’s equations yields \[ \delta m'(r)=4\pi r^{2}\,\Delta\rho(r), \qquad \delta\Phi'(r)=\frac{\delta m(r)+4\pi r^{3}\,\Delta p_r(r)}{r^{2}A_0(r)}. \] After integration by parts, the metric response takes the integral-kernel form \[ \delta f(r)=\int_{0}^{\infty}\!\Big[\mathcal{K}_{\rho}(r,\bar r)\,\Delta\rho(\bar r)+\mathcal{K}_{p}(r,\bar r)\,\Delta p_r(\bar r)\Big]\,d\bar r. \] Propagation to the shadow shift then gives the Shadow Accrual Map \[ \delta R_{\rm sh}=\int_{0}^{\infty}\!\Big[K_{\rho}^{\rm sh}(\bar r)\,\Delta\rho(\bar r)+K_{p}^{\rm sh}(\bar r)\,\Delta p_r(\bar r)\Big]\,d\bar r. \] The kernels are written in a "local + tail" decomposition. Local contributions arise from mass loading and from the pressure term entering through \(\partial_r\mathcal{K}_p\) for \(\bar r