This work presents the second-order forward and adjoint sensitivity analysis procedures (SO-FSAP and SO-ASAP) for computing exactly and efficiently the second-order functional derivatives of physical (engineering, biological, etc.) system responses (i.e., "system performance parameters") to the system's model parameters. The definition of "system parameters" used in this work includes all computational input data, correlations, initial and/or boundary conditions, etc. For a physical system comprising N α parameters and r N responses, we note that the SO-FSAP requires a total of ( )computations for obtaining all of the first-and second-order sensitivities, for all r N system responses. On the other hand, the SO-ASAP requires a total of ( )large-scale computations for obtaining all of the first-and second-order sensitivities, for one functionaltype system responses. Therefore, the SO-FSAP should be used when r N N α , while the SO-ASAP should be used when r N N α . The original SO-ASAP presented in this work should enable the hitherto very difficult, if not intractable, exact computation of all of the second-order response sensitivities (i.e., functional Gateaux-derivatives) for large-systems involving many parameters, as usually encountered in practice. Very importantly, the implementation of the SO-ASAP requires very little additional effort beyond the construction of the adjoint sensitivity system needed for computing the first-order sensitivities.KEYWORDS: second-order adjoint sensitivity analysis procedure (SO-ASAP); exact and efficient computation of first-and second-order functional derivatives; large-scale systems.2