THE NEW HORIZONS MISSION HAS RETURNED A TREASURE TROVE OF HIGH-RESOLUTION IMAGE DATA FOR THE PLUTO-CHARON SYSTEM REVEALING COMPLEX WORLDS WITH BROAD TECTONIC SYSTEMS DIVERSE RESURFACING MECHANISMS AND SUGGESTIONS OF CRYOVOLCANIC ACTIVITY. WE AIM TO EXPLOIT THE NEW HORIZONS IMAGE DATASET AND PRODUCTS DERIVED FROM IT SUCH AS TOPOGRAPHY TO CONSTRAIN THE TECTONIC AND CRYOVOLCANIC HISTORY OF PLUTO AND CHARON AT REGIONAL TO GLOBAL SCALES THROUGH AN INTEGRATED PROGRAM OF DETAILED GEOLOGIC/STRUCTURAL MAPPING AND QUANTITATIVE ANALYSIS OF ICY SHELL LITHOSPHERE STRESS AND DEFORMATION STATES VIA FINITE ELEMENT METHOD MODELS. WE WILL USE THE SOPHISTICATED AND POWERFUL COMSOL MULTIPHYSICS FINITE ELEMENT CODE TO CREATE MODELS OF ICY SHELL LITHOSPHERES SUBJECT TO SEVERAL MODES OF LOADING. WE WILL CREATE GLOBAL-SCALE MODELS OF THE ICY SHELL LITHOSPHERES OF PLUTO AND CHARON TO EVALUATE THE EFFECTS OF PHENOMENA SUCH AS PLANETARY DE- OR RE-SPINNING TRUE POLAR WANDER AND THERMAL EXPANSION AND RELATED EFFECTS FROM A GROWING LITHOSPHERE OR FREEZING OCEAN. SUCH MODELS EXHIBIT SPHERICAL GEOMETRY USING APPROPRIATE SYMMETRY BOUNDARY CONDITIONS TO REDUCE MODEL SIZE WHERE ALLOWABLE. THESE MODELS WILL INCORPORATE EFFECTS THAT ARE MISSING FROM CONVENTIONAL ANALYTIC CALCULATIONS SUCH AS VARIATIONS IN SHELL THICKNESS (INCLUDING TOPOGRAPHIC BASINS) TIME-DEPENDENT BEHAVIOR (INCLUDING STRESS RELIEF FROM BRITTLE YIELDING AND VISCOELASTIC RELAXATION AND GROWING BASIN INFILL LOADS) AND BIASES TOWARD EXTENSIONAL STRESS STATES. MODELS WILL BE EVALUATED ON THEIR CONSISTENCY WITH THE OBSERVED STATES OF STRESS AND DEFORMATION ON PLUTO AND CHARON AS INDICATED BY IMAGE-BASED MAPPING OF TECTONIC FEATURES WITH THE GOAL OF CONSTRAINING THE CONDITIONS THAT CONTROLLED THE EVOLUTION OF THOSE BODIES. THESE MIGHT INCLUDE SHELL AND OCEAN THICKNESSES SPIN RATE CHANGES POLAR WANDER MAGNITUDES AND RATES AND BASIN FILL AND LOADING RATES. WE WILL ALSO USE SMALLER-SCALE MODELS TO ADDRESS THE FORMATION MECHANISMS OF PARTICULAR REGIONS AND FEATURES ON PLUTO AND CHARON. THESE MODELS MAY HAVE PLANAR GEOMETRY (FOR SMALL FEATURES) OR CURVED PARTIAL SPHERE GEOMETRIES (FOR REGIONAL-SCALE FEATURES). FOR EXAMPLE A RADIATING FRACTURE SYSTEM NORTH OF TARTARUS DORSA BEARS RESEMBLANCE TO RADIATING FRACTURE SYSTEMS ON VENUS THAT ARE RELATED TO INTRUSIVE MAGMATIC ACTIVITY IN THE PRESENCE OF A REGIONAL STRESS FIELD. MODELS CAN CONSTRAIN THE VALUES OF REGIONAL STRESS COMPONENTS AND PRESSURIZATIONS OF POTENTIAL CRYOMAGMA RESERVOIRS THAT PRODUCE STRESS STATES CONSISTENT WITH THE OBSERVED FRACTURE PATTERNS. THE STRESS STATES CALCULATED FOR MODELS OF THE VARIOUS SCALES STATED ABOVE HAVE IMPLICATIONS FOR THE TRANSPORT OF CRYOMAGMAS THROUGH THE LITHOSPHERIC SHELLS OF PLUTO AND CHARON. WE WILL USE PREVIOUSLY DEVELOPED EXPRESSIONS FOR FORCE BALANCE ON DIKE WALLS TO CALCULATE ASCENT-ENHANCING OR -INHIBITING EFFECTS OF LITHOSPHERIC STRESS STATES FOR CRYOMAGMAS ALLOWING US TO IDENTIFY LOCATIONS OR REGIONS WHERE ERUPTION OF CRYOLAVAS IS FACILITATED. SUCH CALCULATIONS IN CONJUNCTION WITH MAPPING OF POTENTIAL CRYOVOLCANIC FEATURES SUCH AS WRIGHT AND PICCARD MONTES MAY ALLOW INSIGHTS INTO ERUPTIONS IN FLOOD (WIDESPREAD) OR EDIFICE-BUILDING (LOCALIZED) MODES. THIS PROPOSAL ADDRESSES OBJECTIVES OF NFDAP TO ENHANCE THE SCIENTIFIC RETURN FROM NEW FRONTIERS MISSIONS BY BROADENING SCIENTIFIC PARTICIPATION IN THE ANALYSIS AND INTERPRETATION OF DATA RETURNED BY THESE MISSIONS. THE PROPOSAL TEAM BLENDS EXPERIENCE AND FACILITY WITH NEW HORIZONS DATA WITH EXPERTISE IN STRUCTURAL MODELING OF LARGE PLANETARY VOLCANO-TECTONIC SYSTEMS AND ANALYSES OF STRESS MODULATED MAGMA ASCENT. WE PLAN TO EXPLOIT DECADES OF EXPERIENCE IN MODELING AND INTERPRETATION OF TECTONIC AND VOLCANIC FEATURES AND LITHOSPHERIC DYNAMICS OF PLANETS SPANNING THE DISTANCE FROM VENUS TO RHEA TO ALLOW US TO EXTRACT NEW INSIGHTS INTO THE EVOLUTION OF PLUTO AND CHARON.