THE FORESHOCK REGION UPSTREAM FROM A PLANETARY BOW SHOCK (BS) CONTAINS A VARIETY OF WAVE-PARTICLE INTERACTION PHENOMENA AND IS CHARACTERIZED BY THE PRESENCE OF BACKSTREAMING IONS AND ELECTRONS. THESE BACKSTREAMING PARTICLES PROVIDE A SOURCE OF FREE ENERGY FOR PLASMA INSTABILITIES WHICH GENERATES PLASMA WAVES PROPAGATING ALONG THE INTERPLANETARY MAGNETIC FIELD (IMF) CONNECTED TO THE BS. THE WAVES CAN ALTER THE DISTRIBUTION OF THE BACKSTREAMING PARTICLES MODIFY THE BS STRUCTURE AND AFFECT MAGNETOSHEATH DYNAMICS AS THEY ARE CARRIED DOWNSTREAM BY THE SOLAR WIND. INITIAL MESSENGER STUDY HAD IDENTIFIED THE EXISTENCE OF AT LEAST THREE TYPES OF UPSTREAM ULF WAVES IN MERCURY S FORESHOCK WITH THE WHISTLER WAVES AT FREQUENCY NEAR 2HZ BEING MOST COMMONLY OBSERVED. RECENTLY A STATISTICAL STUDY OF 0.05 0.41HZ ULF WAVES (SIMILAR TO THE 30-S WAVES AT EARTH) FOUND THAT THESE MAGNETOSONIC WAVES HAVE SMALLER WAVE AMPLITUDES AND LOWER OCCURRENCE RATE THAN THEIR EARTH S COUNTERPART. HOWEVER ULF WAVES SIMILAR TO THE 3-S WAVES AT EARTH AND MAGNETOSONIC WAVES WITH STEEPENING EDGES (I.E. SHOCKLETS) HAVE NOT YET BEEN OBSERVED IN MERCURY S FORESHOCK. TO DATE MANY QUESTIONS REGARDING THE PROPERTIES AND GENERATION MECHANISMS OF FORESHOCK WAVES AND THEIR INTERACTION WITH MERCURY'S BS AND MAGNETOSHEATH REMAIN UNANSWERED. WE PROPOSE TO USE THE ENTIRE MESSENGER DATASET AND UNIVERSITY OF MICHIGAN GLOBAL MHD WITH EMBEDDED PARTICLE-IN-CELL (MHD-EPIC) MODEL TO ANSWER KEY SCIENCE QUESTIONS ON THE CHARACTERISTICS OF UPSTREAM ULF WAVES THEIR INTERACTION WITH THE PLASMA ENVIRONMENT SURROUNDING THE BS AND THE FUNDAMENTAL PHYSICS BEHIND THEIR GENERATION MECHANISMS IN MERCURY S FORESHOCK. CORRELATION ANALYSIS WILL BE CONDUCTED TO DETERMINE THE RELATIONSHIP BETWEEN THE WAVE PROPERTIES SOLAR WIND AND IMF CONDITIONS AND SPACECRAFT LOCATION IN THE FORESHOCK. THE MHDEPIC MODEL WILL BE USED FOR THE FIRST TIME IN THIS STUDY TO INVESTIGATE THE ION AND ELECTRON KINETIC PHYSICS OF THE WAVE-PARTICLE INTERACTIONS AND WAVE GENERATION MECHANISMS IN MERCURY S FORESHOCK BY SIMULATING THE BS AND ITS FORESHOCK WITH A FULLY KINETIC CODE. THE PROPOSED WORK WILL GREATLY ENHANCE THE SCIENCE RETURN FROM MESSENGER IN THESE IMPORTANT AREAS AND IT IS HIGHLY RELEVANT TO THE ROSES-2020 DDAP.