In this Problem Set we will modify the 2D problem we studied in Problem Set #2 and explore transport solutions in three dimensions. We will also use three different source terms and perform a remediation simulation.PART 2- Chemical ReactionsGEOLOGY 727 Problem Set #3—Source Terms& Chemical ReactionsIn this Problem Set we will modify the 2D problem we studied in Problem Set #2 and explore transport solutions in three dimensions. We will also use three different source terms and perform a remediation simulation.Problem DesignAdd a layer 10 m thick to the system you modeled in Problem Set #2 so that we can study transport in a three-dimensional system. Keep the constant head boundaries and the injection well in the top layer but move the pumping well (and the observation point) from the top layer to the bottom layer (layer 2, row 11,column 7). Assume all no flow boundary conditions in the bottom layer for the flow problem and zero mass flux boundary conditions for the transport problem. Hydraulic conductivity and porosity in the bottom layer (layer 2) are the same as in the top layer (i.e., K = 1.474 x 10-4 m/s and  = 0.3), except thereis no low K zone in layer 2. Dispersivities are the same as in Problem Set #2 (L = 20 m; Th = 4 m) but now use a transverse vertical dispersivity equal to 0.04 m. Dispersivities are the same in both layers of the model. Designate both layers of the model as “confined” in the BCF package in MODFLOW (and be sure that the “compute Leakance” box is checked). Designating the upper layer as confined means that MODFLOW will use a constant value of transmissivity for the layer (where transmissivity is equal to hydraulic conductivity times the thickness of the layer and thickness is 10 m). Use the TVD solution option with the GCG solver for all simulations in this problem set.PART 1 – Source Termsa) Injection Well. Change the concentration of the solute being injected via the injection well to 200 ppm (gm/m3). Run the model and produce a breakthrough curve showing concentrations at the observation point (i.e., the pumping well in layer 2, row 11, column 7) for 1 year of transport time.b) Constant Concentration Cell. Replace the injection well by a constant concentration cell equal to 200 ppm. Allow the cell to be an active head cell so that the model will calculate the head in the constant concentration cell (i.e, in the “BC>constant head/conc” window, you should check the box labeled “concentration value is constant” but do not check the box labeled “head value is constant”). Run the model and produce a breakthrough curve showing concentrations at the observation point (the pumping well) for 1 year of transport time.c) Recharge. Modify the model to allow the upper layer to be an unconfined layer (set in the BCF package in MODFLOW) and check the box at the bottom of the BCF window “update VCONT when transmissivity is updated for unconfined layers”. Assume that there is a uniform recharge rate of 0.25 m/yr over the entire top layer (entered in the Recharge properties window). Delete theconstant concentration cell and assume that there is an additional recharge equal to 3.1536 m/yr added to this cell. Also assume that the composite recharge to this cell has a concentration of 200 ppm (also entered in the Recharge properties window). It will be helpful to change the top elevation of layer 1 from 20 m to 250 m in the “top elevation” properties window. (Layer thickness is used in the leakance calculation and is also used to configure the display of the cross section on the main screen in GWV, where you will see the location of the water table.) Run the model and produce a breakthrough curve showing concentrations at the observation point (the pumping well) for 1 year of transport time.d) Plot these curves and discuss and explain the differences in the three breakthrough curves .PART 2- Chemical ReactionsStart with the model considered in Part 1.b above (i.e., the constant concentration cell). We will now add chemical reactions. Enter the appropriate parameters in the chemical reaction package window and the chemical reactions properties menu in GWV.a) Re-run the model using a linear adsorption isotherm with R = 2.b) Re-run the model without adsorption but assume a half-life of 0.1 year.c) Re-run the model with R=2 and a half-life of 0.1 year.Produce breakthrough curves at the observation point (i.e., the pumping well) for all three cases. The relation among the curves should be consistent with the curves from Problem Set #1. Plot these three curves and the curve from Part 1.b (no chemical reactions) on one graph.PART 3- RemediationStart with the model considered in Part 1.b above (constant concentration cell without chemical reactions). Suppose that after 1 year, we remove the source. To simulate this we will set the constant concentration cell to zero after 1 year and allow the mass present in the aquifer to attenuate. The aquifer isremediated by natural attenuation and by continued pumping from the well in layer 2, row 11, column 7. The flow field remains the same (at steady state) since we continue to run the pumping well throughout thesimulation.Even though the flow field is at steady state throughout the simulation, we need to set up two stress periods to handle the change in concentration of the source cell. In stress period 1 we will re-create the solution generated in Part 1.b; in stress period 2 we will turn off the source and perform the remediation. Stress period one is one year long and stress period 2 is two years long, for a total simulation time of 3 years. In MODFLOW, use the “basic package” and the “stress periods setup” windows to enter information for the two stress periods. Designate two stress periods for a steady state simulation. In MT3D, use the BC>constant head/conc window to enter transient data for the constant concentration cell. The concentration in the cell is equal to 200 ppm for stress period 1 and equal to 0 for stress period 2. Add a second observation point in layer 1, row 7, column 7 and run the model. Import the *.obs file into Excel (or some other spreadsheet) and plot concentration vs. time at the two observation points or use the Plot>hydrograph option in GWV and then right click on the mouse and use the pop up menu to print the breakthrough curve. You should plot the concentration history from the start of injection through the two year