Recent Advances in Adaptive Sampling and Reconstruction for Monte Carlo Rendering

Monte Carlo integration is firmly established as the basis for most practical realistic image synthesis algorithms because of its flexibility and generality. However, the visual quality of rendered images often suffers from estimator variance, which appears as visually distracting noise. Adaptive sampling and reconstruction algorithms reduce variance by controlling the sampling density and aggregating samples in a reconstruction step, possibly over large image regions. In this paper we survey recent advances in this area. We distinguish between “a priori” methods that analyze the light transport equations and derive sampling rates and reconstruction filters from this analysis, and “a posteriori” methods that apply statistical techniques to sets of samples to drive the adaptive sampling and reconstruction process. They typically estimate the errors of several reconstruction filters, and select the best filter locally to minimize error. We discuss advantages and disadvantages of recent state-of-the-art techniques, and provide visual and quantitative comparisons. Some of these techniques are proving useful in real-world applications, and we aim to provide an overview for practitioners and researchers to assess these approaches. In addition, we discuss directions for potential further improvements.

Regularizing Image Reconstruction for Gradient-Domain Rendering with Feature Patches

We present a novel algorithm to reconstruct high-quality images from sampled pixels and gradients in gradient-domain rendering. Our approach extends screened Poisson reconstruction by adding additional regularization constraints. Our key idea is to exploit local patches in feature images, which contain per-pixels normals, textures, position, etc., to formulate these constraints. We describe a GPU implementation of our approach that runs on the order of seconds on megapixel images. We demonstrate a significant improvement in image quality over screened Poisson reconstruction under the L1 norm. Because we adapt the regularization constraints to the noise level in the input, our algorithm is consistent and converges to the ground truth.

Denoising Your Monte Carlo Renders: Recent Advances in Image Space-Adaptive Sampling and Reconstruction

With the ongoing shift in the computer graphics industry toward Monte Carlo rendering, there is a need for effective, practical noise-reduction techniques that are applicable to a wide range of rendering effects and easily integrated into existing production pipelines. This course surveys recent advances in image-space adaptive sampling and reconstruction algorithms for noise reduction, which have proven very effective at reducing the computational cost of Monte Carlo techniques in practice. These approaches leverage advanced image-filtering techniques with statistical methods for error estimation. They are attractive because they can be integrated easily into conventional Monte Carlo rendering frameworks, they are applicable to most rendering effects, and their computational overhead is modest.

A first chironomid-based summer temperature reconstruction (13–5 ka BP) around 49°N in inland Europe compared with local lake development

Temperature reconstructions for the end of the Pleistocene and the first half of the Holocene based on biotic proxies are rare for inland Europe around 49°N. We analysed a 7 m long sequence of lake deposits in the Vihorlat Mts in eastern Slovakia (820 m a.s.l.). Chironomid head capsules were used to reconstruct mean July temperature (TJuly), other proxies (diatoms, green algae, pollen, geochemistry) were used to reconstruct local environmental changes that might have affected the climate reconstruction, such as epilimnetic total phosphorus concentrations (TP), lake level changes and development of surrounding vegetation. During the Younger Dryas (YD), temperature fluctuated between 7 and 11 °C, with distinct, decadal to centennial scale variations, that agree with other palaeoclimate records in Europe such as δ18O content in stalagmites or Greenland ice cores. The results indicate that the site was somewhat colder than expected from the general south-to-north YD temperature gradient within Europe, possibly because of north-facing exposition. The warmer phases of the YD were characterised by low water level or even complete desiccation of the lake (12,200-12,400 cal yr BP). At the Late-Glacial/Holocene transition TJuly steeply increased from from 11 to 15.5 °C (11,700-11,400 cal yr BP) - the highest TJuly for entire sequence. This rapid climate change was reflected by all proxies as a compositional change and increasing species diversity. The open woodlands of Pinus, Betula, Larix and Picea were replaced by broad-leaved temperate forests dominated by Betula, later by Ulmus and finally by Corylus (ca 9700 cal yr BP). At the same time, input of eroded coarse-grained material into the lake decreased and organic matter (LOI) and biogenic silica increased. The Early-Holocene climate was rather stable till 8700 cal yr BP, with temporary decrease in TJuly around 11,200 cal yr BP. The lake was productive with a well-developed littoral, as indicated by both diatoms and chironomids. A distinct decline of TJuly to 10 °C between 8700 and 8000 cal yr BP was associated with decreasing chironomid diversity and increasing climate moistening indicated by pollen. Tychoplanktonic and phosphorus-demanding diatoms increased which might be explained by hydrological and land-cover changes. Later, a gradual warming started after 7000 cal yr BP and representation of macrophytes, periphytic diatoms and littoral chironomids increased. Our results suggest that the Holocene thermal maximum was taking place unusually early in the Holocene at our study site, but its timing might be affected by topography and mesoclimate. We further demonstrated that temperature changes had coincided with variations in local hydrology

Modern diatom, cladocera, chironomid, and chrysophyte cyst assemblages as quantitative indicators for the reconstruction of past environmental conditions in the Alps. II. Nutrients

Surface sediments from 68 small lakes in the Alps and 9 well-dated sediment core samples that cover a gradient of total phosphorus (TP) concentrations of 6 to 520 μg TP l-1 were studied for diatom, chrysophyte cyst, cladocera, and chironomid assemblages. Inference models for mean circulation log10 TP were developed for diatoms, chironomids, and benthic cladocera using weighted-averaging partial least squares. After screening for outliers, the final transfer functions have coefficients of determination (r2, as assessed by cross-validation, of 0.79 (diatoms), 0.68 (chironomids), and 0.49 (benthic cladocera). Planktonic cladocera and chrysophytes show very weak relationships to TP and no TP inference models were developed for these biota. Diatoms showed the best relationship with TP, whereas the other biota all have large secondary gradients, suggesting that variables other than TP have a strong influence on their composition and abundance. Comparison with other diatom – TP inference models shows that our model has high predictive power and a low root mean squared error of prediction, as assessed by cross-validation.

New algorithms for automated phylogenetic reconstruction using artificial intelligence and data mining techniques

Academic and industrial research in the late 90s have brought about an exponential explosion of DNA sequence data. Automated expert systems are being created to help biologists to extract patterns, trends and links from this ever-deepening ocean of information. Two such systems aimed on retrieving and subsequently utilizing phylogenetically relevant information have been developed in this dissertation, the major objective of which was to automate the often difficult and confusing phylogenetic reconstruction process. ^ Popular phylogenetic reconstruction methods, such as distance-based methods, attempt to find an optimal tree topology (that reflects the relationships among related sequences and their evolutionary history) by searching through the topology space. Various compromises between the fast (but incomplete) and exhaustive (but computationally prohibitive) search heuristics have been suggested. An intelligent compromise algorithm that relies on a flexible “beam” search principle from the Artificial Intelligence domain and uses the pre-computed local topology reliability information to adjust the beam search space continuously is described in the second chapter of this dissertation. ^ However, sometimes even a (virtually) complete distance-based method is inferior to the significantly more elaborate (and computationally expensive) maximum likelihood (ML) method. In fact, depending on the nature of the sequence data in question either method might prove to be superior. Therefore, it is difficult (even for an expert) to tell a priori which phylogenetic reconstruction method—distance-based, ML or maybe maximum parsimony (MP)—should be chosen for any particular data set. ^ A number of factors, often hidden, influence the performance of a method. For example, it is generally understood that for a phylogenetically “difficult” data set more sophisticated methods (e.g., ML) tend to be more effective and thus should be chosen. However, it is the interplay of many factors that one needs to consider in order to avoid choosing an inferior method (potentially a costly mistake, both in terms of computational expenses and in terms of reconstruction accuracy.) ^ Chapter III of this dissertation details a phylogenetic reconstruction expert system that selects a superior proper method automatically. It uses a classifier (a Decision Tree-inducing algorithm) to map a new data set to the proper phylogenetic reconstruction method. ^

Structural conservations and diversities of dsRNA viruses: Structural studies of the cytoplasmic polyhedrosis virus by electron cryomicroscopy and 3D reconstruction

The Reoviridae virus family is a group of economically and pathologically important viruses that have either single-, double-, or triple-shelled protein layers enclosing a segmented double stranded RNA genome. Each virus particle in this family has its own viral RNA dependent RNA polymerase and the enzymatic activities necessary for the mature RNA synthesis. Based on the structure of the inner most cores of the viruses, the Reoviridae viruses can be divided into two major groups. One group of viruses has a smooth surfaced inner core, surrounded by complete outer shells of one or two protein layers. The other group has an inner core decorated with turrets on the five-fold vertices, and could either completely lack or have incomplete outer protein layers. The structural difference is one of the determinant factors for their biological differences during the infection. ^ Cytoplasmic polyhedrosis virus (CPV) is a single-shelled, turreted virus and the structurally simplest member in Reoviridae. It causes specific chronic infections in the insect gut epithelial cells. Due to its wide range of insect hosts, CPV has been engineered as a potential insecticide for use in fruit and vegetable farming. Its unique structural simplicity, unparalleled capsid stability and ease of purification make CPV an ideal model system for studying the structural basis of dsRNA virus assembly at the highest possible resolution by electron cryomicroscopy (cryoEM) and three-dimensional (3D) reconstruction. ^ In this thesis work, I determined the first 3D structure of CPV capsids using 100 kV cryoEM. At an effective resolution of 17 Å, the full capsid reveals a 600-Å diameter, T = 1 icosahedral shell decorated with A and B spikes at the 5-fold vertices. The internal space of the empty CPV is unoccupied except for 12 mushroom-shaped densities that are attributed to the transcriptional enzyme complexes. The inside of the full capsid is packed with icosahedrally-ordered viral genomic RNA. The interactions of viral RNA with the transcriptional enzyme complexes and other capsid proteins suggest a mechanism for RNA transcription and subsequent release. ^ Second, the interactions between the turret proteins (TPs) and the major capsid shell protein (CSPs) have been identified through 3D structural comparisons of the intact CPV capsids with the spikeless CPV capsids, which were generated by chemical treatments. The differential effects of these chemical treatment experiments also indicated that CPV has a significantly stronger structural integrity than other dsRNA viruses, such as the orthoreovirus subcores, which are normally enclosed within outer protein shells. ^ Finally, we have reconstructed the intact CPV to an unprecendented 8 Å resolution from several thousand of 400kV cryoEM images. The 8 Å structure reveals interactions among the 120 molecules of each of the capsid shell protein (CSP), the large protrusion protein (LPP), and 60 molecules of the turret protein (TP). A total of 1980 α-helices and 720 β-sheets have been identified in these capsid proteins. The CSP structure is largely conserved, with the majority of the secondary structures homologous to those observed in the x-ray structures of corresponding proteins of other reoviruses, such as orthoreovirus and bluetongue virus. The three domains of TP are well positioned to play multifunctional roles during viral transcription. The completely non-equivalent interactions between LPP and CSP and those between the anchoring domain of TP and CSP account for the unparalleled stability of this structurally simplest member of the Reoviridae. ^

Late Pliocene and Early Pleistocene dinoflagellate cysts and sea surface temperture reconstruction for several IODP and DSDP sites in the North Atlantic

In an attempt to document the palaeoecological affinities of individual extant and extinct dinoflagellate cysts, Late Pliocene and Early Pleistocene dinoflagellate cyst assemblages have been compared with geochemical data from the same samples. Mg/Ca ratios of Globigerina bulloides were measured to estimate the spring-summer sea-surface temperatures from four North Atlantic IODP/DSDP sites. Currently, our Pliocene-Pleistocene database contains 204 dinoflagellate cyst samples calibrated to geochemical data. This palaeo-database is compared with modern North Atlantic and global datasets. The focus lies in the quantitative relationship between Mg/Ca-based (i.e. spring-summer) sea-surface temperature (SSTMg/Ca) and dinoflagellate cyst distributions. In general, extant species are shown to have comparable spring-summer SST ranges in the past and today, demonstrating that our new approach is valid for inferring spring-summer SST ranges for extinct species. For example, Habibacysta tectata represents SSTMg/Ca values between 10° and 15°C when it exceeds 30% of the assemblage, and Invertocysta lacrymosa exceeds 15% when SSTMg/Ca values are between 18.6° and 23.5°C. However, comparing Pliocene and Pleistocene SSTMg/Ca values with present day summer values for the extant Impagidinium pallidum suggests a greater tolerance of higher temperatures in the past. This species occupies more than 5% of the assemblage at SSTMg/Ca values of 11.6-17.9°C in the Pliocene and Pleistocene, whereas present day summer SSTs are around -1.7 to 6.9°C. This observation questions the value of Impagidinium pallidum as reliable indicator of cold waters in older deposits, and may explain its bipolar distribution.

Stable isotope analysis and temperature reconstruction data from DSDP Hole 94-609B and ODP Hole 162-984B

The Greenland ice sheet is accepted as a key factor controlling the Quaternary glacial scenario. However, the origin and mechanisms of major Arctic glaciation starting at 3.15 Ma and culminating at 2.74 Ma are still controversial. For this phase of intense cooling Ravelo et al. proposed a complex gradual forcing mechanism. In contrast, our new submillennial-scale paleoceanographic records from the Pliocene North Atlantic suggest a far more precise timing and forcing for the initiation of northern hemisphere glaciation (NHG), since it was linked to a 2-3 °C surface water warming during warm stages from 2.95 to 2.82 Ma. These records support previous models, claiming that the final closure of the Panama Isthmus (3.0- ~2.5 Ma induced an increased poleward salt and heat transport. Associated strengthening of North Atlantic Thermohaline Circulation and in turn, an intensified moisture supply to northern high latitudes resulted in the build-up of NHG, finally culminating in the great, irreversible climate crash at marine isotope stage G6 (2.74 Ma). In summary, there was a two-step threshold mechanism that marked the onset of NHG with glacial-to-interglacial cycles quasi-persistent until today.

Geochemistry and sea-surface temperature reconstruction for sediment cores of the Iberian Margin

Biogenic records of the marine palaeoproductivity (carbonates, organic carbon, and C37 alkenones) and the molecular stratigraphy of past sea surface temperatures (SSTs; UK'37) were studied at high resolution in two cores of the Iberian Margin. The comparison of these records indicates that the oceanographic conditions switched abruptly during the past 160 kyr between three kinds of regimes. A first regime with high (17-22°C) SST and low productivity typifies the interglacial periods, marine isotopic stages (MIS) 5 and 1. Several periods during MIS 6, 2, and the terminations II and I are characterised by about 4-5°C colder SST and a higher organic matter accumulation, both of which define the second regime. This anticorrelation between SST and marine productivity suggests that these variations are related to the intensity of the coastal upwelling. By contrast with this upwelling behaviour, extremely low biological productivity and very cold SST (6-12°C) occurred during short phases of glacial MIS 6, 4, and 2, and as abrupt events (~1 kyr or less) during MIS 3. The three oceanographic regimes are consistent with micropalaeontological changes in the same cores based on foraminifera and diatoms. The general trend of these hydrologic changes follows the long-term glacial/interglacial cycle, but the millennium scale variability is clearly related to Heinrich events and Dansgaard-Oeschger cycles. Strengthening of the upwelling corresponds probably to an intensification of the subtropical atmospheric circulation over the North Atlantic which was influenced by the presence of continental ice sheets. However, extreme glacial conditions due to massive discharges of icebergs interrupted the upwelling. Interestingly, both terminations II and I coincided with strong but transient intensification of the upwelling.

Age determination and sea-surface reconstruction of sediment cores from the Norwegian Sea

Polycystine radiolarians are used to reconstruct summer sea surface temperatures (SSSTs) for the Late Pleistocene-Holocene (600-13,400 14C years BP) in the Norwegian Sea. At 13,200 14C years BP, the SSST was close to the average Holocene SSST (~12°C). It then gradually dropped to 7.1°C in the Younger Dryas. Near the Younger Dryas-Holocene transition (~10,000 14C years BP), the SSST increased 5°C in about 530 years. Four abrupt cooling events, with temperature drops of up to 2.1°C, are recognized during the Holocene: at 9340, 7100 ("8200 calendar years event"), 6400 and 1650 14C years BP. Radiolarian SSSTs and the isotopic signal from the GISP2 ice core are strongly coupled, stressing the importance of the Norwegian Sea as a mediator of heat/precipitation exchange between the North Atlantic, the atmosphere, and the Greenland ice sheet. Radiolarian and diatom-derived SSSTs display similarities, with the former not showing the recently reported Holocene cooling trend.

Long chain 1,13- and 1,15-diols in surface sediments ans SST reconstruction for sediment core GeoB6518-1

Although commonly reported in marine and freshwater environments, little is known about the biological sources of long chain alkyl 1,13- and 1,15-diols, and factors controlling their distributions. Here we analyzed the occurrence and distribution of these lipids in a comprehensive set of marine surface sediments and compare their distributions with environmental conditions like sea surface temperature (SST), salinity and nutrient concentrations. Fractional abundances of the C28 1,13-, C30 1,13- and C30 1,15-diols show a strong correlation with SST and based on these results, we propose the Long chain Diol Index (LDI), which expresses the C30 1,15-diol abundance relative to those of C28 1,13-, C30 1,13- and C30 1,15-diols. The LDI shows a strong linear correlation with SST (LDI = 0.033 × SST + 0.095; R2 = 0.969, n = 162) over a temperature range of -3 to 27 °C. Long chain diol distributions in sediments from the South Atlantic close to the Congo River outflow (West Africa) provided a 43 kyr LDI SST record. This record reflects several known climatic events and shows similarities with an alkenone-derived SST record obtained using the same suite of sediments, both in trend and in terms of absolute SST. This confirms the potential of the LDI as a proxy for palaeo-SST reconstruction.

Planktonic foraminifera and sea surface temperture reconstruction for Red Sea sediment cores

In order to assess how insolation-driven climate change superimposed on sea level rise and millennial events influenced the Red Sea during the Holocene, we present new paleoceanographic records from two sediment cores to develop a comprehensive reconstruction of Holocene circulation dynamics in the basin. We show that the recovery of the planktonic foraminiferal fauna after the Younger Dryas was completed earlier in the northern than in the central Red Sea, implying significant changes in the hydrological balance of the northern Red Sea region during the deglaciation. In the early part of the Holocene, the environment of the Red Sea closely followed the development of the Indian summer monsoon and was dominated by a circulation mode similar to the current summer circulation, with low productivity throughout the central and northern Red Sea. The climatic signal during the late Holocene is dominated by a faunal transient event centered around 2.4 ka BP. Its timing corresponds to that of North Atlantic Bond event 2 and to a widespread regionally recorded dry period. This faunal transient is characterized by a more productive foraminiferal fauna and can be explained by an intensification of the winter circulation mode and high evaporation. The modern distribution pattern of planktonic foraminifera, reflecting the prevailing circulation system, was established after 1.7 ka BP.