The Maldives archipelago acts for over 25 myrs as a giant natural sediment trap in the eastern Arabian Sea. Drifts and periplatform deposits bear the record of environmental changes such as sea-level fluctuations but also of monsoon-driven changes of the surface and intermediate water mass current regime, and of wind-driven dust influx. Carbonate drifts in the Inner Sea indicate the establishment of a strong wind-driven current regime in the Maldives at 12.9 – 13 Ma. Ten unconformities, dissecting the Miocene to Recent drift sequences, attest to changes in current strength or direction. A major shift in the drift packages is dated at 3.8 Ma that coincides with the end of stepwise platform drowning and a reduction of the OMZ in the Inner Sea. The lithogenic fraction of the Maldives carbonate drifts provides a unique record of atmospheric dust transport during the past 4 myrs as grain size provides proxies for dust flux as well as wind transport capacity. Entrainment and long-range transport of dust in the medium to coarse silt size range is linked to the strength of the Arabian Shamal winds and the occurrence of convective storms which prolong dust transport. Dust flux and the size of dust particles increased between 4.0 and 3.3 Ma, corresponding to the closure of the Indonesian seaway and the intensification of the South Asian Monsoon. Between 1.6 Ma and the Recent, dust flux again increased and shows higher variability, especially during the last 500 kyr. Transport capacity increased between 1.2 and 0.5 Ma but slightly decreased since then. Dust transport varies on orbital timescales, with eccentricity control being the most prominent (400 kyr throughout the record, 100 kyr between 2.0 and 1.3 Ma, and since 1.0 Ma). Higher frequency cycles (obliquity and precession) are most pronounced in wind transport capacity. The published and ongoing studies of IODP Expedition 359 cores show that deposits surrounding carbonate platforms, i.e. carbonate drifts, bear a previously underestimated potential to add substantial knowledge for the understanding of the monsoon evolution on million-year, but also on shorter time scales. Potential targets for further research and drilling are for example the Laccadives, the Mascarene Plateau or the South China Sea platforms.

The earliest human migration from Sunda (South-East Asian archipelago) to Sahul (Australia and New Guinea) is still heavily debated with proposed timings between c.65-45kaBP depending on the evidence base and interpretation of the data. As part of the EU funded ACROSS project, focused on the mode and route of early migration in to SAHUL, we are undertaking an integrated interpretative study of the evolving submerged landscapes for the Late Pleistocene of the NW Australian Shelf. Oil and gas industry 3D and 2D seismic data, with some core/borehole data, are being used to determine lowstand palaeo-environments and shoreline positions. This information is informing modelling of ocean tide and current patterns that may have been influenced. The seismic is being interpreted supplemented by using time-slices on relative impedance inverted post-stack data. Layer stripping, seismic geomorphology, sequence boundary and depth analysis are being applied to datasets in the Bonaparte Basin, Kimberley Shelf and Arafura regions of Australia’s North-West Shelf area. Interpretation of the seismic data is constrained by dated stratigraphy in shallow cores with lower bounds determined from oil/gas well bores. MIS stages 1-4 are identified, however, the seismic response is a composite of time periods due to varying sedimentation rates, non-depositional hiatuses and minimal vertical seismic travel time covering this interval which limits the analysis to the top 50ms TWT (c. 40-45 m) of events below the seabed. This paper reviews the workflows that have been developed to maximise the fine scale detail that can be recovered for a range of terrestrial and marine environments. Procedures include inverse-Q, impedance inversion, spectral decomposition and time-slicing relative to seabed. High resolution 2D seismic data is also being used to augment and inform the interpretation of the conventional oil/gas 3D seismic data. Data examples will be presented showing the geomorphological characteristics (river channels, avulsions, levees, drainage channels, dunes and near shore carbonate reefs) of the lowstand and transgressive landscapes during this period. The palaeo-reconstructions are now being developed from the interpreted seismic geomorphology for the specific consideration of human seaborne travel.

Tropical modes of variability, including the Madden‐Julian Oscillation (MJO) and the El Niño‐Southern Oscillation (ENSO), are challenging to represent in climate models. Previous studies suggest their fundamental dependence on zonal asymmetry, but such dependence is rarely addressed with fully coupled ocean dynamics. This study fills the gap by using fully coupled, idealized Community Earth System Model (CESM) and comparing two nominally ocean-covered configurations with and without a meridional boundary. For the MJO-like intraseasonal mode, its separation from equatorial Kelvin waves and the eastward propagation of its convective and dynamic signals depend on the zonal gradient of the mean state. For the ENSO-like interannual mode, in the absence of the ocean’s meridional boundary, a circum-equatorial dominant mode emerges with distinct ocean dynamics. The interpretation of the dependence of these modes on zonal asymmetry is relevant to their representation in realistic climate models.

In Feb. 2017, a five-day sequence of atmospheric river storms in California, USA, resulted in extreme inflows to Lake Oroville, the state’s second-largest reservoir. Damage to the reservoir’s spillway infrastructure necessitated evacuation of 188,000 people; subsequent infrastructure repairs cost $1 billion. We assess the atmospheric conditions, snowmelt, and runoff against major historical events. The event generated exceptional runoff volumes (second-largest in a 30 year record) partially at odds with the event precipitation totals (ninth-largest). We explain the discrepancy with observed record melt of deep antecedent snowpack, heavy rainfall extending to unusually high elevations, and high water vapor transport during the atmospheric river storms. An analysis of distributed snow water equivalent indicates that snowmelt increased water available for runoff watershed-wide by 37% (25-52% at 90% confidence). The results highlight an acute flood risk to public safety and infrastructure projected to increase in severity in a warmer and more variable climate.

Shallow cloud fields over the subtropical ocean exhibit many spatial patterns. The frequency of occurrence of these patterns can change under global warming. Hence, they may influence subtropical marine clouds’ climate feedback. While numerous metrics have been proposed to quantify cloud patterns, a systematic, widely accepted description is still missing. Therefore, this paper suggests one. We compute 21 metrics for 5000 satellite scenes of shallow clouds over the subtropical Atlantic Ocean and translate the resulting dataset to its principal components (PCs). This yields a unimodal, continuous distribution without distinct classes, whose first four PCs explain 82% of all 21 metrics’ variance. The PCs correspond to four interpretable dimensions: Characteristic length, void size, directional alignment and horizontal cloud-top height variance. These dimensions span a space in which an effective pattern description can be given, which may be used to better understand the patterns’ underlying physics and feedback on climate.

Mosquitoes in recent years have increased greatly in numbers due to the rapidly changing climate and rising temperatures. With this change comes suitable habitats for mosquitoes which are the most efficient killers in all of the animal kingdom due to the number of death from mosquito-borne diseases. If we were able to pinpoint the certain areas that mosquitoes are most attracted to we could in theory slow or even prevent the spread of mosquitoes. In our project, the research was conducted to find a correlation between the color and size of the traps to the amount of mosquitos that are present. With our findings, we were able to conclude that the bigger the traps, the faster the mosquitoes would be attracted to that area. We also found that the different traps would hold similar densities of mosquito larvae per square inch.

Globally researchers have unraveled unique locations that helped to understand the chronology of the critical events concerning the Earth’s past. Among such geological events, the time-shot of the Permian-Triassic (P-Tr) extinction event is one of the significant revelations concerning the end and start of life on Earth. Among various geological sites in the world that contain the critical information regarding the P-Tr extinction event, Guryul ravine in Kashmir India is geologically a treasure. It bears specimens of primordial corals, small invertebrates, plants, and a group of mammal resembling reptiles, called therapsids. Due to its immense importance, the Government of India had decided to accredit the site of Guryul ravine section as an international fossil park. However, due to political turmoil in the region and unabated mining and industrial activities within the vicinity of it, has threatened the very existence of this scientific wealth. This paper reviews the importance of the Guryul Ravine Section, paleoclimatic conditions of that time, and the current threats it is facing to stimulate the stakeholders for the conservation of this site in the global scientific interest.

To clarify the impact of basal melting of the Antarctic ice sheet and biological productivity on biogeochemical processes in Antarctic coastal waters, concentrations of dissolved inorganic carbon (DIC), total alkalinity (TA), inorganic nutrients, chlorophyll a, and stable oxygen isotopic ratios (δ18O) were measured from the offshore slope to the ice front of the Totten Ice Shelf (TIS) during the spring/summer of 2018, 2019, and 2020. Off the TIS, modified Circumpolar Deep Water (mCDW) intruded onto the continental shelf and flowed along bathymetric troughs into the TIS cavity, where it met the ice shelf base and formed a buoyant mixture with glacial meltwater. Physical oceanographic processes mostly determined the distributions of DIC, TA, and nutrient concentrations. However, DIC, TA, and nutrient concentrations on the surface of the ice front were decreased by photosynthesis and the dilution effect of meltwater from sea ice and the base of the ice shelf. The partial pressure of CO2 (pCO2) in surface water was reduced by photosynthesis and dilution, and the surface water became a strong CO2 sink for the atmosphere. The DIC and TA (normalized to salinity of 34.3 to correct for dilution effects) changed in a molar ratio of 106:16 because of phytoplankton photosynthesis. The decrease of pCO2 by more than 100 μatm with respect to mCDW was thus the result of photosynthesis. The nutrient consumption ratio suggested that enough iron was present in the water column to supply the surface layer via buoyancy-driven upwelling and basal melting of the TIS.

The Delaware River is a major freshwater supplier of New York City (NYC). Nearly half of NYC drinking water is supplied by inter-basin transfer of surface water stored in reservoirs within the upper reaches of the Delaware River. In its lower reaches, the Delaware River is a tidal estuary, and upstream freshwater discharge provides a critical control on estuary salinity. During the record 1950–1960’s drought, NYC water withdrawals exacerbated low flows. Estuary salinity reached levels that threatened freshwater intakes and groundwater recharge, resulting in legal action and Supreme Court decrees. We revisit this classic case study in coupled human and natural systems using the Energy Exascale Earth System Model (E3SM). The E3SM water management sub-model is updated to include inter-basin water transfer and reservoir-specific operating rules. Model simulations are developed to investigate competition between NYC water demand and in-stream flow targets needed to maintain estuary salinity within regulatory guidelines under historic and future climate. To our knowledge, this is a first demonstration of an Earth System Model simulation with inter-basin water transfer, which, in this study area, provides water for nearly five million people living outside the Delaware River basin in New York City and New Jersey.

Particulate inorganic carbon (PIC) indicates CaCO3, which is produced biologically and constitutes about a quarter of total marine sediments. It is produced by coral reefs, foraminifera, coccolithophores and some macrophytic algae (Balch et al., 2005). The dominating calcifiers in the oceans are the coccolithophores, producing roughly half of the global ocean carbonates. They also cover 20 % of the blooming phytoplankton (Krumhardt et al., 2017). An algorithm has been developed to estimate the PIC concentrations from satellites. The data is collected by MODIS Aqua and SeaWiFS. Observational studies show mixed results for coccolithophore abundance as both increase/decrease in various places. However, it can be agreed that anthropogenic climate change has a significant impact on the PIC, mainly due to global warming and ocean acidification.

Permafrost active layer thickness (ALT) is a sensitive indicator of permafrost response to climate change. In recent decades, ALT has increased at sites across the Arctic, concurrent with observed increases in annual minimum streamflow (baseflow). The trends in ALT and baseflow are thought to be linked via: 1) increased soil water storage capacity due to an increased active layer, and 2) enhanced soil water mobility within a more continuous active layer, both of which support higher baseflow in Arctic rivers. One approach to analyzing these changes in ALT and baseflow is to use baseflow recession analysis, which is a classical method in hydrology that relates groundwater storage S to baseflow Q with a power law-like relationship Q = aSb. For the special case of a linear reservoir (b=1.0), the baseflow recession method has been extended to quantify changes in ALT from streamflow measurements alone. We test this approach at sites across the North American Arctic and find that catchments underlain by permafrost behave as nonlinear reservoirs, with scaling exponents b~1.5–3.0, undermining the key assumption of linearity that is commonly applied in this method. Despite this limitation, trends in a provide insight into the relationship between changing ALT and changing Arctic baseflow. Although care should be taken to ensure the theoretical assumptions are met, baseflow recession analysis shows promise as an empirical approach to constrain modeled permafrost change at the river basin scale.

The Canadian Arctic Archipelago (CAA) is vulnerable to climate warming, and with over 300 tidewater glaciers, is a hotspot for enhanced glacial retreat and meltwater runoff to the ocean. In contrast to Greenlandic and Antarctic systems, CAA glaciers and their impact on the marine environment remain largely unexplored. Here we investigate how CAA glaciers impact nutrient delivery to surface waters. We compare water column properties in the nearshore coastal zone along a continuum of locations, spanning those with glaciers (glacierized) to those without (non-glacierized), in Jones Sound, eastern CAA. We find that surface waters of glacierized regions contain significantly more macronutrients (nitrogen, silica, phosphorus) and micronutrients (iron, manganese) than their non-glacierized counterparts. Water column structure and chemical composition suggest that macronutrient enrichments are a result of upwelling induced by rising submarine discharge plumes, while micronutrient enrichments are driven directly by glacial discharge. Generally, the strength of upwelling and associated macronutrient delivery scales with tidewater discharge volume. Glacier-driven delivery of the limiting macronutrient, nitrate, is of particular importance for local productivity, while metal delivery may have consequences for regional micronutrient cycling given Jones Sound’s important role in modifying water masses flowing into the North Atlantic. Finally, we use the natural variability in glacier characteristics observed in Jones Sound to consider how nutrient delivery may be affected as glaciers retreat. The impacts of melting glaciers on marine ecosystems through both these mechanisms will likely be amplified with increased meltwater fluxes in the short-term, but eventually muted as CAA ice masses diminish.

Mixed-phase clouds (MPCs), which consist of both supercooled cloud droplets and ice crystals, play an important role in the Earth’s radiative energy budget and hydrological cycle. In particular, the fraction of ice crystals in MPCs determines their radiative effects, precipitation formation and lifetime. In order for ice crystals to form in MPCs, ice-nucleating particles (INPs) are required. However, a large-scale relationship between INPs and ice initiation in clouds has yet to be observed. By analyzing satellite observations of the typical transition temperature (T*) where MPCs become more frequent than liquid clouds, we constrain the importance of INPs in MPC formation. We find that over the Arctic and Southern Ocean, snow and sea ice cover significantly reduces T*. This indicates that the availability of INPs is essential in controlling cloud phase evolution and that local sources of INPs in the high-latitudes play a key role in the formation of MPCs.

We develop a hierarchy of simplified ocean models for coupled ocean, atmosphere, and sea ice climate simulations using the Community Earth System Model version 1 (CESM1). The hierarchy has four members: a slab ocean model, a mixed-layer model with entrainment and detrainment, an Ekman mixed-layer model, and an ocean general circulation model (OGCM). Flux corrections of heat and salt are applied to the simplified models ensuring that all hierarchy members have the same climatology. We diagnose the needed flux corrections from auxiliary simulations in which we restore the temperature and salinity to the daily climatology obtained from a target CESM1 simulation. The resulting 3-dimensional corrections contain the interannual variability fluxes that maintain the correct vertical gradients of temperature and salinity in the tropics. We find that the inclusion of mixed-layer entrainment and Ekman flow produces sea surface temperature and surface air temperature fields whose means and variances are progressively more similar to those produced by the target CESM1 simulation. We illustrate the application of the hierarchy to the problem of understanding the response of the climate system to the loss of Arctic sea ice. We find that the shifts in the positions of the mid-latitude westerly jet and of the Inter-tropical Convergence Zone (ITCZ) in response to sea-ice loss depend critically on upper ocean processes. Specifically, heat uptake associated with the mixed-layer entrainment influences the shift in the westerly jet and ITCZ. Moreover, the shift of ITCZ is sensitive to the form of Ekman flow parameterization.

Climate data from Earth System Models are increasingly being used to study the impacts of climate change on a broad range of biogeophysical (forest fires, fisheries, etc.) and human systems (reservoir operations, urban heat waves, etc.). Before this data can be used to study many of these systems, post-processing steps commonly referred to as bias correction and statistical downscaling must be performed. “Bias correction” is used to correct persistent biases in climate model output and “statistical downscaling” is used to increase the spatiotemporal resolution of the model output (i.e. 1 deg to 1/16th deg grid boxes). For our purposes, we’ll refer to both parts as “downscaling”. In the past few decades, the applications community has developed a plethora of downscaling methods. Many of these methods are ad-hoc collections of post processing routines while others target very specific applications. The proliferation of downscaling methods has left the climate applications community with an overwhelming body of research to sort through without much in the form of synthesis guiding method selection or applicability. Motivated by the pressing socio-environmental challenges of climate change – and with the learnings from previous downscaling efforts in mind – we have begun working on a community-centered open framework for climate downscaling: scikit-downscale. We believe that the community will benefit from the presence of a well-designed open source downscaling toolbox with standard interfaces alongside a repository of benchmark data to test and evaluate new and existing downscaling methods. In this notebook, we provide an overview of the scikit-downscale project, detailing how it can be used to downscale a range of surface climate variables such as air temperature and precipitation. We also highlight how scikit-downscale framework is being used to compare existing methods and how it can be extended to support the development of new downscaling methods.

Solar-Induced chlorophyll Fluorescence (SIF) provides a powerful proxy for determining forest gross primary production (GPP), particularly in evergreen ecosystems where traditional measures of greenness fail. The dynamics of the SIF/GPP relationship, however, are poorly understood under varying viewing directions and light conditions. This is, in large part, due to challenges in measuring SIF at the spatiotemporal scale that is necessary to understand these effects. Therefore, the aim of this work is to utilize high-temporal and spatial resolution SIF measurements to better constrain the response of SIF to ambient canopy illumination and viewing geometry. We use a PhotoSpec instrument and eddy covariance measurements to explore the SIF/GPP relationship under various viewing directions and light conditions during the 2019 and 2020 growing seasons at the Old Black Spruce site in Saskatchewan, Canada. PhotoSpec is a tower-based 2-D scanning spectrometer system capable of taking Fraunhofer-line based SIF retrievals in the red and far-red wavelength ranges with a 0.7 degree field of view at a ~30 second time resolution. Measured SIF and GPP are combined with SCOPE modelling results to provide a mechanistic understanding of the physical and ecophysiological drivers for the SIF/GPP relationship in the Boreal Forest. Our results show that viewing direction and solar zenith/azimuth angles are important for the SIF signal under direct light conditions, but not under diffuse. Furthermore, the SIF/GPP relationship changes under direct and diffuse light conditions at a 30 minute, daily, and monthly resolution. Our ability to use SIF as a proxy for GPP depends on a quantitative understanding of radiative transfer within the canopy and how scanning geometry impacts SIF measurements. These results provide an important insight into these relationships in the Boreal forest, a region where GPP has been traditionally difficult to track using remote sensing.

The increasing role of women in leadership roles may be both a cause and effect of public attitudes. Recent politicization of science amidst increasing polarization of American politics juxtaposed with examples of female leadership (Fleadership) throughout the U.S (and abroad) begs the question of how gender impacts crisis response decision-making. Here, we investigate at the state level, female gubernatorial leadership, recent presidential (and VP) voting patterns, climate denial, and COVID/100k. These reveal correlations between climate denial, high COVID/100k, and presidential votes for exclusively male candidates. Americans who embraced federal Fleadership in 2016 and again in 2020 were less likely to deny climate change and spread COVID-19. The results suggest that rather than Fleadership inspiring responsible COVID preventative behavior at state or national levels, or leading to better understanding of climate change, the populations who tend to elect females, also understand science more than those who do not.

The Taurus Mountains form the southern margin of the Central Anatolian Plateau of Turkey and form an orographic barrier separating the cold, semi-arid interior to the north from the mild Mediterranean coast to the south. When and how they formed, and the extent which they have influenced the regional climate remains poorly constrained. The Attepe iron deposits sit on the northern part of the Eastern Taurus mountains at altitude of 1.5-2 km and consequently are ideally located to record interactions between climate and tectonics. (U-Th)/He ages of iron-oxide-oxyhydroxides from four mines within the Attepe iron deposits record ages of 1-5 Ma consistent with the persistence of hot humid climate conditions throughout the Pliocene and Pleistocene. In mines where samples are measured from different depths the age data are consistent with water table lowering rate of between 12.3 to 6.4 m/Myr. Translating these to rock uplift rates they are close to uplift/incision recorded within the Central Anatolian Plateau over the past 2 Ma, suggesting that the region was already at or close to its current elevation by the late Miocene. The latest goethite precipitation constrains the cessation of hot-humid climate to sometime in the last million years and implies that regional climate cooling, rather than surface uplift, was the main driver of aridification.

As the largest hydroelectric projects worldwide, the Three Gorges Dam (TGD) affects local precipitation because of the changes of hydrological cycle caused by the impounding and draining of the TGD. However, the influencing characteristics of the TGD on local precipitation remain elusive. In this study, we used precipitation anomaly data derived from long time-series grid precipitation datasets between 1988 and 2017 to understand the changes of precipitation caused by the TGD between 2 epochs, before and after the construction of the TGD (i.e., 1988–2002 and 2003–2017), in the Three Gorges Reservoir Area (TGRA). Results showed that the annual and dry season precipitation anomaly in the TGRA showed an increasing trend, and the flood season precipitation anomaly showed a slight decrease. After the impoundment of the TGD, the precipitation concentration degree in the TGRA was decreased, indicating that the precipitation became increasingly uniform, and the precipitation concentration period was insignificantly increased. An obvious resonance phenomenon between the monthly average water level and precipitation anomaly occurred in the TGRA after 2011 and showed a positive correlation. Our findings excavated the change of local precipitation characteristics before and after the impoundment of the TGRA and proved that this change had a close relationship with the water level.

The diatom oxygen isotope composition (δ18Odiatom) from lacustrine sediments helps tracing the hydrological and climate dynamics in individual lake catchments, and is generally linked to changes in temperature and δ18Olake. Lake Bolshoye Shchuchye (67°53’N; 66°19’ E; 186 m a.s.l) is the largest and deepest freshwater reservoir in the Polar Urals, Arctic Russia. Its δ18Odiatom record generally follows a decrease in summer insolation and the northern hemisphere (NH) temperature history. However, it displays exceptional, short-term variations exceeding 5‰, especially in Mid and Late Holocene. This centennial-scale variability occurs contemporaneously with and similarly to Holocene NH glacier advances. However, larger Holocene glacier advances in the Lake Bolshoye Shchuchye catchment are unknown and have not left any significant imprint on the lake sediment record. As Lake Bolshoye Shchuchye is deep and voluminous, about 30−50% of its volume needs to be exchanged with isotopically different water within decades to account for these shifts in the δ18Odiatom record. A plausible source of water with light isotope composition inflow is snow, known to be transported in surplus by snow redistribution from the windward to the leeward side of the Polar Urals. Here, we propose snow melt and influx changes being the dominant mechanism responsible for the observed short-term changes in the δ18Odiatom record. This is the first time such drastic, centennial-scale hydrological changes in a catchment have been identified in Holocene lacustrine diatom oxygen isotopes, which, for Lake Bolshoye Shchuchye, are interpreted as proxy for summer temperatures and palaeo precipitation.

In operational analyses of the surface moisture imbalance that defines drought, the supply aspect has generally been well characterized by precipitation; however, the same count be said of the demand side—a function of evaporative demand (E0) and surface moisture availability. In drought monitoring, E0 is often poorly parameterized by a climatological mean, by non-physically based estimates, or is neglected entirely. One problem has been a paucity of driver data—on temperature, humidity, solar radiation, and wind speed—required to fully characterize E0. This deficient E0 modeling is particularly troublesome over data-sparse regions that are also home to drought-vulnerable populations, such as across much of Africa. There is thus urgent need for global E0 estimates for physically accurate drought analyses and food security assessments; further we need an improved understanding of how E0 and drought interact and to exploit these interactions in drought monitoring. In this presentation we explore ways to meet these needs. From MERRA-2—an accurate, fine-resolution land-surface/atmosphere reanalysis—we have developed a >38-year, daily, global Penman-Monteith reference ET dataset as a fully physical metric of E0. This dataset is valuable for examining hydroclimatic changes and extremes. A novel drought index based on this dataset—the Evaporative Demand Drought Index (EDDI)—represents drought’s demand perspective, and permits early warning and ongoing monitoring of agricultural flash drought and hydrologic drought. We highlight the findings of our examination of E0-drought interactions and using EDDI in Africa. Using reference ET as an E0 metric has permitted explicit attribution of the variability of E0 across Africa, and of E0 anomalies associated with canonical droughts in the Sahel region. This analysis determines where, when, and to what relative degree each of the individual drivers of E0 affects the demand side of drought. Using independent estimates of drought across space and time—CHIRPS precipitation and the Normalized Difference Vegetation Index for 1982-2015—we examine the differences between drought and non-drought periods, and between precipitation-forced droughts and droughts forced by a combination of precipitation and E0.

The role of individual and collective human action is increasingly recognized as a prominent and arguably paramount determinant in shaping the behavior, trajectory, and vulnerability of multisector systems. This human influence operates at multiple scales: from short-term (hourly to daily) to long-term (annually to centennial) timescales, and from the local to the global, pushing systems towards either desirable or undesirable outcomes. However, the effort to represent human systems in multisector models has been fragmented across philosophical, methodological, and disciplinary lines. To cohere insights across diverse modeling approaches, we present a new typology for classifying how human actors are represented in the broad suite of coupled human-natural system models that are applied in MultiSector Dynamics (MSD) research. The typology conceptualizes a “sector” as a system-of-systems that includes a diverse group of human actors, defined across individual to collective social levels, involved in governing, provisioning, and utilizing products, goods, or services towards some human end. We trace the salient features of modeled representations of human systems by organizing the typology around three key questions: 1) Who are the actors in MSD systems? 2) What are their actions? 3) How and for what purpose are these actors and actions operationalized in a computational model? We use this typology to critically examine existing models and chart the frontier of human systems modeling for MSD research.

In recent years, many researchers and advocates have noted the potential of religious groups and institutions to leverage their significant influence in favor of addressing environmental challenges. However, in the United States, many scientists struggle to communicate the implications of their work on climate change with faith communities who may be skeptical of both climate science and scientists. Recent polls from the Pew Research Center show that white evangelical Protestants are the least likely to believe climate change is caused by human activity and the most likely to assert that there is no solid scientific support for a changing climate. However, the full picture is more nuanced than can be captured in a news headline or polling survey, and evangelical Christiantiy is a diverse movement that is also found at the forefront of enviornmental and climate science and action. Drawing on more than six years of experience working on climate science communication and climate action solutions among fellow evangelicals in the United States, this presentation highlights best practices for communicating climate science to faith communities. Showcasing examples of work advanced through the Evangelical Environmental Network, Young Evangelicals for Climate Action, and PBS Global Weirding series with Dr. Katharine Hayhoe, I present a hopeful view of efforts to communicate climate change in a way that intentionally and genuinely connects with people's values, and ultimately motivates action. Additionally, this presentation discusses the challenges of and opportunities for engaging communities of faith as scientists with a different or no faith affiliation.

Tropical peatlands are among the most carbon-dense ecosystems on Earth, and their water storage dynamics strongly control these carbon stocks. The hydrological functioning of tropical peatlands differs from that of northern peatlands, which has not yet been accounted for in global land surface models (LSMs). Here, we integrated tropical peat-specific hydrology modules into a global LSM for the first time, by utilizing the peatland-specific model structure adaptation (PEATCLSM) of the NASA Catchment Land Surface Model (CLSM). We developed literature-based parameter sets for natural (PEATCLSMTrop,Nat) and drained (PEATCLSMTrop,Drain) tropical peatlands. The operational CLSM version (which includes peat as a soil class) and PEATCLSMTrop,Nat were forced with global meteorological input data and evaluated over the major tropical peatland regions in Central and South America, the Congo Basin, and Southeast Asia. Evaluation against a unique and extensive data set of in situ water level and eddy covariance-derived evapotranspiration showed an overall improvement in bias and correlation over all three study regions. Over Southeast Asia, an additional simulation with PEATCLSMTrop,Drain was run to address the large fraction of drained tropical peatlands in this region. PEATCLSMTrop,Drain outperformed both CLSM and PEATCLSMTrop,Nat over drained sites. Despite the overall improvements of both tropical PEATCLSM modules, there are strong differences in performance between the three study regions. We attribute these performance differences to regional differences in accuracy of meteorological forcing data, and differences in peatland hydrologic response that are not yet captured by our model.