COMET 103P/HARTLEY 2 IS A FASCINATING OBJECT. IT IS THE SMALLEST OF THE COMETS VISITED BY SPACECRAFT BUT HAD VERY WELL DEVELOPED AND PRODUCTIVE COMA. IT IS ONE OF THE SO-CALLED "HYPERACTIVE" COMETS THAT INCLUDE SHORT-PERIOD COMETS 46P/WIRTANEN AND 73P/SCHWASSMAN- WACHMANN 3 AS WELL AS SOME LONG-PERIOD (OORT CLOUD) COMETS. THE GAS PRODUCTION RATES OF THESE COMETS EXCEED THE RATE THAT WOULD BE EXPECTED IF 100% OF THE SURFACE WERE COVERED IN PURE WATER ICE. HARTLEY 2 WAS VISITED BY THE DEEP IMPACT SPACECRAFT AS PART OF THE EXTENDED EPOXI MISSION ON 4 NOVEMBER 2010. THIS WAS A PARTICULARLY FAVORABLE APPARITION OF THIS COMET WHICH PASSED WITHIN 0.15 AU FROM THE EARTH AND WAS WIDELY OBSERVED FROM GROUND BASED AND SPACE BASED FACILITIES. ONE OF THE SURPRISING ASPECTS OF THE EPOXI MISSION IN COMBINATION WITH REMOTE OBSERVATIONS WAS THAT ITS HYPERACTIVITY WAS DUE TO DIRECT SUBLIMATION OF MAJOR REGIONS OF THE NUCLEUS HAVING BEEN DOMINATED BY CO2 WHICH LIBERATED MANY WATER ICE GRAINS AND CHUNKS THAT SUBSEQUENTLY SUBLIMATED IN THE COMA. A STUDY BY OUR GROUP (FOUGERE ET AL. 2013 ICARUS 225 688) SUGGESTED THAT SIMILAR AMOUNTS OF H2O AND CO2 WERE SUBLIMATING FROM THE NUCLEUS DIRECTLY BUT MOST OF THE WATER PRODUCTION RESULTED FROM THE SUBLIMATION OF THE EXTENDED ICY GRAIN/CHUNK HALO RELEASED BY CO2 ACTIVITY. THIS PICTURE WAS QUANTITATIVELY SPECIFIED ONLY FOR THE SNAPSHOT THAT REPRESENTED THE TIME PERIOD OF THE FLYBY AND USED A 2D-AXISYMMETRIC APPROXIMATION FOR THE NUCLEUS SHAPE AND ACTIVITY. WE PROPOSE HERE TO EXTEND OUR ANALYSIS TO THE LARGER SET OF EPOXI SPECTRAL IMAGES TAKEN DURING THE ENTIRE PERIOD AROUND CLOSE APPROACH USING OUR 3D ANALYSIS AND MODELING TOOLS DEVELOPED SINCE 2013 IN CONNECTION WITH THE ROSETTA MISSION (FOUGERE ET AL. 2016 A&A 588 A134) AND CASSINI DATA ANALYSIS (TENISHEV ET AL. 2010 JGR 115 A092302; AND 2014 JGR 119 2658). EPOXI SPECTRA (A'HEARN ET AL. 2011 SCIENCE 332 1396) TAKEN OVER SEVERAL DAYS INDICATED A VARIATION IN THE RELATIVE PRODUCTION OF H2O AND CO2 AS THE BI-LOBED NUCLEUS ROTATED EVERY 18 HOURS CHANGING ILLUMINATION BY THE SUN FROM ONE LOBE TO THE OTHER. WE WILL MODIFY OUR SPHERICAL HARMONIC BASED SURFACE ACTIVITY INVERSION TECHNIQUE TO USE IMAGE DATA. WE HAVE ALREADY DONE THIS FOR A SIMILAR INVERSION METHOD APPLIED TO BOTH CASSINI INMS IN SITU DENSITIES AND CASSINI UVS STELLAR OCCULTATION COLUMN DENSITIES. WE WILL THEN USE THE SPHERICAL HARMONICS INVERSION OF THE H2O AND CO2 POTENTIAL ACTIVITY TOGETHER WITH OUR COMET ADAPTIVE MESH SIMULATION (C-AMPS) DIRECT SIMULATION MONTE CARLO MODEL THAT NOW CAN RUN FULLY IN 3D WITH AN IRREGULAR NUCLEUS AND ANY SOLAR ILLUMINATION GEOMETRY TO DETERMINE THE FULL ACTIVITY DISTRIBUTION OF THE NUCLEUS AS IT CHANGES IN TIME. IN ADDITION TO C-AMPS WE HAVE A NEWLY DEVELOPED MULTINEUTRAL- GAS-DUST-FLUID MODEL THAT IS BASED ON THE U. MICHIGAN BATSRUS FLUID MODEL FOR FORWARD MODELING OF EXPLICIT TIME DEPENDENCE OF THE ROTATING NUCLEUS TO COMPARE EPOXI SPECTRAL IMAGES WITH OTHER REMOTE EARTH AND NEAR-EARTH BASED OBSERVATIONS. THE PROPOSED INVESTIGATION COMBINES EXPERTISE WITH THE EPOXI MEASUREMENTS AND DATA AT THE UNIVERSITY OF MARYLAND AND SOPHISTICATED PHYSICAL MODEL ANALYSIS AT THE UNIVERSITY OF MICHIGAN AND APPLIES THEM DIRECTLY TO SEVERAL OF THE FUNDAMENTAL GOALS OF THE DISCOVERY DATA ANALYSIS PROGRAM. DATA ARE AVAILABLE IN THE PDS SMALL BODIES NODE (SBN) AND IN FACT SCIENTISTS THAT ARE PART OF THE SBN AND WERE MEMBERS OF THE ORIGINAL EPOXI TEAM ARE INVESTIGATORS ON THIS PROPOSAL. THE RESULTS TO BE OBTAINED BY THIS PROPOSED WORK WOULD REPRESENT A SIGNIFICANT INTEREST TO THE DDAP AND TO NASA AND ITS PLANETARY PROGRAMS.